Genome-wide association studies (GWAS) of diseases and traits have found associations to gene regions but not the functional SNP or the gene mediating the effect. Difference in gene regulatory signals can be detected using chromatin immunoprecipitation and next-gen sequencing (ChIP-seq) of transcription factors or histone modifications by aligning reads to known polymorphisms in individual genomes. The aim was to identify such regulatory elements in the human liver to understand the genetics behind type 2 diabetes and metabolic diseases. The genome of liver tissue was sequenced using 10X Genomics technology to call polymorphic positions. Using ChIP-seq for two histone modifications, H3K4me3 and H3K27ac, and the transcription factor CTCF, and our established bioinformatics pipeline, we detected sites with significant difference in signal between the alleles. We detected 2329 allele-specific SNPs (AS-SNPs) including 25 associated to GWAS SNPs linked to liver biology, e.g., 4 AS-SNPs at two type 2 diabetes loci. Two hundred ninety-two AS-SNPs were associated to liver gene expression in GTEx, and 134 AS-SNPs were located on 166 candidate functional motifs and most of them in EGR1-binding sites. This study provides a valuable collection of candidate liver regulatory elements for further experimental validation.

Members of the Mycobacterium chelonae-abscessus complex (MCAC) are close to the mycobacterial ancestor and includes both human, animal and fish pathogens. We present the genomes of 14 members of this complex: the complete genomes of Mycobacterium salmoniphilum and Mycobacterium chelonae type strains, seven M. salmoniphilum isolates, and five M. salmoniphilum-like strains including strains isolated during an outbreak in an animal facility at Uppsala University. Average nucleotide identity (ANI) analysis and core gene phylogeny revealed that the M. salmoniphilum-like strains are variants of the human pathogen Mycobacterium franklinii and phylogenetically close to Mycobacterium abscessus. Our data further suggested that M. salmoniphilum separates into three branches named group I, II and III with the M. salmoniphilum type strain belonging to group II. Among predicted virulence factors, the presence of phospholipase C (plcC), which is a major virulence factor that makes M. abscessus highly cytotoxic to mouse macrophages, and that M. franklinii originally was isolated from infected humans make it plausible that the outbreak in the animal facility was caused by a M. salmoniphilum-like strain. Interestingly, M. salmoniphilum-like was isolated from tap water suggesting that it can be present in the environment. Moreover, we predicted the presence of mutational hotspots in the M. salmoniphilum isolates and 26% of these hotspots overlap with genes categorized as having roles in virulence, disease and defense. We also provide data about key genes involved in transcription and translation such as sigma factor, ribosomal protein and tRNA genes.

This study aimed to characterize the population pharmacokinetics and pharmacogenetics of ethambutol in tuberculosis-HIV-coinfected adult patients. Ethambutol plasma concentrations, determined by liquid chromatography-tandem mass spectrometry, in 63 patients receiving ethambutol as part of rifampin-based fixed-dose combination therapy for tuberculosis were analyzed using nonlinear mixed-effects modeling. A one-compartment disposition model with first-order elimination and four transit compartments prior to first-order absorption was found to adequately describe the concentration-time profiles of ethambutol in plasma. Body weight was implemented as an allometric function on the clearance and volume parameters. Estimates of oral clearance and volume of distribution were 77.4 liters/h and 76.2 liters, respectively. A G/A mutation with regard to CYP1A2 2159 G>A was associated with a 50% reduction in relative bioavailability. Simulations revealed that doses of 30 mg/kg of body weight and 50 mg/kg for G/G and G/A carriers, respectively, would result in clinically adequate exposure. The results presented here suggest that CYP1A2 polymorphism affects ethambutol exposure in this population and that current treatment guidelines may result in underexposure in patients coinfected with tuberculosis and HIV. Based on simulations, a dose increase from15 to 20 mg/kg to 30 mg/kg is suggested. However, the 50-mg/kg dose required to reach therapeutic exposure in G/A carriers may be inappropriate due to the dose-dependent toxicity of ethambutol. Additional studies are required to further investigate CYP450 polymorphism effects on ethambutol pharmacokinetics.

There is an important interplay between epigenetic factors and body weight, and previous work has identified ten sites where DNA methylation is robustly associated with body mass index (BMI) cross-sectionally. However, interpretation of the associations is complicated by the substantial changes in BMI often occurring in late-life, and the fact that methylation is often driven by genetic variation. This study therefore investigated the longitudinal association between these ten sites and BMI from midlife to late-life, and whether associations persist after controlling for genetic factors. We used data from 535 individuals (mean age 68) in the Swedish Adoption/Twin Study of Aging (SATSA) with longitudinal measures of both DNA methylation from blood samples and BMI, spanning up to 20 years. Methylation levels were measured with the Infinium Human Methylation 450K or Infinium MethylationEpic array, with seven of the ten sites passing quality control. Latent growth curve models were applied to investigate longitudinal associations between methylation and BMI, and between-within models to study associations within twin pairs, thus adjusting for genetic factors. Baseline DNA methylation levels at five of the seven sites were associated with BMI level at age 65 (cg00574958 [CPT1A]; cg11024682 [SREBF1]), and/or change (cg06192883 [MYO5C]; cg06946797 [RMI2]; cg08857797 [VPS25]). For four of the five sites, the associations remained comparable within twin pairs. However, the effects of cg06192883 were substantially attenuated within pairs. No change in DNA methylation was detected for any of the seven evaluated sites. Five of the seven sites investigated were associated with late-life level and/or change in BMI. The effects for four of the sites remained similar when examined within twin pairs, indicating that the associations are mainly environmentally driven. However, the substantial attenuation in the association between cg06192883 and late-life BMI within pairs points to the importance of genetic factors in this association.

Molecular characterization of genome-wide association study (GWAS) loci can uncover key genes and biological mechanisms underpinning complex traits and diseases. Here we present deep, high-throughput characterization of gene regulatory mechanisms underlying prostate cancer risk loci. Our methodology integrates data from 295 prostate cancer chromatin immunoprecipitation and sequencing experiments with genotype and gene expression data from 602 prostate tumor samples. The analysis identifies new gene regulatory mechanisms affected by risk locus SNPs, including widespread disruption of ternary androgen receptor (AR)-FOXA1 and AR-HOXB13 complexes and competitive binding mechanisms. We identify 57 expression quantitative trait loci at 35 risk loci, which we validate through analysis of allele-specific expression. We further validate predicted regulatory SNPs and target genes in prostate cancer cell line models. Finally, our integrated analysis can be accessed through an interactive visualization tool. This analysis elucidates how genome sequence variation affects disease predisposition via gene regulatory mechanisms and identifies relevant genes for downstream biomarker and drug development.

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool-seq data to generate a de novo genome assembly for mining exons, upon which Pool-seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual-based single nucleotide polymorphisms [SNPs]) and from mapping the Pool-seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F ST) between the two introduced populations exceeds that of the naturally sympatric populations (F ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ ( ≈ 0.002 and π ¯ ≈ 0.001 when mapping to π ¯S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high-quality reference assembly from a divergent species. We conclude that the Pool-seq-only approach can be suitable for detecting and quantifying genome-wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.

The Norwegian lemming (Lemmus lemmus) is a small rodent distributed across the Fennoscandian mountain tundra and the Kola Peninsula. The Norwegian lemming likely evolved during the Late Pleistocene and inhabited Fennoscandia shortly prior to the Last Glacial Maximum. However, the exact timing and origins of the species, and its phylogenetic position relative to the closely related Siberian lemming (Lemmus sibiricus) remain disputed. Moreover, the presence of ancient or contemporary gene flow between both species is largely untested. The Norwegian lemming displays characteristic phenotypic and behavioral adaptations (e.g., coat color, aggression) that are not present in other Lemmus species. We generated a de novo genome assembly for the Norwegian lemming and resequenced nine modern and two ancient Lemmus spp. genomes. We show that all Lemmus species form distinct monophyletic clades, with concordant topology between the mitochondrial and nuclear genome phylogenies. The Siberian lemming is divided into two distinct but paraphyletic clades, one in the east and one in the west, where the western clade represents a sister taxon to the Norwegian lemming. We estimate that the Norwegian and western Siberian lemming diverged shortly before the Last Glacial Maximum, making the Norwegian lemming one of the youngest known mammalian species. We did not find any indication of gene flow between L. lemmus and L. sibiricus, suggesting postglacial isolation of L. lemmus. Furthermore, we identify species-specific genomic differences in genes related to coat color and fat transport, which are likely associated with the distinctive coloration and overwintering behavior observed in the Norwegian lemming.