Background/Aims In primary biliary cirrhosis (PBC), pathogenesis is influenced by genetic factors that remain poorly elucidated up to now. We investigated the impact of sequence diversity in ...candidate genes involved in immunity ( CTLA-4 and TNFα ), in bile formation (10 hepatobiliary transporter genes) and in the adaptative response to cholestasis (three nuclear receptor genes) on the susceptibility and severity of PBC. Methods A total of 42 Ht SNPs were identified and compared in 258 PBC patients and two independent groups of 286 and 269 healthy controls. All participants were white continental individuals with French ancestry. Results Ht SNPs of CTLA-4 and TNFα genes were significantly associated with susceptibility to PBC. The progression rate of liver disease under ursodeoxycholic acid (UDCA) therapy was significantly linked to SNPs of TNFα and SLC4A2 /anion exchanger 2 (AE2) genes. A multivariate Cox regression analysis including clinical and biochemical parameters showed that SLC4A2 /AE2 variant was an independent prognostic factor. Conclusions These data point to a primary role of genes encoding regulators of the immune system in the susceptibility to PBC. They also demonstrate that allelic variations in TNFα and SLC4A2 /AE2 have a significant impact on the evolutive profile of PBC under UDCA therapy.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
We apply a known algorithm for computing exactly inequalities between Beta distributions to assess whether a given position in a genome is differentially methylated across samples. We discuss the ...advantages brought by the adoption of this solution with respect to two approximations (Fisher's test and Z score). The same formalism presented here can be applied in a similar way to variant calling.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Binding of the programmed death-1 (PD-1) receptor to its ligands (PD-L1/2) transduces inhibitory signals that promote exhaustion of activated T cells. Blockade of the PD-1 pathway is widely used for ...cancer treatment, yet the inhibitory signals transduced by PD-1 in T cells remain elusive.
Expression profiles of human CD8
T cells in resting, activated (CD3 + CD28) and PD-1-stimulated cells (CD3 + CD28 + PD-L1-Fc) conditions were evaluated by RNA-seq. Bioinformatic analyses were used to identify signaling pathways differentially regulated in PD-1-stimulated cells. Metabolic analyses were performed with SeaHorse technology, and mitochondrial ultrastructure was determined by transmission electron microscopy. PD-1-regulated mitochondrial genes were silenced using short-hairpin RNA in primary cells. Blue native gel electrophoresis was used to determine respiratory supercomplex assembly.
PD-1 engagement in human CD8
T cells triggers a specific, progressive genetic program different from that found in resting cells. Gene ontology identified metabolic processes, including glycolysis and oxidative phosphorylation (OXPHOS), as the main pathways targeted by PD-1. We observed severe functional and structural alterations in the mitochondria of PD-1-stimulated cells, including a reduction in the number and length of mitochondrial cristae. These cristae alterations were associated with reduced expression of CHCHD3 and CHCHD10, two proteins that form part of the mitochondrial contact site and cristae organizing system (MICOS). Although PD-1-stimulated cells showed severe cristae alterations, assembly of respiratory supercomplexes was unexpectedly greater in these cells than in activated T cells. CHCHD3 silencing in primary CD8
T cells recapitulated some effects induced by PD-1 stimulation, including reduced mitochondrial polarization and interferon-γ production following T cell activation with anti-CD3 and -CD28 activating antibodies.
Our results suggest that mitochondria are the main targets of PD-1 inhibitory activity. PD-1 reprograms CD8
T cell metabolism for efficient use of fatty acid oxidation; this mitochondrial phenotype might explain the long-lived phenotype of PD-1-engaged T cells.
Dairy cattle breeds have been subjected over the last fifty years to intense artificial selection towards improvement of milk production traits. In this study, we performed a whole genome scan for ...differentiation using 42,486 SNPs in the three major French dairy cattle breeds (Holstein, Normande and Montbéliarde) to identify the main physiological pathways and regions which were affected by this selection. After analyzing the population structure, we estimated F(ST) within and across the three breeds for each SNP under a pure drift model. We further considered two different strategies to evaluate the effect of selection at the genome level. First, smoothing F(ST) values over each chromosome with a local variable bandwidth kernel estimator allowed identifying 13 highly significant regions subjected to strong and/or recent positive selection. Some of them contained genes within which causal variants with strong effect on milk production traits (GHR) or coloration (MC1R) have already been reported. To go further in the interpretation of the observed signatures of selection we subsequently concentrated on the annotation of differentiated genes defined according to the F(ST) value of SNPs localized close or within them. To that end we performed a comprehensive network analysis which suggested a central role of somatotropic and gonadotropic axes in the response to selection. Altogether, these observations shed light on the antagonism, at the genome level, between milk production and reproduction traits in highly producing dairy cows.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Genome-wide association studies have identified several common genetic variants associated with breast cancer risk. It is likely, however, that a substantial proportion of such loci have not yet been ...discovered.
We compared 296,114 tagging single-nucleotide polymorphisms in 1694 breast cancer case subjects (92% with two primary cancers or at least two affected first-degree relatives) and 2365 control subjects, with validation in three independent series totaling 11,880 case subjects and 12,487 control subjects. Odds ratios (ORs) and associated 95% confidence intervals (CIs) in each stage and all stages combined were calculated using unconditional logistic regression. Heterogeneity was evaluated with Cochran Q and I(2) statistics. All statistical tests were two-sided.
We identified a novel risk locus for breast cancer at 9q31.2 (rs865686: OR = 0.89, 95% CI = 0.85 to 0.92, P = 1.75 × 10(-10)). This single-nucleotide polymorphism maps to a gene desert, the nearest genes being Kruppel-like factor 4 (KLF4, 636 kb centromeric), RAD23 homolog B (RAD23B, 794 kb centromeric), and actin-like 7A (ACTL7A, 736 kb telomeric). We also identified two variants (rs3734805 and rs9383938) mapping to 6q25.1 estrogen receptor 1 (ESR1), which were associated with breast cancer in subjects of northern European ancestry (rs3734805: OR = 1.19, 95% CI = 1.11 to 1.27, P = 1.35 × 10(-7); rs9383938: OR = 1.18, 95% CI = 1.11 to 1.26, P = 1.41 × 10(-7)). A variant mapping to 10q26.13, approximately 300 kb telomeric to the established risk locus within the second intron of FGFR2, was also associated with breast cancer risk, although not at genome-wide statistical significance (rs10510102: OR = 1.12, 95% CI = 1.07 to 1.17, P = 1.58 × 10(-6)).
These findings provide further evidence on the role of genetic variation in the etiology of breast cancer. Fine mapping will be needed to identify causal variants and to determine their functional effects.
To identify novel susceptibility loci for Crohn disease (CD), we undertook a genome-wide association study with more than 300,000 SNPs characterized in 547 patients and 928 controls. We found three ...chromosome regions that provided evidence of disease association with p-values between 10(-6) and 10(-9). Two of these (IL23R on Chromosome 1 and CARD15 on Chromosome 16) correspond to genes previously reported to be associated with CD. In addition, a 250-kb region of Chromosome 5p13.1 was found to contain multiple markers with strongly suggestive evidence of disease association (including four markers with p < 10(-7)). We replicated the results for 5p13.1 by studying 1,266 additional CD patients, 559 additional controls, and 428 trios. Significant evidence of association (p < 4 x 10(-4)) was found in case/control comparisons with the replication data, while associated alleles were over-transmitted to affected offspring (p < 0.05), thus confirming that the 5p13.1 locus contributes to CD susceptibility. The CD-associated 250-kb region was saturated with 111 SNP markers. Haplotype analysis supports a complex locus architecture with multiple variants contributing to disease susceptibility. The novel 5p13.1 CD locus is contained within a 1.25-Mb gene desert. We present evidence that disease-associated alleles correlate with quantitative expression levels of the prostaglandin receptor EP4, PTGER4, the gene that resides closest to the associated region. Our results identify a major new susceptibility locus for CD, and suggest that genetic variants associated with disease risk at this locus could modulate cis-acting regulatory elements of PTGER4.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Venous thromboembolism (VTE) is a complex disease that has a major genetic component of risk. To identify genetic factors that may modify the risk of VTE, we conducted a genome-wide association study ...by analyzing approximately 317 000 single nucleotide polymorphisms (SNPs) in 453 VTE cases and 1327 controls. Only 3 SNPs located in the FV and ABO blood group genes were found associated with VTE at a genome-wide significant level of 1.7 × 10−7. Detailed analysis of these SNPs in additional cohorts of more than 1700 cases and 1400 controls revealed that the association observed at the FV locus was the result of the increased risk mediated by the FV Leiden mutation, whereas O and A2 blood groups were found to be at lower risk for VTE. Apart from the FV and ABO loci, no other locus was found strongly associated with VTE. However, using this large cohort of subjects, we were able to replicate the mild effects of 2 nonsynonymous SNPs, rs1613662 in GP6 and rs13146272 in CYP4V2, recently suspected to be associated with VTE.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Emery-Dreifuss muscular dystrophy (EDMD) is a rare disorder characterized by early joint contractures, muscular dystrophy, and cardiac involvement with conduction defects and arrhythmias. So far, ...only 35% of EDMD cases are genetically elucidated and associated with
EMD or
LMNA gene mutations, suggesting the existence of additional major genes. By whole-genome scan, we identified linkage to the Xq26.3 locus containing the
FHL1 gene in three informative families belonging to our
EMD- and
LMNA-negative cohort. Analysis of the
FHL1 gene identified seven mutations, in the distal exons of
FHL1 in these families, three additional families, and one isolated case, which differently affect the three FHL1 protein isoforms: two missense mutations affecting highly conserved cysteines, one abolishing the termination codon, and four out-of-frame insertions or deletions. The predominant phenotype was characterized by myopathy with scapulo-peroneal and/or axial distribution, as well as joint contractures, and associated with a peculiar cardiac disease characterized by conduction defects, arrhythmias, and hypertrophic cardiomyopathy in all index cases of the seven families. Heterozygous female carriers were either asymptomatic or had cardiac disease and/or mild myopathy. Interestingly, four of the
FHL1-mutated male relatives had isolated cardiac disease, and an overt hypertrophic cardiomyopathy was present in two. Expression and functional studies demonstrated that the FHL1 proteins were severely reduced in all tested patients and that this was associated with a severe delay in myotube formation in the two patients for whom myoblasts were available. In conclusion,
FHL1 should be considered as a gene associated with the X-linked EDMD phenotype, as well as with hypertrophic cardiomyopathy.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
SNP calling by sequencing pooled samples Raineri, Emanuele; Ferretti, Luca; Esteve-Codina, Anna ...
BMC bioinformatics,
09/2012, Volume:
13, Issue:
1
Journal Article
Peer reviewed
Open access
Performing high throughput sequencing on samples pooled from different individuals is a strategy to characterize genetic variability at a small fraction of the cost required for individual ...sequencing. In certain circumstances some variability estimators have even lower variance than those obtained with individual sequencing. SNP calling and estimating the frequency of the minor allele from pooled samples, though, is a subtle exercise for at least three reasons. First, sequencing errors may have a much larger relevance than in individual SNP calling: while their impact in individual sequencing can be reduced by setting a restriction on a minimum number of reads per allele, this would have a strong and undesired effect in pools because it is unlikely that alleles at low frequency in the pool will be read many times. Second, the prior allele frequency for heterozygous sites in individuals is usually 0.5 (assuming one is not analyzing sequences coming from, e.g. cancer tissues), but this is not true in pools: in fact, under the standard neutral model, singletons (i.e. alleles of minimum frequency) are the most common class of variants because P(f) ∝ 1/f and they occur more often as the sample size increases. Third, an allele appearing only once in the reads from a pool does not necessarily correspond to a singleton in the set of individuals making up the pool, and vice versa, there can be more than one read - or, more likely, none - from a true singleton.
To improve upon existing theory and software packages, we have developed a Bayesian approach for minor allele frequency (MAF) computation and SNP calling in pools (and implemented it in a program called snape): the approach takes into account sequencing errors and allows users to choose different priors. We also set up a pipeline which can simulate the coalescence process giving rise to the SNPs, the pooling procedure and the sequencing. We used it to compare the performance of snape to that of other packages.
We present a software which helps in calling SNPs in pooled samples: it has good power while retaining a low false discovery rate (FDR). The method also provides the posterior probability that a SNP is segregating and the full posterior distribution of f for every SNP. In order to test the behaviour of our software, we generated (through simulated coalescence) artificial genomes and computed the effect of a pooled sequencing protocol, followed by SNP calling. In this setting, snape has better power and False Discovery Rate (FDR) than the comparable packages samtools, PoPoolation, Varscan : for N = 50 chromosomes, snape has power ≈ 35%and FDR ≈ 2.5%. snape is available at http://code.google.com/p/snape-pooled/ (source code and precompiled binaries).
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK