Search

Abstract 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.

Abstract Background Brucella is an intracellular pathogen capable of infecting animals and humans. There are six recognized species of Brucella that differ in their host preference. The genomes of the three Brucella species have been recently sequenced. Comparison of the three revealed over 98% sequence similarity at the protein level and enabled computational identification of common and differentiating genes. We validated these computational predictions and examined the expression patterns of the putative unique and differentiating genes, using genomic and reverse transcription PCR. We then screened a set of differentiating genes against classical Brucella biovars and showed the applicability of these regions in the design of diagnostic tests. Results We have identified and tested set of molecular targets that are associated in unique patterns with each of the sequenced Brucella spp. A comprehensive comparison was made among the published genome sequences of B. abortus, B. melitensis and B. suis. The comparison confirmed published differences between the three Brucella genomes, and identified subsets of features that were predicted to be of interest in a functional comparison of B. melitensis and B. suis to B. abortus. Differentiating sequence regions from B. abortus, B. melitensis and B. suis were used to develop PCR primers to test for the existence and in vitro transcription of these genes in these species. Only B. suis is found to have a significant number of unique genes, but combinations of genes and regions that exist in only two out of three genomes and are therefore useful for diagnostics were identified and confirmed. Conclusion Although not all of the differentiating genes identified were transcribed under steady state conditions, a group of genes sufficient to discriminate unambiguously between B. suis, B. melitensis, and B. abortus was identified. We present an overview of these genomic differences and the use of these features to discriminate among a number of Brucella biovars.

Abstract 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.

Genetic and epigenetic variations in regulatory enhancer elements increase susceptibility to a range of pathologies. Despite recent advances, linking enhancer elements to target genes and predicting transcriptional outcomes of enhancer dysfunction remain significant challenges. Using 3D chromatin conformation assays, we generated an extensive enhancer interaction dataset for the human pancreas, encompassing more than 20 donors and five major cell types, including both exocrine and endocrine compartments. We employed a network approach to parse chromatin interactions into enhancer-promoter tree models, facilitating a quantitative, genome-wide analysis of enhancer connectivity. With these tree models, we developed a machine learning algorithm to estimate the impact of enhancer perturbations on cell type-specific gene expression in the human pancreas. Orthogonal to our computational approach, we perturbed enhancer function in primary human pancreas cells using CRISPR interference and quantified the effects at the single-cell level through RNA FISH coupled with high-throughput imaging. Our enhancer tree models enabled the annotation of common germline risk variants associated with pancreas diseases, linking them to putative target genes in specific cell types. For pancreatic ductal adenocarcinoma, we found a stronger enrichment of disease susceptibility variants within acinar cell regulatory elements, despite ductal cells historically being assumed as the primary cell-of-origin. Our integrative approach-combining cell type-specific enhancer-promoter interaction mapping, computational models, and single-cell enhancer perturbation assays-produced a robust resource for studying the genetic basis of pancreas disorders.

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., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

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 4 -induced reduction of ac4C and conversion to thymidine followed by sequencing (RedaC:T-seq). Total RNA is ribodepleted and then treated with NaBH 4 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). 1 ., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Inc.)

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.

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., Competing Interests: The authors declare that they have no competing interests., (Copyright: © 2021 Faculty Opinions Ltd.)

Background: Coronavirus disease 2019 (COVID-19) is often compared with seasonal influenza and the two diseases have similarities, including the risk of systemic manifestations such as AKI. The aim of this study was to perform a comparative analysis of the prevalence, risk factors, and outcomes of AKI in patients who were hospitalized with COVID-19 and influenza., Methods: Retrospective cohort study of patients who were hospitalized with COVID-19 ( n =325) or seasonal influenza ( n =433). AKI was defined by Kidney Disease: Improving Global Outcomes (KDIGO) criteria. Baseline characteristics and hospitalization data were collected, and multivariable analysis was performed to determine the independent predictors for AKI., Results: AKI occurred in 33% of COVID-19 hospitalizations (COV-AKI) and 33% of influenza hospitalizations (FLU-AKI). After adjusting for age, sex, and comorbidity count, the risk of stage 3 AKI was significantly higher in COV-AKI (OR, 3.46; 95% CI, 1.63 to 7.37). Pre-existing CKD was associated with a six- to seven-fold increased likelihood for FLU-AKI and COV-AKI. Mechanical ventilation was associated with a higher likelihood of developing AKI in the COVID-19 cohort (OR, 5.85; 95% CI, 2.30 to 15.63). Black race, after adjustment for comorbidities, was an independent risk for COV-AKI., Conclusions: Pre-existing CKD was a major risk factor for AKI in both cohorts. Black race (independent of comorbidities) and mechanical ventilation were associated with a higher risk of developing COV-AKI, which is characterized by a higher burden of stage 3 AKI and overall poorer prognosis., Competing Interests: M. Ozieh is supported by the National Center for Advancing Translational Sciences, National Institutes of Health award KL2TR001438. K.R. Regner reports having consultancy agreements with Mallinckrodt Pharma and reports receiving honoraria from Mallinckrodt. All remaining authors have nothing to disclose., (Copyright © 2021 by the American Society of Nephrology.)

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)., Competing Interests: Competing interestsThe authors declare no competing interests., (Copyright © 2019 The Authors; exclusive licensee Bio-protocol LLC.)

Patients with advanced chronic kidney disease (CKD), including end-stage renal disease (ESRD), have a life-threatening illness complicated by high morbidity and mortality and, therefore, should be suitable candidates for early intervention by palliative care specialists. However, the average patient with CKD does not have an advanced care plan, has multiple debilitating symptoms, and does not utilise hospice care at the end of life. In this review, we outline the scope of the problem of unmet palliative care needs for patients with advanced CKD and ESRD, barriers to improving palliative care for patients with renal failure, and possible future directions for palliative nephrology., (© Royal College of Physicians 2019. All rights reserved.)

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., (Published by Elsevier Inc.)

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., Competing Interests: The authors have declared that no competing interests exist.

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., (Published by Elsevier Inc.)

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., (Published by Oxford University Press on behalf of Nucleic Acids Research 2017.)

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.

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., (Published by Oxford University Press on behalf of Nucleic Acids Research 2015. This work is written by (a) US Government employee(s) and is in the public domain in the US.)

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.

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., (Published 2015. This article is a U.S. Government work and is in the public domain in the USA.)

Somatic sexual determination and behavior in Drosophila melanogaster are under the control of a genetic cascade initiated by Sex lethal (Sxl). In the female soma, SXL RNA-binding protein regulates the splicing of transformer (tra) transcripts into a female-specific form. The RNA-binding protein TRA and its cofactor TRA2 function in concert in females, whereas SXL, TRA, and TRA2 are thought to not function in males. To better understand sex-specific regulation of gene expression, we analyzed male and female head transcriptome datasets for expression levels and splicing, quantifying sex-biased gene expression via RNA-Seq and qPCR. Our data uncouple the effects of Sxl and tra/tra2 in females in the-sex-biased alternative splicing of head transcripts from the X-linked locus found in neurons (fne), encoding a pan-neuronal RNA-binding protein of the ELAV family. We show that FNE protein levels are downregulated by Sxl in female heads, also independently of tra/tra2. We argue that this regulation may have important sexually dimorphic consequences for the regulation of nervous system development or function., (Copyright © 2015 Sun et al.)

Primary sex-determination "switches" evolve rapidly, but Doublesex (DSX)-related transcription factors (DMRTs) act downstream of these switches to control sexual development in most animal species. Drosophila dsx encodes female- and male-specific isoforms (DSX(F) and DSX(M)), but little is known about how dsx controls sexual development, whether DSX(F) and DSX(M) bind different targets, or how DSX proteins direct different outcomes in diverse tissues. We undertook genome-wide analyses to identify DSX targets using in vivo occupancy, binding site prediction, and evolutionary conservation. We find that DSX(F) and DSX(M) bind thousands of the same targets in multiple tissues in both sexes, yet these targets have sex- and tissue-specific functions. Interestingly, DSX targets show considerable overlap with targets identified for mouse DMRT1. DSX targets include transcription factors and signaling pathway components providing for direct and indirect regulation of sex-biased expression., (Copyright © 2014 Elsevier Inc. All rights reserved.)

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.

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., (© 2014 Chen et al.; Published by Cold Spring Harbor Laboratory Press.)

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.

Drosophila melanogaster is one of the most well studied genetic model organisms; nonetheless, its genome still contains unannotated coding and non-coding genes, transcripts, exons and RNA editing sites. Full discovery and annotation are pre-requisites for understanding how the regulation of transcription, splicing and RNA editing directs the development of this complex organism. Here we used RNA-Seq, tiling microarrays and cDNA sequencing to explore the transcriptome in 30 distinct developmental stages. We identified 111,195 new elements, including thousands of genes, coding and non-coding transcripts, exons, splicing and editing events, and inferred protein isoforms that previously eluded discovery using established experimental, prediction and conservation-based approaches. These data substantially expand the number of known transcribed elements in the Drosophila genome and provide a high-resolution view of transcriptome dynamics throughout development.

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.

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.

Comparative analysis of multiple genomes in a phylogenetic framework dramatically improves the precision and sensitivity of evolutionary inference, producing more robust results than single-genome analyses can provide. The genomes of 12 Drosophila species, ten of which are presented here for the first time (sechellia, simulans, yakuba, erecta, ananassae, persimilis, willistoni, mojavensis, virilis and grimshawi), illustrate how rates and patterns of sequence divergence across taxa can illuminate evolutionary processes on a genomic scale. These genome sequences augment the formidable genetic tools that have made Drosophila melanogaster a pre-eminent model for animal genetics, and will further catalyse fundamental research on mechanisms of development, cell biology, genetics, disease, neurobiology, behaviour, physiology and evolution. Despite remarkable similarities among these Drosophila species, we identified many putatively non-neutral changes in protein-coding genes, non-coding RNA genes, and cis-regulatory regions. These may prove to underlie differences in the ecology and behaviour of these diverse species.

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.

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.

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.

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.

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 < or = 1 x 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.

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.

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.

West Virginia leads the nation in mortality from coronary heart disease among both men and women aged 35 to 74. Although there has been some research with respect to behavioral risk factors, little is known about the prevalence of high serum cholesterol among West Virginia residents. The present paper begins this examination by reporting the results of a medical chart audit conducted recently in a rural medicine clinic in Matewan, West Virginia. Results revealed that only 17 per cent of the 501 charts reviewed reported serum cholesterol. Significant differences were noted between the local sample and a national comparison for two groups. The implications of these findings are discussed.