Detection of ac4C in human mRNA is preserved upon data reassessment

We recently reported the distribution of N4-acetylcytidine (ac4C) in HeLa mRNA at base resolution through chemical reduction and the induction of C:T mismatches in sequencing (RedaC:T-seq). Our results contradicted an earlier report from Schwartz and colleagues utilizing a similar method termed ac4C-seq. Here, we revisit both datasets and reaffirm our findings. Through RedaC:T-seq reanalysis, we establish a low basal error rate at unmodified nucleotides that is not skewed to any specific mismatch type and a prominent increase in C:T substitutions as the dominant mismatch type in both treated wild-type replicates, with a high degree of reproducibility across replicates. In contrast, through ac4C-seq reanalysis, we uncover significant data quality issues including insufficient depth, with one wild-type replicate yielding 2.7 million reads, inconsistencies in reduction efficiencies between replicates, and an overall increase in mismatches involving thymine that could obscure ac4C detection. These analyses bolster the detection of ac4C in HeLa mRNA through RedaC:T-seq.

PancrESS - a meta-analysis resource for understanding cell-type specific expression in the human pancreas

BACKGROUND

The human pancreas is composed of specialized cell types producing hormones and enzymes critical to human health. These specialized functions are the result of cell type-specific transcriptional programs which manifest in cell-specific gene expression. Understanding these programs is essential to developing therapies for pancreatic disorders. Transcription in the human pancreas has been widely studied by single-cell RNA technologies, however the diversity of protocols and analysis methods hinders their interpretability in the aggregate.

RESULTS

In this work, we perform a meta-analysis of pancreatic single-cell RNA sequencing data. We present a database for reference transcriptome abundances and cell-type specificity metrics. This database facilitates the identification and definition of marker genes within the pancreas. Additionally, we introduce a versatile tool which is freely available as an R package, and should permit integration into existing workflows. Our tool accepts count data files generated by widely-used single-cell gene expression platforms in their original format, eliminating an additional pre-formatting step. Although we designed it to calculate expression specificity of pancreas cell types, our tool is agnostic to the biological source of count data, extending its applicability to other biological systems.

CONCLUSIONS

Our findings enhance the current understanding of expression specificity within the pancreas, surpassing previous work in terms of scope and detail. Furthermore, our database and tool enable researchers to perform similar calculations in diverse biological systems, expanding the applicability of marker gene identification and facilitating comparative analyses.

Protocol for base resolution mapping of ac4C using RedaC:T-seq

N4-acetylcytidine (ac4C) is an mRNA modification catalyzed by the enzyme N-acetyltransferase 10 (NAT10), with position-dependent effects on mRNA translation. This protocol details a procedure to map ac4C at base resolution using NaBH₄-induced reduction of ac4C and conversion to thymidine followed by sequencing (RedaC:T-seq). Total RNA is ribodepleted and then treated with NaBH₄ to reduce ac4C to tetrahydro-ac4C, which specifically alters base pairing during cDNA synthesis, allowing the detection of ac4C at positions called as thymidine following Illumina sequencing. For complete details on the use and execution of this protocol, please refer to Arango et al. (2022).¹.

Direct epitranscriptomic regulation of mammalian translation initiation through N4-acetylcytidine

mRNA function is influenced by modifications that modulate canonical nucleobase behavior. We show that a single modification mediates distinct impacts on mRNA translation in a position-dependent manner. Although cytidine acetylation (ac4C) within protein-coding sequences stimulates translation, ac4C within 5' UTRs impacts protein synthesis at the level of initiation. 5' UTR acetylation promotes initiation at upstream sequences, competitively inhibiting annotated start codons. Acetylation further directly impedes initiation at optimal AUG contexts: ac4C within AUG-flanking Kozak sequences reduced initiation in base-resolved transcriptome-wide HeLa results and in vitro utilizing substrates with site-specific ac4C incorporation. Cryo-EM of mammalian 80S initiation complexes revealed that ac4C in the -1 position adjacent to an AUG start codon disrupts an interaction between C and hypermodified t6A at nucleotide 37 of the initiator tRNA. These findings demonstrate the impact of RNA modifications on nucleobase function at a molecular level and introduce mRNA acetylation as a factor regulating translation in a location-specific manner.

Oncogenic lncRNAs alter epigenetic memory at a fragile chromosomal site in human cancer cells

Chromosome instability is a critical event in cancer progression. Histone H3 variant CENP-A plays a fundamental role in defining centromere identity, structure, and function but is innately overexpressed in several types of solid cancers. In the cancer background, excess CENP-A is deposited ectopically on chromosome arms, including 8q24/cMYC locus, by invading transcription-coupled H3.3 chaperone pathways. Up-regulation of lncRNAs in many cancers correlates with poor prognosis and recurrence in patients. We report that transcription of 8q24-derived oncogenic lncRNAs plays an unanticipated role in altering the 8q24 chromatin landscape by H3.3 chaperone-mediated deposition of CENP-A-associated complexes. Furthermore, a transgene cassette carrying specific 8q24-derived lncRNA integrated into a naïve chromosome locus recruits CENP-A to the new location in a cis-acting manner. These data provide a plausible mechanistic link between locus-specific oncogenic lncRNAs, aberrant local chromatin structure, and the generation of new epigenetic memory at a fragile site in human cancer cells.

The X chromosome from telomere to telomere: key achievements and future opportunities

While the human genome represents the most accurate vertebrate reference assembly to date, it still contains numerous gaps, including centromeric and other large repeat-containing regions - often termed the "dark side" of the genome - many of which are of fundamental biological importance. Miga et al.^1,2 present the first gapless assembly of the human X chromosome, with the help of ultra-long-read nanopore reads generated for the haploid complete hydatidiform mole (CHM13) genome. They reconstruct the ~3.1 megabase centromeric satellite DNA array and map DNA methylation patterns across complex tandem repeats and satellite arrays. This Telomere-to-Telomere assembly provides a superior human X chromosome reference enabling future sex-determination and X-linked disease research, and provides a path towards finishing the entire human genome sequence.

Activation of Crtc2/Creb1 in skeletal muscle enhances weight loss during intermittent fasting

The Creb-Regulated Transcriptional Coactivator (Crtc) family of transcriptional coregulators drive Creb1-mediated transcription effects on metabolism in many tissues, but the in vivo effects of Crtc2/Creb1 transcription on skeletal muscle metabolism are not known. Skeletal muscle-specific overexpression of Crtc2 (Crtc2 mice) induced greater mitochondrial activity, metabolic flux capacity for both carbohydrates and fats, improved glucose tolerance and insulin sensitivity, and increased oxidative capacity, supported by upregulation of key metabolic genes. Crtc2 overexpression led to greater weight loss during alternate day fasting (ADF), selective loss of fat rather than lean mass, maintenance of higher energy expenditure during the fast and reduced binge-eating during the feeding period. ADF downregulated most of the mitochondrial electron transport genes, and other regulators of mitochondrial function, that were substantially reversed by Crtc2-driven transcription. Glucocorticoids acted with AMPK to drive atrophy and mitophagy, which was reversed by Crtc2/Creb1 signaling. Crtc2/Creb1-mediated signaling coordinates metabolic adaptations in skeletal muscle that explain how Crtc2/Creb1 regulates metabolism and weight loss.

Immunoprecipitation and Sequencing of Acetylated RNA

Generation of the epitranscriptome through chemical modifications of protein-coding messenger RNAs (mRNAs) has emerged as a new mechanism of post-transcriptional gene regulation. While most mRNA modifications are methylation events, a single acetylated ribonucleoside has been described in eukaryotes, occurring at the N4-position of cytidine (N4-acetylcytidine or ac4C). Using a combination of antibody-based enrichment of acetylated regions and deep sequencing, we recently reported ac4C as a novel mRNA modification that is catalyzed by the N-acetyltransferase enzyme NAT10. In this protocol, we describe in detail the procedures to identify acetylated mRNA regions transcriptome-wide using acetylated RNA immunoprecipitation and sequencing (acRIP-seq).

Acetylation of Cytidine in mRNA Promotes Translation Efficiency

Generation of the "epitranscriptome" through post-transcriptional ribonucleoside modification embeds a layer of regulatory complexity into RNA structure and function. Here, we describe N4-acetylcytidine (ac4C) as an mRNA modification that is catalyzed by the acetyltransferase NAT10. Transcriptome-wide mapping of ac4C revealed discretely acetylated regions that were enriched within coding sequences. Ablation of NAT10 reduced ac4C detection at the mapped mRNA sites and was globally associated with target mRNA downregulation. Analysis of mRNA half-lives revealed a NAT10-dependent increase in stability in the cohort of acetylated mRNAs. mRNA acetylation was further demonstrated to enhance substrate translation in vitro and in vivo. Codon content analysis within ac4C peaks uncovered a biased representation of cytidine within wobble sites that was empirically determined to influence mRNA decoding efficiency. These findings expand the repertoire of mRNA modifications to include an acetylated residue and establish a role for ac4C in the regulation of mRNA translation.

HJURP antagonizes CENP-A mislocalization driven by the H3.3 chaperones HIRA and DAXX

The centromere specific histone H3 variant CENP-A/CENH3 specifies where the kinetochore is formed in most eukaryotes. Despite tight regulation of CENP-A levels in normal cells, overexpression of CENP-A is a feature shared by various types of solid tumors and results in its mislocalization to non-centromeric DNA. How CENP-A is assembled ectopically and the consequences of this mislocalization remain topics of high interest. Here, we report that in human colon cancer cells, the H3.3 chaperones HIRA and DAXX promote ectopic CENP-A deposition. Moreover, the correct balance between levels of the centromeric chaperone HJURP and CENP-A is essential to preclude ectopic assembly by H3.3 chaperones. In addition, we find that ectopic localization can recruit kinetochore components, and correlates with mitotic defects and DNA damage in G1 phase. Finally, CENP-A occupancy at the 8q24 locus is also correlated with amplification and overexpression of the MYC gene within that locus. Overall, these data provide insights into the causes and consequences of histone variant mislocalization in human cancer cells.

Independence between pre-mRNA splicing and DNA methylation in an isogenic minigene resource

Actively transcribed genes adopt a unique chromatin environment with characteristic patterns of enrichment. Within gene bodies, H3K36me3 and cytosine DNA methylation are elevated at exons of spliced genes and have been implicated in the regulation of pre-mRNA splicing. H3K36me3 is further responsive to splicing, wherein splicing inhibition led to a redistribution and general reduction over gene bodies. In contrast, little is known of the mechanisms supporting elevated DNA methylation at actively spliced genic locations. Recent evidence associating the de novo DNA methyltransferase Dnmt3b with H3K36me3-rich chromatin raises the possibility that genic DNA methylation is influenced by splicing-associated H3K36me3. Here, we report the generation of an isogenic resource to test the direct impact of splicing on chromatin. A panel of minigenes of varying splicing potential were integrated into a single FRT site for inducible expression. Profiling of H3K36me3 confirmed the established relationship to splicing, wherein levels were directly correlated with splicing efficiency. In contrast, DNA methylation was equivalently detected across the minigene panel, irrespective of splicing and H3K36me3 status. In addition to revealing a degree of independence between genic H3K36me3 and DNA methylation, these findings highlight the generated minigene panel as a flexible platform for the query of splicing-dependent chromatin modifications.

Replication Stress Shapes a Protective Chromatin Environment across Fragile Genomic Regions

Recent integrative epigenome analyses highlight the importance of functionally distinct chromatin states for accurate cell function. How these states are established and maintained is a matter of intense investigation. Here, we present evidence for DNA damage as an unexpected means to shape a protective chromatin environment at regions of recurrent replication stress (RS). Upon aberrant fork stalling, DNA damage signaling and concomitant H2AX phosphorylation coordinate the FACT-dependent deposition of macroH2A1.2, a histone variant that promotes DNA repair by homologous recombination (HR). MacroH2A1.2, in turn, facilitates the accumulation of the tumor suppressor and HR effector BRCA1 at replication forks to protect from RS-induced DNA damage. Consequently, replicating primary cells steadily accrue macroH2A1.2 at fragile regions, whereas macroH2A1.2 loss in these cells triggers DNA damage signaling-dependent senescence, a hallmark of RS. Altogether, our findings demonstrate that recurrent DNA damage contributes to the chromatin landscape to ensure the epigenomic integrity of dividing cells.

Cajal body function in genome organization and transcriptome diversity

Nuclear bodies contribute to non-random organization of the human genome and nuclear function. Using a major prototypical nuclear body, the Cajal body, as an example, we suggest that these structures assemble at specific gene loci located across the genome as a result of high transcriptional activity. Subsequently, target genes are physically clustered in close proximity in Cajal body-containing cells. However, Cajal bodies are observed in only a limited number of human cell types, including neuronal and cancer cells. Ultimately, Cajal body depletion perturbs splicing kinetics by reducing target small nuclear RNA (snRNA) transcription and limiting the levels of spliceosomal snRNPs, including their modification and turnover following each round of RNA splicing. As such, Cajal bodies are capable of shaping the chromatin interaction landscape and the transcriptome by influencing spliceosome kinetics. Future studies should concentrate on characterizing the direct influence of Cajal bodies upon snRNA gene transcriptional dynamics. Also see the video abstract here.

Linking Genes and Brain Development of Honeybee Workers: A Whole-Transcriptome Approach

Honeybees live in complex societies whose capabilities far exceed those of the sum of their single members. This social synergism is achieved mainly by the worker bees, which form a female caste. The worker bees display diverse collaborative behaviors and engage in different behavioral tasks, which are controlled by the central nervous system (CNS). The development of the worker brain is determined by the female sex and the worker caste determination signal. Here, we report on genes that are controlled by sex or by caste during differentiation of the worker's pupal brain. We sequenced and compared transcriptomes from the pupal brains of honeybee workers, queens and drones. We detected 333 genes that are differently expressed and 519 genes that are differentially spliced between the sexes, and 1760 genes that are differentially expressed and 692 genes that are differentially spliced between castes. We further found that 403 genes are differentially regulated by both the sex and caste signals, providing evidence of the integration of both signals through differential gene regulation. In this gene set, we found that the molecular processes of restructuring the cell shape and cell-to-cell signaling are overrepresented. Our approach identified candidate genes that may be involved in brain differentiation that ensures the various social worker behaviors.

Identifying Centromeric RNAs Involved in Histone Dynamics In Vivo

Over the last decade, the long accepted dogma that heterochromatin is silent has been challenged by increasing evidence of active transcription in these apocryphally annotated quiescent regions of the genome. The recent discovery of noncoding RNAs (ncRNAs) originating from, or localizing to, centromeres, pericentromeres, and telomeres (ie, constitutive heterochromatin) suggest a potential role for ncRNAs in genome integrity. This new paradigm suggests that ncRNAs may recruit chromatin-binding factors, stabilize the higher order folded state of the chromatin fiber, and participate in regulation of processes such as transcription-mediated nucleosome assembly. Thus, identifying, purifying, and elucidating the function of ncRNAs has the potential to provide key insights into genome organization and is currently a topic of intense experimental investigation.

Cajal bodies are linked to genome conformation

The mechanisms underlying nuclear body (NB) formation and their contribution to genome function are unknown. Here we examined the non-random positioning of Cajal bodies (CBs), major NBs involved in spliceosomal snRNP assembly and their role in genome organization. CBs are predominantly located at the periphery of chromosome territories at a multi-chromosome interface. Genome-wide chromosome conformation capture analysis (4C-seq) using CB-interacting loci revealed that CB-associated regions are enriched with highly expressed histone genes and U small nuclear or nucleolar RNA (sn/snoRNA) loci that form intra- and inter-chromosomal clusters. In particular, we observed a number of CB-dependent gene-positioning events on chromosome 1. RNAi-mediated disassembly of CBs disrupts the CB-targeting gene clusters and suppresses the expression of U sn/snoRNA and histone genes. This loss of spliceosomal snRNP production results in increased splicing noise, even in CB-distal regions. Therefore, we conclude that CBs contribute to genome organization with global effects on gene expression and RNA splicing fidelity.

Nuclear bodies can nucleate at sites of active transcription and are beneficial for efficient gene expression. Here, the authors show that Cajal bodies, a prominent type of nuclear body, contribute to genome organization with global effects on gene expression and RNA splicing fidelity.

TET-catalyzed oxidation of intragenic 5-methylcytosine regulates CTCF-dependent alternative splicing

Intragenic 5-methylcytosine and CTCF mediate opposing effects on pre-mRNA splicing: CTCF promotes inclusion of weak upstream exons through RNA polymerase II pausing, whereas 5-methylcytosine evicts CTCF, leading to exon exclusion. However, the mechanisms governing dynamic DNA methylation at CTCF-binding sites were unclear. Here, we reveal the methylcytosine dioxygenases TET1 and TET2 as active regulators of CTCF-mediated alternative splicing through conversion of 5-methylcytosine to its oxidation derivatives. 5-hydroxymethylcytosine and 5-carboxylcytosine are enriched at an intragenic CTCF-binding sites in the CD45 model gene and are associated with alternative exon inclusion. Reduced TET levels culminate in increased 5-methylcytosine, resulting in CTCF eviction and exon exclusion. In vitro analyses establish the oxidation derivatives are not sufficient to stimulate splicing, but efficiently promote CTCF association. We further show genomewide that reciprocal exchange of 5-hydroxymethylcytosine and 5-methylcytosine at downstream CTCF-binding sites is a general feature of alternative splicing in naïve and activated CD4(+) T cells. These findings significantly expand our current concept of the pre-mRNA "splicing code" to include dynamic intragenic DNA methylation catalyzed by the TET proteins.

Controlled DNA double-strand break induction in mice reveals post-damage transcriptome stability

DNA double-strand breaks (DSBs) and their repair can cause extensive epigenetic changes. As a result, DSBs have been proposed to promote transcriptional and, ultimately, physiological dysfunction via both cell-intrinsic and cell-non-autonomous pathways. Studying the consequences of DSBs in higher organisms has, however, been hindered by a scarcity of tools for controlled DSB induction. Here, we describe a mouse model that allows for both tissue-specific and temporally controlled DSB formation at ∼140 defined genomic loci. Using this model, we show that DSBs promote a DNA damage signaling-dependent decrease in gene expression in primary cells specifically at break-bearing genes, which is reversed upon DSB repair. Importantly, we demonstrate that restoration of gene expression can occur independently of cell cycle progression, underlining its relevance for normal tissue maintenance. Consistent with this, we observe no evidence for persistent transcriptional repression in response to a multi-day course of continuous DSB formation and repair in mouse lymphocytes in vivo Together, our findings reveal an unexpected capacity of primary cells to maintain transcriptome integrity in response to DSBs, pointing to a limited role for DNA damage as a mediator of cell-autonomous epigenetic dysfunction.

Comparative validation of the D. melanogaster modENCODE transcriptome annotation

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community.

Diversity and dynamics of the Drosophila transcriptome

Animal transcriptomes are dynamic, with each cell type, tissue and organ system expressing an ensemble of transcript isoforms that give rise to substantial diversity. Here we have identified new genes, transcripts and proteins using poly(A)+ RNA sequencing from Drosophila melanogaster in cultured cell lines, dissected organ systems and under environmental perturbations. We found that a small set of mostly neural-specific genes has the potential to encode thousands of transcripts each through extensive alternative promoter usage and RNA splicing. The magnitudes of splicing changes are larger between tissues than between developmental stages, and most sex-specific splicing is gonad-specific. Gonads express hundreds of previously unknown coding and long non-coding RNAs (lncRNAs), some of which are antisense to protein-coding genes and produce short regulatory RNAs. Furthermore, previously identified pervasive intergenic transcription occurs primarily within newly identified introns. The fly transcriptome is substantially more complex than previously recognized, with this complexity arising from combinatorial usage of promoters, splice sites and polyadenylation sites.

Design of RNA splicing analysis null models for post hoc filtering of Drosophila head RNA-Seq data with the splicing analysis kit (Spanki)

BACKGROUND

The production of multiple transcript isoforms from one gene is a major source of transcriptome complexity. RNA-Seq experiments, in which transcripts are converted to cDNA and sequenced, allow the resolution and quantification of alternative transcript isoforms. However, methods to analyze splicing are underdeveloped and errors resulting in incorrect splicing calls occur in every experiment.

RESULTS

We used RNA-Seq data to develop sequencing and aligner error models. By applying these error models to known input from simulations, we found that errors result from false alignment to minor splice motifs and antisense stands, shifted junction positions, paralog joining, and repeat induced gaps. By using a series of quantitative and qualitative filters, we eliminated diagnosed errors in the simulation, and applied this to RNA-Seq data from Drosophila melanogaster heads. We used high-confidence junction detections to specifically interrogate local splicing differences between transcripts. This method out-performed commonly used RNA-seq methods to identify known alternative splicing events in the Drosophila sex determination pathway. We describe a flexible software package to perform these tasks called Splicing Analysis Kit (Spanki), available at http://www.cbcb.umd.edu/software/spanki.

CONCLUSIONS

Splice-junction centric analysis of RNA-Seq data provides advantages in specificity for detection of alternative splicing. Our software provides tools to better understand error profiles in RNA-Seq data and improve inference from this new technology. The splice-junction centric approach that this software enables will provide more accurate estimates of differentially regulated splicing than current tools.

Evidence for compensatory upregulation of expressed X-linked genes in mammals, Caenorhabditis elegans and Drosophila melanogaster

Many animal species use a chromosome-based mechanism of sex determination, which has led to the coordinate evolution of dosage-compensation systems. Dosage compensation not only corrects the imbalance in the number of X chromosomes between the sexes but also is hypothesized to correct dosage imbalance within cells that is due to monoallelic X-linked expression and biallelic autosomal expression, by upregulating X-linked genes twofold (termed 'Ohno's hypothesis'). Although this hypothesis is well supported by expression analyses of individual X-linked genes and by microarray-based transcriptome analyses, it was challenged by a recent study using RNA sequencing and proteomics. We obtained new, independent RNA-seq data, measured RNA polymerase distribution and reanalyzed published expression data in mammals, C. elegans and Drosophila. Our analyses, which take into account the skewed gene content of the X chromosome, support the hypothesis of upregulation of expressed X-linked genes to balance expression of the genome.

Identification of functional elements and regulatory circuits by Drosophila modENCODE

To gain insight into how genomic information is translated into cellular and developmental programs, the Drosophila model organism Encyclopedia of DNA Elements (modENCODE) project is comprehensively mapping transcripts, histone modifications, chromosomal proteins, transcription factors, replication proteins and intermediates, and nucleosome properties across a developmental time course and in multiple cell lines. We have generated more than 700 data sets and discovered protein-coding, noncoding, RNA regulatory, replication, and chromatin elements, more than tripling the annotated portion of the Drosophila genome. Correlated activity patterns of these elements reveal a functional regulatory network, which predicts putative new functions for genes, reveals stage- and tissue-specific regulators, and enables gene-expression prediction. Our results provide a foundation for directed experimental and computational studies in Drosophila and related species and also a model for systematic data integration toward comprehensive genomic and functional annotation.

Germline-dependent gene expression in distant non-gonadal somatic tissues of Drosophila

Background

Drosophila females commit tremendous resources to egg production and males produce some of the longest sperm in the animal kingdom. We know little about the coordinated regulation of gene expression patterns in distant somatic tissues that support the developmental cost of gamete production.

Results

We determined the non-gonadal gene expression patterns of Drosophila females and males with or without a germline. Our results show that germline-dependent expression in the non-gonadal soma is extensive. Interestingly, gene expression patterns and hormone titers are consistent with a hormone axis between the gonads and non-gonadal soma.

Conclusions

The germline has a long-range influence on gene expression in the Drosophila sexes. We suggest that this is the result of a germline/soma hormonal axis.

Constraint and turnover in sex-biased gene expression in the genus Drosophila

Both genome content and deployment contribute to phenotypic differences between species. Sex is the most important difference between individuals in a species and has long been posited to be rapidly evolving. Indeed, in the Drosophila genus, traits such as sperm length, genitalia, and gonad size are the most obvious differences between species. Comparative analysis of sex-biased expression should deepen our understanding of the relationship between genome content and deployment during evolution. Using existing and newly assembled genomes, we designed species-specific microarrays to examine sex-biased expression of orthologues and species-restricted genes in D. melanogaster, D. simulans, D. yakuba, D. ananassae, D. pseudoobscura, D. virilis and D. mojavensis. We show that averaged sex-biased expression changes accumulate monotonically over time within the genus. However, different genes contribute to expression variance within species groups compared to between groups. We observed greater turnover of species-restricted genes with male-biased expression, indicating that gene formation and extinction may play a significant part in species differences. Genes with male-biased expression also show the greatest expression and DNA sequence divergence. This higher divergence and turnover of genes with male-biased expression may be due to high transcription rates in the male germline, greater functional pleiotropy of genes expressed in females, and/or sexual competition.

Demasculinization of X chromosomes in the Drosophila genus

X chromosomes evolve differently from autosomes, but general governing principles have not emerged. For example, genes with male-biased expression are under-represented on the X chromosome of D. melanogaster, but are randomly distributed in the genome of Anopheles gambiae. In direct global profiling experiments using species-specific microarrays, we find a nearly identical paucity of genes with male-biased expression on D. melanogaster, D. simulans, D. yakuba, D. ananassae, D. virilis and D. mojavensis X chromosomes. We observe the same under-representation on the neo-X of D. pseudoobscura. It has been suggested that precocious meiotic silencing of the X chromosome accounts for reduced X chromosome male-biased expression in nematodes, mammals and Drosophila. We show that X chromosome genes with male-biased expression are under-represented in somatic cells and in mitotic male germ cells. These data are incompatible with simple X chromosome inactivation models. Using expression profiling and comparative sequence analysis, we show that selective gene extinction on the X chromosome, creation of new genes on autosomes and changed genomic location of existing genes contribute to the unusual X chromosome gene content.

Comparative genomics of Drosophila and human core promoters

Comparison of DNA sequence distributions in Drosophila and human promoters suggests that different motifs have distinct functional roles.

Background

The core promoter region plays a critical role in the regulation of eukaryotic gene expression. We have determined the non-random distribution of DNA sequences relative to the transcriptional start site in Drosophila melanogaster promoters to identify sequences that may be biologically significant. We compare these results with those obtained for human promoters.

Results

We determined the distribution of all 65,536 octamer (8-mers) DNA sequences in 10,914 Drosophila promoters and two sets of human promoters aligned relative to the transcriptional start site. In Drosophila, 298 8-mers have highly significant (p ≤ 1 × 10^-16) non-random distributions peaking within 100 base-pairs of the transcriptional start site. These sequences were grouped into 15 DNA motifs. Ten motifs, termed directional motifs, occur only on the positive strand while the remaining five motifs, termed non-directional motifs, occur on both strands. The only directional motifs to localize in human promoters are TATA, INR, and DPE. The directional motifs were further subdivided into those precisely positioned relative to the transcriptional start site and those that are positioned more loosely relative to the transcriptional start site. Similar numbers of non-directional motifs were identified in both species and most are different. The genes associated with all 15 DNA motifs, when they occur in the peak, are enriched in specific Gene Ontology categories and show a distinct mRNA expression pattern, suggesting that there is a core promoter code in Drosophila.

Conclusion

Drosophila and human promoters use different DNA sequences to regulate gene expression, supporting the idea that evolution occurs by the modulation of gene regulation.

Drosophila mojoless, a retroposed GSK-3, has functionally diverged to acquire an essential role in male fertility

Retroposition is increasingly recognized as an important mechanism for the acquisition of new genes. We show that a glycogen synthase kinase-3 gene, shaggy (sgg), retroposed at least 50 MYA in the Drosophila genus to generate a new gene, mojoless (mjl). We have extensively analyzed the function of mjl and examined its functional divergence from the parental gene sgg in Drosophila melanogaster. Unlike Sgg, which is expressed in many tissues of both sexes, Mjl is expressed specifically in the male germ line, where it is required for male germ line survival. Our analysis indicates that mjl has acquired a specific function in the maintenance of male germ line viability. However, it has not completely lost its ancestral biochemical function and can partially compensate for loss of the parental gene sgg when ectopically expressed in somatic cells. We postulate that mjl has undergone functional diversification and is now under stabilizing selection in the Drosophila genus.

Genetic interactions between Drosophila melanogaster menin and Jun/Fos

Menin is a tumor suppressor required to prevent multiple endocrine neoplasia in humans. Mammalian menin protein is associated with chromatin modifying complexes and has been shown to bind a number of nuclear proteins, including the transcription factor JunD. Menin shows bidirectional effects acting positively on c-Jun and negatively on JunD. We have produced protein null alleles of Drosophila menin (mnn1) and have over expressed the Mnn1 protein. Flies homozygous for protein-null mnn1 alleles are viable and fertile. Localized over-expression of Mnn1 causes defects in thoracic closure, a phenotype that sometimes results from insufficient Jun activity. We observed complex genetic interactions between mnn1 and jun in different developmental settings. Our data support the idea that one function of menin is to modulate Jun activity in a manner dependent on the cellular context.

Global analysis of X-chromosome dosage compensation

Background

Drosophila melanogaster females have two X chromosomes and two autosome sets (XX;AA), while males have a single X chromosome and two autosome sets (X;AA). Drosophila male somatic cells compensate for a single copy of the X chromosome by deploying male-specific-lethal (MSL) complexes that increase transcription from the X chromosome. Male germ cells lack MSL complexes, indicating that either germline X-chromosome dosage compensation is MSL-independent, or that germ cells do not carry out dosage compensation.

Results

To investigate whether dosage compensation occurs in germ cells, we directly assayed X-chromosome transcripts using DNA microarrays and show equivalent expression in XX;AA and X;AA germline tissues. In X;AA germ cells, expression from the single X chromosome is about twice that of a single autosome. This mechanism ensures balanced X-chromosome expression between the sexes and, more importantly, it ensures balanced expression between the single X chromosome and the autosome set. Oddly, the inactivation of an X chromosome in mammalian females reduces the effective X-chromosome dose and means that females face the same X-chromosome transcript deficiency as males. Contrary to most current dosage-compensation models, we also show increased X-chromosome expression in X;AA and XX;AA somatic cells of Caenorhabditis elegans and mice.

Conclusion

Drosophila germ cells compensate for X-chromosome dose. This occurs by equilibrating X-chromosome and autosome expression in X;AA cells. Increased expression of the X chromosome in X;AA individuals appears to be phylogenetically conserved.

Core promoter sequences contribute to ovo-B regulation in the Drosophila melanogaster germline

Utilization of tightly linked ovo-A vs. ovo-B germline promoters results in the expression of OVO-A and OVO-B, C(2)H(2) transcription factors with different N -termini, and different effects on target gene transcription and on female germline development. We show that two sex-determination signals, the X chromosome number within the germ cells and a female soma, differentially regulate ovo-B and ovo-A. We have previously shown that OVO regulates ovarian tumor transcription by binding the transcription start site. We have explored the regulation of the ovo-B promoter using an extensive series of transgenic reporter gene constructs to delimit cis-regulatory sequences as assayed in wild-type and sex-transformed flies and flies with altered ovo dose. Minimum regulated expression of ovo-B requires a short region flanking the transcription start site, suggesting that the ovo-B core promoter bears regulatory information in addition to a "basal" activity. In support of this idea, the core promoter region binds distinct factors in ovary and testis extracts, but not in soma extracts, suggesting that regulatory complexes form at the start site. This idea is further supported by the evolutionarily conserved organization of OVO binding sites at or near the start sites of ovo loci in other flies.

Lexical effects and lexical properties associated with National Adult Reading Test (NART) stimuli in healthy young adults and healthy elderly adults

Fifty healthy younger adults and 45 healthy elderly adults performed a speeded lexical decision task (LDT). Stimuli consisted of 57 National Adult Reading Test (NART) words (the NART consists of "irregular" words that violate standard spelling-to-sound correspondence rules) and 57 pronounceable pseudowords (e.g., blant). Both groups displayed statistically equivalent lexicality (PW-W) effects, and error rates were lower in elderly adults. With groups equated on vocabulary ability, lexicality effects remained the same and error rates did not differ across either group. Correlational analyses confirmed the role played by word frequency and word familiarity in latencies to NART stimuli. Results are discussed regarding the importance of obtaining speeded latency measures for age-related word recognition comparisons.

Cardiac performance after reduction of myocardial hypertrophy

PURPOSE

The current study was performed to assess the functional sequelae of reducing left ventricular hypertrophy in patients with essential hypertension.

PATIENTS AND METHODS

To analyze left ventricular function and contractility in patients with essential hypertension after reduction of left ventricular hypertrophy, 14 patients with essential hypertension and left ventricular hypertrophy were studied prospectively by echocardiogram (1) before, (2) during, and (3) after left ventricular mass had been reduced by antihypertensive therapy of 19 +/- 3 months' duration. All drugs were discontinued four weeks before the first and the third study.

RESULTS

At the time of the third study, arterial pressure had returned to pretreatment values, and mean, peak, and isovolumetric (but not end-systolic) wall stress increased, whereas left ventricular mass remained diminished. Despite the increased pressure load to the heart, myocardial contractility was maintained or improved after reduction of left ventricular hypertrophy, as indicated by the ratio of end-systolic wall stress to end-systolic volume index (p less than 0.02) and by the relation of fractional shortening to end-systolic wall stress (p less than 0.06). End-diastolic volume, an indicator of preload, remained reduced after therapy (p less than 0.05). As a result, pump function of the left ventricle improved as shown by an increase in the ejection fraction (p less than 0.05), fractional fiber shortening (p less than 0.05), and velocity of circumferential fiber shortening (p less than 0.01).

CONCLUSION

Thus, in patients with essential hypertension, reduction of myocardial hypertrophy by antihypertensive therapy appears to be beneficial rather than detrimental to cardiac pump performance.