Candidate Causal Regulatory Effects by Integration of Expression QTLs with Complex Trait Genetic Associations

The recent success of genome-wide association studies (GWAS) is now followed by the challenge to determine how the reported susceptibility variants mediate complex traits and diseases. Expression quantitative trait loci (eQTLs) have been implicated in disease associations through overlaps between eQTLs and GWAS signals. However, the abundance of eQTLs and the strong correlation structure (LD) in the genome make it likely that some of these overlaps are coincidental and not driven by the same functional variants. In the present study, we propose an empirical methodology, which we call Regulatory Trait Concordance (RTC) that accounts for local LD structure and integrates eQTLs and GWAS results in order to reveal the subset of association signals that are due to cis eQTLs. We simulate genomic regions of various LD patterns with both a single or two causal variants and show that our score outperforms SNP correlation metrics, be they statistical (r2) or historical (D'). Following the observation of a significant abundance of regulatory signals among currently published GWAS loci, we apply our method with the goal to prioritize relevant genes for each of the respective complex traits. We detect several potential disease-causing regulatory effects, with a strong enrichment for immunity-related conditions, consistent with the nature of the cell line tested (LCLs). Furthermore, we present an extension of the method in trans, where interrogating the whole genome for downstream effects of the disease variant can be informative regarding its unknown primary biological effect. We conclude that integrating cellular phenotype associations with organismal complex traits will facilitate the biological interpretation of the genetic effects on these traits

Coevolution of Interacting Fertilization Proteins

Reproductive proteins are among the fastest evolving in the proteome, often due to the consequences of positive selection, and their rapid evolution is frequently attributed to a coevolutionary process between interacting female and male proteins. Such a process could leave characteristic signatures at coevolving genes. One signature of coevolution, predicted by sexual selection theory, is an association of alleles between the two genes. Another predicted signature is a correlation of evolutionary rates during divergence due to compensatory evolution. We studied female–male coevolution in the abalone by resequencing sperm lysin and its interacting egg coat protein, VERL, in populations of two species. As predicted, we found intergenic linkage disequilibrium between lysin and VERL, despite our demonstration that they are not physically linked. This finding supports a central prediction of sexual selection using actual genotypes, that of an association between a male trait and its female preference locus. We also created a novel likelihood method to show that lysin and VERL have experienced correlated rates of evolution. These two signatures of coevolution can provide statistical rigor to hypotheses of coevolution and could be exploited for identifying coevolving proteins a priori. We also present polymorphism-based evidence for positive selection and implicate recent selective events at the specific structural regions of lysin and VERL responsible for their species-specific interaction. Finally, we observed deep subdivision between VERL alleles in one species, which matches a theoretical prediction of sexual conflict. Thus, abalone fertilization proteins illustrate how coevolution can lead to reproductive barriers and potentially drive speciation

Interaction of Crohn's Disease Susceptibility Genes in an Australian Paediatric Cohort

Genetic susceptibility is an important contributor to the pathogenesis of Crohn's disease (CD). We investigated multiple CD susceptibility genes in an Australian paediatric onset CD cohort. Newly diagnosed paediatric onset CD patients (n = 72) and controls (n = 98) were genotyped for 34 single nucleotide polymorphisms (SNPs) in 18 genetic loci. Gene-gene interaction analysis, gene-disease phenotype analysis and genetic risk profiling were performed for all SNPs and all genes. Of the 34 SNPs analysed, four polymorphisms on three genes (NOD2, IL23R, and region 3p21) were significantly associated with CD status (p<0.05). All three CD specific paediatric polymorphisms on PSMG1 and TNFRSF6B showed a trend of association with p<0.1. An additive gene-gene interaction involving TLR4, PSMG1, TNFRSF6B and IRGM was identified with CD. Genes involved in microbial processing (TLR4, PSMG1, NOD2) were significantly associated either at the individual level or in gene-gene interactive roles. Colonic disease was significantly associated with disease SNP rs7517847 (IL23R) (p<0.05) and colonic and ileal/colonic disease was significantly associated with disease SNP rs125221868 (IBD5) and SLC22A4 & SLC22A4/5 variants (p<0.05). We were able to demonstrate genetic association of several genes to CD in a paediatric onset cohort. Several of the observed associations have not been reported previously in association with paediatric CD patients. Our findings demonstrate that CD genetic susceptibility in paediatric patients presents as a complex interaction between numerous genes

A Survey of Bayesian Statistical Approaches for Big Data

The modern era is characterised as an era of information or Big Data. This has motivated a huge literature on new methods for extracting information and insights from these data. A natural question is how these approaches differ from those that were available prior to the advent of Big Data. We present a review of published studies that present Bayesian statistical approaches specifically for Big Data and discuss the reported and perceived benefits of these approaches. We conclude by addressing the question of whether focusing only on improving computational algorithms and infrastructure will be enough to face the challenges of Big Data

Exploration of Shared Genetic Architecture Between Subcortical Brain Volumes and Anorexia Nervosa

In MRI scans of patients with anorexia nervosa (AN), reductions in brain volume are often apparent. However, it is unknown whether such brain abnormalities are influenced by genetic determinants that partially overlap with those underlying AN. Here, we used a battery of methods (LD score regression, genetic risk scores, sign test, SNP effect concordance analysis, and Mendelian randomization) to investigate the genetic covariation between subcortical brain volumes and risk for AN based on summary measures retrieved from genome-wide association studies of regional brain volumes (ENIGMA consortium, n = 13,170) and genetic risk for AN (PGC-ED consortium, n = 14,477). Genetic correlations ranged from − 0.10 to 0.23 (all p > 0.05). There were some signs of an inverse concordance between greater thalamus volume and risk for AN (permuted p = 0.009, 95% CI: [0.005, 0.017]). A genetic variant in the vicinity of ZW10, a gene involved in cell division, and neurotransmitter and immune system relevant genes, in particular DRD2, was significantly associated with AN only after conditioning on its association with caudate volume (pFDR = 0.025). Another genetic variant linked to LRRC4C, important in axonal and synaptic development, reached significance after conditioning on hippocampal volume (pFDR = 0.021). In this comprehensive set of analyses and based on the largest available sample sizes to date, there was weak evidence for associations between risk for AN and risk for abnormal subcortical brain volumes at a global level (that is, common variant genetic architecture), but suggestive evidence for effects of single genetic markers. Highly powered multimodal brain- and disorder-related genome-wide studies are needed to further dissect the shared genetic influences on brain structure and risk for AN