Rat B-13 progenitor cells are readily converted into functional hepatocyte-like B-13/H cells capable of phase I cytochrome P450-dependent activation of pro-carcinogens and induction of DNA damage. The aim of the present study was to investigate whether the cells are also capable of Phase II sulphotransferase (SULT)-dependent activation of a pro-carcinogen to an ultimate carcinogen. To this end we therefore examined the bioactivation of the model hepatic (hepato- and cholangio-) carcinogen estragole and its proximate SULT1A1-activated genotoxic metabolite 1'-hydroxyestragole. Exposing B-13 or B-13/H cells to estragole (at concentrations up to 1mM) resulted in the production of low levels of 1'-hydroxyestragole, but did not result in detectable DNA damage. Exposing B-13/H cells - but not B-13 cells - to 1'-hydroxyestragole resulted in a dose-dependent increase in DNA damage in comet assays, confirmed by detection of N(2)-(trans-isoestragol-3'-yl)-2'-deoxyguanosine adducts. Genotoxicity was inhibited by general SULT inhibitors, supporting a role for SULTS in the activation of 1-hydroxyestragole in B-13/H cells. However, B-13 and B-13/H cells did not express biologically significant levels of SULT1A1 as determined by qRT-PCR, Western blotting and its associated 7-hydroxycoumarin sulphation activity. B-13 and B-13/H cells expressed - relative to intact rat liver - high levels of SULT2B1 (primarily the b isoform) and SULT4A1 mRNAs and proteins. B-13 and B-13/H cells also expressed the 3'-phosphoadenosine 5'-phosphosulphate synthase 1 required for the generation of activated sulphate cofactor 3'-phosphoadenosine 5'-phosphosulphate. However, only B-13/H cells expressed functional SULT activities towards SULT2B1 substrates DHEA, pregnenolone and 4 methylumbelliferone. Since liver progenitor cells are bi-potential and also form cholangiocytes, we therefore hypothesised that B-13 cells express a cholangiocyte-like SULT profile. To test this hypothesis, the expression of SULTs was examined in liver by RT-PCR and immunohistochemistry. SULT2B1 - but not SULT1A1 - was determined to be expressed in both rat and human cholangiocytes. Since 1'-hydroxyestragole exposure readily produced DNA injury in B-13/H cells, these data suggest that cholangiocarcinomas generated in rats fed estragole may be dependent, in part, on SULT2B1 activation of the 1'-hydroxyestragole metabolite.

Some chemicals are ligands to efflux transporters which may result in high concentrations in milk. Limited knowledge is available on the influence of maternal exposure to chemicals on the expression and function of transporters in the lactating mammary gland. We determined gene expression of ABC and SLC transporters in murine mammary tissue of different gestation and lactation stages, in murine mammary cells (HC11) featuring resting and secreting phenotypes and in bovine mammary tissue and cells (BME-UV). Effects on transporter expression and function of the imidazole fungicide prochloraz, previously reported to influence BCRP in mammary cells, was investigated on transporter expression and function in the two cell lines. Transporters studied were BCRP, MDR1, MRP1, OATP1A5/OATP1A2, OCTN1 and OCT1. Gene expressions of BCRP and OCT1 in murine mammary glands were increased during gestation and lactation, whereas MDR1, MRP1, OATP1A5 and OCTN1 were decreased, compared to expressions in virgins. All transporters measured in mammary glands of mice were detected in bovine mammary tissue and in HC11 cells, while only MDR1 and MRP1 were detected in BME-UV cells. Prochloraz treatment induced MDR1 gene and protein expression in both differentiated HC11 and BME-UV cells and increased protein function in HC11 cells, resulting in decreased accumulation of the MDR1 substrate digoxin. In conclusion, our results demonstrate that murine (HC11) and bovine (BME-UV) mammary epithelial cells can be applied to characterize expression and function of transporters as well as effects of contaminants on the mammary transporters. An altered expression, induced by a drug or toxic chemical, on any of the transporters expressed in the mammary epithelial cells during lactation may modulate the well-balanced composition of nutrients and/or secretion of contaminants in milk with potential adverse effects on breast-fed infants and dairy consumers.

Post exertional muscle fatigue is a key feature in Chronic Fatigue Syndrome (CFS). Abnormalities of skeletal muscle function have been identified in some but not all patients with CFS. To try to limit potential confounders that might contribute to this clinical heterogeneity, we developed a novel in vitro system that allows comparison of AMP kinase (AMPK) activation and metabolic responses to exercise in cultured skeletal muscle cells from CFS patients and control subjects.

Thymic stromal lymphopoietin (TSLP), an IL7-like cytokine produced by bronchial epithelial cells is upregulated in asthma and induces dendritic cell maturation supporting a Th2 response. Environmental pollutants, including tobacco smoke and diesel exhaust particles upregulate TSLP suggesting that TSLP may be an interface between environmental pollution and immune responses in asthma. Since asthma is prevalent in urban communities, variants in the TSLP gene may be important in asthma susceptibility in these populations.

Primary biliary cirrhosis (PBC) is a classical autoimmune liver disease for which effective immunomodulatory therapy is lacking. Here we perform meta-analyses of discovery data sets from genome-wide association studies of European subjects (n=2,764 cases and 10,475 controls) followed by validation genotyping in an independent cohort (n=3,716 cases and 4,261 controls). We discover and validate six previously unknown risk loci for PBC (Pcombined<5 × 10(-8)) and used pathway analysis to identify JAK-STAT/IL12/IL27 signalling and cytokine-cytokine pathways, for which relevant therapies exist.

Patients with primary sclerosing cholangitis (PSC) frequently highlight the impact of fatigue on their life quality. The study aims were to evaluate fatigue and its associations in PSC and investigate whether overt autonomic dysfunction contributes to the expression of fatigue.

In this study, we demonstrate the use of natural language processing methods to extract, from nanomedicine literature, numeric values of biomedical property terms of poly(amidoamine) dendrimers. We have developed a method for extracting these values for properties taken from the NanoParticle Ontology, using the General Architecture for Text Engineering and a Nearly-New Information Extraction System. We also created a method for associating the identified numeric values with their corresponding dendrimer properties, called NanoSifter. We demonstrate that our system can correctly extract numeric values of dendrimer properties reported in the cancer treatment literature with high recall, precision, and f-measure. The micro-averaged recall was 0.99, precision was 0.84, and f-measure was 0.91. Similarly, the macro-averaged recall was 0.99, precision was 0.87, and f-measure was 0.92. To our knowledge, these results are the first application of text mining to extract and associate dendrimer property terms and their corresponding numeric values.

Diabetic kidney disease (DKD) is the most common etiology of chronic kidney disease (CKD) in the industrialized world and accounts for much of the excess mortality in patients with diabetes mellitus. Approximately 45% of U.S. patients with incident end-stage kidney disease (ESKD) have DKD. Independent of glycemic control, DKD aggregates in families and has higher incidence rates in African, Mexican, and American Indian ancestral groups relative to European populations. The Family Investigation of Nephropathy and Diabetes (FIND) performed a genome-wide association study (GWAS) contrasting 6,197 unrelated individuals with advanced DKD with healthy and diabetic individuals lacking nephropathy of European American, African American, Mexican American, or American Indian ancestry. A large-scale replication and trans-ethnic meta-analysis included 7,539 additional European American, African American and American Indian DKD cases and non-nephropathy controls. Within ethnic group meta-analysis of discovery GWAS and replication set results identified genome-wide significant evidence for association between DKD and rs12523822 on chromosome 6q25.2 in American Indians (P = 5.74x10-9). The strongest signal of association in the trans-ethnic meta-analysis was with a SNP in strong linkage disequilibrium with rs12523822 (rs955333; P = 1.31x10-8), with directionally consistent results across ethnic groups. These 6q25.2 SNPs are located between the SCAF8 and CNKSR3 genes, a region with DKD relevant changes in gene expression and an eQTL with IPCEF1, a gene co-translated with CNKSR3. Several other SNPs demonstrated suggestive evidence of association with DKD, within and across populations. These data identify a novel DKD susceptibility locus with consistent directions of effect across diverse ancestral groups and provide insight into the genetic architecture of DKD.

The minor allele of the R620W missense single-nucleotide polymorphism (SNP) (rs2476601) in the hematopoietic-specific protein tyrosine phosphatase gene, PTPN22, has been associated with multiple autoimmune diseases, including rheumatoid arthritis (RA). These genetic data, combined with biochemical evidence that this SNP affects PTPN22 function, suggest that this phosphatase is a key regulator of autoimmunity. To determine whether other genetic variants in PTPN22 contribute to the development of RA, we sequenced the coding regions of this gene in 48 white North American patients with RA and identified 15 previously unreported SNPs, including 2 coding SNPs in the catalytic domain. We then genotyped 37 SNPs in or near PTPN22 in 475 patients with RA and 475 individually matched controls (sample set 1) and selected a subset of markers for replication in an additional 661 patients with RA and 1,322 individually matched controls (sample set 2). Analyses of these results predict 10 common (frequency >1%) PTPN22 haplotypes in white North Americans. The sole haplotype found to carry the previously identified W620 risk allele was strongly associated with disease in both sample sets, whereas another haplotype, identical at all other SNPs but carrying the R620 allele, showed no association. R620W, however, does not fully explain the association between PTPN22 and RA, since significant differences between cases and controls persisted in both sample sets after the haplotype data were stratified by R620W. Additional analyses identified two SNPs on a single common haplotype that are associated with RA independent of R620W, suggesting that R620W and at least one additional variant in the PTPN22 gene region influence RA susceptibility.

Many lines of evidence implicate mitochondria in phenotypic variation: (a) rare mutations in mitochondrial proteins cause metabolic, neurological, and muscular disorders; (b) alterations in oxidative phosphorylation are characteristic of type 2 diabetes, Parkinson disease, Huntington disease, and other diseases; and (c) common missense variants in the mitochondrial genome (mtDNA) have been implicated as having been subject to natural selection for adaptation to cold climates and contributing to "energy deficiency" diseases today. To test the hypothesis that common mtDNA variation influences human physiology and disease, we identified all 144 variants with frequency >1% in Europeans from >900 publicly available European mtDNA sequences and selected 64 tagging single-nucleotide polymorphisms that efficiently capture all common variation (except the hypervariable D-loop). Next, we evaluated the complete set of common mtDNA variants for association with type 2 diabetes in a sample of 3,304 diabetics and 3,304 matched nondiabetic individuals. Association of mtDNA variants with other metabolic traits (body mass index, measures of insulin secretion and action, blood pressure, and cholesterol) was also tested in subsets of this sample. We did not find a significant association of common mtDNA variants with these metabolic phenotypes. Moreover, we failed to identify any physiological effect of alleles that were previously proposed to have been adaptive for energy metabolism in human evolution. More generally, this comprehensive association-testing framework can readily be applied to other diseases for which mitochondrial dysfunction has been implicated.

To discover new rheumatoid arthritis (RA) risk loci, we systematically examined 370 SNPs from 179 independent loci with P < 0.001 in a published meta-analysis of RA genome-wide association studies (GWAS) of 3,393 cases and 12,462 controls. We used Gene Relationships Across Implicated Loci (GRAIL), a computational method that applies statistical text mining to PubMed abstracts, to score these 179 loci for functional relationships to genes in 16 established RA disease loci. We identified 22 loci with a significant degree of functional connectivity. We genotyped 22 representative SNPs in an independent set of 7,957 cases and 11,958 matched controls. Three were convincingly validated: CD2-CD58 (rs11586238, P = 1 x 10(-6) replication, P = 1 x 10(-9) overall), CD28 (rs1980422, P = 5 x 10(-6) replication, P = 1 x 10(-9) overall) and PRDM1 (rs548234, P = 1 x 10(-5) replication, P = 2 x 10(-8) overall). An additional four were replicated (P < 0.0023): TAGAP (rs394581, P = 0.0002 replication, P = 4 x 10(-7) overall), PTPRC (rs10919563, P = 0.0003 replication, P = 7 x 10(-7) overall), TRAF6-RAG1 (rs540386, P = 0.0008 replication, P = 4 x 10(-6) overall) and FCGR2A (rs12746613, P = 0.0022 replication, P = 2 x 10(-5) overall). Many of these loci are also associated to other immunologic diseases.

Accounting for population genetic substructure is important in reducing type 1 errors in genetic studies of complex disease. As efforts to understand complex genetic disease are expanded to different continental populations the understanding of genetic substructure within these continents will be useful in design and execution of association tests. In this study, population differentiation (Fst) and Principal Components Analyses (PCA) are examined using >200 K genotypes from multiple populations of East Asian ancestry. The population groups included those from the Human Genome Diversity Panel [Cambodian, Yi, Daur, Mongolian, Lahu, Dai, Hezhen, Miaozu, Naxi, Oroqen, She, Tu, Tujia, Naxi, Xibo, and Yakut], HapMap [ Han Chinese (CHB) and Japanese (JPT)], and East Asian or East Asian American subjects of Vietnamese, Korean, Filipino and Chinese ancestry. Paired Fst (Wei and Cockerham) showed close relationships between CHB and several large East Asian population groups (CHB/Korean, 0.0019; CHB/JPT, 00651; CHB/Vietnamese, 0.0065) with larger separation with Filipino (CHB/Filipino, 0.014). Low levels of differentiation were also observed between Dai and Vietnamese (0.0045) and between Vietnamese and Cambodian (0.0062). Similarly, small Fst's were observed among different presumed Han Chinese populations originating in different regions of mainland of China and Taiwan (Fst's <0.0025 with CHB). For PCA, the first two PC's showed a pattern of relationships that closely followed the geographic distribution of the different East Asian populations. PCA showed substructure both between different East Asian groups and within the Han Chinese population. These studies have also identified a subset of East Asian substructure ancestry informative markers (EASTASAIMS) that may be useful for future complex genetic disease association studies in reducing type 1 errors and in identifying homogeneous groups that may increase the power of such studies.

Case-control genetic studies of complex human diseases can be confounded by population stratification. This issue can be addressed using panels of ancestry informative markers (AIMs) that can provide substantial population substructure information. Previously, we described a panel of 128 SNP AIMs that were designed as a tool for ascertaining the origins of subjects from Europe, Sub-Saharan Africa, Americas, and East Asia.

For admixture mapping studies in Mexican Americans (MAM), we define a genomewide single-nucleotide-polymorphism (SNP) panel that can distinguish between chromosomal segments of Amerindian (AMI) or European (EUR) ancestry. These studies used genotypes for >400,000 SNPs, defined in EUR and both Pima and Mayan AMI, to define a set of ancestry-informative markers (AIMs). The use of two AMI populations was necessary to remove a subset of SNPs that distinguished genotypes of only one AMI subgroup from EUR genotypes. The AIMs set contained 8,144 SNPs separated by a minimum of 50 kb with only three intermarker intervals >1 Mb and had EUR/AMI FST values >0.30 (mean FST = 0.48) and Mayan/Pima FST values <0.05 (mean FST < 0.01). Analysis of a subset of these SNP AIMs suggested that this panel may also distinguish ancestry between EUR and other disparate AMI groups, including Quechuan from South America. We show, using realistic simulation parameters that are based on our analyses of MAM genotyping results, that this panel of SNP AIMs provides good power for detecting disease-associated chromosomal segments for genes with modest ethnicity risk ratios. A reduced set of 5,287 SNP AIMs captured almost the same admixture mapping information, but smaller SNP sets showed substantial drop-off in admixture mapping information and power. The results will enable studies of type 2 diabetes, rheumatoid arthritis, and other diseases among which epidemiological studies suggest differences in the distribution of ancestry-associated susceptibility.

Primary biliary cholangitis (PBC) is an autoimmune liver disease with a strong hereditary component. Here, we report a genome-wide association study that included 1,122 PBC cases and 4,036 controls of Han Chinese descent, with subsequent replication in a separate cohort of 907 PBC cases and 2,127 controls. Our results show genome-wide association of 14 PBC risk loci including previously identified 6p21 (HLA-DRA and DPB1), 17q12 (ORMDL3), 3q13.33 (CD80), 2q32.3 (STAT1/STAT4), 3q25.33 (IL12A), 4q24 (NF-κB) and 22q13.1 (RPL3/SYNGR1). We also identified variants in IL21, IL21R, CD28/CTLA4/ICOS, CD58, ARID3A and IL16 as novel PBC risk loci. These new findings and histochemical studies showing enhanced expression of IL21 and IL21R in PBC livers (particularly in the hepatic portal tracks) support a disease mechanism in which the deregulation of the IL21 signalling pathway, in addition to CD4 T-cell activation and T-cell co-stimulation are critical components in the development of PBC.

Allele expression (AE) analysis robustly measures cis-regulatory effects. Here, we present and demonstrate the utility of a vast AE resource generated from the GTEx v8 release, containing 15,253 samples spanning 54 human tissues for a total of 431 million measurements of AE at the SNP level and 153 million measurements at the haplotype level. In addition, we develop an extension of our tool phASER that allows effect sizes of cis-regulatory variants to be estimated using haplotype-level AE data. This AE resource is the largest to date, and we are able to make haplotype-level data publicly available. We anticipate that the availability of this resource will enable future studies of regulatory variation across human tissues.

The European Commission requested EFSA to carry out a scientific evaluation on four studies on the potential toxicity of titanium dioxide (TiO2) used as a food additive (E 171) and to indicate whether they would merit re-opening the existing opinion of EFSA on the safety of TiO2 (E 171) as a food additive. The results of the Bettini et al. (2017) study did not provide enough justification for a new carcinogenicity study, but, should additional useful mechanistic information become available, this could be reconsidered in future. The new in vitro findings in the Proquin et al. (2017) study did not modify the conclusion on the genotoxicity of TiO2 as stated in the previous EFSA opinion of 2016 on the safety of TiO2 (E 171) as a food additive. The effects of engineered TiO2 nanoparticles reported by the Guo et al. (2017) study were of uncertain biological significance and therefore of limited relevance for the risk assessment of the food additive TiO2 (E 171). There was significant uncertainty in the risk assessment performed by Heringa et al. (2016), which did not include a weight of evidence analysis of the whole database. The Panel considered that the four studies evaluated, highlighted some concerns but with uncertainties, therefore their relevance for the risk assessment was considered limited and further research would be needed to decrease the level of uncertainties. Overall, three of the studies, reporting that TiO2 induced various effects in in vitro and in vivo models, may be useful for hazard identification of TiO2. In the fourth study by Heringa et al. (2016), numerous assumptions were made, which resulted in large uncertainty in their conclusion. Altogether, the Panel concluded that the outcome of the four studies did not merit re-opening the existing opinion of EFSA related to the safety of TiO2 (E 171) as a food additive.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Despite the growing number of genome-wide association studies (GWASs), it remains unclear to what extent gene-by-gene and gene-by-environment interactions influence complex traits in humans. The magnitude of genetic interactions in complex traits has been difficult to quantify because GWASs are generally underpowered to detect individual interactions of small effect. Here, we develop a method to test for genetic interactions that aggregates information across all trait-associated loci. Specifically, we test whether SNPs in regions of European ancestry shared between European American and admixed African American individuals have the same causal effect sizes. We hypothesize that in African Americans, the presence of genetic interactions will drive the causal effect sizes of SNPs in regions of European ancestry to be more similar to those of SNPs in regions of African ancestry. We apply our method to two traits: gene expression in 296 African Americans and 482 European Americans in the Multi-Ethnic Study of Atherosclerosis (MESA) and low-density lipoprotein cholesterol (LDL-C) in 74K African Americans and 296K European Americans in the Million Veteran Program (MVP). We find significant evidence for genetic interactions in our analysis of gene expression; for LDL-C, we observe a similar point estimate, although this is not significant, most likely due to lower statistical power. These results suggest that gene-by-gene or gene-by-environment interactions modify the effect sizes of causal variants in human complex traits.

Multi-omics datasets are becoming more common, necessitating better integration methods to realize their revolutionary potential. Here, we introduce multi-set correlation and factor analysis (MCFA), an unsupervised integration method tailored to the unique challenges of high-dimensional genomics data that enables fast inference of shared and private factors. We used MCFA to integrate methylation markers, protein expression, RNA expression, and metabolite levels in 614 diverse samples from the Trans-Omics for Precision Medicine/Multi-Ethnic Study of Atherosclerosis multi-omics pilot. Samples cluster strongly by ancestry in the shared space, even in the absence of genetic information, while private spaces frequently capture dataset-specific technical variation. Finally, we integrated genetic data by conducting a genome-wide association study (GWAS) of our inferred factors, observing that several factors are enriched for GWAS hits and trans-expression quantitative trait loci. Two of these factors appear to be related to metabolic disease. Our study provides a foundation and framework for further integrative analysis of ever larger multi-modal genomic datasets.