Family-based Designs Research Articles

Quantifying Bias from Measurable and Unmeasurable Confounders Across Three Domains of Individual Determinants of Political Preferences

AbstractA core part of political research is to identify how political preferences are shaped. The nature of these questions is such that robust causal identification is often difficult to achieve, and we are not seldom stuck with observational methods that we know have limited causal validity. The purpose of this paper is to measure the magnitude of bias stemming from both measurable and unmeasurable confounders across three broad domains of individual determinants of political preferences: socio-economic factors, moral values, and psychological constructs. We leverage a unique combination of rich Swedish registry data for a large sample of identical twins, with a comprehensive battery of 34 political preference measures, and build a meta-analytical model comparing our most conservative observational (naive) estimates with discordant twin estimates. This allows us to infer the amount of bias from unobserved genetic and shared environmental factors that remains in the naive models for our predictors, while avoiding precision issues common in family-based designs. The results are sobering: in most cases, substantial bias remains in naive models. A rough heuristic is that about half of the effect size even in conservative observational estimates is composed of confounding.

Associations of ARHGAP26 Polymorphisms with Alzheimer's Disease and Cardiovascular Disease.

The Rho GTPase activating protein 26 (ARHGAP26) gene has been reported to be associated with neuropsychiatric diseases and neurodegenerative diseases including Parkinson's disease. We examined whether the ARHGAP26 gene is associated with Alzheimer's disease (AD) and/or cardiovascular disease (CVD). Multivariable logistic regression model was used to examine the associations of 154 single nucleotide polymorphisms (SNPs) within the ARHGAP26 gene with AD and CVD using the Alzheimer's Disease Neuroimaging Initiative 1 (ADNI-1) cohort. Fourteen SNPs were associated with AD (top SNP rs3776362 with p = 3.43 × 10-3), while 37 SNPs revealed associations with CVD (top SNP rs415235 with p = 2.06 × 10-4). Interestingly, 13 SNPs were associated with both AD and CVD. SNP rs3776362 was associated with CVD, Functional Activities Questionnaire (FAQ), and Clinical Dementia Rating Sum of Boxes (CDR-SB). A replication study using a Caribbean Hispanics sample showed that 17 SNPs revealed associations with AD, and 12 SNPs were associated with CVD. The third sample using a family-based study design showed that 9 SNPs were associated with AD, and 3 SNPs were associated with CVD. SNP rs6836509 within the ARHGAP10 gene (an important paralogon of ARHGAP26) was associated with AD and cerebrospinal fluid total tau (t-tau) level in the ADNI sample. Several SNPs were functionally important using the RegulomeDB, while a number of SNPs were associated with significant expression quantitative trait loci (eQTLs) using Genotype-Tissue Expression (GTEx) databases. In conclusion, genetic variants within ARHGAP26 were associated with AD and CVD. These findings add important new insights into the potentially shared pathogenesis of AD and CVD.

Indirect Genetic Effects of ADIPOQ Variants on Lipid Levels in a Sibling Study of a Rural Chinese Population

Variations in lipid levels are the result of combinations of genetic and environmental factors. We aim to investigate the indirect effect between siblings of the three polymorphisms of ADIPOQ on serum lipid levels in rural Chinese populations. A total of 2571 sibling pairs were enrolled as study participants. A generalized estimating equation was used to accommodate a family-based design. We used stratified analysis to detect sex combination differences in the indirect genetic effect. We found a significant association between the number of altered risk alleles of rs182052 and ego lipid levels of TG (β = 0.177, P = 0.003), TC (β = 0.140, P = 0.004) and LDL-C (β = 0.098, P = 0.014). Ego and altered genotypes of rs182052 demonstrated a joint effect on ego lipid levels of TC (β = 0.212, P = 0.019), HDL-C (β = 0.099, P = 0.002) and LDL-C (β = 0.177, P = 0.013) in recessive inheritance mode. In opposite-sex siblings, the altered GG genotype of rs182052 increased the ego lipid levels. Thus, our findings demonstrate that ADIPOQ has an indirect genetic effect on lipid levels in sibling pairs, and there are sex-combination differences in the indirect genetic effect in siblings.

Candidate gene family-based and case-control studies of susceptibility to high Schistosoma mansoni worm burden in African children: a protocol.

Background: Approximately 25% of the risk of Schistosoma mansoni is associated with host genetic variation. We will test 24 candidate genes, mainly in the T h2 and T h17 pathways, for association with S. mansoni infection intensity in four African countries, using family based and case-control approaches. Methods: Children aged 5-15 years will be recruited in S. mansoni endemic areas of Ivory Coast, Cameroon, Uganda and the Democratic Republic of Congo (DRC). We will use family based (study 1) and case-control (study 2) designs. Study 1 will take place in Ivory Coast, Cameroon, Uganda and the DRC. We aim to recruit 100 high worm burden families from each country except Uganda, where a previous study recruited at least 40 families. For phenotyping, cases will be defined as the 20% of children in each community with heaviest worm burdens as measured by the circulating cathodic antigen (CCA) assay. Study 2 will take place in Uganda. We will recruit 500 children in a highly endemic community. For phenotyping, cases will be defined as the 20% of children with heaviest worm burdens as measured by the CAA assay, while controls will be the 20% of infected children with the lightest worm burdens. Deoxyribonucleic acid (DNA) will be genotyped on the Illumina H3Africa SNP (single nucleotide polymorphisms) chip and genotypes will be converted to sets of haplotypes that span the gene region for analysis. We have selected 24 genes for genotyping that are mainly in the Th2 and Th17 pathways and that have variants that have been demonstrated to be or could be associated with Schistosoma infection intensity. Analysis: In the family-based design, we will identify SNP haplotypes disproportionately transmitted to children with high worm burden. Case-control analysis will detect overrepresentation of haplotypes in extreme phenotypes with correction for relatedness by using whole genome principal components.

Candidate gene family-based and case-control studies of susceptibility to high Schistosoma mansoni worm burden in African children: a protocol

Background: Approximately 25% of the risk of Schistosoma mansoni is associated with host genetic variation. We will test 24 candidate genes, mainly in the T h2 and T h17 pathways, for association with S. mansoni infection intensity in four African countries, using family based and case-control approaches. Methods: Children aged 5-15 years will be recruited in S. mansoni endemic areas of Ivory Coast, Cameroon, Uganda and the Democratic Republic of Congo (DRC). We will use family based (study 1) and case-control (study 2) designs. Study 1 will take place in Ivory Coast, Cameroon, Uganda and the DRC. We aim to recruit 100 high worm burden families from each country except Uganda, where a previous study recruited at least 40 families. For phenotyping, cases will be defined as the 20% of children in each community with heaviest worm burdens as measured by the circulating cathodic antigen (CCA) assay. Study 2 will take place in Uganda. We will recruit 500 children in a highly endemic community. For phenotyping, cases will be defined as the 20% of children with heaviest worm burdens as measured by the CAA assay, while controls will be the 20% of infected children with the lightest worm burdens. Deoxyribonucleic acid (DNA) will be genotyped on the Illumina H3Africa SNP (single nucleotide polymorphisms) chip and genotypes will be converted to sets of haplotypes that span the gene region for analysis. We have selected 24 genes for genotyping that are mainly in the Th2 and Th17 pathways and that have variants that have been demonstrated to be or could be associated with Schistosoma infection intensity. Analysis: In the family-based design, we will identify SNP haplotypes disproportionately transmitted to children with high worm burden. Case-control analysis will detect overrepresentation of haplotypes in extreme phenotypes with correction for relatedness by using whole genome principal components.

Investigation of discordant sibling pairs from hereditary breast cancer families and analysis of a rare PMS1 variant

BackgroundIt is likely that additional genes for hereditary breast cancer can be identified using a discordant sib pair design. Using this design we identified individuals harboring a rare PMS1 c.605G>A variant previously predicted to result in loss of function. ObjectivesA family-based design and predictive algorithms were used to prioritize candidate variants possibly associated with an increased risk of hereditary breast cancer. Functional analyses were performed for one of the candidate variants, PMS1 c.605G>A. Methods1) 14 discordant sister-pairs from hereditary breast cancer families were identified. 2) Whole exome sequencing was performed and candidate risk variants identified. 3) A rare PMS variant was identified in 2 unrelated affected sisters but no unaffected siblings. 4) Functional analysis of this variant was carried out using targeted mRNA sequencing. ResultsGenotype-phenotype correlation did not demonstrate tracking of the variant with cancer in the family. Functional analysis revealed no difference in exon 6 incorporation, which was validated by analyzing PMS1 allele specific expression. ConclusionsThe PMS1 c.605G>A variant did not segregate with disease, and there was no variant-dependent impact on PMS1 exon 6 splicing, supporting this variant is likely benign. Functional analyses are imperative to understanding the clinical significance of predictive algorithms.

Quantitative Trait of Neutrophil Count Is Influenced By Variants in the Region of Gsdma and PSMD3-CSF3, aging, Cardiometabolic Comorbidities but Not By Chip

IntroductionHematopoiesis ensures lifelong oxygenation, immune and hemostatic functions through the production of billions of blood cells on a daily basis. However, little is known about the factors regulating and influencing this highly coordinated process and the resulting peripheral blood cell quantitative traits during the aging process. In this study, we investigated germline and acquired factors associated with blood cell counts in a cohort of normal aging individuals.MethodsDeterminants underlying blood cell trait variability were assessed in a cohort of 2996 related and unrelated women of French-Canadian ancestry aged 55 to 101 years. All participants answered a medical questionnaire and provided a blood sample for complete blood count with differential (GenS, Beckman Coulter) and DNA analysis. Potential hereditary variants of significance were assessed using a genome-wide association study (GWAS). Genotyping was performed according to the manufacturer's specifications on the Illumina Infinium Global Array v3-MD (Illumina, San Diego, CA). Variants with a completion rate of ≥98%, minor allele frequency >1% and imputation probability of ≥0.80 were retained. Genome-wide association testing was conducted with SAIGE_0.43.3 to account for the family-based design of our study. Genome-wide significance threshold was at 5x10 -8. Potential acquired factors of importance such as chronic comorbidities and clonal hematopoiesis of indeterminate potential (CHIP) were also sought. CHIP status was assessed with targeted next-generation sequencing of polymorphonuclear cells using the Ampliseq AML panel. Germline and acquired factors were then integrated in generalized linear mixed models to identify factors associated with blood cell indices in multivariate analysis and to characterize their relative contribution. Statistical analyses were conducted with R.ResultsMean age of study participants was 69.2 years (standard deviation 9.1 years). Mean values (and standard deviation) of blood cell indices were as follows: total white blood count (WBC) 6.43x10 9/L (1.63), absolute lymphocyte count 1.86x10 9/L (0.57), absolute monocyte count 0.46x10 9/L (0.16), absolute neutrophil count 3.93x10 9/L (1.29), hemoglobin 131.2 g/L (9.67), platelet count 249.2x10 9/L (57.6). Neutrophil and monocyte counts increased with age (b1 coefficient 0.02 (p-value 6.06E-08) and 0.004 (p-value 2.2E-16) respectively), while lymphocyte and platelet counts decreased with advancing age (b1 coefficient -0.01 (p-value 3.65E-07) and -0.29 (p-value 0.023) respectively). GWAS identified 13 variants in the region of GSDMA and PSMD3-CSF3 (chromosome 17) that met genome-wide requirements for WBC and neutrophil counts. Platelet count was significantly associated with a variant intronic to ARHGEF3 (chromosome 3). Among acquired factors, smoking was positively associated with neutrophil, monocyte, lymphocyte counts and hemoglobin levels. Distinctly, we document that cardiometabolic comorbidities (diabetes, coronary heart disease, hypertension and dyslipidemia) are associated with a statistically higher count of myeloid-derived cells (neutrophil count 4.1 vs 3.83 (95%CI 0.28-0.37, p-value <0.001), monocyte count 0.50 vs 0.45 (95%CI 0.02-0.05, p-value <0.001), and platelet count 259 vs 243 (95%CI 10-20, p-value <0.001). These results remained significant in multivariate analysis. In accordance with previous reports, CHIP, which was documented in 14% of the cohort, had no influence on blood counts.ConclusionWe document that individual variation in blood counts in individuals is influenced by several factors: (i) germline variants related to GSDMA and PSMD3-CSF3 contribute to white blood cell and neutrophil counts and a variant intronic to ARHGEF3 is influential in determining platelet counts; (ii) aging is associated with increasing levels of neutrophils and monocytes, and reduced lymphocyte and platelet counts, indicating a shift towards myelopoiesis; (iii) this myeloid-biased skewing is further increased among individuals with cardiometabolic comorbidities; (iv) CHIP does not contribute to the age-associated myeloid shift. These findings support that chronic age-related diseases may promote myelopoiesis and contribute to population variability in peripheral blood traits, possibly through a state of low-grade inflammation. [Display omitted] DisclosuresBusque: Novartis: Consultancy.

Data Resource Profile: The Manitoba Multigenerational Cohort.

Families are important components of society; given their shared genetic and social environments, studying families can provide critical insight into health and social outcomes within family members and across generations. Few places in the world have high-quality linkages within families at a population level; one such location is Manitoba, Canada. The Manitoba Multigenerational Cohort (MMC) has been developed to facilitate health and social research using family-based designs. The MMC is derived from the Manitoba Health Insurance Registry (‘the Registry’)—a population-based registry of all individuals registered with Manitoba Health and Seniors Care (MHSC).1 The Registry is updated at the Manitoba Centre for Health Policy (MCHP) twice a year and is integrated with historical registry data to create a longitudinal population-based registry. Over 80 administrative and survey-based data sets are linkable at the individual level; these data sets make up the Manitoba Population Research Data Repository (‘ the Repository’) and are housed in a secure environment at the MCHP.

Candidate gene family-based and case-control studies of susceptibility to high Schistosoma mansoni worm burden in African children: a protocol

Background: Approximately 25% of the risk of Schistosoma mansoni is associated with host genetic variation. We will test 24 candidate genes, mainly in the Th2 and Th17 pathways, for association with S. mansoni infection intensity in four African countries, using family based and case-control approaches. Methods: Children aged 5-15 years will be recruited in S. mansoni endemic areas of Ivory Coast, Cameroon, Uganda and the Democratic Republic of Congo (DRC). We will use family based (study 1) and case-control (study 2) designs. Study 1 will take place in Ivory Coast, Cameroon, Uganda and the DRC. We aim to recruit 100 high worm burden families from each country except Uganda, where a previous study recruited at least 40 families. For phenotyping, cases will be defined as the 20% of children in each community with heaviest worm burdens as measured by the circulating cathodic antigen (CCA) assay. Study 2 will take place in Uganda. We will recruit 500 children in a highly endemic community. For phenotyping, cases will be defined as the 20% of children with heaviest worm burdens as measured by the CAA assay, while controls will be the 20% of infected children with the lightest worm burdens. Deoxyribonucleic acid (DNA) will be genotyped on the Illumina H3Africa SNP (single nucleotide polymorphisms) chip and genotypes will be converted to sets of haplotypes that span the gene region for analysis. We have selected 24 genes for genotyping that are mainly in the Th2 and Th17 pathways and that have variants that have been demonstrated to be or could be associated with Schistosoma infection intensity. Analysis: In the family-based design, we will identify SNP haplotypes disproportionately transmitted to children with high worm burden. Case-control analysis will detect overrepresentation of haplotypes in extreme phenotypes with correction for relatedness by using whole genome principal components.

What Have We Learned From Family-Based Studies About Spondyloarthritis?

Spondyloarthritis (SpA) is a chronic inflammatory disorder with a high familial aggregation, emphasizing the existence of genetic susceptibility factors. In the last decades, family-based studies have contributed to better understand the genetic background of SpA, in particular by showing that the most likely model of transmission is oligogenic with multiplicative effects. Coexistence of different SpA subtypes within families also highlighted the complex interplay between all subtypes. Several whole-genome linkage analyses using sib-pairs or multiplex families were performed in the 1990s to try to identify genetic susceptibility factors besides HLA-B27. Unfortunately, no consistent results were obtained and family-based studies have been progressively set aside in favor of case-control designs. In particular, case-control genome-wide association studies allowed the identification of more than 40 susceptibility regions. However, all these loci explain only a small fraction of disease predisposition. Several hypotheses have been advanced to account for this unexplained heritability, including rare variants involvement, leading to a renewed interest in family-based designs, which are probably more powerful in the detection of such variants. In this review, our purpose is to summarize what has been learned to date regarding SpA genetics from family-based studies, with a special focus on recent identification of rare associated variants through next-generation sequencing studies.

Characterizing sensitivity and coverage of clinical WGS as a diagnostic test for genetic disorders

BackgroundDue to its reduced cost and incomparable advantages, WGS is likely to lead to changes in clinical diagnosis of rare and undiagnosed diseases. However, the sensitivity and breadth of coverage of clinical WGS as a diagnostic test for genetic disorders has not been fully evaluated.MethodsHere, the performance of WGS in NA12878, the YH cell line, and the Chinese trios were measured by assessing their sensitivity, PPV, depth and breadth of coverage using MGISEQ-2000. We also compared the performance of WES and WGS using NA12878. The sensitivity and PPV were tested using the family-based trio design for the Chinese trios. We further developed a systematic WGS pipeline for the analysis of 8 clinical cases.ResultsIn general, the sensitivity and PPV for SNV/indel detection increased with mean depth and reached a plateau at an ~ 40X mean depth using down-sampling samples of NA12878. With a mean depth of 40X, the sensitivity of homozygous and heterozygous SNPs of NA12878 was > 99.25% and > 99.50%, respectively, and the PPV was 99.97% and 98.96%. Homozygous and heterozygous indels showed lower sensitivity and PPV. The sensitivity and PPV were still not 100% even with a mean depth of ~ 150X. We also observed a substantial variation in the sensitivity of CNV detection across different tools, especially in CNVs with a size less than 1 kb. In general, the breadth of coverage for disease-associated genes and CNVs increased with mean depth. The sensitivity and coverage of WGS (~ 40X) was better than WES (~ 120X). Among the Chinese trios with an ~ 40X mean depth, the sensitivity among offspring was > 99.48% and > 96.36% for SNP and indel detection, and the PPVs were 99.86% and 97.93%. All 12 previously validated variants in the 8 clinical cases were successfully detected using our WGS pipeline.ConclusionsThe current standard of a mean depth of 40X may be sufficient for SNV/indel detection and identification of most CNVs. It would be advisable for clinical scientists to determine the range of sensitivity and PPV for different classes of variants for a particular WGS pipeline, which would be useful when interpreting and delivering clinical reports.

Plasma metabolomics of autism spectrum disorder and influence of shared components in proband families.

Prevalence of autism spectrum disorder (ASD) has been increasing in the United States in the past decades. The exact mechanisms remain enigmatic, and diagnosis of the disease still relies primarily on assessment of behavior. We first used a case–control design (75 idiopathic cases and 29 controls, enrolled at Boston Children’s Hospital from 2007-2012) to identify plasma biomarkers of ASD through a metabolome-wide association study approach. Then we leveraged a family-based design (31 families) to investigate the influence of shared genetic and environmental components on the autism-associated features. Using untargeted high-resolution mass spectrometry metabolomics platforms, we detected 19 184 features. Of these, 191 were associated with ASD (false discovery rate < 0.05). We putatively annotated 30 features that had an odds ratio (OR) between <0.01 and 5.84. An identified endogenous metabolite, O-phosphotyrosine, was associated with an extremely low autism odds (OR 0.17; 95% confidence interval 0.06-0.39). We also found that glutathione metabolism was associated with ASD (P = 0.048). Correlations of the significant features between proband and parents were low (median = 0.09). Of the 30 annotated features, the median correlations within families (proband–parents) were −0.15 and 0.24 for the endogenous and exogenous metabolites, respectively. We hypothesize that, without feature identification, family-based correlation analysis of autism-associated features can be an alternative way to assist the prioritization of potentially diagnostic features. A panel of ASD diagnostic metabolic markers with high specificity could be derived upon further studies.

Linkage of Alzheimer disease families with Puerto Rican ancestry identifies a chromosome 9 locus

The genetic admixture of Caribbean Hispanics provides an opportunity to discover novel genetic factors in Alzheimer disease (AD). We sought to identify genetic variants for AD through a family-based design using the Puerto Rican (PR) Alzheimer Disease Initiative (PRADI). Whole-genome sequencing (WGS) and parametric linkage analysis were performed for 100 individuals from 23 multiplex PRADI families. Variants were prioritized by minor allele frequency (<0.01), functional potential [combined annotation dependent depletion score (CADD) >10], and co-segregation with AD. Variants were further ranked using an independent PR case-control WGS dataset (PR10/66). A genome-wide significant linkage peak was found in 9p21 with a heterogeneity logarithm of the odds score (HLOD) >5.1, which overlaps with an AD linkage region from two published independent studies. The region harbors C9orf72, but no expanded repeats were observed in the families. Seven variants prioritized by the PRADI families also displayed evidence for association in the PR10/66 (p < 0.05), including a missense variant in UNC13B. Our study demonstrated the importance of family-based design and WGS in genetic study of AD.

Investigating the impact of Down syndrome on methylation and glycomics with two-stage PO2PLS.

Down syndrome (DS) is a condition that leads to precocious and accelerated aging in affected subjects. Several alterations in DS cases have been reported at a molecular level, particularly in methylation and glycosylation. Investigating the relation between methylation, glycomics and DS can lead to new insights underlying the atypical aging. We consider a data integration approach, where we investigate how DS affects the parts of glycomics and methylation which are correlated, and which CpG sites and glycans are relevant. Our motivating datasets consist of methylation and glycomics data, measured on 29 DS patients and their unaffected siblings and mothers. The family-based case-control design needs to be taken into account when studying the relationship between methylation, glycomics and DS. We propose a two-stage approach to first integrate methylation and glycomics data, and then link the joint information to Down syndrome. For the data integration step, we consider probabilistic two-way orthogonal partial least squares (PO2PLS). PO2PLS models two omics datasets in terms of low-dimensional joint and omic-specific latent components, and takes into account heterogeneity across the omics data. The relationship between the omics data can be statistically tested. The joint components represent the joint information in methylation and glycomics. In the second stage, we apply a linear mixed model to the relationship between DS and the joint methylation and glycomics components. For the components that are significantly as sociated with DS, we identify the most important CpG sites and glycans. A simulation study is conducted to evaluate the performance of our approach. The results showed that the effects of DS on the omics data can be detected in a large sample size, and the accuracy of the feature selection was high in both small and large sample sizes. Our approach is applied to the DS datasets, a significant effect of DS on the joint components is found. The identified CpG sites and glycans appeared to be related to DS. Our proposed method that jointly analyzes multiple omics data with an outcome variable may provide new insight into the molecular implications of DS at different omics levels.

A unifying framework for rare variant association testing in family-based designs, including higher criticism approaches, SKATs, and burden tests.

Analysis of rare variants in family-based studies remains a challenge. Transmission-based approaches provide robustness against population stratification, but the evaluation of the significance of test statistics based on asymptotic theory can be imprecise. Also, power will depend heavily on the choice of the test statistic and on the underlying genetic architecture of the locus, which will be generally unknown. In our proposed framework, we utilize the FBAT haplotype algorithm to obtain the conditional offspring genotype distribution under the null hypothesis given the sufficient statistic. Based on this conditional offspring genotype distribution, the significance of virtually any association test statistic can be evaluated based on simulations or exact computations, without the need for asymptotic approximations. Besides standard linear burden-type statistics, this enables our approach to also evaluate other test statistics such as variance components statistics, higher criticism approaches, and maximum-single-variant-statistics, where asymptotic theory might be involved or does not provide accurate approximations for rare variant data. Based on these P-values, combined test statistics such as the aggregated Cauchy association test (ACAT) can also be utilized. In simulation studies, we show that our framework outperforms existing approaches for family-based studies in several scenarios. We also applied our methodology to a TOPMed whole-genome sequencing dataset with 897 asthmatic trios from Costa Rica. FBAT software is available at https://sites.google.com/view/fbatwebpage. Simulation code is available at https://github.com/julianhecker/FBAT_rare_variant_test_simulations. Whole-genome sequencing data for 'NHLBI TOPMed: The Genetic Epidemiology of Asthma in Costa Rica' is available at https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000988.v4.p1. Supplementary data are available at Bioinformatics online.

Case-Parent Trio Studies in Cleft Lip and Palate.

Cleft lip with or without cleft palate (CL/P) is one of the most common congenital malformations in humans involving various genetic and environmental risk factors. The prevalence of CL/P varies according to geographical location, ethnicity, race, gender, and socioeconomic status, affecting approximately 1 in 800 live births worldwide. Genetic studies aim to understand the mechanisms contributory to a phenotype by measuring the association between genetic variants and also between genetic variants and phenotype population. Genome-wide association studies are standard tools used to discover genetic loci related to a trait of interest. Genetic association studies are generally divided into two main design types: population-based studies and family-based studies. The epidemiological population-based studies comprise unrelated individuals that directly compare the frequency of genetic variants between (usually independent) cases and controls. The alternative to population-based studies (case–control designs) includes various family-based study designs that comprise related individuals. An example of such a study is a case–parent trio design study, which is commonly employed in genetics to identify the variants underlying complex human disease where transmission of alleles from parents to offspring is studied. This article describes the fundamentals of case–parent trio study, trio design and its significances, statistical methods, and limitations of the trio studies.

Innovative approach to identify multigenomic and environmental interactions associated with birth defects in family-based hybrid designs.

Genes, including those with transgenerational effects, work in concert with behavioral, environmental, and social factors via complex biological networks to determine human health. Understanding complex relationships between causal factors underlying human health is an essential step towards deciphering biological mechanisms. We propose a new analytical framework to investigate the interactions between maternal and offspring genetic variants or their surrogate single nucleotide polymorphisms (SNPs) and environmental factors using family-based hybrid study design. The proposed approach can analyze diverse genetic and environmental factors and accommodate samples from a variety of family units, including case/control-parental triads, and case/control-parental dyads, while minimizing potential bias introduced by population admixture. Comprehensive simulations demonstrated that our innovative approach outperformed the log-linear approach, the best available method for case-control family data. The proposed approach had greater statistical power and was capable to unbiasedly estimate the maternal and child genetic effects and the effects of environmental factors, while controlling the Type I error rate against population stratification. Using our newly developed approach, we analyzed the associations between maternal and fetal SNPs and obstructive and conotruncal heart defects, with adjustment for demographic and lifestyle factors and dietary supplements. Fourteen and 11 fetal SNPs were associated with obstructive and conotruncal heart defects, respectively. Twenty-seven and 17 maternal SNPs were associated with obstructive and conotruncal heart defects, respectively. In addition, maternal body mass index was a significant risk factor for obstructive defects. The proposed approach is a powerful tool for interrogating the etiological mechanism underlying complex traits.

Huntingtin gene CAG repeat size affects autism risk: Family-based and case-control association study.

The Huntingtin (HTT) gene contains a CAG repeat in exon 1, whose expansion beyond 39 repeats consistently leads to Huntington's disease (HD), whereas normal-to-intermediate alleles seemingly modulate brain structure, function and behavior. The role of the CAG repeat in Autism Spectrum Disorder (ASD) was investigated applying both family-based and case-control association designs, with the SCA3 repeat as a negative control. Significant overtransmission of "long" CAG alleles (≥17 repeats) to autistic children and of "short" alleles (≤16 repeats) to their unaffected siblings (all p < 10-5 ) was observed in 612 ASD families (548 simplex and 64 multiplex). Surprisingly, both 193 population controls and 1,188 neurological non-HD controls have significantly lower frequencies of "short" CAG alleles compared to 185 unaffected siblings and higher rates of "long" alleles compared to 548 ASD patients from the same families (p < .05-.001). The SCA3 CAG repeat displays no association. "Short" HTT alleles seemingly exert a protective effect from clinically overt autism in families carrying a genetic predisposition for ASD, while "long" alleles may enhance autism risk. Differential penetrance of autism-inducing genetic/epigenetic variants may imply atypical developmental trajectories linked to HTT functions, including excitation/inhibition imbalance, cortical neurogenesis and apoptosis, neuronal migration, synapse formation, connectivity and homeostasis.

Parental education and the risk of cerebral palsy for children: an evaluation of causality.

To explore whether increasing parental education has a causal effect on risk of cerebral palsy (CP) in the child, or whether unobserved confounding is a more likely explanation. We used data from Norwegian registries on approximately 1.5million children born between 1967 and 2011. We compared results from a traditional cohort design with results from a family-based matched case-control design, in which children with CP were matched to their first cousins without CP. In addition, we performed a simulation study to assess the role of unobserved confounding. In the cohort design, the odds of CP were reduced in children of mothers and fathers with higher education (adjusted odds ratio [OR] 0.67, 95% confidence interval [CI] 0.60-0.75 for maternal education, and adjusted OR 0.75, 95% CI 0.67-0.85 for paternal education). In the family-based case-control design, only an association for maternal education remained (adjusted OR 0.80, 95% CI 0.64-0.99). Results from a simulation study suggested that this association could be explained by unobserved confounding. A causal effect of obtaining higher education on risk of CP in the child is unlikely. Results stress the importance of continued research on the role of genetic and environmental risk factors that vary by parents' educational level. Children of higher-educated parents had significantly lower odds of cerebral palsy (CP). There was no evidence of difference in risk of CP within first cousins whose mothers or fathers had different educational levels. Association between parental education and odds of CP did not reflect a causal effect.

Family Study Designs Informed by Tumor Heterogeneity and Multi-Cancer Pleiotropies: The Power of the Utah Population Database.

Previously, family-based designs and high-risk pedigrees have illustrated value for the discovery of high- and intermediate-risk germline breast cancer susceptibility genes. However, genetic heterogeneity is a major obstacle hindering progress. New strategies and analytic approaches will be necessary to make further advances. One opportunity with the potential to address heterogeneity via improved characterization of disease is the growing availability of multisource databases. Specific to advances involving family-based designs are resources that include family structure, such as the Utah Population Database (UPDB). To illustrate the broad utility and potential power of multisource databases, we describe two different novel family-based approaches to reduce heterogeneity in the UPDB. Our first approach focuses on using pedigree-informed breast tumor phenotypes in gene mapping. Our second approach focuses on the identification of families with similar pleiotropies. We use a novel network-inspired clustering technique to explore multi-cancer signatures for high-risk breast cancer families. Our first approach identifies a genome-wide significant breast cancer locus at 2q13 [P = 1.6 × 10-8, logarithm of the odds (LOD) equivalent 6.64]. In the region, IL1A and IL1B are of particular interest, key cytokine genes involved in inflammation. Our second approach identifies five multi-cancer risk patterns. These clusters include expected coaggregations (such as breast cancer with prostate cancer, ovarian cancer, and melanoma), and also identify novel patterns, including coaggregation with uterine, thyroid, and bladder cancers. Our results suggest pedigree-informed tumor phenotypes can map genes for breast cancer, and that various different cancer pleiotropies exist for high-risk breast cancer pedigrees. Both methods illustrate the potential for decreasing etiologic heterogeneity that large, population-based multisource databases can provide.See all articles in this CEBP Focus section, "Modernizing Population Science."