A genome-wide association study identifies multiple loci for variation in human ear morphology

Here we report a genome-wide association study for non-pathological pinna morphology in over 5,000 Latin Americans. We find genome-wide significant association at seven genomic regions affecting: lobe size and attachment, folding of antihelix, helix rolling, ear protrusion and antitragus size (linear regression P values 2 × 10−8 to 3 × 10−14). Four traits are associated with a functional variant in the Ectodysplasin A receptor (EDAR) gene, a key regulator of embryonic skin appendage development. We confirm expression of Edar in the developing mouse ear and that Edar-deficient mice have an abnormally shaped pinna. Two traits are associated with SNPs in a region overlapping the T-Box Protein 15 (TBX15) gene, a major determinant of mouse skeletal development. Strongest association in this region is observed for SNP rs17023457 located in an evolutionarily conserved binding site for the transcription factor Cartilage paired-class homeoprotein 1 (CART1), and we confirm that rs17023457 alters in vitro binding of CART1

A Genome-Wide Association Study of Neuroticism in a Population-Based Sample

Neuroticism is a moderately heritable personality trait considered to be a risk factor for developing major depression, anxiety disorders and dementia. We performed a genome-wide association study in 2,235 participants drawn from a population-based study of neuroticism, making this the largest association study for neuroticism to date. Neuroticism was measured by the Eysenck Personality Questionnaire. After Quality Control, we analysed 430,000 autosomal SNPs together with an additional 1.2 million SNPs imputed with high quality from the Hap Map CEU samples. We found a very small effect of population stratification, corrected using one principal component, and some cryptic kinship that required no correction. NKAIN2 showed suggestive evidence of association with neuroticism as a main effect (p<10−6) and GPC6 showed suggestive evidence for interaction with age (p≈10−7). We found support for one previously-reported association (PDE4D), but failed to replicate other recent reports. These results suggest common SNP variation does not strongly influence neuroticism. Our study was powered to detect almost all SNPs explaining at least 2% of heritability, and so our results effectively exclude the existence of loci having a major effect on neuroticism

A Whole-Genome SNP Association Study of NCI60 Cell Line Panel Indicates a Role of Ca2+ Signaling in Selenium Resistance

Epidemiological studies have suggested an association between selenium intake and protection from a variety of cancer. Considering this clinical importance of selenium, we aimed to identify the genes associated with resistance to selenium treatment. We have applied a previous methodology developed by our group, which is based on the genetic and pharmacological data publicly available for the NCI60 cancer cell line panel. In short, we have categorized the NCI60 cell lines as selenium resistant and sensitive based on their growth inhibition (GI50) data. Then, we have utilized the Affymetrix 125K SNP chip data available and carried out a genome-wide case-control association study for the selenium sensitive and resistant NCI60 cell lines. Our results showed statistically significant association of four SNPs in 5q33–34, 10q11.2, 10q22.3 and 14q13.1 with selenium resistance. These SNPs were located in introns of the genes encoding for a kinase-scaffolding protein (AKAP6), a membrane protein (SGCD), a channel protein (KCNMA1), and a protein kinase (PRKG1). The knock-down of KCNMA1 by siRNA showed increased sensitivity to selenium in both LNCaP and PC3 cell lines. Furthermore, SNP-SNP interaction (epistasis) analysis indicated the interactions of the SNPs in AKAP6 with SGCD as well as SNPs in AKAP6 with KCNMA1 with each other, assuming additive genetic model. These genes were also all involved in the Ca2+ signaling, which has a direct role in induction of apoptosis and induction of apoptosis in tumor cells is consistent with the chemopreventive action of selenium. Once our findings are further validated, this knowledge can be translated into clinics where individuals who can benefit from the chemopreventive characteristics of the selenium supplementation will be easily identified using a simple DNA analysis

Human Genetics in Rheumatoid Arthritis Guides a High-Throughput Drug Screen of the CD40 Signaling Pathway

Although genetic and non-genetic studies in mouse and human implicate the CD40 pathway in rheumatoid arthritis (RA), there are no approved drugs that inhibit CD40 signaling for clinical care in RA or any other disease. Here, we sought to understand the biological consequences of a CD40 risk variant in RA discovered by a previous genome-wide association study (GWAS) and to perform a high-throughput drug screen for modulators of CD40 signaling based on human genetic findings. First, we fine-map the CD40 risk locus in 7,222 seropositive RA patients and 15,870 controls, together with deep sequencing of CD40 coding exons in 500 RA cases and 650 controls, to identify a single SNP that explains the entire signal of association (rs4810485, P = 1.4×10(−9)). Second, we demonstrate that subjects homozygous for the RA risk allele have ∼33% more CD40 on the surface of primary human CD19+ B lymphocytes than subjects homozygous for the non-risk allele (P = 10(−9)), a finding corroborated by expression quantitative trait loci (eQTL) analysis in peripheral blood mononuclear cells from 1,469 healthy control individuals. Third, we use retroviral shRNA infection to perturb the amount of CD40 on the surface of a human B lymphocyte cell line (BL2) and observe a direct correlation between amount of CD40 protein and phosphorylation of RelA (p65), a subunit of the NF-κB transcription factor. Finally, we develop a high-throughput NF-κB luciferase reporter assay in BL2 cells activated with trimerized CD40 ligand (tCD40L) and conduct an HTS of 1,982 chemical compounds and FDA–approved drugs. After a series of counter-screens and testing in primary human CD19+ B cells, we identify 2 novel chemical inhibitors not previously implicated in inflammation or CD40-mediated NF-κB signaling. Our study demonstrates proof-of-concept that human genetics can be used to guide the development of phenotype-based, high-throughput small-molecule screens to identify potential novel therapies in complex traits such as RA

RNA-Seq Identifies SNP Markers for Growth Traits in Rainbow Trout

Fast growth is an important and highly desired trait, which affects the profitability of food animal production, with feed costs accounting for the largest proportion of production costs. Traditional phenotype-based selection is typically used to select for growth traits; however, genetic improvement is slow over generations. Single nucleotide polymorphisms (SNPs) explain 90% of the genetic differences between individuals; therefore, they are most suitable for genetic evaluation and strategies that employ molecular genetics for selective breeding. SNPs found within or near a coding sequence are of particular interest because they are more likely to alter the biological function of a protein. We aimed to use SNPs to identify markers and genes associated with genetic variation in growth. RNA-Seq whole-transcriptome analysis of pooled cDNA samples from a population of rainbow trout selected for improved growth versus unselected genetic cohorts (10 fish from 1 full-sib family each) identified SNP markers associated with growth-rate. The allelic imbalances (the ratio between the allele frequencies of the fast growing sample and that of the slow growing sample) were considered at scores >5.0 as an amplification and <0.2 as loss of heterozygosity. A subset of SNPs (n = 54) were validated and evaluated for association with growth traits in 778 individuals of a three-generation parent/offspring panel representing 40 families. Twenty-two SNP markers and one mitochondrial haplotype were significantly associated with growth traits. Polymorphism of 48 of the markers was confirmed in other commercially important aquaculture stocks. Many markers were clustered into genes of metabolic energy production pathways and are suitable candidates for genetic selection. The study demonstrates that RNA-Seq at low sequence coverage of divergent populations is a fast and effective means of identifying SNPs, with allelic imbalances between phenotypes. This technique is suitable for marker development in non-model species lacking complete and well-annotated genome reference sequences

Effects of the Cryptochrome CryB from Rhodobacter sphaeroides on Global Gene Expression in the Dark or Blue Light or in the Presence of Singlet Oxygen

Several regulators are controlling the formation of the photosynthetic apparatus in the facultatively photosynthetic bacterium Rhodobacter sphaeroides. Among the proteins affecting photosynthesis gene expression is the blue light photoreceptor cryptochrome CryB. This study addresses the effect of CryB on global gene expression. The data reveal that CryB does not only influence photosynthesis gene expression but also genes for the non-photosynthetic energy metabolism like citric acid cycle and oxidative phosphorylation. In addition several genes involved in RNA processing and in transcriptional regulation are affected by a cryB deletion. Although CryB was shown to undergo a photocycle it does not only affect gene expression in response to blue light illumination but also in response to singlet oxygen stress conditions. While there is a large overlap in these responses, some CryB-dependent effects are specific for blue-light or photooxidative stress. In addition to protein-coding genes some genes for sRNAs show CryB-dependent expression. These findings give new insight into the function of bacterial cryptochromes and demonstrate for the first time a function in the oxidative stress response

Genetic variability of the grey wolf Canis lupus in the Caucasus in comparison with Europe and the Middle East: distinct or intermediary population?

Despite continuous historical distribution of the grey wolf (Canis lupus) throughout Eurasia, the species displays considerable morphological differentiation that resulted in delimitation of a number of subspecies. However, these morphological discontinuities are not always consistent with patterns of genetic differentiation. Here we assess genetic distinctiveness of grey wolves from the Caucasus (a region at the border between Europe and West Asia) that have been classified as a distinct subspecies C. l. cubanensis. We analysed their genetic variability based on mtDNA control region, microsatellite loci and genome-wide SNP genotypes (obtained for a subset of the samples), and found similar or higher levels of genetic diversity at all these types of loci as compared with other Eurasian populations. Although we found no evidence for a recent genetic bottleneck, genome-wide linkage disequilibrium patterns suggest a long-term demographic decline in the Caucasian population – a trend consistent with other Eurasian populations. Caucasian wolves share mtDNA haplotypes with both Eastern European and West Asian wolves, suggesting past or ongoing gene flow. Microsatellite data also suggest gene flow between the Caucasus and Eastern Europe. We found evidence for moderate admixture between the Caucasian wolves and domestic dogs, at a level comparable with other Eurasian populations. Taken together, our results show that Caucasian wolves are not genetically isolated from other Eurasian populations, share with them the same demographic trends, and are affected by similar conservation problems

A High Density SNP Array for the Domestic Horse and Extant Perissodactyla: Utility for Association Mapping, Genetic Diversity, and Phylogeny Studies

An equine SNP genotyping array was developed and evaluated on a panel of samples representing 14 domestic horse breeds and 18 evolutionarily related species. More than 54,000 polymorphic SNPs provided an average inter-SNP spacing of ∼43 kb. The mean minor allele frequency across domestic horse breeds was 0.23, and the number of polymorphic SNPs within breeds ranged from 43,287 to 52,085. Genome-wide linkage disequilibrium (LD) in most breeds declined rapidly over the first 50–100 kb and reached background levels within 1–2 Mb. The extent of LD and the level of inbreeding were highest in the Thoroughbred and lowest in the Mongolian and Quarter Horse. Multidimensional scaling (MDS) analyses demonstrated the tight grouping of individuals within most breeds, close proximity of related breeds, and less tight grouping in admixed breeds. The close relationship between the Przewalski's Horse and the domestic horse was demonstrated by pair-wise genetic distance and MDS. Genotyping of other Perissodactyla (zebras, asses, tapirs, and rhinoceros) was variably successful, with call rates and the number of polymorphic loci varying across taxa. Parsimony analysis placed the modern horse as sister taxa to Equus przewalski. The utility of the SNP array in genome-wide association was confirmed by mapping the known recessive chestnut coat color locus (MC1R) and defining a conserved haplotype of ∼750 kb across all breeds. These results demonstrate the high quality of this SNP genotyping resource, its usefulness in diverse genome analyses of the horse, and potential use in related species