Joseph has written what purports to be a refutation of studies of Twins Reared-Apart (TRAs) with a singular focus on the Minnesota Study of Twins Reared-Apart (MISTRA). I show, in detail, that (a) his criticisms of previous TRA studies depend on sources that were discredited prior to MISTRA, as they all failed the test of replicability, (b) the list of biases he uses to invalidate MISTRA do not support his arguments, (c) the accusations of questionable research practices are unsubstantiated, (d) his claim that MISTRA should be evaluated in the context of psychology’s replication crisis is refuted. The TRA studies are constructive replications. Like many other scholars, past and present, he has been misled by the variation introduced by small samples (sampling error) and the distortion created by walking in the garden of forking paths. His endeavor is a concatenation of elision and erroneous statistical/scientific reasoning.

Estimated heritability of educational attainment (EA) varies widely, from 23% to 80%, with growing evidence suggesting the degree to which genetic variation contributes to individual differences in EA is highly dependent upon situational factors. We aimed to decompose EA into influences attributable to genetic propensity and to environmental context and their interplay, while considering influences of rearing household economic status (HES) and sex. We use the Project Talent Twin and Sibling Study, drawn from the population-representative cohort of high school students assessed in 1960 and followed through 2014, to ages 68−72. Data from 3552 twins and siblings from 1741 families were analyzed using multilevel regression and multiple group structural equation models. Individuals from less-advantaged backgrounds had lower EA and less variation. Genetic variance accounted for 51% of the total variance, but within women and men, 40% and 58% of the total variance respectively. Men had stable genetic variance on EA across all HES strata, whereas high HES women showed the same level of genetic influence as men, and lower HES women had constrained genetic influence on EA. Unexpectedly, middle HES women showed the largest constraints in genetic influence on EA. Shared family environment appears to make an outsized contribution to greater variability for women in this middle stratum and whether they pursue more EA. Implications are that without considering early life opportunity, genetic studies on education may mischaracterize sex differences because education reflects different degrees of genetic and environmental influences for women and men.

Establishing causal relationships in observational studies is an important step in research and policy decision making. The association between an exposure and an outcome can be confounded by multiple factors, often making it hard to draw causal conclusions. The co-twin control design (CTCD) is a powerful approach that allows for the investigation of causal effects while controlling for genetic and shared environmental confounding factors. This article introduces the CTCD and offers an overview of analysis methods for binary and continuous outcome and exposure variables. Tools for data simulation are provided, along with practical guidance and accompanying scripts for implementing the CTCD in R, SPSS, and Stata. While the CTCD offers valuable insights into causal inference, it depends on several assumptions that are important when interpreting CTCD results. By presenting a broad overview of the CTCD, this article aims to equip researchers with actionable recommendations and a comprehensive understanding of the design’s strengths and limitations.

Female fertility is a complex trait with age-specific changes in spontaneous dizygotic (DZ) twinning and fertility. To elucidate factors regulating female fertility and infertility, we conducted a genome-wide association study (GWAS) on mothers of spontaneous DZ twins (MoDZT) versus controls (3273 cases, 24,009 controls). This is a follow-up study to the Australia/New Zealand (ANZ) component of that previously reported (Mbarek et al., 2016), with a sample size almost twice that of the entire discovery sample meta-analysed in the previous article (and five times the ANZ contribution to that), resulting from newly available additional genotyping and representing a significant increase in power. We compare analyses with and without male controls and show unequivocally that it is better to include male controls who have been screened for recent family history, than to use only female controls. Results from the SNP based GWAS identified four genomewide significant signals, including one novel region, ZFPM1 (Zinc Finger Protein, FOG Family Member 1), on chromosome 16. Previous signals near FSHB (Follicle Stimulating Hormone beta subunit) and SMAD3 (SMAD Family Member 3) were also replicated (Mbarek et al., 2016). We also ran the GWAS with a dominance model that identified a further locus ADRB2 on chr 5. These results have been contributed to the International Twinning Genetics Consortium for inclusion in the next GWAS meta-analysis (Mbarek et al., in press).

In studies of singletons, a range of early-life characteristics have been reported to be associated with handedness, but some of these associations have failed to replicate. We examined associations between 23 early life characteristics with handedness in a large sample of 37,495 5-year-old twins. We considered three definitions of handedness: left-handedness (LH), mixed-handedness (MH), and non-right-handedness (NRH). Our main aim was to test whether the associations with sex, birth weight, gestational age, and season of birth — as reported in singletons — replicate in twins, and to examine twin-specific variables, including zygosity, chorionicity, birth order, and intertwin delivery time. Compared to previously published data from adults born as singletons (7.23%), the prevalence of NRH was higher in both twins (16.19%) and their parents (15.09%). In the twins, LH and NRH were associated with parents’ LH. Male sex and lower gestational age were associated with NRH, and LH was associated with not being breastfed. MH was related to neurodevelopmental delays and higher externalizing problems later in childhood. Other previously reported associations were not replicated, and no twin-specific characteristics were related to handedness. These results emphasize the importance of considering multiple definitions of handedness and indicate a small number of replicated associations across studies.

Reading difficulties are prevalent worldwide, including in economically developed countries, and are associated with low academic achievement and unemployment. Longitudinal studies have identified several early childhood predictors of reading ability, but studies frequently lack genotype data that would enable testing of predictors with heritable influences. The National Child Development Study (NCDS) is a UK birth cohort study containing direct reading skill variables at every data collection wave from age 7 years through to adulthood with a subsample (final n = 6431) for whom modern genotype data are available. It is one of the longest running UK cohort studies for which genotyped data are currently available and is a rich dataset with excellent potential for future phenotypic and gene-by-environment interaction studies in reading. Here, we carry out imputation of the genotype data to the Haplotype Reference Panel, an updated reference panel that offers greater imputation quality. Guiding phenotype choice, we report a principal components analysis of nine reading variables, yielding a composite measure of reading ability in the genotyped sample. We include recommendations for use of composite scores and the most reliable variables for use during childhood when conducting longitudinal, genetically sensitive analyses of reading ability.

Self-perception in early childhood and self-esteem in adulthood are related to a variety of aspects of psychological wellbeing. The goal of the present study was to examine genetic and familial influences on self-perception and self-esteem in separate samples of children (153 twin pairs of 5-year-olds) and adults (753 twin pairs between the ages of 25–75 years). Genetic common factor modeling showed that three facets of self-perception (physical competence, peer acceptance, and maternal acceptance) loaded onto a single heritable factor in children. Multilevel modeling showed no effects of self or co-twin sex on self-perception, but authoritative parenting style was negatively related to self-perception in boys. Similarly, in Study 2, with the adult sample, five self-esteem items loaded on a single heritable factor with no effects of co-twin sex on adult self-esteem. Remembered maternal affection, paternal affection, and maternal discipline were positively related to self-esteem in adults; maternal affection was especially significant for women. The reversal in direction of parenting effects between early childhood and adulthood suggests that parents may play different roles in shaping how children and adults think of themselves. These results suggest that self-perception in childhood and self-esteem in adulthood are both influenced by genetic and environmental factors and that parenting is an important environmental factor for both children and adults.

Individual differences in educational attainment (EA) and physical health, as indexed by body mass index (BMI), are correlated within individuals and across generations. The aim of our study was to assess the transmission of these traits from parents to their offspring in childhood and adolescence. We analyzed BMI and EA in 13,916 families from the Netherlands. Data were available for 27,577 parents (mean age 33) and 26,855 of their offspring at 4 and 12 years of age. We employed structural equation modeling to simultaneously estimate the phenotypic transmission of BMI and EA from parents to offspring, the spousal correlations, and the residual child BMI-EA associations after accounting for intergenerational transmission and testing for gender differences therein. We found a significant intergenerational transmission of BMI to BMI in childhood (age 4; standardized regression coefficient β = .10) and adolescence (age 12; β = .20), and of EA to academic achievement in adolescence (β = .19). Cross-trait parent-to-offspring transmission was weak. All transmission effects were independent of parent or offspring gender. We observed within-person EA-BMI correlations that were negative in parents (∼−.09), positive in children (∼.05) and negative in adolescents (∼−.06). Residual EA-BMI were positive in children (∼.05) and insignificant in adolescents. Spousal correlations were .46 for EA, .21 for BMI, and ∼−.09 cross-trait. After accounting for spousal correlations, the intergenerational transmission for BMI and EA is mainly predictive within, but not across, traits. The within-person correlation between BMI and EA can change in direction between childhood and adulthood.

The current study explored the impact of genetic relatedness differences (ΔH) and sample size on the performance of nonclassical ACE models, with a focus on same-sex and opposite-sex twin groups. The ACE model is a statistical model that posits that additive genetic factors (A), common environmental factors (C), and specific (or nonshared) environmental factors plus measurement error (E) account for individual differences in a phenotype. By extending Visscher’s (2004) least squares paradigm and conducting simulations, we illustrated how genetic relatedness of same-sex twins (HSS) influences the statistical power of additive genetic estimates (A), AIC-based model performance, and the frequency of negative estimates. We found that larger HSS and increased sample sizes were positively associated with increased power to detect additive genetic components and improved model performance, and reduction of negative estimates. We also found that the common solution of fixing the common environment correlation for sex-limited effects to .95 caused slightly worse model performance under most circumstances. Further, negative estimates were shown to be possible and were not always indicative of a failed model, but rather, they sometimes pointed to low power or model misspecification. Researchers using kin pairs with ΔH less than .5 should carefully consider performance implications and conduct comprehensive power analyses. Our findings provide valuable insights and practical guidelines for those working with nontwin kin pairs or situations where zygosity is unavailable, as well as areas for future research.

The recruitment of participants for research studies may be subject to bias. The Prospective Imaging Study of Ageing (PISA) aims to characterize the phenotype and natural history of healthy adult Australians at high future risk of Alzheimer’s disease (AD). Participants approached to take part in PISA were selected from existing cohort studies with available genomewide genetic data for both successfully and unsuccessfully recruited participants, allowing us to investigate the genetic contribution to voluntary recruitment, including the genetic predisposition to AD. We use a polygenic risk score (PRS) approach to test to what extent the genetic risk for AD, and related risk factors predict participation in PISA. We did not identify a significant association of genetic risk for AD with study participation, but we did identify significant associations with PRS for key causal risk factors for AD, IQ, household income and years of education. We also found that older and female participants were more likely to take part in the study. Our findings highlight the importance of considering bias in key risk factors for AD in the recruitment of individuals for cohort studies.

Phenotypic variation is the result of gene expression based on complex interaction between genetic and environmental factors. It is well known that genetic and environmental factors influence gene expression, but our understanding of their relative importance remains limited. To obtain a hint for the understanding of their contributions, we took advantage of monozygotic twins, as they share genetic and shared environmental factors but differ in nonshared factors, such as environmental differences and stochastic factors. In this study, we performed cap analysis of gene expression on three pairs of twins and clustered each individual based on their expression profiles of annotated genes. The dendrogram of annotated gene transcripts showed a monophyletic clade for each twin pair. We also analyzed the expression of retrotransposons, such as human endogenous retroviruses (HERVs) and long interspersed nuclear elements (LINEs), given their abundance in the genome. Clustering analyses demonstrated that HERV and LINE expression diverged even within monozygotic twin pairs. Thus, HERVs and LINEs are more susceptible to nonshared factors than annotated genes. Motif analysis of differentially expressed annotated genes suggests that specificity protein/Krüppel-like factor family transcription factors are involved in the expression divergence of annotated gene influenced by nonshared factors. Collectively, our findings suggest that expressions of annotated genes and retrotransposons are differently regulated, and that the expression of retrotransposons is more susceptible to nonshared factors than annotated genes.

Prior research based on Swedish data suggests that collective optimism, as measured by monthly incidence of suicides, correlates inversely with selection in utero against male twins in a population. We test this finding in the US, which reports the highest suicide rate of all high-income countries, and examine whether monthly changes in overall suicides precede changes in the ratio of male twin to male singleton live births. Consistent with prior work, we also examine as a key independent variable, suicides among women aged 15−49 years. We retrieved monthly data on suicides and the ratio of male twin to singleton live births from CDC WONDER, 2003 to 2019, and applied Box-Jenkins iterative time-series routines to detect and remove autocorrelation from both series. Results indicate that a 1% increase in monthly change in overall suicides precedes a 0.005 unit decline in male twin live births ratio 6 months later (coefficient = −.005, p value = .004). Results remain robust to use of suicides among reproductive-aged women as the independent variable (coefficient = −.0012, p value = .014). Our study lends external validity to prior research and supports the notion that a decline in collective optimism corresponds with greater selection in utero.

Monozygotic twins (MZT) are 2.5 times more frequent in ART than in natural conceptions. A number of ART-related mechanisms have been probably linked with MZT. Studies that retrospectively analyze the time-lapse (TL) records resulting in MZT suggest that some morphokinetic traits of the inner cell mass and the trophectoderm could be predictors of MZT, but results are controversial. We present the complete TL record of one case of MZT that split itself at the very moment of the division into two cells, with one of the cells coming out through a hole in the zona pellucida (ZP). Both resulting embryos developed normally, and were vitrified. It is suggested that the hole in the ZP may facilitate the extrusion of some cells of the <day 4 embryo and that this cell development is not constrained by being inside the ZP. Despite the lack of the inhibition of the ZP itself or the influence of the other embryo cells, the totipotent cell was then able to develop correctly from the start. Moreover, the embryo inside the ZP compensated for the loss of this cell apparently without problems. Our findings are discussed in the context of previous literature and ethical problems are addressed.

A model in the form of a Markov chain is constructed to mimic variations in the human sex ratio. It is illustrated by simulation. The equilibrium distribution is shown to be a simple modification of the binomial distribution. This enables an easy calculation of the variation in sex ratio which could be expected in small populations.

Parents’ alcohol use is associated with alcohol use of their adolescent offspring, but does this association extend to the adulthood of the offspring? We examined associations of paternal and maternal problem drinking with lifetime problem drinking of their adult offspring prospectively assessed in a population-based Finnish twin-family cohort (FinnTwin16). Problem drinking (Malmö-modified Michigan Alcoholism Screening Test) was self-reported separately by mothers and fathers when their children were 16. The children reported on an extended lifetime version of the same measure during their mid-twenties (21-28 years) and mid-thirties (31-37 years). 1235 sons and 1461 daughters in mid-twenties and 991 sons and 1278 daughters in mid-thirties had complete data. Correlations between fathers’ and their adult children’s problem drinking ranged from .12 to .18. For mothers and their adult children, these correlations ranged from .09 to .14. In multivariate models, adjustment for potential confounders had little effect on the observed associations. In this study, parental problem drinking was modestly associated with lifetime problem drinking of their adult children. This association could be detected even when the children had reached the fourth decade of life.

Almost all creatinine is excreted by the kidney in individuals. Serum creatinine concentration, a widely used renal function index in clinical practice, can be affected by both genetic and environmental factors, as evidenced by current research exploring the relationship between these factors and kidney function. However, few studies have explored the heritability of serum creatinine in Asian populations. Therefore, we explored the genetic and environmental factors that affect the serum creatinine level in Asian populations. Participants in this study came from the Qingdao Twin Registry in China, and 374 pairs of twins were included, of which 139 pairs were dizygotic twins, whose ages ranged from 40 to 80 years old, and the serum creatinine level ranged from 10 to 126 μmol/L. Structural equation models were constructed using Mx software to calculate heritability, with adjusted covariates being age, sex, and body mass index. The results of heritability analysis showed that ACE was the best fit model. Serum creatinine level is influenced by genetic and environmental factors. The result of heritability was 35.44%, and the influence of shared environmental factors accounted for 52.13%. This study provided the relevant basis for future research on genetic and environmental factors affecting serum creatinine levels in Asian populations.

A simple model by which Hardy-Weinberg proportions are attained in a single generation while maintaining gene frequencies is stated and illustrated. The title ‘Quasi-random mating’ is proposed. Confusion about the Hardy-Weinberg principle can be avoided only if there is clear separation between the basic deterministic model and factors influencing a population’s structure. Eighty years passed before C. C. Li coined the term ‘pseudo-random mating’. The lesson taught by Li has not been taken on board.

Parkinson’s disease (PD) is a complex disorder with a significant genetic component. Genetic variations associated with PD play a crucial role in the disease’s inheritance and prognosis. Currently, 31 genes have been linked to PD in the OMIM database, and the number of genes and genetic variations identified is steadily increasing. To establish a robust correlation between phenotype and genotype, it is essential to compare research findings with existing literature. In this study, we aimed to identify genetic variants associated with PD using a targeted gene panel with next-generation sequencing (NGS) technology. Our objective was also to explore the idea of re-analyzing genetic variants of unknown significance (VUS). We screened 18 genes known to be related to PD using NGS in 43 patients who visited our outpatient clinic between 2018−2019. After 12−24 months, we re-evaluated the detected variants. We found 14 different heterozygous variants classified as pathogenic, likely pathogenic, or VUS in 14 individuals from nonconsanguineous families. We re-evaluated 15 variants and found changes in their interpretation. Targeted gene panel analysis with NGS can help identify genetic variants associated with PD with confidence. Re-analyzing certain variants at specific time intervals can be especially beneficial in selected situations. Our study aims to expand the clinical and genetic understanding of PD and emphasizes the importance of re-analysis.

Dysmenorrhea is associated with epilepsy. Existing evidence is mostly limited to observational studies, which are liable to confounding and bias. This study investigated the causal relevance of dysmenorrhea on epilepsy using Mendelian randomization (MR). We extracted instrumental variants for dysmenorrhea and epilepsy from published genomewide association study data, focusing on individuals of East Asian descent. A comprehensive suite of MR estimations and sensitivity analyses was performed to ensure the robustness of the findings. Each outcome database was analyzed separately in both directions. For dysmenorrhea and epilepsy, 7 and 3 genetic variants respectively were selectively extracted as instrumental variants. The results suggest that dysmenorrhea is causally associated with an elevated risk of epilepsy (inverse variance weighted [IVW]: OR = 1.26; 95% CI [1.07, 1.47]; p = 4.42 × 10−3); conversely, no strong evidence was found to corroborate that epilepsy exerts a causal effect on the incidence of dysmenorrhea (IVW: OR = 1.04; 95% CI [0.82, 1.33]; p = .72). These findings provide novel insights into the causal relationship between dysmenorrhea and epilepsy, which may have implications for clinical decision-making in patients with epilepsy and dysmenorrhea.

An abnormal alanine aminotransferase (ALT) level is predictive of disease and all-cause mortality and may indicate liver injury. Using twin modeling, the genetic and environmental factors that affect human serum ALT levels have been well studied for the populations in the different countries, and the results showed moderate-to-high heritability. However, the heritability of ALT level has not been explored in Chinese population. Thus, we recruited 369 pairs of twins (233 monozygotic and 136 dizygotic) from the Qingdao Twin Registry in China with a median age of 50 years (40−80 years). Correlation analysis and a structural equation model (SEM) were conducted to evaluate the heritability of ALT level. The data for age, gender, body mass index and alcohol consumption were set as covariates. Intrapair correlation in monozygotic twins was 0.64 (95%CI [.56, .71]) and 0.42 (95% CI [.28, .55]) in dizygotic twins. The SEM analysis indicated that 65% (95% CI [57%, 71%]) of the variation in ALT levels can be explained by additive genetics and 35% (95% CI [29%, 44%]) of the variation is attributed to unique environmental factors or residuals. Shared environmental influences were not significant. In conclusion, serum ALT variations exhibited strong genetic effects. The variation could also be explained by unique environmental factors. However, shared environmental factors have a minor impact on the serum ALT level.