Transposon Removal Reveals Their Adaptive Fitness Contribution.

Transposable elements are molecular parasites that persist in their host genome by generating new copies to outpace natural selection. Transposable elements exert a large influence on host genome evolution, in some cases providing adaptive changes. Here we measure the fitness effect of the transposable element insertions in the fission yeast Schizosaccharomyces pombe type strain by removing all insertions of its only native transposable element family, the long terminal repeat retrotransposon Tf2. We show that Tf2 elements provide a positive fitness contribution to its host. Tf2 ablation results in changes to the regulation of a mitochondrial gene and, consistently, the fitness effect are sensitive to growth conditions. We propose that Tf2 influences host fitness in a directed manner by dynamically rewiring the transcriptional response to metabolic stress.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai, Drosophila takahashii, Drosophila bipectinata, and Drosophila ananassae.

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma. However, the assembled portion of the Muller F element (orthologous to the fourth chromosome in Drosophila melanogaster) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to >18 Mb. Here, we present chromosome-level long-read genome assemblies for 4 Drosophila species with expanded F elements ranging in size from 2.3 to 20.5 Mb. Each Muller element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai , Drosophila takahashii , Drosophila bipectinata , and Drosophila ananassae.

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma . However, the assembled portion of the Muller F Element (orthologous to the fourth chromosome in Drosophila melanogaster ) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to > 18 Mb. Here, we present chromosome-level long read genome assemblies for four Drosophila species with expanded F Elements ranging in size from 2.3 Mb to 20.5 Mb. Each Muller Element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

Mode and Tempo of 3D Genome Evolution in Drosophila.

Topologically associating domains (TADs) are thought to play an important role in preventing gene misexpression by spatially constraining enhancer-promoter contacts. The deleterious nature of gene misexpression implies that TADs should, therefore, be conserved among related species. Several early studies comparing chromosome conformation between species reported high levels of TAD conservation; however, more recent studies have questioned these results. Furthermore, recent work suggests that TAD reorganization is not associated with extensive changes in gene expression. Here, we investigate the evolutionary conservation of TADs among 11 species of Drosophila. We use Hi-C data to identify TADs in each species and employ a comparative phylogenetic approach to derive empirical estimates of the rate of TAD evolution. Surprisingly, we find that TADs evolve rapidly. However, we also find that the rate of evolution depends on the chromatin state of the TAD, with TADs enriched for developmentally regulated chromatin evolving significantly slower than TADs enriched for broadly expressed, active chromatin. We also find that, after controlling for differences in chromatin state, highly conserved TADs do not exhibit higher levels of gene expression constraint. These results suggest that, in general, most TADs evolve rapidly and their divergence is not associated with widespread changes in gene expression. However, higher levels of evolutionary conservation and gene expression constraints in TADs enriched for developmentally regulated chromatin suggest that these TAD subtypes may be more important for regulating gene expression, likely due to the larger number of long-distance enhancer-promoter contacts associated with developmental genes.

Sex-specific variation in R-loop formation in Drosophila melanogaster.

R-loops are three-stranded nucleotide structures consisting of a DNA:RNA hybrid and a displaced ssDNA non-template strand. Previous work suggests that R-loop formation is primarily determined by the thermodynamics of DNA:RNA binding, which are governed by base composition (e.g., GC skew) and transcription-induced DNA superhelicity. However, R-loops have been described at genomic locations that lack these properties, suggesting that they may serve other context-specific roles. To better understand the genetic determinants of R-loop formation, we have characterized the Drosophila melanogaster R-loop landscape across strains and between sexes using DNA:RNA immunoprecipitation followed by high-throughput sequencing (DRIP-seq). We find that R-loops are associated with sequence motifs that are G-rich or exhibit G/C skew, as well as highly expressed genes, tRNAs, and small nuclear RNAs, consistent with a role for DNA sequence and torsion in R-loop specification. However, we also find motifs associated with R-loops that are A/T-rich and lack G/C skew as well as a subset of R-loops that are enriched in polycomb-repressed chromatin. Differential enrichment analysis reveals a small number of sex-biased R-loops: while non-differentially enriched and male-enriched R-loops form at similar genetic features and chromatin states and contain similar sequence motifs, female-enriched R-loops form at unique genetic features, chromatin states, and sequence motifs and are associated with genes that show ovary-biased expression. Male-enriched R-loops are most abundant on the dosage-compensated X chromosome, where R-loops appear stronger compared to autosomal R-loops. R-loop-containing genes on the X chromosome are dosage-compensated yet show lower MOF binding and reduced H4K16ac compared to R-loop-absent genes, suggesting that H4K16ac or MOF may attenuate R-loop formation. Collectively, these results suggest that R-loop formation in vivo is not fully explained by DNA sequence and topology and raise the possibility that a distinct subset of these hybrid structures plays an important role in the establishment and maintenance of epigenetic differences between sexes.

A transposon expression burst accompanies the activation of Y-chromosome fertility genes during Drosophila spermatogenesis.

Transposable elements (TEs) must replicate in germline cells to pass novel insertions to offspring. In Drosophila melanogaster ovaries, TEs can exploit specific developmental windows of opportunity to evade host silencing and increase their copy numbers. However, TE activity and host silencing in the distinct cell types of Drosophila testis are not well understood. Here, we reanalyze publicly available single-cell RNA-seq datasets to quantify TE expression in the distinct cell types of the Drosophila testis. We develop a method for identification of TE and host gene expression modules and find that a distinct population of early spermatocytes expresses a large number of TEs at much higher levels than other germline and somatic components of the testes. This burst of TE expression coincides with the activation of Y chromosome fertility factors and spermatocyte-specific transcriptional regulators, as well as downregulation of many components of the piRNA pathway. The TEs expressed by this cell population are specifically enriched on the Y chromosome and depleted on the X chromosome, relative to other active TEs. These data suggest that some TEs may achieve high insertional activity in males by exploiting a window of opportunity for mobilization created by the activation of spermatocyte-specific and Y chromosome-specific transcriptional programs.

The nuclear receptor HNF4 drives a brush border gene program conserved across murine intestine, kidney, and embryonic yolk sac.

The brush border is comprised of microvilli surface protrusions on the apical surface of epithelia. This specialized structure greatly increases absorptive surface area and plays crucial roles in human health. However, transcriptional regulatory networks controlling brush border genes are not fully understood. Here, we identify that hepatocyte nuclear factor 4 (HNF4) transcription factor is a conserved and important regulator of brush border gene program in multiple organs, such as intestine, kidney and yolk sac. Compromised brush border gene signatures and impaired transport were observed in these tissues upon HNF4 loss. By ChIP-seq, we find HNF4 binds and activates brush border genes in the intestine and kidney. H3K4me3 HiChIP-seq identifies that HNF4 loss results in impaired chromatin looping between enhancers and promoters at gene loci of brush border genes, and instead enhanced chromatin looping at gene loci of stress fiber genes in the intestine. This study provides comprehensive transcriptional regulatory mechanisms and a functional demonstration of a critical role for HNF4 in brush border gene regulation across multiple murine epithelial tissues.

Three-dimensional interactions between enhancers and promoters during intestinal differentiation depend upon HNF4.

Cells in renewing tissues exhibit dramatic transcriptional changes as they differentiate. The contribution of chromatin looping to tissue renewal is incompletely understood. Enhancer-promoter interactions could be relatively stable as cells transition from progenitor to differentiated states; alternatively, chromatin looping could be as dynamic as the gene expression from their loci. The intestinal epithelium is the most rapidly renewing mammalian tissue. Proliferative cells in crypts of Lieberkühn sustain a stream of differentiated cells that are continually shed into the lumen. We apply chromosome conformation capture combined with chromatin immunoprecipitation (HiChIP) and sequencing to measure enhancer-promoter interactions in progenitor and differentiated cells of the intestinal epithelium. Despite dynamic gene regulation across the differentiation axis, we find that enhancer-promoter interactions are relatively stable. Functionally, we find HNF4 transcription factors are required for chromatin looping at target genes. Depletion of HNF4 disrupts local chromatin looping, histone modifications, and target gene expression. This study provides insights into transcriptional regulatory mechanisms governing homeostasis in renewing tissues.

Telomeric TART elements target the piRNA machinery in Drosophila.

Coevolution between transposable elements (TEs) and their hosts can be antagonistic, where TEs evolve to avoid silencing and the host responds by reestablishing TE suppression, or mutualistic, where TEs are co-opted to benefit their host. The TART-A TE functions as an important component of Drosophila telomeres but has also reportedly inserted into the Drosophila melanogaster nuclear export factor gene nxf2. We find that, rather than inserting into nxf2, TART-A has actually captured a portion of nxf2 sequence. We show that TART-A produces abundant Piwi-interacting small RNAs (piRNAs), some of which are antisense to the nxf2 transcript, and that the TART-like region of nxf2 is evolving rapidly. Furthermore, in D. melanogaster, TART-A is present at higher copy numbers, and nxf2 shows reduced expression, compared to the closely related species Drosophila simulans. We propose that capturing nxf2 sequence allowed TART-A to target the nxf2 gene for piRNA-mediated repression and that these 2 elements are engaged in antagonistic coevolution despite the fact that TART-A is serving a critical role for its host genome.

3D genome evolution and reorganization in the Drosophila melanogaster species group.

Topologically associating domains, or TADs, are functional units that organize chromosomes into 3D structures of interacting chromatin. TADs play an important role in regulating gene expression by constraining enhancer-promoter contacts and there is evidence that deletion of TAD boundaries leads to aberrant expression of neighboring genes. While the mechanisms of TAD formation have been well-studied, current knowledge on the patterns of TAD evolution across species is limited. Due to the integral role TADs play in gene regulation, their structure and organization is expected to be conserved during evolution. However, more recent research suggests that TAD structures diverge relatively rapidly. We use Hi-C chromosome conformation capture to measure evolutionary conservation of whole TADs and TAD boundary elements between D. melanogaster and D. triauraria, two early-branching species from the melanogaster species group which diverged ∼15 million years ago. We find that the majority of TADs have been reorganized since the common ancestor of D. melanogaster and D. triauraria, via a combination of chromosomal rearrangements and gain/loss of TAD boundaries. TAD reorganization between these two species is associated with a localized effect on gene expression, near the site of disruption. By separating TADs into subtypes based on their chromatin state, we find that different subtypes are evolving under different evolutionary forces. TADs enriched for broadly expressed, transcriptionally active genes are evolving rapidly, potentially due to positive selection, whereas TADs enriched for developmentally-regulated genes remain conserved, presumably due to their importance in restricting gene-regulatory element interactions. These results provide novel insight into the evolutionary dynamics of TADs and help to reconcile contradictory reports related to the evolutionary conservation of TADs and whether changes in TAD structure affect gene expression.

Nanopore sequencing and Hi-C scaffolding provide insight into the evolutionary dynamics of transposable elements and piRNA production in wild strains of Drosophila melanogaster.

Illumina sequencing has allowed for population-level surveys of transposable element (TE) polymorphism via split alignment approaches, which has provided important insight into the population dynamics of TEs. However, such approaches are not able to identify insertions of uncharacterized TEs, nor can they assemble the full sequence of inserted elements. Here, we use nanopore sequencing and Hi-C scaffolding to produce de novo genome assemblies for two wild strains of Drosophila melanogaster from the Drosophila Genetic Reference Panel (DGRP). Ovarian piRNA populations and Illumina split-read TE insertion profiles have been previously produced for both strains. We find that nanopore sequencing with Hi-C scaffolding produces highly contiguous, chromosome-length scaffolds, and we identify hundreds of TE insertions that were missed by Illumina-based methods, including a novel micropia-like element that has recently invaded the DGRP population. We also find hundreds of piRNA-producing loci that are specific to each strain. Some of these loci are created by strain-specific TE insertions, while others appear to be epigenetically controlled. Our results suggest that Illumina approaches reveal only a portion of the repetitive sequence landscape of eukaryotic genomes and that population-level resequencing using long reads is likely to provide novel insight into the evolutionary dynamics of repetitive elements.

Adenine Methylation in Drosophila Is Associated with the Tissue-Specific Expression of Developmental and Regulatory Genes.

N6-methyladenine (6mA or m6dA) is a DNA modification that has long been known to play an important role in a variety of biological functions in prokaryotes. This modification has only recently been described in eukaryotes, where it seems to have evolved species-specific functions ranging from nucleosome positioning to transposon repression. In Drosophila, 6mA has been shown to be important for enforcing the tissue specificity of neuronal genes in the brain and suppressing transposable element expression in the ovaries. In this study, we have analyzed the raw signal data from nanopore sequencing to identify 6mA positions in the D. melanogaster genome at single-base resolution. We find that this modification is enriched upstream from transcription start sites, within the introns and 3' UTRs of genes, as well as in simple repeats. These 6mA positions are enriched for sequence motifs that are recognized by known transcriptional activators involved in development, such as Bicoid and Caudal, and the genes that carry this modification are enriched for functions involved in development, regulation of transcription, and neuronal activity. These genes show high expression specificity in a variety of tissues besides the brain, suggesting that this modification may play a more general role in enforcing the specificity of gene expression across many tissues, throughout development, and between the sexes.

The SR protein B52/SRp55 regulates splicing of the period thermosensitive intron and mid-day siesta in Drosophila.

Similar to many diurnal animals, Drosophila melanogaster exhibits a mid-day siesta that is more robust as temperature increases, an adaptive response that aims to minimize the deleterious effects from exposure to heat. This temperature-dependent plasticity in mid-day sleep levels is partly based on the thermal sensitive splicing of an intron in the 3' untranslated region (UTR) of the circadian clock gene termed period (per). In this study, we evaluated a possible role for the serine/arginine-rich (SR) splicing factors in the regulation of dmpi8 splicing efficiency and mid-day siesta. Using a Drosophila cell culture assay we show that B52/SRp55 increases dmpi8 splicing efficiency, whereas other SR proteins have little to no effect. The magnitude of the stimulatory effect of B52 on dmpi8 splicing efficiency is modulated by natural variation in single nucleotide polymorphisms (SNPs) in the per 3' UTR that correlate with B52 binding levels. Down-regulating B52 expression in clock neurons increases mid-day siesta and reduces dmpi8 splicing efficiency. Our results establish a novel role for SR proteins in sleep and suggest that polymorphisms in the per 3' UTR contribute to natural variation in sleep behavior by modulating the binding efficiencies of SR proteins.

Mid-day siesta in natural populations of D. melanogaster from Africa exhibits an altitudinal cline and is regulated by splicing of a thermosensitive intron in the period clock gene.

BACKGROUND: Many diurnal animals exhibit a mid-day 'siesta', generally thought to be an adaptive response aimed at minimizing exposure to heat on warm days, suggesting that in regions with cooler climates mid-day siestas might be a less prominent feature of animal behavior. Drosophila melanogaster exhibits thermal plasticity in its mid-day siesta that is partly governed by the thermosensitive splicing of the 3'-terminal intron (termed dmpi8) from the key circadian clock gene period (per). For example, decreases in temperature lead to progressively more efficient splicing, which increasingly favors activity over sleep during the mid-day. In this study we sought to determine if the adaptation of D. melanogaster from its ancestral range in the lowlands of tropical Africa to the cooler temperatures found at high altitudes involved changes in mid-day sleep behavior and/or dmpi8 splicing efficiency. RESULTS: Using natural populations of Drosophila melanogaster from different altitudes in tropical Africa we show that flies from high elevations have a reduced mid-day siesta and less consolidated sleep. We identified a single nucleotide polymorphism (SNP) in the per 3' untranslated region that has strong effects on dmpi8 splicing and mid-day sleep levels in both low and high altitude flies. Intriguingly, high altitude flies with a particular variant of this SNP exhibit increased dmpi8 splicing efficiency compared to their low altitude counterparts, consistent with reduced mid-day siesta. Thus, a boost in dmpi8 splicing efficiency appears to have played a prominent but not universal role in how African flies adapted to the cooler temperatures at high altitude. CONCLUSIONS: Our findings point towards mid-day sleep behavior as a key evolutionary target in the thermal adaptation of animals, and provide a genetic framework for investigating daytime sleep in diurnal animals which appears to be driven by mechanisms distinct from those underlying nighttime sleep.

A novel pathway for sensory-mediated arousal involves splicing of an intron in the period clock gene.

STUDY OBJECTIVES: D. melanogaster is an excellent animal model to study how the circadian (≅24-h) timing system and sleep regulate daily wake-sleep cycles. Splicing of a temperature-sensitive 3'-terminal intron (termed dmpi8) from the circadian clock gene period (per) regulates the distribution of daily activity in Drosophila. The role of dmpi8 splicing on daily behavior was further evaluated by analyzing sleep. DESIGN: Transgenic flies of the same genetic background but expressing either a wild-type recombinant per gene or one where the efficiency of dmpi8 splicing was increased were exposed to different temperatures in daily light-dark cycles and sleep parameters measured. In addition, transgenic flies were briefly exposed to a variety of sensory-mediated stimuli to measure arousal responses. RESULTS: Surprisingly, we show that the effect of dmpi8 splicing on daytime activity levels does not involve a circadian role for per but is linked to adjustments in sensory-dependent arousal and sleep behavior. Genetically altered flies with high dmpi8 splicing efficiency remain aroused longer following short treatments with light and non-photic cues such as mechanical stimulation. CONCLUSIONS: We propose that the thermal regulation of dmpi8 splicing acts as a temperature-calibrated rheostat in a novel arousal mechanism, so that on warm days the inefficient splicing of the dmpi8 intron triggers an increase in quiescence by decreasing sensory-mediated arousal, thus ensuring flies minimize being active during the hot midday sun despite the presence of light in the environment, which is usually a strong arousal cue for diurnal animals.

FKBP immunophilins and Alzheimer's disease: a chaperoned affair.

The FK506-binding protein (FKBP) family of immunophilins consists of proteins with a variety of protein-protein interaction domains and versatile cellular functions. Analysis of the functions of immunophilins has been the focus of studies in recent years and has led to the identification of various molecular pathways in which FKBPs play an active role. All FKBPs contain a domain with prolyl cis/trans isomerase (PPIase) activity. Binding of the immunosuppressant molecule FK506 to this domain inhibits their PPIase activity while mediating immune suppression through inhibition of calcineurin. The larger members, FKBP51 and FKBP52, interact with Hsp90 and exhibit chaperone activity that is shown to regulate steroid hormone signalling. From these studies it is clear that FKBP proteins are expressed ubiquitously but show relatively high levels of expression in the nervous system. Consistent with this expression, FKBPs have been implicated with both neuroprotection and neurodegeneration. This review will focus on recent studies involving FKBP immunophilins in Alzheimer's-disease-related pathways.

Control of Alzheimer's amyloid beta toxicity by the high molecular weight immunophilin FKBP52 and copper homeostasis in Drosophila.

FK506 binding proteins (FKBPs), also called immunophilins, are prolyl-isomerases (PPIases) that participate in a wide variety of cellular functions including hormone signaling and protein folding. Recent studies indicate that proteins that contain PPIase activity can also alter the processing of Alzheimer's Amyloid Precursor Protein (APP). Originally identified in hematopoietic cells, FKBP52 is much more abundantly expressed in neurons, including the hippocampus, frontal cortex, and basal ganglia. Given the fact that the high molecular weight immunophilin FKBP52 is highly expressed in CNS regions susceptible to Alzheimer's, we investigated its role in Abeta toxicity. Towards this goal, we generated Abeta transgenic Drosophila that harbor gain of function or loss of function mutations of FKBP52. FKBP52 overexpression reduced the toxicity of Abeta and increased lifespan in Abeta flies, whereas loss of function of FKBP52 exacerbated these Abeta phenotypes. Interestingly, the Abeta pathology was enhanced by mutations in the copper transporters Atox1, which interacts with FKBP52, and Ctr1A and was suppressed in FKBP52 mutant flies raised on a copper chelator diet. Using mammalian cultures, we show that FKBP52 (-/-) cells have increased intracellular copper and higher levels of Abeta. This effect is reversed by reconstitution of FKBP52. Finally, we also found that FKBP52 formed stable complexes with APP through its FK506 interacting domain. Taken together, these studies identify a novel role for FKBP52 in modulating toxicity of Abeta peptides.

Identification of novel genes that modify phenotypes induced by Alzheimer's beta-amyloid overexpression in Drosophila.

Sustained increases in life expectancy have underscored the importance of managing diseases with a high incidence in late life, such as various neurodegenerative conditions. Alzheimer's disease (AD) is the most common among these, and consequently significant research effort is spent on studying it. Although a lot is known about the pathology of AD and the role of beta-amyloid (Abeta) peptides, the complete network of interactions regulating Abeta metabolism and toxicity still eludes us. To address this, we have conducted genetic interaction screens using transgenic Drosophila expressing Abeta and we have identified mutations that affect Abeta metabolism and toxicity. These analyses highlight the involvement of various biochemical processes such as secretion, cholesterol homeostasis, and regulation of chromatin structure and function, among others, in mediating toxic Abeta effects. Several of the mutations that we identified have not been linked to Abeta toxicity before and thus constitute novel potential targets for AD intervention. We additionally tested these mutations for interactions with tau and expanded-polyglutamine overexpression and found a few candidate mutations that may mediate common mechanisms of neurodegeneration. Our data offer insight into the toxicity of Abeta and open new areas for further study into AD pathogenesis.

Replication and inheritance of Nocardia plasmid pC1.

Nocardia sp. C-14-1, isolated from acrylic fiber wastewater, can degrade long-chain alkanes and succinonitrile efficiently. Here we report the characterization of an indigenous plasmid pC1. The overall nucleotide sequence of pC1 consisted of 5841 bp. The five ORFs, encoding a DNA recombinase, replication protein (Rep(pC1)) and three proteins of unknown function, were predicted on pC1. The Rep(pC1) displayed its homology with the Rep of Rhodococcus large plasmid p33701, suggesting a theta type of replication. An Escherichia coli plasmid (containing the single rep(pC1) gene) propagated autonomously in low copy number in Nocardia or Rhodococcus, suggesting that rep(pC1) was an essential gene for plasmid replication. The plasmid (containing the single rep(pC1) gene) presented as inheritance unstable hints that other pC1 loci were required for the stable inheritance of plasmids. By comparison of the plasmid-borne Rep proteins, we classify Rhodococcus or Nocardia plasmids into four groups.

[Screening of high efficiency succinonitrile degrading bacterium strains and their degradability].

Succinonitrile used as a sole source of carbon and nitrogen, two bacterium strains, named as J-1-3 and J-13-1, which had high degrading capacity for succinonitrile were isolated and screened out from acrylic fiber wastewater and biofilm in its treatment structure. Through morphological and biochemical method, the two strains were primarily identified as Pseudomonas spp. By tests in shaking flasks, it was determined that the strains can be optimum growth at 30 degrees C, with shaker rotary speed which indirectly reflected aeration capacity at 250 r/min, inoculum amount of 0.1%, and initial pH6. On the optimum conditions for growth, the degradation experiments on different initial concentrations of succinonitrile were carried out. The results indicate that the two strains, especially J-13-1 had high degrading efficiency for succinonitrile. With the initial concentration of succinonitrile at ca. 6000, 8000 and 10,000 mg/L, the degrading rates of succinonitrile by strain J-13-1 reached to 100% after 12.5 h, 14 h and 16 h, respectively.