Description: Background The high prevalence of tobacco use in some developing nations, including Bangladesh, poses several public health challenges for these populations. Smoking behaviour is determined by genetic and environmental factors; however, the genetic determinants of smoking behaviour have not been previously examined in a Bangladeshi or South Asian population. We performed a genome-wide association study (GWAS) of tobacco smoking behaviour among a population-based sample of 5354 (2035 ever smokers and 3319 never smokers) men and women in Bangladesh. Methods Genome-wide association analyses were conducted for smoking initiation (ever vs never smokers), smoking quantity (cigarettes per day), age of smoking initiation, and smoking cessation (former vs current smokers). Sex-stratified associations were performed for smoking initiation. Results We observed associations for smoking initiation in the SLC39A11 region at 17q21.31 (rs2567519, p=1.33 x 10 –7 ) among men and in the SLCO3A1 region at 15q26 (rs12912184, p=9.32 x 10 –8 ) among women. Conclusions These findings suggest possible underlying mechanisms related to solute carrier transporter genes, which transport neurotransmitters, nutrients, heavy metals and other substrates into cells, for smoking initiation in a South Asian population in a sex-specific pattern. Genetic markers could have potential translational implications for the prevention or treatment of tobacco use and addiction in South Asian populations and warrant further exploration.

Description: Changes in chromatin state play important roles in cell fate transitions. Current computational approaches to analyze chromatin modifications across multiple cell types do not model how the cell types are related on a lineage or over time . To overcome this limitation, we developed a method called Chromatin Module INference on Trees (CMINT), a probabilistic clustering approach to systematically capture chromatin state dynamics across multiple cell types. Compared to existing approaches, CMINT can handle complex lineage topologies, capture higher quality clusters, and reliably detect chromatin transitions between cell types. We applied CMINT to gain novel insights in two complex processes: reprogramming to induced pluripotent stem cells (iPSCs) and hematopoiesis. In reprogramming, chromatin changes could occur without large gene expression changes, different combinations of activating marks were associated with specific reprogramming factors, there was an order of acquisition of chromatin marks at pluripotency loci, and multivalent states (comprising previously undetermined combinations of activating and repressive histone modifications) were enriched for CTCF. In the hematopoietic system, we defined critical decision points in the lineage tree, identified regulatory elements that were enriched in cell-type–specific regions, and found that the underlying chromatin state was achieved by specific erasure of preexisting chromatin marks in the precursor cell or by de novo assembly. Our method provides a systematic approach to model the dynamics of chromatin state to provide novel insights into the relationships among cell types in diverse cell-fate specification processes.

Description: Comparative functional genomics studies the evolution of biological processes by analyzing functional data, such as gene expression profiles, across species. A major challenge is to compare profiles collected in a complex phylogeny. Here, we present Arboretum, a novel scalable computational algorithm that integrates expression data from multiple species with species and gene phylogenies to infer modules of coexpressed genes in extant species and their evolutionary histories. We also develop new, generally applicable measures of conservation and divergence in gene regulatory modules to assess the impact of changes in gene content and expression on module evolution. We used Arboretum to study the evolution of the transcriptional response to heat shock in eight species of Ascomycota fungi and to reconstruct modules of the ancestral environmental stress response (ESR). We found substantial conservation in the stress response across species and in the reconstructed components of the ancestral ESR modules. The greatest divergence was in the most induced stress, primarily through module expansion. The divergence of the heat stress response exceeds that observed in the response to glucose depletion in the same species. Arboretum and its associated analyses provide a comprehensive framework to systematically study regulatory evolution of condition-specific responses.

Description: Background Leucocyte telomere length (TL) is a potential biomarker of ageing and risk for age-related disease. Leucocyte TL is heritable and shows substantial differences by race/ethnicity. Recent genome-wide association studies (GWAS) report ~10 loci harbouring SNPs associated with leucocyte TL, but these studies focus primarily on populations of European ancestry. Objective This study aims to enhance our understanding of genetic determinants of TL across populations. Methods We performed a GWAS of TL using data on 5075 Bangladeshi adults. We measured TL using one of two technologies (qPCR or a Luminex-based method) and used standardised variables as TL phenotypes. Results Our results replicate previously reported associations in the TERC and TERT regions (P=2.2 x 10 –8 and P=6.4 x 10 –6 , respectively). We observed a novel association signal in the RTEL1 gene (intronic SNP rs2297439; P=2.82 x 10 –7 ) that is independent of previously reported TL-associated SNPs in this region. The minor allele for rs2297439 is common in South Asian populations (≥0.25) but at lower frequencies in other populations (eg, 0.07 in Northern Europeans). Among the eight other previously reported association signals, all were directionally consistent with our study, but only rs8105767 ( ZNF208 ) was nominally significant (P=0.003). SNP-based heritability estimates were as high as 44% when analysing close relatives but much lower when analysing distant relatives only. Conclusions In this first GWAS of TL in a South Asian population, we replicate some, but not all, of the loci reported in prior GWAS of individuals of European ancestry, and we identify a novel second association signal at the RTEL1 locus.

Description: Background Opioids such as morphine are widely used for the management of pain associated with acute pancreatitis. Interestingly, opioids are also known to affect the immune system and modulate inflammatory pathways in non-pancreatic diseases. However, the impact of morphine on the progression of acute pancreatitis has never been evaluated. In the current study, we evaluated the impact of morphine on the progression and severity of acute pancreatitis. Methods Effect of morphine treatment on acute pancreatitis in caerulein, L-arginine and ethanol–palmitoleic acid models was evaluated after induction of the disease. Inflammatory response, gut permeability and bacterial translocation were compared. Experiments were repeated in mu (µ) opioid receptor knockout mice (MORKO) and in wild-type mice in the presence of opioid receptor antagonist naltrexone to evaluate the role of µ-opioid receptors in morphine’s effect on acute pancreatitis. Effect of morphine treatment on pathways activated during pancreatic regeneration like sonic Hedgehog and activation of embryonic transcription factors like pdx-1 and ptf-1 were measured by immunofluorescence and quantitative PCR. Results Histological data show that treatment with morphine after induction of acute pancreatitis exacerbates the disease with increased pancreatic neutrophilic infiltration and necrosis in all three models of acute pancreatitis. Morphine also exacerbated acute pancreatitis-induced gut permeabilisation and bacteraemia. These effects were antagonised in the MORKO mice or in the presence of naltrexone suggesting that morphine’s effect on severity of acute pancreatitis are mediated through the µ-opioid receptors. Morphine treatment delayed macrophage infiltration, sonic Hedgehog pathway activation and expression of pdx-1 and ptf-1. Conclusion Morphine treatment worsens the severity of acute pancreatitis and delays resolution and regeneration. Considering our results, the safety of morphine for analgesia during acute pancreatitis should be re-evaluated in future human studies.

Description: Gaining insights on gene regulation from large-scale functional data sets is a grand challenge in systems biology. In this article, we develop and apply methods for transcriptional regulatory network inference from diverse functional genomics data sets and demonstrate their value for gene function and gene expression prediction. We formulate the network inference problem in a machine-learning framework and use both supervised and unsupervised methods to predict regulatory edges by integrating transcription factor (TF) binding, evolutionarily conserved sequence motifs, gene expression, and chromatin modification data sets as input features. Applying these methods to Drosophila melanogaster, we predict ~300,000 regulatory edges in a network of ~600 TFs and 12,000 target genes. We validate our predictions using known regulatory interactions, gene functional annotations, tissue-specific expression, protein–protein interactions, and three-dimensional maps of chromosome conformation. We use the inferred network to identify putative functions for hundreds of previously uncharacterized genes, including many in nervous system development, which are independently confirmed based on their tissue-specific expression patterns. Last, we use the regulatory network to predict target gene expression levels as a function of TF expression, and find significantly higher predictive power for integrative networks than for motif or ChIP-based networks. Our work reveals the complementarity between physical evidence of regulatory interactions (TF binding, motif conservation) and functional evidence (coordinated expression or chromatin patterns) and demonstrates the power of data integration for network inference and studies of gene regulation at the systems level.

Description: The order of genes in eukaryotic genomes has generally been assumed to be neutral, since gene order is largely scrambled over evolutionary time. Only a handful of exceptional examples are known, typically involving deeply conserved clusters of tandemly duplicated genes (e.g., Hox genes and histones). Here we report the first systematic survey of microsynteny conservation across metazoans, utilizing 17 genome sequences. We identified nearly 600 pairs of unrelated genes that have remained tightly physically linked in diverse lineages across over 600 million years of evolution. Integrating sequence conservation, gene expression data, gene function, epigenetic marks, and other genomic features, we provide extensive evidence that many conserved ancient linkages involve (1) the coordinated transcription of neighboring genes, or (2) genomic regulatory blocks (GRBs) in which transcriptional enhancers controlling developmental genes are contained within nearby bystander genes. In addition, we generated ChIP-seq data for key histone modifications in zebrafish embryos, which provided further evidence of putative GRBs in embryonic development. Finally, using chromosome conformation capture (3C) assays and stable transgenic experiments, we demonstrate that enhancers within bystander genes drive the expression of genes such as Otx and Islet , critical regulators of central nervous system development across bilaterians. These results suggest that ancient genomic functional associations are far more common than previously thought—involving ~12% of the ancestral bilaterian genome—and that cis -regulatory constraints are crucial in determining metazoan genome architecture.

Description: Background Creatine transporter deficiency is a monogenic cause of X-linked intellectual disability. Since its first description in 2001 several case reports have been published but an overview of phenotype, genotype and phenotype–genotype correlation has been lacking. Methods We performed a retrospective study of clinical, biochemical and molecular genetic data of 101 males with X-linked creatine transporter deficiency from 85 families with a pathogenic mutation in the creatine transporter gene ( SLC6A8 ). Results and conclusions Most patients developed moderate to severe intellectual disability; mild intellectual disability was rare in adult patients. Speech language development was especially delayed but almost a third of the patients were able to speak in sentences. Besides behavioural problems and seizures, mild to moderate motor dysfunction, including extrapyramidal movement abnormalities, and gastrointestinal problems were frequent clinical features. Urinary creatine to creatinine ratio proved to be a reliable screening method besides MR spectroscopy, molecular genetic testing and creatine uptake studies, allowing definition of diagnostic guidelines. A third of patients had a de novo mutation in the SLC6A8 gene. Mothers with an affected son with a de novo mutation should be counselled about a recurrence risk in further pregnancies due to the possibility of low level somatic or germline mosaicism. Missense mutations with residual activity might be associated with a milder phenotype and large deletions extending beyond the 3' end of the SLC6A8 gene with a more severe phenotype. Evaluation of the biochemical phenotype revealed unexpected high creatine levels in cerebrospinal fluid suggesting that the brain is able to synthesise creatine and that the cerebral creatine deficiency is caused by a defect in the reuptake of creatine within the neurones.