Genomic and transcriptomic analyses reveal a tandem amplification unit of 11 genes and mutations in mismatch repair genes in methotrexate-resistant HT-29 cells

DHFR gene amplification is commonly present in methotrexate (MTX)-resistant colon cancer cells and acute lymphoblastic leukemia. In this study, we proposed an integrative framework to characterize the amplified region by using a combination of single-molecule real-time sequencing, next-generation optical mapping, and chromosome conformation capture (Hi-C). We identified an amplification unit spanning 11 genes, from the DHFR gene to the ATP6AP1L gene position, with high adjusted interaction frequencies on chromosome 5 (~2.2 Mbp) and a twenty-fold tandemly amplified region, and novel inversions at the start and end positions of the amplified region as well as frameshift insertions in most of the MSH and MLH genes were detected. These mutations might stimulate chromosomal breakage and cause the dysregulation of mismatch repair. Characterizing the tandem gene-amplified unit may be critical for identifying the mechanisms that trigger genomic rearrangements. These findings may provide new insight into the mechanisms underlying the amplification process and the evolution of drug resistance.

Reversible gene silencing through frameshift indels and frameshift scars provide adaptive plasticity for Mycobacterium tuberculosis

Mycobacterium tuberculosis can adapt to changing environments by non-heritable mechanisms. Frame-shifting insertions and deletions (indels) may also participate in adaptation through gene disruption, which could be reversed by secondary introduction of a frame-restoring indel. We present ScarTrek, a program that scans genomic data for indels, including those that together disrupt and restore a gene’s reading frame, producing “frame-shift scars” suggestive of reversible gene inactivation. We use ScarTrek to analyze 5977 clinical M. tuberculosis isolates. We show that indel frequency inversely correlates with genomic linguistic complexity and varies with gene-position and gene-essentiality. Using ScarTrek, we detect 74 unique frame-shift scars in 48 genes, with a 3.74% population-level incidence of unique scar events. We find multiple scars in the ESX-1 gene cluster. Six scars show evidence of convergent evolution while the rest shared a common ancestor. Our results suggest that sequential indels are a mechanism for reversible gene silencing and adaptation in M. tuberculosis.

Consequences of producing DNA gyrase from a synthetic gyrBA operon in Salmonella enterica serovar Typhimurium

DNA gyrase is an essential type II topoisomerase that is composed of two subunits, GyrA and GyrB and has an A2B2 structure. Although both subunits are required in equal proportions to form DNA gyrase, the gyrA and gyrB genes that encode them in Salmonella (and in many other bacteria) are at widely separated locations on the chromosome, are under separate transcriptional control and are present in different copy numbers in rapidly growing bacteria (gyrA is near the terminus of chromosome replication while gyrB is near the origin). We generated a synthetic gyrBA operon at the oriC-proximal location of gyrB to test the significance of the gyrase gene position for Salmonella physiology. Producing gyrase from an operon did not alter growth kinetics, cell morphology, competitive fitness index, or sensitivity to some gyrase-inhibiting antibiotics. However, the operon strain had altered DNA supercoiling set points, its SPI-2 virulence genes were expressed at a reduced level and its survival was reduced in macrophage. The gyrB gene could not be deleted from its oriC-proximal location, even in a gyrB merodiploid strain. We discuss the physiological significance of the different gyrA and gyrB gene arrangements found naturally in Salmonella and other bacteria.

Genes on Different Strands Mark Boundaries Associated with Co-regulation Domains

ABSTRACTGene regulation is influenced by chromatin conformation. Current models suggest that topologically associating domains (TADs) act as regulatory units, which could also include distinct co-expression domains (CODs) favouring correlated gene expression. We integrated publicly available RNA-seq, ChIP-seq and Hi-C data from A549 cells stimulated with the glucocorticoid dexamethasone to explore how differentially expressed genes are co-regulated among TADs and CODs. Interestingly, we found that gene position and orientation also impact co-regulation. Indeed, divergent and convergent pairs of genes we enriched at sub-TAD boundaries, forming distinct CODs. We also found that genes at COD boundaries were less likely to be separated by structural proteins such as Cohesin and CTCF. A complementary analysis of lung expression quantitative trait loci (eQTL) demonstrated that genes affected by the same variant were more likely to be found on the same strand while lacking a TAD boundary. Taken together, these results suggest a model where gene orientation can provide a boundary between CODs, at the sub-TAD level, thus affecting their likelihood of co-regulation.