Phase analysis of circadian-related genes in two tissues

BACKGROUND: Recent circadian clock studies using gene expression microarray in two different tissues of mouse have revealed not all circadian-related genes are synchronized in phase or peak expression times across tissues in vivo. Instead, some circadian-related genes may be delayed by 4–8 hrs in peak expression in one tissue relative to the other. These interesting biological observations prompt a statistical question regarding how to distinguish the synchronized genes from genes that are systematically lagged in phase/peak expression time across two tissues. RESULTS: We propose a set of techniques from circular statistics to analyze phase angles of circadian-related genes in two tissues. We first estimate the phases of a cycling gene separately in each tissue, which are then used to estimate the paired angular difference of the phase angles of the gene in the two tissues. These differences are modeled as a mixture of two von Mises distributions which enables us to cluster genes into two groups; one group having synchronized transcripts with the same phase in the two tissues, the other containing transcripts with a discrepancy in phase between the two tissues. For each cluster of genes we assess the association of phases across the tissue types using circular-circular regression. We also develop a bootstrap methodology based on a circular-circular regression model to evaluate the improvement in fit provided by allowing two components versus a one-component von-Mises model. CONCLUSION: We applied our proposed methodologies to the circadian-related genes common to heart and liver tissues in Storch et al. [2], and found that an estimated 80% of circadian-related transcripts common to heart and liver tissues were synchronized in phase, and the other 20% of transcripts were lagged about 8 hours in liver relative to heart. The bootstrap p-value for being one cluster is 0.063, which suggests the possibility of two clusters. Our methodologies can be extended to analyze peak expression times of circadian-related genes across more than two tissues, for example, kidney, heart, liver, and the suprachiasmatic nuclei (SCN) of the hypothalamus

Learning about the X from our parents

The X chromosome is generally understudied in association studies, in part because the analyst has had limited methodological options. For nuclear-family-based association studies, most current methods extend the transmission disequilibrium test (TDT) to the X chromosome. We present a new method to study association in case-parent triads: the parent-informed likelihood ratio test for the X chromosome (PIX-LRT). Our method enables estimation of relative risks and takes advantage of parental genotype information and the sex of the affected offspring to increase statistical power to detect an effect. Under a parental exchangeability assumption for the X, if case-parent triads are complete, the parents of affected offspring provide an independent replication sample for estimates based on transmission distortion to their affected offspring. For each offspring sex we combine the parent-level and the offspring-level information to form a likelihood ratio test statistic; we then combine the two to form a combined test statistic. Our method can estimate relative risks under different modes of inheritance or a more general co-dominant model. In triads with missing parental genotypes, the method accounts for missingness with the Expectation-Maximization algorithm. We calculate non-centrality parameters to assess the power gain and robustness of our method compared to alternative methods. We apply PIX-LRT to publically available data from an international consortium of genotyped families affected by the birth defect oral cleft and find a strong, internally-replicated signal for a SNP marker related to cleft lip with or without cleft palate

Learning about the X from our parents

The X chromosome is generally understudied in association studies, in part because the analyst has had limited methodological options. For nuclear-family-based association studies, most current methods extend the transmission disequilibrium test (TDT) to the X chromosome. We present a new method to study association in case-parent triads: the parent-informed likelihood ratio test for the X chromosome (PIX-LRT). Our method enables estimation of relative risks and takes advantage of parental genotype information and the sex of the affected offspring to increase statistical power to detect an effect. Under a parental exchangeability assumption for the X, if case-parent triads are complete, the parents of affected offspring provide an independent replication sample for estimates based on transmission distortion to their affected offspring. For each offspring sex we combine the parent-level and the offspring-level information to form a likelihood ratio test statistic; we then combine the two to form a combined test statistic. Our method can estimate relative risks under different modes of inheritance or a more general co-dominant model. In triads with missing parental genotypes, the method accounts for missingness with the Expectation-Maximization algorithm. We calculate non-centrality parameters to assess the power gain and robustness of our method compared to alternative methods. We apply PIX-LRT to publically available data from an international consortium of genotyped families affected by the birth defect oral cleft and find a strong, internally-replicated signal for a SNP marker related to cleft lip with or without cleft palate

"Toward a Clearer Definition of Confounding" Revisited With Directed Acyclic Graphs

In a 1993 paper (Am J Epidemiol. 1993;137(1):1–8), Weinberg considered whether a variable that is associated with the outcome and is affected by exposure but is not an intermediate variable between exposure and outcome should be considered a confounder in etiologic studies. As an example, she examined the common practice of adjusting for history of spontaneous abortion when estimating the effect of an exposure on the risk of spontaneous abortion. She showed algebraically that such an adjustment could substantially bias the results even though history of spontaneous abortion would meet some definitions of a confounder. Directed acyclic graphs (DAGs) were introduced into epidemiology several years later as a tool with which to identify confounders. The authors now revisit Weinberg's paper using DAGs to represent scenarios that arise from her original assumptions. DAG theory is consistent with Weinberg's finding that adjusting for history of spontaneous abortion introduces bias in her original scenario. In the authors' examples, treating history of spontaneous abortion as a confounder introduces bias if it is a descendant of the exposure and is associated with the outcome conditional on exposure or is a child of a collider on a relevant undirected path. Thoughtful DAG analyses require clear research questions but are easily modified for examining different causal assumptions that may affect confounder assessment

Hormone Therapy and Young-Onset Breast Cancer

Estrogen plus progestin hormone therapy (HT) is associated with an increased risk of postmenopausal breast cancer, but few studies have examined the impact of HT use on the risk of breast cancer in younger women. We assessed the association between estrogen plus progestin HT or unopposed estrogen HT and young-onset breast cancer using data from the Two Sister Study (2008–2010), a sister-matched study of 1,419 cases diagnosed with breast cancer before the age of 50 years and 1,665 controls. We assessed exposures up to a family-specific index age to ensure comparable opportunities for exposures and used propensity scores to control for birth cohort effects on HT use. Ever HT use was uncommon (7% and 11% in cases and controls, respectively). Use of estrogen plus progestin was not associated with an increased risk of young-onset breast cancer (odds ratio = 0.80, 95% confidence interval: 0.41, 1.59). Unopposed estrogen use was inversely associated with the risk of young-onset breast cancer (odds ratio = 0.58, 95% confidence interval: 0.34, 0.99). Duration of use, age at first use, and recency of use did not modify these associations

Genetic Determinants of Facial Clefting: Analysis of 357 Candidate Genes Using Two National Cleft Studies from Scandinavia

Facial clefts are common birth defects with a strong genetic component. To identify fetal genetic risk factors for clefting, 1536 SNPs in 357 candidate genes were genotyped in two population-based samples from Scandinavia (Norway: 562 case-parent and 592 control-parent triads; Denmark: 235 case-parent triads).We used two complementary statistical methods, TRIMM and HAPLIN, to look for associations across these two national samples. TRIMM tests for association in each gene by using multi-SNP genotypes from case-parent triads directly without the need to infer haplotypes. HAPLIN on the other hand estimates the full haplotype distribution over a set of SNPs and estimates relative risks associated with each haplotype. For isolated cleft lip with or without cleft palate (I-CL/P), TRIMM and HAPLIN both identified significant associations with IRF6 and ADH1C in both populations, but only HAPLIN found an association with FGF12. For isolated cleft palate (I-CP), TRIMM found associations with ALX3, MKX, and PDGFC in both populations, but only the association with PDGFC was identified by HAPLIN. In addition, HAPLIN identified an association with ETV5 that was not detected by TRIMM.Strong associations with seven genes were replicated in the Scandinavian samples and our approach effectively replicated the strongest previously known association in clefting--with IRF6. Based on two national cleft cohorts of similar ancestry, two robust statistical methods and a large panel of SNPs in the most promising cleft candidate genes to date, this study identified a previously unknown association with clefting for ADH1C and provides additional candidates and analytic approaches to advance the field

Do Genetic Variants Modify the Effect of Smoking on Risk of Preeclampsia in Pregnancy?

Under embargo until: 2022-11-28Objective Maternal smoking is associated with as much as a 50% reduced risk of preeclampsia, despite increasing risk of other poor pregnancy outcomes that often co-occur with preeclampsia, such as preterm birth and fetal growth restriction. Researchers have long sought to understand whether this perplexing association is biologically based, or a result of noncausal mechanisms. We examined whether smoking-response genes modify the smoking-preeclampsia association to investigate potential biological explanations. Study Design We conducted a nested case–control study within the Norwegian Mother, Father and Child Birth Cohort (1999–2008) of 2,596 mother–child dyads. We used family-based log-linear Poisson regression to examine modification of the maternal smoking-preeclampsia relationship by maternal and fetal single nucleotide polymorphisms involved in cellular processes related to components of cigarette smoke (n = 1,915 with minor allele frequency ≥10%). We further investigated the influence of smoking cessation during pregnancy. Results Three polymorphisms showed overall (p < 0.001) multiplicative interaction between smoking and maternal genotype. For rs3765692 (TP73) and rs10770343 (PIK3C2G), protection associated with smoking was reduced with two maternal copies of the risk allele and was stronger in continuers than quitters (interaction p = 0.02 for both loci, based on testing 3-level smoking by 3-level genotype). For rs2278361 (APAF1) the inverse smoking-preeclampsia association was eliminated by the presence of a single risk allele, and again the trend was stronger in continuers than in quitters (interaction p = 0.01). Conclusion Evidence for gene–smoking interaction was limited, but differences by smoking cessation warrant further investigation. We demonstrate the potential utility of expanded dyad methods and gene–environment interaction analyses for outcomes with complex relationships between maternal and fetal genotypes and exposures.acceptedVersio

A family-based study of gene variants and maternal folate and choline in neuroblastoma: a report from the Children’s Oncology Group

Neuroblastoma is a childhood cancer of the sympathetic nervous system with embryonic origins. Previous epidemiologic studies suggest maternal vitamin supplementation during pregnancy reduces the risk of neuroblastoma. We hypothesized offspring and maternal genetic variants in folate-related and choline-related genes are associated with neuroblastoma and modify the effects of maternal intake of folate, choline and folic acid

Anti-dense fine speckled 70 (DFS70) autoantibodies: correlates and increasing prevalence in the United States

ObjectiveRecent studies report high-titer anti-dense fine speckled 70 (DFS70) autoantibodies in persons with inflammatory conditions, but the clinical significance remains unclear. Our goals were to estimate anti-DFS70 autoantibody prevalence, identify correlates, and assess time trends.MethodsSerum antinuclear antibodies (ANA) were measured by indirect immunofluorescence assay on HEp-2 cells in 13,519 participants ≥12 years old from three time periods (1988–1991, 1999–2004, 2011–2012) of the National Health and Nutrition Examination Survey. ANA-positive participants with dense fine speckled staining were evaluated for anti-DFS70 antibodies by enzyme-linked immunosorbent assay. We used logistic models adjusted for survey-design variables to estimate period-specific anti-DFS70 antibody prevalence in the US, and we further adjusted for sex, age, and race/ethnicity to identify correlates and assess time trends.ResultsWomen were more likely than men (odds ratio (OR)=2.97), black persons were less likely than white persons (OR=0.60), and active smokers were less likely than nonsmokers (OR=0.28) to have anti-DFS70 antibodies. The prevalence of anti-DFS70 antibodies increased from 1.6% in 1988-1991 to 2.5% in 1999-2004 to 4.0% in 2011-2012, which corresponds to 3.2 million, 5.8 million, and 10.4 million seropositive individuals, respectively. This increasing time trend in the US population (P&lt;0.0001) was modified in some subgroups and was not explained by concurrent changes in tobacco smoke exposure. Some, but not all, anti-DFS70 antibody correlates and time trends resembled those reported for total ANA.ConclusionMore research is needed to elucidate anti-DFS70 antibody triggers, their pathologic or potentially protective influences on disease, and their possible clinical implications