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.

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.

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.

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.

Schizophrenia (SZ) is a severe, complex, and common mental disorder with high heritability (80%), an adult age of onset, and high discordance (∼50%) in monozygotic twins (MZ). Extensive studies on familial and non-familial cases have implicated a number of segregating mutations and de novo changes in SZ that may include changes to the mitochondrial genome. Yet, no single universally causal variant has been identified, highlighting its extensive genetic heterogeneity. This report specifically focuses on the assessment of changes in the mitochondrial genome in a unique set of monozygotic twins discordant (MZD) for SZ using blood. Genomic DNA from six pairs of MZD twins and two sets of parents (N = 16) was hybridized to the Affymetrix Human SNP Array 6.0 to assess mitochondrial DNA copy number (mtDNA-CN). Whole genome sequencing (WGS) and quantitative polymerase chain reaction (qPCR) was performed for a subset of MZD pairs and their parents and was also used to derive mtDNA-CN estimates. The WGS data were further analyzed to generate heteroplasmy (HP) estimates. Our results show that mtDNA-CN estimates for within-pair and mother-child differences were smaller than comparisons involving unrelated individuals, as expected. MZD twins showed discordance in mtDNA-CN estimates and displayed concordance in directionality of differences for mtDNA-CN across all technologies. Further, qPCR performed better than Affymetrix in estimating mtDNA-CN based on relatedness. No reliable differences in HP were detected between MZD twins. The within-MZD differences in mtDNA-CN observed represent postzygotic somatic changes that may contribute to discordance of MZ twins for diseases, including SZ.

Astigmatism and myopia are two common ocular refractive errors that can impact daily life, including learning and productivity. Current knowledge suggests that the etiology of these conditions is the result of a complex interplay between genetic and environmental factors. Studies in populations of European ancestry have demonstrated a higher concordance of refractive errors in monozygotic (MZ) twins compared to dizygotic (DZ) twins. However, there is a lack of studies on genetically informative samples of multi-ethnic ancestry. This study aimed to estimate the genetic contribution to astigmatism and myopia in the Mexican population. A sample of 1399 families, including 243 twin pairs and 1156 single twins, completed a medical questionnaire about their own and their co-twin’s diagnosis of astigmatism and myopia. Concordance rates for astigmatism and myopia were estimated, and heritability and genetic correlations were determined using a bivariate ACE Cholesky decomposition method, decomposed into A (additive genetic), C (shared environmental) and E (unique environmental) components. The results showed a higher concordance rate for astigmatism and myopia for MZ twins (.74 and .74, respectively) than for DZ twins (.50 and .55). The AE model, instead of the ACE model, best fitted the data. Based on this, heritability estimates were .81 for astigmatism and .81 for myopia, with a cross-trait genetic correlation of rA = .80, nonshared environmental correlation rE = .89, and a phenotypic correlation of rP = .80. These results are consistent with previous findings in other populations, providing evidence for a similar genetic architecture of these conditions in the multi-ethnic Mexican population.

The purpose of this research was to determine the frequency of mutation of the cytochrome CYP3A5 genes and transport proteins SLCO1B1 and MDR1 in patients with coronary heart disease in the Kazakh nation. A prospective cohort clinical and genetic study was conducted. The study was conducted in 2017−2019. Medical records containing information about drug prescription conducted in hospitals and outpatient departments were carefully analyzed. In the examined group of 178 patients treated with statins, a significant frequency of genetic variants that determine the increased risk of complications of statin use was revealed. There was a tendency toward an increase in the activity of creatine phosphokinase (CPK) in the blood upon detection of the A6986G mutation of the cytochrome gene and SLCO1B1 (c.521T>C) gene of the transport protein OATP1B1. In the studied Kazakh population, the presence of a homozygous mutant SLCO1B1 gene of the transport protein can be recommended as a genetic marker for the undesirability of using antihypercholesterolemic therapy with statins, which simultaneously leads to a decrease in the effectiveness of treatment and an increase in the risk of side effects.

Fetal sex contributes to the determination of obstetric outcome, as pregnancies carrying male babies seem to have an increased risk of maternal-fetal complications. Most studies have been conducted on singleton pregnancies, whereas less evidence is available for twins. A 10-year retrospective observational study was conducted on a cohort of 1180 women with twin pregnancy delivered at a single tertiary hospital. Clinical data on maternal characteristics, and obstetric and neonatal outcomes were collected, and the analysis was performed on monochorionic (MC) and dichorionic (DC) diamniotic twins separately. The group of DC twins included 837 cases, and those conceived by assisted reproductive technologies (ART) were more likely to have one or both female fetuses rather than males. The incidence of hypertensive disorders of pregnancy (HDP) was higher in same-sex pairs than in opposite-sex pairs. No differences were found regarding other obstetric and neonatal outcomes among the three sex-pairing groups. The MC twins group included 228 cases, and in female-carrying pregnancies a higher incidence of gestational diabetes (GDM) was observed compared to the male group. Furthermore, male pairs had significantly lower Apgar scores than females. Fetal sex seems to have a mild effect in twins compared to singleton pregnancies, suggesting a more complex set of factors contributing to pregnancy outcome in multiple pregnancies. However, we observed a higher incidence of HDP among same-sex DC pairs, a higher rate of GDM among MC female-female pairs, and a worse adaptation to extrauterine life among male-male pairs in MC twins.

The aim of this study was to investigate the prevalence of anemia in twin pregnancies and the influence of anemia on maternal and neonatal outcomes. This retrospective study included twin pregnant women who delivered in a tertiary hospital in China from January 2018 to December 2018. Patients were divided by WHO criteria (hemoglobin <11.0 g/dL): the anemic and nonanemic groups. Patients with anemia were further classified as recovered or unrecovered subgroup after oral iron therapy. Maternal and neonatal outcomes in women carrying twins were compared using Student’s t test and the chi-squared test or the Fisher exact test. Univariable and multivariable logistic regression models were used to determine the association of maternal and neonatal characteristics with anemia. Linear regression analysis was used to estimate mean birth weight and gestational week. The prevalence of anemia was 42.6% (182/427) in twin pregnancies. The anemic group had higher rates of low 1-minute Apgar score (4.4% vs. 1.8%, p = .028), perinatal death (1.9% vs. 0.2%, p = .012) and neonatal intensive care unit (NICU) admission (27.2% vs. 20.2%, p = .017; adjusted OR, 1.478; 95% CI [1.07, 2.044]). The recovered subgroup had lower NICU admission rate (13.5% vs. 30.3%, p = .006; OR, 0.388; 95% CI [0.186, 0.809]), higher gestational week and birth weight (β, 0.954 week; 95% CI [0.114, 1.794] and β, 171.01 g; 95% CI [9.894, 332.126] respectively). The prevalence of anemia in twin gestation is high. Anemia is associated with adverse neonatal outcomes, and correction of anemia significantly improved the pregnancy outcomes.

The 18th International Society for Twin Studies convention took place between June 15−17, 2023, in Budapest, Hungary. A selective sampling of highlights from the meeting are presented. This is followed by a brief overview of ongoing twin research at the University of Aberystwyth in Wales and some twin treasures in the National Library. Next, reviews of timely research that examine temperamental similarities and twin relations, the maximization of twin research designs, a rare ectopic twin pregnancy, and twins’ academic self-concept formation are presented. The final section covers stories appearing in various media sources that are informative, enlightening, and interesting, namely a reunion of identical infant twins, the birth of identical-fraternal quadruplets, identical twin comedians, meeting an unrelated look-alike, and twins in the Tour de France.