Cis and Trans Effects of Human Genomic Variants on Gene Expression

This work was funded by the Louis-Jeantet Foundation (http://www.jeantet.ch/), the European Research Council (Grant ID: 260927 http://erc.europa.eu/), the Swiss National Foundation (Grant ID: 130342 http://www.snf.ch), NCCR Frontiers In Genetics (http://www.frontiers-in-genetics.org), the UK Medical Research Council (http://www.mrc.ac.uk) and the Wellcome Trust (Grant ID: 092731).

Characterization and Comparison of 2 Distinct Epidemic Community-Associated Methicillin-Resistant Staphylococcus aureus Clones of ST59 Lineage.

Sequence type (ST) 59 is an epidemic lineage of community-associated (CA) methicillin-resistant Staphylococcus aureus (MRSA) isolates. Taiwanese CA-MRSA isolates belong to ST59 and can be grouped into 2 distinct clones, a virulent Taiwan clone and a commensal Asian-Pacific clone. The Taiwan clone carries the Panton-Valentine leukocidin (PVL) genes and the staphylococcal chromosomal cassette mec (SCCmec) VT, and is frequently isolated from patients with severe disease. The Asian-Pacific clone is PVL-negative, carries SCCmec IV, and a frequent colonizer of healthy children. Isolates of both clones were characterized by their ability to adhere to respiratory A549 cells, cytotoxicity to human neutrophils, and nasal colonization of a murine and murine sepsis models. Genome variation was determined by polymerase chain reaction of selected virulence factors and by multi-strain whole genome microarray. Additionally, the expression of selected factors was compared between the 2 clones. The Taiwan clone showed a much higher cytotoxicity to the human neutrophils and caused more severe septic infections with a high mortality rate in the murine model. The clones were indistinguishable in their adhesion to A549 cells and persistence of murine nasal colonization. The microarray data revealed that the Taiwan clone had lost the ø3-prophage that integrates into the β-hemolysin gene and includes staphylokinase- and enterotoxin P-encoding genes, but had retained the genes for human immune evasion, scn and chps. Production of the virulence factors did not differ significantly in the 2 clonal groups, although more α-toxin was expressed in Taiwan clone isolates from pneumonia patients. In conclusion, the Taiwan CA-MRSA clone was distinguished by enhanced virulence in both humans and an animal infection model. The evolutionary acquisition of PVL, the higher expression of α-toxin, and possibly the loss of a large portion of the β-hemolysin-converting prophage likely contribute to its higher pathogenic potential than the Asian-Pacific clone

Patterns of Cis Regulatory Variation in Diverse Human Populations

The genetic basis of gene expression variation has long been studied with the aim to understand the landscape of regulatory variants, but also more recently to assist in the interpretation and elucidation of disease signals. To date, many studies have looked in specific tissues and population-based samples, but there has been limited assessment of the degree of inter-population variability in regulatory variation. We analyzed genome-wide gene expression in lymphoblastoid cell lines from a total of 726 individuals from 8 global populations from the HapMap3 project and correlated gene expression levels with HapMap3 SNPs located in cis to the genes. We describe the influence of ancestry on gene expression levels within and between these diverse human populations and uncover a non-negligible impact on global patterns of gene expression. We further dissect the specific functional pathways differentiated between populations. We also identify 5,691 expression quantitative trait loci (eQTLs) after controlling for both non-genetic factors and population admixture and observe that half of the cis-eQTLs are replicated in one or more of the populations. We highlight patterns of eQTL-sharing between populations, which are partially determined by population genetic relatedness, and discover significant sharing of eQTL effects between Asians, European-admixed, and African subpopulations. Specifically, we observe that both the effect size and the direction of effect for eQTLs are highly conserved across populations. We observe an increasing proximity of eQTLs toward the transcription start site as sharing of eQTLs among populations increases, highlighting that variants close to TSS have stronger effects and therefore are more likely to be detected across a wider panel of populations. Together these results offer a unique picture and resource of the degree of differentiation among human populations in functional regulatory variation and provide an estimate for the transferability of complex trait variants across populations

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

The Architecture of Gene Regulatory Variation across Multiple Human Tissues: The MuTHER Study

While there have been studies exploring regulatory variation in one or more tissues, the complexity of tissue-specificity in multiple primary tissues is not yet well understood. We explore in depth the role of cis-regulatory variation in three human tissues: lymphoblastoid cell lines (LCL), skin, and fat. The samples (156 LCL, 160 skin, 166 fat) were derived simultaneously from a subset of well-phenotyped healthy female twins of the MuTHER resource. We discover an abundance of cis-eQTLs in each tissue similar to previous estimates (858 or 4.7% of genes). In addition, we apply factor analysis (FA) to remove effects of latent variables, thus more than doubling the number of our discoveries (1,822 eQTL genes). The unique study design (Matched Co-Twin Analysis—MCTA) permits immediate replication of eQTLs using co-twins (93%–98%) and validation of the considerable gain in eQTL discovery after FA correction. We highlight the challenges of comparing eQTLs between tissues. After verifying previous significance threshold-based estimates of tissue-specificity, we show their limitations given their dependency on statistical power. We propose that continuous estimates of the proportion of tissue-shared signals and direct comparison of the magnitude of effect on the fold change in expression are essential properties that jointly provide a biologically realistic view of tissue-specificity. Under this framework we demonstrate that 30% of eQTLs are shared among the three tissues studied, while another 29% appear exclusively tissue-specific. However, even among the shared eQTLs, a substantial proportion (10%–20%) have significant differences in the magnitude of fold change between genotypic classes across tissues. Our results underline the need to account for the complexity of eQTL tissue-specificity in an effort to assess consequences of such variants for complex traits

Alpha-Toxin Induces Programmed Cell Death of Human T cells, B cells, and Monocytes during USA300 Infection

This investigation examines the influence of alpha-toxin (Hla) during USA300 infection of human leukocytes. Survival of an USA300 isogenic deletion mutant of hla (USA300Δhla) in human blood was comparable to the parental wild-type strain and polymorphonuclear leukocyte (PMN) plasma membrane permeability caused by USA300 did not require Hla. Flow cytometry analysis of peripheral blood mononuclear cells (PBMCs) following infection by USA300, USA300Δhla, and USA300Δhla transformed with a plasmid over-expressing Hla (USA300Δhla Comp) demonstrated this toxin plays a significant role inducing plasma membrane permeability of CD14+, CD3+, and CD19+ PBMCs. Rapid plasma membrane permeability independent of Hla was observed for PMNs, CD14+ and CD19+ PBMCs following intoxication with USA300 supernatant while the majority of CD3+ PBMC plasma membrane permeability induced by USA300 required Hla. Addition of recombinant Hla to USA300Δhla supernatant rescued CD3+ and CD19+ PBMC plasma membrane permeability generated by USA300 supernatant. An observed delay in plasma membrane permeability caused by Hla in conjunction with Annexin V binding and ApoBrdU Tunel assays examining PBMCs intoxicated with recombinant Hla or infected with USA300, USA300Δhla, USA300Δhla Comp, and USA300ΔsaeR/S suggest Hla induces programmed cell death of monocytes, B cells, and T cells that results in plasma membrane permeability. Together these findings underscore the importance of Hla during S. aureus infection of human tissue and specifically demonstrate Hla activity during USA300 infection triggers programmed cell death of human monocytes, T cells and B cells that leads to plasma membrane permeability

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition

Recurrent Zoonotic Transmission of Nipah Virus into Humans, Bangladesh, 2001–2007

More than half of identified cases result from person-to-person transmission