Accounting for Population Structure in Gene-by-Environment Interactions in Genome-Wide Association Studies Using Mixed Models

Although genome-wide association studies (GWASs) have discovered numerous novel genetic variants associated with many complex traits and diseases, those genetic variants typically explain only a small fraction of phenotypic variance. Factors that account for phenotypic variance include environmental factors and gene-by-environment interactions (GEIs). Recently, several studies have conducted genome-wide gene-by-environment association analyses and demonstrated important roles of GEIs in complex traits. One of the main challenges in these association studies is to control effects of population structure that may cause spurious associations. Many studies have analyzed how population structure influences statistics of genetic variants and developed several statistical approaches to correct for population structure. However, the impact of population structure on GEI statistics in GWASs has not been extensively studied and nor have there been methods designed to correct for population structure on GEI statistics. In this paper, we show both analytically and empirically that population structure may cause spurious GEIs and use both simulation and two GWAS datasets to support our finding. We propose a statistical approach based on mixed models to account for population structure on GEI statistics. We find that our approach effectively controls population structure on statistics for GEIs as well as for genetic variants.     

High-Field fMRI Reveals Brain Activation Patterns Underlying Saccade Execution in the Human Superior Colliculus

Background

The superior colliculus (SC) has been shown to play a crucial role in the initiation and coordination of eye- and head-movements. The knowledge about the function of this structure is mainly based on single-unit recordings in animals with relatively few neuroimaging studies investigating eye-movement related brain activity in humans.

Methodology/Principal Findings

The present study employed high-field (7 Tesla) functional magnetic resonance imaging (fMRI) to investigate SC responses during endogenously cued saccades in humans. In response to centrally presented instructional cues, subjects either performed saccades away from (centrifugal) or towards (centripetal) the center of straight gaze or maintained fixation at the center position. Compared to central fixation, the execution of saccades elicited hemodynamic activity within a network of cortical and subcortical areas that included the SC, lateral geniculate nucleus (LGN), occipital cortex, striatum, and the pulvinar.

Conclusions/Significance

Activity in the SC was enhanced contralateral to the direction of the saccade (i.e., greater activity in the right as compared to left SC during leftward saccades and vice versa) during both centrifugal and centripetal saccades, thereby demonstrating that the contralateral predominance for saccade execution that has been shown to exist in animals is also present in the human SC. In addition, centrifugal saccades elicited greater activity in the SC than did centripetal saccades, while also being accompanied by an enhanced deactivation within the prefrontal default-mode network. This pattern of brain activity might reflect the reduced processing effort required to move the eyes toward as compared to away from the center of straight gaze, a position that might serve as a spatial baseline in which the retinotopic and craniotopic reference frames are aligned.     

Phylogeography of the large Myotis bats (Chiroptera: Vespertilionidae) in Europe,Asia Minor,and Transcaucasia

The large Myotis complex in continental Europe, Asia Minor, and Transcaucasia comprises two sibling bat species, the greater mouse‐eared bat, Myotis myotis, and the lesser mouse‐eared bat, Myotis blythii, also referred to as Myotis oxygnathus. Here, we investigate the phylogeography of these bats using two mitochondrial markers: the second hypervariable domain of the control region (HVII) and a fragment of the cytochrome b gene (cyt b). The HVII haplotypes formed six distinct haplogroups associated with different geographical regions. Most of the European HVII haplotypes were exclusive to M. myotis, whereas the majority of HVII haplotypes found in Asia Minor were exclusive to M. blythii/M. oxygnathus. The phylogenetic reconstruction based on the concatenated cyt b and HVII fragments recovered two major lineages. The first lineage comprised samples from Europe (western lineage), and the second lineage included samples from Asia Minor, Transcaucasia, Crimea, Western Ukraine, Thrace, the Balkans, and Eastern Europe (eastern lineage). The mitochondrial lineage of M. blythii, reported from Kyrgyzstan, was not present in Asia Minor and Transcaucasia. Therefore, we consider the possibility that the M. blythii/M. oxygnathus found in Europe, Asia Minor, and Transcaucasia are not recent descendants of the Central Asian M. blythii. Instead, we suggest that M. blythii/M. oxygnathus and M. myotis diverged through allopatric speciation in Asia Minor and Europe, and that they are represented by the eastern and western mitochondrial lineages. We also examine an alternative hypothesis: that the large Myotis complex consists of more than two species that diverged independently in Asia Minor and Europe through ecological speciation. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ??, ??–??.     

Mouse genome-wide association and systems genetics identify Asxl2 as a regulator of bone mineral density and osteoclastogenesis

Significant advances have been made in the discovery of genes affecting bone mineral density (BMD); however, our understanding of its genetic basis remains incomplete. In the current study, genome-wide association (GWA) and co-expression network analysis were used in the recently described Hybrid Mouse Diversity Panel (HMDP) to identify and functionally characterize novel BMD genes. In the HMDP, a GWA of total body, spinal, and femoral BMD revealed four significant associations (-log10P>5.39) affecting at least one BMD trait on chromosomes (Chrs.) 7, 11, 12, and 17. The associations implicated a total of 163 genes with each association harboring between 14 and 112 genes. This list was reduced to 26 functional candidates by identifying those genes that were regulated by local eQTL in bone or harbored potentially functional non-synonymous (NS) SNPs. This analysis revealed that the most significant BMD SNP on Chr. 12 was a NS SNP in the additional sex combs like-2 (Asxl2) gene that was predicted to be functional. The involvement of Asxl2 in the regulation of bone mass was confirmed by the observation that Asxl2 knockout mice had reduced BMD. To begin to unravel the mechanism through which Asxl2 influenced BMD, a gene co-expression network was created using cortical bone gene expression microarray data from the HMDP strains. Asxl2 was identified as a member of a co-expression module enriched for genes involved in the differentiation of myeloid cells. In bone, osteoclasts are bone-resorbing cells of myeloid origin, suggesting that Asxl2 may play a role in osteoclast differentiation. In agreement, the knockdown of Asxl2 in bone marrow macrophages impaired their ability to form osteoclasts. This study identifies a new regulator of BMD and osteoclastogenesis and highlights the power of GWA and systems genetics in the mouse for dissecting complex genetic traits.     

Purification and Characterization of Carbonic Anhydrase from Ağrı Balık Lake Trout Gill (Salmo trutta labrax) and Effects of Sulfonamides on Enzyme Activity

Carbonic anhydrase (CA) was purified from A?r? Bal?k Lake trout gill (fCA) by affinity chromatography on a sepharose 4B‐tyrosine‐sulfanilamide column. The fCA enzyme was purified with about a 303.9 purification factor, a specific activity 4130.4 EU (mg‐protein)^–1, and a yield of 79.3 by using sepharose‐4B‐l tyrosine‐sulfanilamide affinity gel chromatography. The molecular weight determined by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS–PAGE) was found to be about 29.9 kDa. The kinetic parameters, K_M and V_max were determined for the 4‐nitrophenyl acetate hydrolysis reaction. Some sulfonamides were tested as inhibitors against the purified CA enzymes. The K_i constants for mafenide ( 1 ), p‐toluenesulfonamide ( 2 ), 2‐bromo‐benzene sulfonamide ( 3 ), 4‐chlorobenzene sulfonamide ( 4 ), 4‐amino‐6‐chloro‐1–3 benzenedisulfonamide ( 5 ), sulfamethazine ( 6 ), sulfaguanidine ( 7 ), sulfadiazine ( 8 ), and acetozazolamide ( 9 ) were in the range of 7.5–108.75 μM.