Genetik Ufku

by Fengzhi Jin, Hong Liu, Ping Li, Hong-Guo Yu, Yanchang Wang

The attachment of sister kinetochores by microtubules emanating from opposite spindle poles establishes chromosome bipolar attachment, which generates tension on chromosomes and is essential for sister-chromatid segregation. Syntelic attachment occurs when both sister kinetochores are attached by microtubules from the same spindle pole and this attachment is unable to generate tension on chromosomes, but a reliable method to induce syntelic attachments is not available in budding yeast. The spindle checkpoint can sense the lack of tension on chromosomes as well as detached kinetochores to prevent anaphase onset. In budding yeast Saccharomyces cerevisiae, tension checkpoint proteins Aurora/Ipl1 kinase and centromere-localized Sgo1 are required to sense the absence of tension but are dispensable for the checkpoint response to detached kinetochores. We have found that the loss of function of a motor protein complex Cik1/Kar3 in budding yeast leads to syntelic attachments. Inactivation of either the spindle or tension checkpoint enables premature anaphase entry in cells with dysfunctional Cik1/Kar3, resulting in co-segregation of sister chromatids. Moreover, the abolished Kar3-kinetochore interaction in cik1 mutants suggests that the Cik1/Kar3 complex mediates chromosome movement along microtubules, which could facilitate bipolar attachment. Therefore, we can induce syntelic attachments in budding yeast by inactivating the Cik1/Kar3 complex, and this approach will be very useful to study the checkpoint response to syntelic attachments.

by Martin Ladouceur, Zari Dastani, Yurii S. Aulchenko, Celia M. T. Greenwood, J. Brent Richards

The role of rare genetic variation in the etiology of complex disease remains unclear. However, the development of next-generation sequencing technologies offers the experimental opportunity to address this question. Several novel statistical methodologies have been recently proposed to assess the contribution of rare variation to complex disease etiology. Nevertheless, no empirical estimates comparing their relative power are available. We therefore assessed the parameters that influence their statistical power in 1,998 individuals Sanger-sequenced at seven genes by modeling different distributions of effect, proportions of causal variants, and direction of the associations (deleterious, protective, or both) in simulated continuous trait and case/control phenotypes. Our results demonstrate that the power of recently proposed statistical methods depend strongly on the underlying hypotheses concerning the relationship of phenotypes with each of these three factors. No method demonstrates consistently acceptable power despite this large sample size, and the performance of each method depends upon the underlying assumption of the relationship between rare variants and complex traits. Sensitivity analyses are therefore recommended to compare the stability of the results arising from different methods, and promising results should be replicated using the same method in an independent sample. These findings provide guidance in the analysis and interpretation of the role of rare base-pair variation in the etiology of complex traits and diseases.

by Victoria Campuzano, Maria Segura-Puimedon, Verena Terrado, Carolina Sánchez-Rodríguez, Mathilde Coustets, Mauricio Menacho-Márquez, Julián Nevado, Xosé R. Bustelo, Uta Francke, Luis A. Pérez-Jurado

A hallmark feature of Williams-Beuren Syndrome (WBS) is a generalized arteriopathy due to elastin deficiency, presenting as stenoses of medium and large arteries and leading to hypertension and other cardiovascular complications. Deletion of a functional NCF1 gene copy has been shown to protect a proportion of WBS patients against hypertension, likely through reduced NADPH-oxidase (NOX)–mediated oxidative stress. DD mice, carrying a 0.67 Mb heterozygous deletion including the Eln gene, presented with a generalized arteriopathy, hypertension, and cardiac hypertrophy, associated with elevated angiotensin II (angII), oxidative stress parameters, and Ncf1 expression. Genetic (by crossing with Ncf1 mutant) and/or pharmacological (with ang II type 1 receptor blocker, losartan, or NOX inhibitor apocynin) reduction of NOX activity controlled hormonal and biochemical parameters in DD mice, resulting in normalized blood pressure and improved cardiovascular histology. We provide strong evidence for implication of the redox system in the pathophysiology of the cardiovascular disease in a mouse model of WBS. The phenotype of these mice can be ameliorated by either genetic or pharmacological intervention reducing NOX activity, likely through reduced angII–mediated oxidative stress. Therefore, anti-NOX therapy merits evaluation to prevent the potentially serious cardiovascular complications of WBS, as well as in other cardiovascular disorders mediated by similar pathogenic mechanism.

by Erin K. Butler, Aaron Voigt, A. Kathrin Lutz, Jane P. Toegel, Ellen Gerhardt, Peter Karsten, Björn Falkenburger, Andrea Reinartz, Konstanze F. Winklhofer, Jörg B. Schulz

Overexpression or mutation of α-Synuclein is associated with protein aggregation and interferes with a number of cellular processes, including mitochondrial integrity and function. We used a whole-genome screen in the fruit fly Drosophila melanogaster to search for novel genetic modifiers of human [A53T]α-Synuclein–induced neurotoxicity. Decreased expression of the mitochondrial chaperone protein tumor necrosis factor receptor associated protein-1 (TRAP1) was found to enhance age-dependent loss of fly head dopamine (DA) and DA neuron number resulting from [A53T]α-Synuclein expression. In addition, decreased TRAP1 expression in [A53T]α-Synuclein–expressing flies resulted in enhanced loss of climbing ability and sensitivity to oxidative stress. Overexpression of human TRAP1 was able to rescue these phenotypes. Similarly, human TRAP1 overexpression in rat primary cortical neurons rescued [A53T]α-Synuclein–induced sensitivity to rotenone treatment. In human (non)neuronal cell lines, small interfering RNA directed against TRAP1 enhanced [A53T]α-Synuclein–induced sensitivity to oxidative stress treatment. [A53T]α-Synuclein directly interfered with mitochondrial function, as its expression reduced Complex I activity in HEK293 cells. These effects were blocked by TRAP1 overexpression. Moreover, TRAP1 was able to prevent alteration in mitochondrial morphology caused by [A53T]α-Synuclein overexpression in human SH-SY5Y cells. These results indicate that [A53T]α-Synuclein toxicity is intimately connected to mitochondrial dysfunction and that toxicity reduction in fly and rat primary neurons and human cell lines can be achieved using overexpression of the mitochondrial chaperone TRAP1. Interestingly, TRAP1 has previously been shown to be phosphorylated by the serine/threonine kinase PINK1, thus providing a potential link of PINK1 via TRAP1 to α-Synuclein.

by Maheen Ferdous, James D. Higgins, Kim Osman, Christophe Lambing, Elisabeth Roitinger, Karl Mechtler, Susan J. Armstrong, Ruth Perry, Mónica Pradillo, Nieves Cuñado, F. Chris H. Franklin

In this study we have analysed AtASY3, a coiled-coil domain protein that is required for normal meiosis in Arabidopsis. Analysis of an Atasy3-1 mutant reveals that loss of the protein compromises chromosome axis formation and results in reduced numbers of meiotic crossovers (COs). Although the frequency of DNA double-strand breaks (DSBs) appears moderately reduced in Atasy3-1, the main recombination defect is a reduction in the formation of COs. Immunolocalization studies in wild-type meiocytes indicate that the HORMA protein AtASY1, which is related to Hop1 in budding yeast, forms hyper-abundant domains along the chromosomes that are spatially associated with DSBs and early recombination pathway proteins. Loss of AtASY3 disrupts the axial organization of AtASY1. Furthermore we show that the AtASY3 and AtASY1 homologs BoASY3 and BoASY1, from the closely related species Brassica oleracea, are co-immunoprecipitated from meiocyte extracts and that AtASY3 interacts with AtASY1 via residues in its predicted coiled-coil domain. Together our results suggest that AtASY3 is a functional homolog of Red1. Since studies in budding yeast indicate that Red1 and Hop1 play a key role in establishing a bias to favor inter-homolog recombination (IHR), we propose that AtASY3 and AtASY1 may have a similar role in Arabidopsis. Loss of AtASY3 also disrupts synaptonemal complex (SC) formation. In Atasy3-1 the transverse filament protein AtZYP1 forms small patches rather than a continuous SC. The few AtMLH1 foci that remain in Atasy3-1 are found in association with the AtZYP1 patches. This is sufficient to prevent the ectopic recombination observed in the absence of AtZYP1, thus emphasizing that in addition to its structural role the protein is important for CO formation.

by Xinhe Huang, Jun Liu, Robert C. Dickson

Knowledge of the mechanisms for regulating lifespan is advancing rapidly, but lifespan is a complex phenotype and new features are likely to be identified. Here we reveal a novel approach for regulating lifespan. Using a genetic or a pharmacological strategy to lower the rate of sphingolipid synthesis, we show that Saccharomyces cerevisiae cells live longer. The longer lifespan is due in part to a reduction in Sch9 protein kinase activity and a consequent reduction in chromosomal mutations and rearrangements and increased stress resistance. Longer lifespan also arises in ways that are independent of Sch9 or caloric restriction, and we speculate on ways that sphingolipids might mediate these aspects of increased lifespan. Sch9 and its mammalian homolog S6 kinase work downstream of the target of rapamycin, TOR1, protein kinase, and play evolutionarily conserved roles in regulating lifespan. Our data establish Sch9 as a focal point for regulating lifespan by integrating nutrient signals from TOR1 with growth and stress signals from sphingolipids. Sphingolipids are found in all eukaryotes and our results suggest that pharmacological down-regulation of one or more sphingolipids may provide a means to reduce age-related diseases and increase lifespan in other eukaryotes.

by Carlos Díaz-Castillo, Xiao-Qin Xia, José M. Ranz

Why gene order is conserved over long evolutionary timespans remains elusive. A common interpretation is that gene order conservation might reflect the existence of functional constraints that are important for organismal performance. Alteration of the integrity of genomic regions, and therefore of those constraints, would result in detrimental effects. This notion seems especially plausible in those genomes that can easily accommodate gene reshuffling via chromosomal inversions since genomic regions free of constraints are likely to have been disrupted in one or more lineages. Nevertheless, no empirical test has been performed to this notion. Here, we disrupt one of the largest conserved genomic regions of the Drosophila genome by chromosome engineering and examine the phenotypic consequences derived from such disruption. The targeted region exhibits multiple patterns of functional enrichment suggestive of the presence of constraints. The carriers of the disrupted collinear block show no defects in their viability, fertility, and parameters of general homeostasis, although their odorant perception is altered. This change in odorant perception does not correlate with modifications of the level of expression and sex bias of the genes within the genomic region disrupted. Our results indicate that even in highly rearranged genomes, like those of Diptera, unusually high levels of gene order conservation cannot be systematically attributed to functional constraints, which raises the possibility that other mechanisms can be in place and therefore the underpinnings of the maintenance of gene organization might be more diverse than previously thought.

by Gregory N. Adams, Alvin H. Schmaier

by Alessia Buscaino, Sharon A. White, Douglas R. Houston, Erwan Lejeune, Femke Simmer, Flavia de Lima Alves, Piyush T. Diyora, Takeshi Urano, Elizabeth H. Bayne, Juri Rappsilber, Robin C. Allshire

Non-coding transcription can trigger histone post-translational modifications forming specialized chromatin. In fission yeast, heterochromatin formation requires RNAi and the histone H3K9 methyltransferase complex CLRC, composed of Clr4, Raf1, Raf2, Cul4, and Rik1. CLRC mediates H3K9 methylation and siRNA production; it also displays E3-ubiquitin ligase activity in vitro. DCAFs act as substrate receptors for E3 ligases and may couple ubiquitination with histone methylation. Here, structural alignment and mutation of signature WDxR motifs in Raf1 indicate that it is a DCAF for CLRC. We demonstrate that Raf1 promotes H3K9 methylation and siRNA amplification via two distinct, separable functions. The association of the DCAF Raf1 with Cul4-Rik1 is critical for H3K9 methylation, but dispensable for processing of centromeric transcripts into siRNAs. Thus the association of a DCAF, Raf1, with its adaptor, Rik1, is required for histone methylation and to allow RNAi to signal to chromatin.

by Brian N. Pasley, Stephen V. David, Nima Mesgarani, Adeen Flinker, Shihab A. Shamma, Nathan E. Crone, Robert T. Knight, Edward F. Chang

How the human auditory system extracts perceptually relevant acoustic features of speech is unknown. To address this question, we used intracranial recordings from nonprimary auditory cortex in the human superior temporal gyrus to determine what acoustic information in speech sounds can be reconstructed from population neural activity. We found that slow and intermediate temporal fluctuations, such as those corresponding to syllable rate, were accurately reconstructed using a linear model based on the auditory spectrogram. However, reconstruction of fast temporal fluctuations, such as syllable onsets and offsets, required a nonlinear sound representation based on temporal modulation energy. Reconstruction accuracy was highest within the range of spectro-temporal fluctuations that have been found to be critical for speech intelligibility. The decoded speech representations allowed readout and identification of individual words directly from brain activity during single trial sound presentations. These findings reveal neural encoding mechanisms of speech acoustic parameters in higher order human auditory cortex.

by Paul Beier, Andrew J. Gregory

by Joe H. Levine, Michelle E. Fontes, Jonathan Dworkin, Michael B. Elowitz

Environmental signals induce diverse cellular differentiation programs. In certain systems, cells defer differentiation for extended time periods after the signal appears, proliferating through multiple rounds of cell division before committing to a new fate. How can cells set a deferral time much longer than the cell cycle? Here we study Bacillus subtilis cells that respond to sudden nutrient limitation with multiple rounds of growth and division before differentiating into spores. A well-characterized genetic circuit controls the concentration and phosphorylation of the master regulator Spo0A, which rises to a critical concentration to initiate sporulation. However, it remains unclear how this circuit enables cells to defer sporulation for multiple cell cycles. Using quantitative time-lapse fluorescence microscopy of Spo0A dynamics in individual cells, we observed pulses of Spo0A phosphorylation at a characteristic cell cycle phase. Pulse amplitudes grew systematically and cell-autonomously over multiple cell cycles leading up to sporulation. This pulse growth required a key positive feedback loop involving the sporulation kinases, without which the deferral of sporulation became ultrasensitive to kinase expression. Thus, deferral is controlled by a pulsed positive feedback loop in which kinase expression is activated by pulses of Spo0A phosphorylation. This pulsed positive feedback architecture provides a more robust mechanism for setting deferral times than constitutive kinase expression. Finally, using mathematical modeling, we show how pulsing and time delays together enable “polyphasic” positive feedback, in which different parts of a feedback loop are active at different times. Polyphasic feedback can enable more accurate tuning of long deferral times. Together, these results suggest that Bacillus subtilis uses a pulsed positive feedback loop to implement a “timer” that operates over timescales much longer than a cell cycle.

by Liza Gross

by Yorann Baron, Sébastien Corre, Nicolas Mouchet, Sophie Vaulont, Sharon Prince, Marie-Dominique Galibert

An important function of all organisms is to ensure that their genetic material remains intact and unaltered through generations. This is an extremely challenging task since the cell's DNA is constantly under assault by endogenous and environmental agents. To protect against this, cells have evolved effective mechanisms to recognize DNA damage, signal its presence, and mediate its repair. While these responses are expected to be highly regulated because they are critical to avoid human diseases, very little is known about the regulation of the expression of genes involved in mediating their effects. The Nucleotide Excision Repair (NER) is the major DNA–repair process involved in the recognition and removal of UV-mediated DNA damage. Here we use a combination of in vitro and in vivo assays with an intermittent UV-irradiation protocol to investigate the regulation of key players in the DNA–damage recognition step of NER sub-pathways (TCR and GGR). We show an up-regulation in gene expression of CSA and HR23A, which are involved in TCR and GGR, respectively. Importantly, we show that this occurs through a p53 independent mechanism and that it is coordinated by the stress-responsive transcription factor USF-1. Furthermore, using a mouse model we show that the loss of USF-1 compromises DNA repair, which suggests that USF-1 plays an important role in maintaining genomic stability.

by Michael P. Heaton, Michael L. Clawson, Carol G. Chitko-Mckown, Kreg A. Leymaster, Timothy P. L. Smith, Gregory P. Harhay, Stephen N. White, Lynn M. Herrmann-Hoesing, Michelle R. Mousel, Gregory S. Lewis, Theodore S. Kalbfleisch, James E. Keen, William W. Laegreid

Visna/Maedi, or ovine progressive pneumonia (OPP) as it is known in the United States, is an incurable slow-acting disease of sheep caused by persistent lentivirus infection. This disease affects multiple tissues, including those of the respiratory and central nervous systems. Our aim was to identify ovine genetic risk factors for lentivirus infection. Sixty-nine matched pairs of infected cases and uninfected controls were identified among 736 naturally exposed sheep older than five years of age. These pairs were used in a genome-wide association study with 50,614 markers. A single SNP was identified in the ovine transmembrane protein (TMEM154) that exceeded genome-wide significance (unadjusted p-value 3×10−9). Sanger sequencing of the ovine TMEM154 coding region identified six missense and two frameshift deletion mutations in the predicted signal peptide and extracellular domain. Two TMEM154 haplotypes encoding glutamate (E) at position 35 were associated with infection while a third haplotype with lysine (K) at position 35 was not. Haplotypes encoding full-length E35 isoforms were analyzed together as genetic risk factors in a multi-breed, matched case-control design, with 61 pairs of 4-year-old ewes. The odds of infection for ewes with one copy of a full-length TMEM154 E35 allele were 28 times greater than the odds for those without (p-value<0.0001, 95% CI 5–1,100). In a combined analysis of nine cohorts with 2,705 sheep from Nebraska, Idaho, and Iowa, the relative risk of infection was 2.85 times greater for sheep with a full-length TMEM154 E35 allele (p-value<0.0001, 95% CI 2.36–3.43). Although rare, some sheep were homozygous for TMEM154 deletion mutations and remained uninfected despite a lifetime of significant exposure. Together, these findings indicate that TMEM154 may play a central role in ovine lentivirus infection and removing sheep with the most susceptible genotypes may help eradicate OPP and protect flocks from reinfection.

by Jérôme Albre, Marjorie A. Liénard, Tamara M. Sirey, Silvia Schmidt, Leah K. Tooman, Colm Carraher, David R. Greenwood, Christer Löfstedt, Richard D. Newcomb

Chemical signals are prevalent in sexual communication systems. Mate recognition has been extensively studied within the Lepidoptera, where the production and recognition of species-specific sex pheromone signals are typically the defining character. While the specific blend of compounds that makes up the sex pheromones of many species has been characterized, the molecular mechanisms underpinning the evolution of pheromone-based mate recognition systems remain largely unknown. We have focused on two sets of sibling species within the leafroller moth genera Ctenopseustis and Planotortrix that have rapidly evolved the use of distinct sex pheromone blends. The compounds within these blends differ almost exclusively in the relative position of double bonds that are introduced by desaturase enzymes. Of the six desaturase orthologs isolated from all four species, functional analyses in yeast and gene expression in pheromone glands implicate three in pheromone biosynthesis, two Δ9-desaturases, and a Δ10-desaturase, while the remaining three desaturases include a Δ6-desaturase, a terminal desaturase, and a non-functional desaturase. Comparative quantitative real-time PCR reveals that the Δ10-desaturase is differentially expressed in the pheromone glands of the two sets of sibling species, consistent with differences in the pheromone blend in both species pairs. In the pheromone glands of species that utilize (Z)-8-tetradecenyl acetate as sex pheromone component (Ctenopseustis obliquana and Planotortrix octo), the expression levels of the Δ10-desaturase are significantly higher than in the pheromone glands of their respective sibling species (C. herana and P. excessana). Our results demonstrate that interspecific sex pheromone differences are associated with differential regulation of the same desaturase gene in two genera of moths. We suggest that differential gene regulation among members of a multigene family may be an important mechanism of molecular innovation in sex pheromone evolution and speciation.

by Silvia Naitza, Eleonora Porcu, Maristella Steri, Dennis D. Taub, Antonella Mulas, Xiang Xiao, James Strait, Mariano Dei, Sandra Lai, Fabio Busonero, Andrea Maschio, Gianluca Usala, Magdalena Zoledziewska, Carlo Sidore, Ilenia Zara, Maristella Pitzalis, Alessia Loi, Francesca Virdis, Roberta Piras, Francesca Deidda, Michael B. Whalen, Laura Crisponi, Antonio Concas, Carlo Podda, Sergio Uzzau, Paul Scheet, Dan L. Longo, Edward Lakatta, Gonçalo R. Abecasis, Antonio Cao, David Schlessinger, Manuela Uda, Serena Sanna, Francesco Cucca

Identifying the genes that influence levels of pro-inflammatory molecules can help to elucidate the mechanisms underlying this process. We first conducted a two-stage genome-wide association scan (GWAS) for the key inflammatory biomarkers Interleukin-6 (IL-6), the general measure of inflammation erythrocyte sedimentation rate (ESR), monocyte chemotactic protein-1 (MCP-1), and high-sensitivity C-reactive protein (hsCRP) in a large cohort of individuals from the founder population of Sardinia. By analysing 731,213 autosomal or X chromosome SNPs and an additional ∼1.9 million imputed variants in 4,694 individuals, we identified several SNPs associated with the selected quantitative trait loci (QTLs) and replicated all the top signals in an independent sample of 1,392 individuals from the same population. Next, to increase power to detect and resolve associations, we further genotyped the whole cohort (6,145 individuals) for 293,875 variants included on the ImmunoChip and MetaboChip custom arrays. Overall, our combined approach led to the identification of 9 genome-wide significant novel independent signals—5 of which were identified only with the custom arrays—and provided confirmatory evidence for an additional 7. Novel signals include: for IL-6, in the ABO gene (rs657152, p = 2.13×10−29); for ESR, at the HBB (rs4910472, p = 2.31×10−11) and UCN119B/SPPL3 (rs11829037, p = 8.91×10−10) loci; for MCP-1, near its receptor CCR2 (rs17141006, p = 7.53×10−13) and in CADM3 (rs3026968, p = 7.63×10−13); for hsCRP, within the CRP gene (rs3093077, p = 5.73×10−21), near DARC (rs3845624, p = 1.43×10−10), UNC119B/SPPL3 (rs11829037, p = 1.50×10−14), and ICOSLG/AIRE (rs113459440, p = 1.54×10−08) loci. Confirmatory evidence was found for IL-6 in the IL-6R gene (rs4129267); for ESR at CR1 (rs12567990) and TMEM57 (rs10903129); for MCP-1 at DARC (rs12075); and for hsCRP at CRP (rs1205), HNF1A (rs225918), and APOC-I (rs4420638). Our results improve the current knowledge of genetic variants underlying inflammation and provide novel clues for the understanding of the molecular mechanisms regulating this complex process.

by Kai Yu, Sholom Wacholder, William Wheeler, Zhaoming Wang, Neil Caporaso, Maria Teresa Landi, Faming Liang

An important follow-up step after genetic markers are found to be associated with a disease outcome is a more detailed analysis investigating how the implicated gene or chromosomal region and an established environment risk factor interact to influence the disease risk. The standard approach to this study of gene–environment interaction considers one genetic marker at a time and therefore could misrepresent and underestimate the genetic contribution to the joint effect when one or more functional loci, some of which might not be genotyped, exist in the region and interact with the environment risk factor in a complex way. We develop a more global approach based on a Bayesian model that uses a latent genetic profile variable to capture all of the genetic variation in the entire targeted region and allows the environment effect to vary across different genetic profile categories. We also propose a resampling-based test derived from the developed Bayesian model for the detection of gene–environment interaction. Using data collected in the Environment and Genetics in Lung Cancer Etiology (EAGLE) study, we apply the Bayesian model to evaluate the joint effect of smoking intensity and genetic variants in the 15q25.1 region, which contains a cluster of nicotinic acetylcholine receptor genes and has been shown to be associated with both lung cancer and smoking behavior. We find evidence for gene–environment interaction (P-value = 0.016), with the smoking effect appearing to be stronger in subjects with a genetic profile associated with a higher lung cancer risk; the conventional test of gene–environment interaction based on the single-marker approach is far from significant.

by Yukinori Okada, Kenichi Shimane, Yuta Kochi, Tomoko Tahira, Akari Suzuki, Koichiro Higasa, Atsushi Takahashi, Tetsuya Horita, Tatsuya Atsumi, Tomonori Ishii, Akiko Okamoto, Keishi Fujio, Michito Hirakata, Hirofumi Amano, Yuya Kondo, Satoshi Ito, Kazuki Takada, Akio Mimori, Kazuyoshi Saito, Makoto Kamachi, Yasushi Kawaguchi, Katsunori Ikari, Osman Wael Mohammed, Koichi Matsuda, Chikashi Terao, Koichiro Ohmura, Keiko Myouzen, Naoya Hosono, Tatsuhiko Tsunoda, Norihiro Nishimoto, Tsuneyo Mimori, Fumihiko Matsuda, Yoshiya Tanaka, Takayuki Sumida, Hisashi Yamanaka, Yoshinari Takasaki, Takao Koike, Takahiko Horiuchi, Kenshi Hayashi, Michiaki Kubo, Naoyuki Kamatani, Ryo Yamada, Yusuke Nakamura, Kazuhiko Yamamoto

Systemic lupus erythematosus (SLE) is an autoimmune disease that causes multiple organ damage. Although recent genome-wide association studies (GWAS) have contributed to discovery of SLE susceptibility genes, few studies has been performed in Asian populations. Here, we report a GWAS for SLE examining 891 SLE cases and 3,384 controls and multi-stage replication studies examining 1,387 SLE cases and 28,564 controls in Japanese subjects. Considering that expression quantitative trait loci (eQTLs) have been implicated in genetic risks for autoimmune diseases, we integrated an eQTL study into the results of the GWAS. We observed enrichments of cis-eQTL positive loci among the known SLE susceptibility loci (30.8%) compared to the genome-wide SNPs (6.9%). In addition, we identified a novel association of a variant in the AF4/FMR2 family, member 1 (AFF1) gene at 4q21 with SLE susceptibility (rs340630; P = 8.3×10−9, odds ratio = 1.21). The risk A allele of rs340630 demonstrated a cis-eQTL effect on the AFF1 transcript with enhanced expression levels (P<0.05). As AFF1 transcripts were prominently expressed in CD4+ and CD19+ peripheral blood lymphocytes, up-regulation of AFF1 may cause the abnormality in these lymphocytes, leading to disease onset.

by Daniel John Lawson, Garrett Hellenthal, Simon Myers, Daniel Falush

The advent of genome-wide dense variation data provides an opportunity to investigate ancestry in unprecedented detail, but presents new statistical challenges. We propose a novel inference framework that aims to efficiently capture information on population structure provided by patterns of haplotype similarity. Each individual in a sample is considered in turn as a recipient, whose chromosomes are reconstructed using chunks of DNA donated by the other individuals. Results of this “chromosome painting” can be summarized as a “coancestry matrix,” which directly reveals key information about ancestral relationships among individuals. If markers are viewed as independent, we show that this matrix almost completely captures the information used by both standard Principal Components Analysis (PCA) and model-based approaches such as STRUCTURE in a unified manner. Furthermore, when markers are in linkage disequilibrium, the matrix combines information across successive markers to increase the ability to discern fine-scale population structure using PCA. In parallel, we have developed an efficient model-based approach to identify discrete populations using this matrix, which offers advantages over PCA in terms of interpretability and over existing clustering algorithms in terms of speed, number of separable populations, and sensitivity to subtle population structure. We analyse Human Genome Diversity Panel data for 938 individuals and 641,000 markers, and we identify 226 populations reflecting differences on continental, regional, local, and family scales. We present multiple lines of evidence that, while many methods capture similar information among strongly differentiated groups, more subtle population structure in human populations is consistently present at a much finer level than currently available geographic labels and is only captured by the haplotype-based approach. The software used for this article, ChromoPainter and fineSTRUCTURE, is available from http://www.paintmychromosomes.com/.