Meta-analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for meta-analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harboring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found that genes implicated in upstream developmental processes had greater cis-regulatory proportions than downstream cellular function genes, and that ligands were associated with higher cis-regulatory proportions than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

A novel disease mechanism leading to the expression of a disallowed gene in the pancreatic beta-cell identified by non-coding, regulatory mutations controlling HK1

Gene expression is tightly regulated with many genes exhibiting cell-specific silencing when their protein product would disrupt normal cellular function. This silencing is largely controlled by non-coding elements and their disruption might cause human disease. We performed gene-agnostic screening of the non-coding regions to discover new molecular causes of congenital hyperinsulinism. This identified 14 non-coding de novo mutations affecting a 42bp conserved region encompassed by a regulatory element in intron 2 of Hexokinase 1 (HK1), a pancreatic beta-cell disallowed gene. We demonstrated that these mutations resulted in expression of HK1 in the pancreatic beta-cells causing inappropriate insulin secretion and congenital hyperinsulinism. These mutations identify a regulatory region critical for cell-specific silencing. Importantly, this has revealed a new disease mechanism for non-coding mutations that cause inappropriate expression of a disallowed gene.

Mutational sources of trans-regulatory variation affecting gene expression in Saccharomyces cerevisiae

Heritable variation in a gene’s expression arises from mutations impacting cis- and trans-acting components of its regulatory network. Here, we investigate how trans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations with trans-regulatory effects on expression of the same focal gene in Saccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that these trans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of the trans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism, and signal transduction. These data show how genetic changes in diverse types of genes can impact a gene’s expression in trans, revealing properties of trans-regulatory mutations that provide the raw material for trans-regulatory variation segregating within natural populations.

OmpF Downregulation Mediated by Sigma E or OmpR Activation Confers Cefalexin Resistance in Escherichia coli in the Absence of Acquired β-Lactamases.

Cefalexin is a widely used 1 st generation cephalosporin, and resistance in Escherichia coli is caused by Extended-Spectrum (e.g., CTX-M) and AmpC β-lactamase production and therefore frequently coincides with 3 rd generation cephalosporin resistance. However, we have recently identified large numbers of E. coli isolates from human infections, and from cattle, where cefalexin resistance is not β-lactamase mediated. Here we show, by studying laboratory selected mutants, clinical isolates, and isolates from cattle, that OmpF porin disruption or downregulation is a major cause of cefalexin resistance in E. coli . Importantly, we identify multiple regulatory mutations that cause OmpF downregulation. In addition to mutation of ompR , already known to downregulate OmpF and OmpC porin production, we find that rseA mutation, which strongly activates the Sigma E regulon, greatly increasing DegP production, which degrades OmpF, OmpC and OmpA. Furthermore, we reveal that mutations affecting lipopolysaccharide structure, exemplified by the loss of GmhB, essential for lipopolysaccharide heptosylation, also modestly activate DegP production, resulting in OmpF degradation. Remarkably, given the critical importance attached to such systems for normal E. coli physiology, we find evidence for DegP-mediated OmpF downregulation, gmhB and rseA loss of function mutation in E. coli isolates derived from human infections. Finally, we show that these regulatory mutations enhance the ability of group 1 CTX-M β-lactamase to confer reduced carbapenem susceptibility, particularly those mutations that cause OmpC in addition to OmpF downregulation.

Functional annotation of noncoding mutations in cancer

In a cancer genome, the noncoding sequence contains the vast majority of somatic mutations. While very few are expected to be cancer drivers, those affecting regulatory elements have the potential to have downstream effects on gene regulation that may contribute to cancer progression. To prioritize regulatory mutations, we screened somatic mutations in the Pan-Cancer Analysis of Whole Genomes cohort of 2,515 cancer genomes on individual bases to assess their potential regulatory roles in their respective cancer types. We found a highly significant enrichment of regulatory mutations associated with the deamination signature overlapping a CpG site in the CCAAT/Enhancer Binding Protein β recognition sites in many cancer types. Overall, 5,749 mutated regulatory elements were identified in 1,844 tumor samples from 39 cohorts containing 11,962 candidate regulatory mutations. Our analysis indicated 20 or more regulatory mutations in 5.5% of the samples, and an overall average of six per tumor. Several recurrent elements were identified, and major cancer-related pathways were significantly enriched for genes nearby the mutated regulatory elements. Our results provide a detailed view of the role of regulatory elements in cancer genomes.

Microsatellites as Agents of Adaptive Change: An RNA-Seq-Based Comparative Study of Transcriptomes from Five Helianthus Species

Mutations that provide environment-dependent selective advantages drive adaptive divergence among species. Many phenotypic differences among related species are more likely to result from gene expression divergence rather than from non-synonymous mutations. In this regard, cis-regulatory mutations play an important part in generating functionally significant variation. Some proposed mechanisms that explore the role of cis-regulatory mutations in gene expression divergence involve microsatellites. Microsatellites exhibit high mutation rates achieved through symmetric or asymmetric mutation processes and are abundant in both coding and non-coding regions in positions that could influence gene function and products. Here we tested the hypothesis that microsatellites contribute to gene expression divergence among species with 50 individuals from five closely related Helianthus species using an RNA-seq approach. Differential expression analyses of the transcriptomes revealed that genes containing microsatellites in non-coding regions (UTRs and introns) are more likely to be differentially expressed among species when compared to genes with microsatellites in the coding regions and transcripts lacking microsatellites. We detected a greater proportion of shared microsatellites in 5′UTRs and coding regions compared to 3′UTRs and non-coding transcripts among Helianthus spp. Furthermore, allele frequency differences measured by pairwise FST at single nucleotide polymorphisms (SNPs), indicate greater genetic divergence in transcripts containing microsatellites compared to those lacking microsatellites. A gene ontology (GO) analysis revealed that microsatellite-containing differentially expressed genes are significantly enriched for GO terms associated with regulation of transcription and transcription factor activity. Collectively, our study provides compelling evidence to support the role of microsatellites in gene expression divergence.

OmpF Downregulation Mediated by Sigma E or OmpR Activation Confers Cefalexin Resistance in Escherichia coli in the Absence of Acquired β-Lactamases.

Cefalexin is a widely used 1st generation cephalosporin, and resistance in Escherichia coli is caused by Extended-Spectrum (e.g. CTX-M) and AmpC β-lactamase production and therefore frequently coincides with 3rd generation cephalosporin resistance. However, we have recently identified large numbers of E. coli isolates from human infections, and from cattle, where cefalexin resistance is not β-lactamase mediated. Here we show, by studying laboratory selected mutants, clinical isolates, and isolates from cattle, that OmpF porin disruption or downregulation is a major cause of cefalexin resistance in E. coli. Importantly, we identify multiple regulatory mutations that cause OmpF downregulation. In addition to mutation of ompR, already known to downregulate OmpF and OmpC porin production, we find that rseA mutation, which strongly activates the Sigma E regulon, greatly increasing DegP production, which degrades OmpF, OmpC and OmpA porins. Furthermore, we reveal that mutations affecting lipopolysaccharide structure, exemplified by the loss of GmhB, essential for lipopolysaccharide heptosylation, also modestly activate DegP production, resulting in OmpF degradation. Remarkably, given the critical importance attached to such systems for normal E. coli physiology, we find evidence for DegP-mediated OmpF downregulation, gmhB and rseA loss of function mutation in E. coli isolates derived from human infections. Finally, we show that these regulatory mutations enhance the ability of group 1 CTX-M β-lactamase to confer reduced carbapenem susceptibility, particularly those mutations that cause OmpC in addition to OmpF downregulation.

Microsatellites as Agents of Adaptive Change: An RNA-seq based Comparative Study of Transcriptomes from Five Helianthus Species

Mutations that provide environment dependent selective advantages drive adaptive divergence among species. Many phenotypic differences among related species are more likely to result from gene expression divergence rather than from non-synonymous mutations. In this regard, cis-regulatory mutations play an important part in generating functionally significant variation. Some proposed mechanisms that explore the role of cis-regulatory mutations in gene expression divergence involve microsatellites. Microsatellites exhibit high mutation rates and are abundant in both coding and non-coding regions and could influence gene function and products. Here we tested the hypothesis that microsatellites contribute to gene expression divergence among species with 50 individuals from nine closely related Helianthus species using an RNA-seq approach. Differential expression analyses of the transcriptomes revealed that genes containing microsatellites in non-coding regions (UTRs and introns) are more likely to be differentially expressed among species when compared to genes with microsatellites in the coding regions and transcripts lacking microsatellites. We detected a greater proportion of shared microsatellites in 5’UTRs and coding regions compared to 3’UTRs and non-coding transcripts among Helianthus spp. Further, allele frequency differences measured by pairwise FST at single nucleotide polymorphisms (SNPs), indicate greater genetic divergence in transcripts containing microsatellites compared to those lacking microsatellites. A gene ontology (GO) analysis revealed that microsatellite-containing differentially expressed genes are significantly enriched for GO terms associated with regulation of transcription and transcription factor activity. Collectively, our study provides compelling evidence to support the role of microsatellites in gene expression divergence.

Rapid adaptation of a complex trait during experimental evolution of Mycobacterium tuberculosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis (M. tb), is a leading cause of death due to infectious disease. TB is not traditionally associated with biofilms, but M. tb biofilms are linked with drug and immune tolerance and there is increasing recognition of their potential role in the recalcitrance of TB infections. Here we used M. tb experimental evolution to investigate this complex phenotype and identify candidate loci controlling biofilm formation. We identified novel candidate loci, adding to our understanding of the genetic architecture underlying M. tb biofilm development. Under selective pressure to grow as a biofilm, regulatory mutations rapidly swept to fixation and were associated with changes in multiple traits including extracellular matrix production, cell size, and growth rate. Genetic and phenotypic paths to enhanced biofilm growth varied according to the genetic background of the parent strain, suggesting that epistatic interactions are important in M. tb adaptation to changing environments.