Statistical analysis of the associations of hereditary factors with the risk of Type 1 diabetes and Diabetic Nephropathy based on the familial data collected through ascertaiment from population-based registers

In genetic epidemiology, population-based disease registries are commonly used to collect genotype or other risk factor information concerning affected subjects and their relatives. This work presents two new approaches for the statistical inference of ascertained data: a conditional and full likelihood approaches for the disease with variable age at onset phenotype using familial data obtained from population-based registry of incident cases. The aim is to obtain statistically reliable estimates of the general population parameters. The statistical analysis of familial data with variable age at onset becomes more complicated when some of the study subjects are non-susceptible, that is to say these subjects never get the disease. A statistical model for a variable age at onset with long-term survivors is proposed for studies of familial aggregation, using latent variable approach, as well as for prospective studies of genetic association studies with candidate genes. In addition, we explore the possibility of a genetic explanation of the observed increase in the incidence of Type 1 diabetes (T1D) in Finland in recent decades and the hypothesis of non-Mendelian transmission of T1D associated genes. Both classical and Bayesian statistical inference were used in the modelling and estimation. Despite the fact that this work contains five studies with different statistical models, they all concern data obtained from nationwide registries of T1D and genetics of T1D. In the analyses of T1D data, non-Mendelian transmission of T1D susceptibility alleles was not observed. In addition, non-Mendelian transmission of T1D susceptibility genes did not make a plausible explanation for the increase in T1D incidence in Finland. Instead, the Human Leucocyte Antigen associations with T1D were confirmed in the population-based analysis, which combines T1D registry information, reference sample of healthy subjects and birth cohort information of the Finnish population. Finally, a substantial familial variation in the susceptibility of T1D nephropathy was observed. The presented studies show the benefits of sophisticated statistical modelling to explore risk factors for complex diseases.

Evaluation of polymorphisms in predicted target sites for micro RNAs differentially expressed in endometriosis

Previous microarray analyses identified 22 microRNAs (miRNAs) differentially expressed in paired ectopic and eutopic endometrium of women with and without endometriosis. To investigate further the role of these miRNAs in women with endometriosis, we conducted an association study aiming to explore the relationship between endometriosis risk and single-nucleotide polymorphisms (SNPs) in miRNA target sites for these differentially expressed miRNAs. A panel of 102 SNPs in the predicted miRNA binding sites were evaluated for an endometriosis association study and an ingenuity pathway analysis was performed. Fourteen rare variants were identified in this study. We found SNP rs14647 in the Wolf-Hirschhorn syndrome candidate gene1 (WHSC1) 3'UTR (untranslated region) was associated with endometriosis-related infertility presenting an odds ratio of 12.2 (95% confidence interval = 2.4-60.7, P = 9.03 x 10(-5)). SNP haplotype AGG in the solute carrier family 22, member 23 (SLC22A23) 3'UTR was associated with endometriosis-related infertility and more severe disease. With the individual genotyping data, ingenuity pathways analysis identified the tumour necrosis factor and cyclin-dependant kinase inhibitor as major factors in the molecular pathways. Significant associations between WHSC1 alleles and endometriosis-related infertility and SLC22A23 haplotypes and the disease severe stage were identified. These findings may help focus future research on subphenotypes of this disease. Replication studies in independent large sample sets to confirm and characterize the involvement of the gene variation in the pathogenesis of endometriosis are needed.

Variation in BMPR1B, TGFRB1 and BMPR2 and control of dizygotic twinning

Genes in the TGF9 signaling pathway play important roles in the regulation of ovarian follicle growth and ovulation rate. Mutations in three genes in this pathway, growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15) and the bone morphogenetic protein receptor B 1 (BMPRB1), influence dizygotic (DZ) twinning rates in sheep. To date, only variants in GDF9 and BMP15, but not their receptors transforming growth factor ss receptor 1 (TGFBR1), bone morphogenetic protein receptor 2 (BMPR2) and BMPR1B, have been investigated with respect to their roles in human DZ twinning. We screened for rare and novel variants in TGFBR1, BMPR2 and BMPR1B in mothers of dizygotic twins (MODZT) from twin-dense families, and assessed association between genotyped and imputed variants and DZ twinning in another large sample of MODZT. Three novel variants were found: a deep intronic variant in BMPR2, and one intronic and one non-synonymous exonic variant in BMPRB1 which would result in the replacement of glutamine by glutamic acid at amino acid position 294 (p.Gln294Glu). None of these variants were predicted to have major impacts on gene function. However, the p.Gln294Glu variant changes the same amino acid as a sheep BMPR1B functional variant and may have functional consequences. Six BMPR1B variants were marginally associated with DZ twinning in the larger case-control sample, but these were no longer significant once multiple testing was taken into account. Our results suggest that variation in the TGF9 signaling pathway type II receptors has limited effects on DZ twinning rates in humans.

A rare P2X7 variant Arg307Gln with absent pore formation function protects against neuroinflammation in multiple sclerosis

Multiple sclerosis (MS) is a chronic relapsing-remitting inflammatory disease of the central nervous system characterized by oligodendrocyte damage, demyelination and neuronal death. Genetic association studies have shown a 2-fold or greater prevalence of the HLA-DRB1*1501 allele in the MS population compared with normal Caucasians. In discovery cohorts of Australasian patients with MS (total 2941 patients and 3008 controls), we examined the associations of 12 functional polymorphisms of P2X7, a microglial/macrophage receptor with proinflammatory effects when activated by extracellular adenosine triphosphate (ATP). In discovery cohorts, rs28360457, coding for Arg307Gln was associated with MS and combined analysis showed a 2-fold lower minor allele frequency compared with controls (1.11% for MS and 2.15% for controls, P = 0.0000071). Replication analysis of four independent European MS case–control cohorts (total 2140 cases and 2634 controls) confirmed this association [odds ratio (OR) = 0.69, P = 0.026]. A meta-analysis of all Australasian and European cohorts indicated that Arg307Gln confers a 1.8-fold protective effect on MS risk (OR = 0.57, P = 0.0000024). Fresh human monocytes heterozygous for Arg307Gln have >85% loss of ‘pore’ function of the P2X7 receptor measured by ATP-induced ethidium uptake. Analysis shows Arg307Gln always occurred with 270His suggesting a single 307Gln–270His haplotype that confers dominant negative effects on P2X7 function and protection against MS. Modeling based on the homologous zP2X4 receptor showed Arg307 is located in a region rich in basic residues located only 12 Å from the ligand binding site. Our data show the protective effect against MS of a rare genetic variant of P2RX7 with heterozygotes showing near absent proinflammatory ‘pore’ function.

Heritability of the network architecture of intrinsic brain functional connectivity

The brain's functional network exhibits many features facilitating functional specialization, integration, and robustness to attack. Using graph theory to characterize brain networks, studies demonstrate their small-world, modular, and "rich-club" properties, with deviations reported in many common neuropathological conditions. Here we estimate the heritability of five widely used graph theoretical metrics (mean clustering coefficient (γ), modularity (Q), rich-club coefficient (ϕnorm), global efficiency (λ), small-worldness (σ)) over a range of connection densities (k=5-25%) in a large cohort of twins (N=592, 84 MZ and 89 DZ twin pairs, 246 single twins, age 23±2.5). We also considered the effects of global signal regression (GSR). We found that the graph metrics were moderately influenced by genetic factors h2 (γ=47-59%, Q=38-59%, ϕnorm=0-29%, λ=52-64%, σ=51-59%) at lower connection densities (≤15%), and when global signal regression was implemented, heritability estimates decreased substantially h2 (γ=0-26%, Q=0-28%, ϕnorm=0%, λ=23-30%, σ=0-27%). Distinct network features were phenotypically correlated (|r|=0.15-0.81), and γ, Q, and λ were found to be influenced by overlapping genetic factors. Our findings suggest that these metrics may be potential endophenotypes for psychiatric disease and suitable for genetic association studies, but that genetic effects must be interpreted with respect to methodological choices.

Novel associations of UDP-glucuronosyltransferase 2B gene variants with prostate cancer risk in a multiethnic study.

BACKGROUND: We have previously shown that a functional polymorphism of the UGT2B15 gene (rs1902023) was associated with increased risk of prostate cancer (PC). Novel functional polymorphisms of the UGT2B17 and UGT2B15 genes have been recently characterized by in vitro assays but have not been evaluated in epidemiologic studies. METHODS: Fifteen functional SNPs of the UGT2B17 and UGT2B15 genes, including cis-acting UGT2B gene SNPs, were genotyped in African American and Caucasian men (233 PC cases and 342 controls). Regression models were used to analyze the association between SNPs and PC risk. RESULTS: After adjusting for race, age and BMI, we found that six UGT2B15 SNPs (rs4148269, rs3100, rs9994887, rs13112099, rs7686914 and rs7696472) were associated with an increased risk of PC in log-additive models (p < 0.05). A SNP cis-acting on UGT2B17 and UGT2B15 expression (rs17147338) was also associated with increased risk of prostate cancer (OR = 1.65, 95% CI = 1.00-2.70); while a stronger association among men with high Gleason sum was observed for SNPs rs4148269 and rs3100. CONCLUSIONS: Although small sample size limits inference, we report novel associations between UGT2B15 and UGT2B17 variants and PC risk. These associations with PC risk in men with high Gleason sum, more frequently found in African American men, support the relevance of genetic differences in the androgen metabolism pathway, which could explain, in part, the high incidence of PC among African American men. Larger studies are required.

Polymorphisms in the kinesin-like factor 1 B gene and risk of epithelial ovarian cancer in Eastern Chinese women.

The kinesin-like factor 1 B (KIF1B) gene plays an important role in the process of apoptosis and the transformation and progression of malignant cells. Genetic variations in KIF1B may contribute to risk of epithelial ovarian cancer (EOC). In this study of 1,324 EOC patients and 1,386 cancer-free female controls, we investigated associations between two potentially functional single nucleotide polymorphisms in KIF1B and EOC risk by the conditional logistic regression analysis. General linear regression model was used to evaluate the correlation between the number of variant alleles and KIF1B mRNA expression levels. We found that the rs17401966 variant AG/GG genotypes were significantly associated with a decreased risk of EOC (adjusted odds ratio (OR) = 0.81, 95 % confidence interval (CI) = 0.68-0.97), compared with the AA genotype, but no associations were observed for rs1002076. Women who carried both rs17401966 AG/GG and rs1002076 AG/AA genotypes of KIF1B had a 0.82-fold decreased risk (adjusted 95 % CI = 0.69-0.97), compared with others. Additionally, there was no evidence of possible interactions between about-mentioned co-variants. Further genotype-phenotype correlation analysis indicated that the number of rs17401966 variant G allele was significantly associated with KIF1B mRNA expression levels (P for GLM = 0.003 and 0.001 in all and Chinese subjects, respectively), with GG carriers having the lowest level of KIF1B mRNA expression. Taken together, the rs17401966 polymorphism likely regulates KIF1B mRNA expression and thus may be associated with EOC risk in Eastern Chinese women. Larger, independent studies are warranted to validate our findings.

Apoptotic variants as predictors of risk of oropharyngeal cancer recurrence after definitive radiotherapy.

Single nucleotide polymorphisms (SNPs) in the promoter region of FAS and FASLG may alter their transcriptional activity. Thus, we determined the associations between four FAS and FASLG promoter variants (FAS1377G>A, rs2234767; 670A>G, rs1800682; FASLG844T>C, rs763110 and 124A>G, rs5030772) and the risk of recurrence of squamous cell carcinoma of the oropharynx (SCCOP). We evaluated the associations between FAS and FASLG genetic variants and the risk of recurrence in a cohort of 1,008 patients. The log-rank test and multivariate Cox models were used to evaluate the associations. Compared with patients with common homozygous genotypes of FAS670 and FASLG844 polymorphisms, patients with variant genotypes had lower disease-free survival rates (log-rank p < 0.0001 and p < 0.0001, respectively) and an approximately threefold higher risk of SCCOP recurrence (HR, 3.2;95% CI, 2.2-4.6; and HR, 3.1; 95% CI, 2.2-4.4, respectively) after multivariate adjustment. Furthermore, among patients with HPV16-positive tumors, those with variant genotypes of these two polymorphisms had lower disease-free survival rates (log-rank, p < 0.0001 and p < 0.0001, respectively) and a higher recurrence risk than did patients with common homozygous genotypes (HR, 12.9; 95% CI, 3.8-43.6; and HR, 8.1; 95% CI, 3.6-18.6, respectively), whereas no significant associations were found for FAS1377 and FASLG124 polymorphisms. Our findings suggest that FAS670 and FASLG844 polymorphisms modulate the risk of recurrence of SCCOP, particularly in patients with HPV16-positive tumors. Larger studies are needed to validate these results.

SNP discovery using Next Generation Transcriptomic Sequencing in Atlantic herring (Clupea harengus)

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of non-model species with unprecedented genomic coverage, allowing evolutionary biologists to address questions previously far beyond the reach of available resources. Furthermore, the simple mutation model of Single Nucleotide Polymorphisms (SNPs) permits cost-effective high-throughput genotyping in thousands of individuals simultaneously. Genomic resources are scarce for the Atlantic herring (Clupea harengus), a small pelagic species that sustains high revenue fisheries. This paper details the development of 578 SNPs using a combined NGS and high-throughput genotyping approach. Eight individuals covering the species distribution in the eastern Atlantic were bar-coded and multiplexed into a single cDNA library and sequenced using the 454 GS FLX platform. SNP discovery was performed by de novo sequence clustering and contig assembly, followed by the mapping of reads against consensus contig sequences. Selection of candidate SNPs for genotyping was conducted using an in silico approach. SNP validation and genotyping were performed simultaneously using an Illumina 1,536 GoldenGate assay. Although the conversion rate of candidate SNPs in the genotyping assay cannot be predicted in advance, this approach has the potential to maximise cost and time efficiencies by avoiding expensive and time-consuming laboratory stages of SNP validation. Additionally, the in silico approach leads to lower ascertainment bias in the resulting SNP panel as marker selection is based only on the ability to design primers and the predicted presence of intron-exon boundaries. Consequently SNPs with a wider spectrum of minor allele frequencies (MAFs) will be genotyped in the final panel. The genomic resources presented here represent a valuable multi-purpose resource for developing informative marker panels for population discrimination, microarray development and for population genomic studies in the wild.

Bases génétiques de la sténose valvulaire aortique calcifiée

La sténose valvulaire aortique (SVA) est une valvulopathie résultant en l'ouverture incomplète de la valve aortique. La calcification des feuillets associée au vieillissement est la cause la plus importante de la SVA. Sa pathogénèse implique des dépôts de lipoprotéines, de l'inflammation et de la calcification des feuillets. Notre étude vise à identifier les gènes associés à une prédisposition à la SVA afin de mieux comprendre les mécanismes sous-jacents à cette maladie et potentiellement identifier de nouvelles cibles thérapeutiques. Pour ce faire, nous avons recruté 190 patients avec SVA dégénérative et 192 témoins, appariés pour l'âge et le sexe, puis effectué une étude d’association par gènes candidats en utilisant des marqueurs génétiques polymorphiques (SNP). Les gènes candidats choisis incluent (1) ceux dont les polymorphismes ont été présumés associés à la SVA dans des études antérieures (APOB, APOE, ESR1, PTH et VDR) (2) des gènes dont les polymorphismes ont été significativement associés et validés pour quelques maladies inflammatoires (IL-10, TNFAIP3) ou pour le métabolisme lipidique (PCSK9, LDLR) dans des études d’association pangénomiques, et (3) des gènes impliqués dans la pathogénie de la SVA à partir d’études faites sur des modèles animaux en lien avec la calcification (BMP2, CCR5, CTGF, LRP5, MSX2, WNT3), le remodelage tissulaire (CTSS, MMP9) ou le métabolisme lipidique (SMPD1). Pour les gènes des groupes (1) et (2), nous avons utilisé les SNPs rapportés dans la littérature comme étant significativement associés. Pour le groupe (3), nous avons effectué une approche par «tagSNP» qui consiste à sélectionner un groupe de SNP capturant la variabilité génétique dans la région ciblée. Au total, 81 SNPs dans 18 gènes ont été testés. Nous avons trouvé une association nominale avec les gènes BMP2 (OR = 1.55, IC95%: 1.14-2.10, p = 0.004) et LRP5 (OR = 1.47, IC95%: 1.06-2.03, p = 0.023) après ajustement pour la maladie coronarienne. Les gènes BMP2 et LRP5, impliqués dans la calcification selon certains modèles expérimentaux, sont donc associés à la SVA. Ce travail devrait être validé dans une cohorte indépendante plus large dans un avenir rapproché et il pourrait être étendu à d’autres gènes.

Why are autism spectrum conditions more prevalent in males?

Autism Spectrum Conditions (ASC) are much more common in males, a bias that may offer clues to the etiology of this condition. Although the cause of this bias remains a mystery, we argue that it occurs because ASC is an extreme manifestation of the male brain. The extreme male brain (EMB) theory, first proposed in 1997, is an extension of the Empathizing-Systemizing (E-S) theory of typical sex differences that proposes that females on average have a stronger drive to empathize while males on average have a stronger drive to systemize. In this first major update since 2005, we describe some of the evidence relating to the EMB theory of ASC and consider how typical sex differences in brain structure may be relevant to ASC. One possible biological mechanism to account for the male bias is the effect of fetal testosterone (fT). We also consider alternative biological theories, the X and Y chromosome theories, and the reduced autosomal penetrance theory. None of these theories has yet been fully confirmed or refuted, though the weight of evidence in favor of the fT theory is growing from converging sources (longitudinal amniocentesis studies from pregnancy to age 10 years old, current hormone studies, and genetic association studies of SNPs in the sex steroid pathways). Ultimately, as these theories are not mutually exclusive and ASC is multi-factorial, they may help explain the male prevalence of ASC.

Variation in the autism candidate gene GABRB3 modulates tactile sensitivity in typically developing children

Background: Autism spectrum conditions have a strong genetic component. Atypical sensory sensitivities are one of the core but neglected features of autism spectrum conditions. GABRB3 is a well-characterised candidate gene for autism spectrum conditions. In mice, heterozygous Gabrb3 deletion is associated with increased tactile sensitivity. However, no study has examined if tactile sensitivity is associated with GABRB3 genetic variation in humans. To test this, we conducted two pilot genetic association studies in the general population, analysing two phenotypic measures of tactile sensitivity (a parent-report and a behavioural measure) for association with 43 SNPs in GABRB3. Findings: Across both tactile sensitivity measures, three SNPs (rs11636966, rs8023959 and rs2162241) were nominally associated with both phenotypes, providing a measure of internal validation. Parent-report scores were nominally associated with six SNPs (P <0.05). Behaviourally measured tactile sensitivity was nominally associated with 10 SNPs (three after Bonferroni correction). Conclusions: This is the first human study to show an association between GABRB3 variation and tactile sensitivity. This provides support for the evidence from animal models implicating the role of GABRB3 variation in the atypical sensory sensitivity in autism spectrum conditions. Future research is underway to directly test this association in cases of autism spectrum conditions.

TRIO LOGIC REGRESSION - DETECTION OF SNP - SNP INTERACTIONS IN CASE-PARENT TRIOS

Statistical approaches to evaluate higher order SNP-SNP and SNP-environment interactions are critical in genetic association studies, as susceptibility to complex disease is likely to be related to the interaction of multiple SNPs and environmental factors. Logic regression (Kooperberg et al., 2001; Ruczinski et al., 2003) is one such approach, where interactions between SNPs and environmental variables are assessed in a regression framework, and interactions become part of the model search space. In this manuscript we extend the logic regression methodology, originally developed for cohort and case-control studies, for studies of trios with affected probands. Trio logic regression accounts for the linkage disequilibrium (LD) structure in the genotype data, and accommodates missing genotypes via haplotype-based imputation. We also derive an efficient algorithm to simulate case-parent trios where genetic risk is determined via epistatic interactions.

Characteristics of a spina bifida population including North American Caucasian and Hispanic individuals.

BACKGROUND: Meningomyelocele (MM) is a common human birth defect. MM is a disorder of neural development caused by contributions from genes and environmental factors that result in the NTD and lead to a spectrum of physical and neurocognitive phenotypes. METHODS: A multidisciplinary approach has been taken to develop a comprehensive understanding of MM through collaborative efforts from investigators specializing in genetics, development, brain imaging, and neurocognitive outcome. Patients have been recruited from five different sites: Houston and the Texas-Mexico border area; Toronto, Canada; Los Angeles, California; and Lexington, Kentucky. Genetic risk factors for MM have been assessed by genotyping and association testing using the transmission disequilibrium test. RESULTS: A total of 509 affected child/parent trios and 309 affected child/parent duos have been enrolled to date for genetic association studies. Subsets of the patients have also been enrolled for studies assessing development, brain imaging, and neurocognitive outcomes. The study recruited two major ethnic groups, with 45.9% Hispanics of Mexican descent and 36.2% North American Caucasians of European descent. The remaining patients are African-American, South and Central American, Native American, and Asian. Studies of this group of patients have already discovered distinct corpus callosum morphology and neurocognitive deficits that associate with MM. We have identified maternal MTHFR 667T allele as a risk factor for MM. In addition, we also found that several genes for glucose transport and metabolism are potential risk factors for MM. CONCLUSIONS: The enrolled patient population provides a valuable resource for elucidating the disease characteristics and mechanisms for MM development.

Further evidence for a parent-of-origin effect at the NOP9 locus on language-related phenotypes

Background - Specific language impairment (SLI) is a common neurodevelopmental disorder, observed in 5–10 % of children. Family and twin studies suggest a strong genetic component, but relatively few candidate genes have been reported to date. A recent genome-wide association study (GWAS) described the first statistically significant association specifically for a SLI cohort between a missense variant (rs4280164) in the NOP9 gene and language-related phenotypes under a parent-of-origin model. Replications of these findings are particularly challenging because the availability of parental DNA is required. Methods - We used two independent family-based cohorts characterised with reading- and language-related traits: a longitudinal cohort (n = 106 informative families) including children with language and reading difficulties and a nuclear family cohort (n = 264 families) selected for dyslexia. Results - We observed association with language-related measures when modelling for parent-of-origin effects at the NOP9 locus in both cohorts: minimum P = 0.001 for phonological awareness with a paternal effect in the first cohort and minimum P = 0.0004 for irregular word reading with a maternal effect in the second cohort. Allelic and parental trends were not consistent when compared to the original study. Conclusions - A parent-of-origin effect at this locus was detected in both cohorts, albeit with different trends. These findings contribute in interpreting the original GWAS report and support further investigations of the NOP9 locus and its role in language-related traits. A systematic evaluation of parent-of-origin effects in genetic association studies has the potential to reveal novel mechanisms underlying complex traits.