Chromatin-Associated Proteins HP1 and Mod(mdg4) Modify Y-Linked Regulatory Variation in the Drosophila Testis

Chromatin remodeling is crucial for gene regulation. Remodeling is often mediated through chemical modifications of the DNA template, DNA-associated proteins, and RNA-mediated processes. Y-linked regulatory variation (YRV) refers to the quantitative effects that polymorphic tracts of Y-linked chromatin exert on gene expression of X-linked and autosomal genes. Here we show that naturally occurring polymorphisms in the Drosophila melanogaster Y chromosome contribute disproportionally to gene expression variation in the testis. The variation is dependent on wild-type expression levels of mod(mdg4) as well as Su(var)205; the latter gene codes for heterochromatin protein 1 (HP1) in Drosophila. Testis-specific YRV is abolished in genotypes with heterozygous loss-of-function mutations for mod(mdg4) and Su(var)205 but not in similar experiments with JIL-1. Furthermore, the Y chromosome differentially regulates several ubiquitously expressed genes. The results highlight the requirement for wild-type dosage of Su(var)205 and mod(mdg4) in enabling naturally occurring Y-linked regulatory variation in the testis. The phenotypes that emerge in the context of wild-type levels of the HP1 and Mod(mdg4) proteins might be part of an adaptive response to the environment.     

Ensemble-Based Network Aggregation Improves the Accuracy of Gene Network Reconstruction

Reverse engineering approaches to constructing gene regulatory networks (GRNs) based on genome-wide mRNA expression data have led to significant biological findings, such as the discovery of novel drug targets. However, the reliability of the reconstructed GRNs needs to be improved. Here, we propose an ensemble-based network aggregation approach to improving the accuracy of network topologies constructed from mRNA expression data. To evaluate the performances of different approaches, we created dozens of simulated networks from combinations of gene-set sizes and sample sizes and also tested our methods on three Escherichia coli datasets. We demonstrate that the ensemble-based network aggregation approach can be used to effectively integrate GRNs constructed from different studies – producing more accurate networks. We also apply this approach to building a network from epithelial mesenchymal transition (EMT) signature microarray data and identify hub genes that might be potential drug targets. The R code used to perform all of the analyses is available in an R package entitled “ENA”, accessible on CRAN (http://cran.r-project.org/web/packages/ENA/).     

Draft genome sequences and description of Lactobacillus rhamnosus strains L31, L34, and L35

Lactobacillus rhamnosus is a facultative, lactic acid bacterium in the phylum Firmicutes. Lactobacillus spp. are generally considered beneficial, and specific strains of L. rhamnosus are validated probiotics. We describe the draft genomes of three L. rhamnosus strains (L31, L34, and L35) isolated from the feces of Thai breastfed infants, which exhibit anti-inflammatory properties in vitro. The three genomes range between 2.8 – 2.9 Mb, and contain approximately 2,700 protein coding genes.     

Functional Genomic Identification of Genes Required for Male Gonadal Differentiation in Caenorhabditis elegans

The Caenorhabditis elegans somatic gonad develops from a four-cell primordium into a mature organ that differs dramatically between the sexes in overall morphology (two arms in hermaphrodites and one in males) and in the cell types comprising it. Gonadal development in C. elegans is well studied, but regulation of sexual differentiation, especially later in gonadal development, remains poorly elucidated. To identify genes involved in this process, we performed a genome-wide RNAi screen using sex-specifically expressed gonadal GFP reporters. This screen identified several phenotypic classes, including ∼70 genes whose depletion feminized male gonadal cells. Among the genes required for male cell fate specification are Wnt/β-catenin pathway members, cell cycle regulators, and genes required for mitotic spindle function and cytokinesis. We find that a Wnt/β-catenin pathway independent of extracellular Wnt ligand is essential for asymmetric cell divisions and male differentiation during gonadal development in larvae. We also find that the cell cycle regulators cdk-1 and cyb-3 and the spindle/cytokinesis regulator zen-4 are required for Wnt/β-catenin pathway activity in the developing gonad. After sex is determined in the gonadal primordium the global sex determination pathway is dispensable for gonadal sexual fate, suggesting that male cell fates are promoted and maintained independently of the global pathway during this period.THE Caenorhabditis elegans gonad derives from a simple primordium of four cells that coalesces during embryogenesis and contains two somatic gonad precursors (SGPs), Z1 and Z4, flanking two germline precursors, Z2 and Z3 (Kimble and Hirsh 1979). The SGPs undergo very different developmental programs in each sex, involving sexually dimorphic cell lineages and migrations and sex-specific cellular differentiation. The result is a two-armed bilaterally symmetrical gonad in the adult hermaphrodite or a single-armed asymmetric gonad in the adult male. The high degree of sexual dimorphism of the mature organ and variety of cellular events that occur sex specifically during its development make the C. elegans gonad an outstanding model for sex-specific organogenesis.Development of the somatic gonad occurs in two phases. The early phase defines the gonadal axes and establishes the precursors of the major gonadal cell types. This takes place during the first larval stage (L1), beginning shortly after hatching with the first division of the SGPs. In both sexes SGP division is asymmetric in terms of both the sizes and the fates of the daughter cells, and establishes the proximal/distal axis of the gonad (Hirsh et al. 1976; Kimble and Hirsh 1979). The global sex determination pathway establishes the future sex of the gonad around the time of hatching (Klass et al. 1976; Nelson et al. 1978), and sexual dimorphism is already apparent when the SGPs divide: the size asymmetry of the SGP daughters is much more pronounced in males than hermaphrodites. In both sexes the asymmetry of the first SGP division requires a Wnt/β-catenin pathway. Mutations compromising this pathway cause a “symmetrical sisters” phenotype in which both daughters adopt the same fate (Miskowski et al. 2001; Siegfried and Kimble 2002; Phillips and Kimble 2009). Sex specificity is imposed on the SGPs by the global sex determining gene tra-1 (Hodgkin 1987) and the gonad-specific sex determining gene fkh-6 (Chang et al. 2004). These genes play opposing roles in SGP sex determination, with tra-1 feminizing and fkh-6 masculinizing the somatic gonad, and they also act redundantly to promote mitotic proliferation of the SGP lineage (Chang et al. 2004). SGP sex determination is linked to cell cycle progression by cyclin D, which is required to overcome repression of fkh-6 expression in the SGPs by E2F (Tilmann and Kimble 2005).The later phase of gonadal development involves the elongation of the gonad, together with cellular proliferation and differentiation, and lasts from L2 to adulthood. During L2 the somatic cells enlarge and leader cells (distal tip cells in the hermaphrodite, linker cell in the male) begin long-range migrations that extend the gonad. During L3, somatic gonad cell division resumes in both sexes, leading to the formation of differentiated somatic cell types by the end of L3 or beginning of L4. Gonadal morphogenesis is completed and gametogenesis begins during L4 (Kimble and Hirsh 1979).Although SGP division and much of hermaphrodite gonadal development have been well studied (Hubbard and Greenstein 2000), sexual cell fate specification in the somatic gonad is more poorly understood, particularly after the L1 stage. Despite the importance of fkh-6 in promoting male differentiation, it is expressed in males only during early L1 and null mutants have incomplete gonadal sex reversal. We have therefore performed a genome-wide RNAi screen to identify additional genes required after hatching for gonadal development in each sex. Among the advantages of this approach is the ability to identify gonadal regulators that also are essential for embryonic development. To our knowledge this is the first functional genomic study of gonadal sex differentiation.The screen identified many genes whose depletion disrupts gonadogenesis in each sex and nearly 70 genes whose depletion causes gonadal feminization in males. Prominent among this latter class were components of a Wnt/β-catenin pathway, cell cycle regulators, and genes involved in mitotic spindle function and cytokinesis. We find that Wnt/β-catenin activity continues in both sexes after SGP division and is required for male cell fate commitment in the gonad. We also find that the cyclin-dependent kinase cdk-1 and its cognate cyclin cyb-3 as well as the mitotic spindle regulator zen-4 are required for gonadal Wnt/β-catenin pathway activity, providing a potential new link between the cell cycle, asymmetric division, and sexual differentiation. The feminization caused by depletion of Wnt/β-catenin pathway components or cdk-1 is independent of the global sex determination pathway, suggesting that sexual fates in the male gonad remain plastic after the primary sex determination decision.     

Non-contiguous finished genome sequence and description of Halopiger goleamassiliensis sp. nov.

Halopiger goleamassiliensis strain IIH3^T sp. nov. is a novel, extremely halophilic archaeon within the genus Halopiger. This strain was isolated from an evaporitic sediment in El Golea Lake, Ghardaïa region (Algeria). The type strain is strain IIH3^T. H. goleamassiliensis is moderately thermophilic, neutrophilic, non-motile and coccus-shaped. Here we describe the features of this organism, together with the complete genome sequence and annotation. The 3,906,923 bp long genome contains 3,854 protein-encoding genes and 49 RNA genes (1 gene is 16S rRNA, 1 gene is 23S rRNA, 3 genes are 5S rRNA, and 44 are tRNA genes).     

The Border Between the Ultrabithorax and abdominal-A Regulatory Domains in the Drosophila Bithorax Complex

The bithorax complex in Drosophila melanogaster includes three homeobox-containing genes—Ultrabithorax (Ubx), abdominal-A (abd-A), and Abdominal-B (Abd-B)—which are required for the proper differentiation of the posterior 10 segments of the body. Each of these genes has multiple distinct regulatory regions; there is one for each segmental unit of the body plan where the genes are expressed. One additional protein- coding gene in the bithorax complex, Glut3, a sugar-transporter homolog, can be deleted without phenotype. We focus here on the upstream regulatory region for Ubx, the bithoraxoid (bxd) domain, and its border with the adjacent infraabdominal-2 (iab-2) domain, which controls abdA. These two domains can be defined by the phenotypes of rearrangement breakpoints, and by the expression patterns of enhancer traps. In D. virilis, the homeotic cluster is split between Ubx and abd-A, and so the border can also be located by a sequence comparison between species. When the border region is deleted in melanogaster, the flies show a dominant phenotype called Front-ultraabdominal (Fub); the first abdominal segment is transformed into a copy of the second abdominal segment. Thus, the border blocks the spread of activation from the bxd domain into the iab-2 domain.     

DNA Replication Error-Induced Extinction of Diploid Yeast

Genetic defects in DNA polymerase accuracy, proofreading, or mismatch repair (MMR) induce mutator phenotypes that accelerate adaptation of microbes and tumor cells. Certain combinations of mutator alleles synergistically increase mutation rates to levels that drive extinction of haploid cells. The maximum tolerated mutation rate of diploid cells is unknown. Here, we define the threshold for replication error-induced extinction (EEX) of diploid Saccharomyces cerevisiae. Double-mutant pol3 alleles that carry mutations for defective DNA polymerase-δ proofreading (pol3-01) and accuracy (pol3-L612M or pol3-L612G) induce strong mutator phenotypes in heterozygous diploids (POL3/pol3-01,L612M or POL3/pol3-01,L612G). Both pol3-01,L612M and pol3-01,L612G alleles are lethal in the homozygous state; cells with pol3-01,L612M divide up to 10 times before arresting at random stages in the cell cycle. Antimutator eex mutations in the pol3 alleles suppress this lethality (pol3-01,L612M,eex or pol3-01,L612G,eex). MMR defects synergize with pol3-01,L612M,eex and pol3-01,L612G,eex alleles, increasing mutation rates and impairing growth. Conversely, inactivation of the Dun1 S-phase checkpoint kinase suppresses strong pol3-01,L612M,eex and pol3-01,L612G,eex mutator phenotypes as well as the lethal pol3-01,L612M phenotype. Our results reveal that the lethal error threshold in diploids is 10 times higher than in haploids and likely determined by homozygous inactivation of essential genes. Pronounced loss of fitness occurs at mutation rates well below the lethal threshold, suggesting that mutator-driven cancers may be susceptible to drugs that exacerbate replication errors.     

From Gigabyte to Kilobyte: A Bioinformatics Protocol for Mining Large RNA-Seq Transcriptomics Data

RNA-Seq techniques generate hundreds of millions of short RNA reads using next-generation sequencing (NGS). These RNA reads can be mapped to reference genomes to investigate changes of gene expression but improved procedures for mining large RNA-Seq datasets to extract valuable biological knowledge are needed. RNAMiner—a multi-level bioinformatics protocol and pipeline—has been developed for such datasets. It includes five steps: Mapping RNA-Seq reads to a reference genome, calculating gene expression values, identifying differentially expressed genes, predicting gene functions, and constructing gene regulatory networks. To demonstrate its utility, we applied RNAMiner to datasets generated from Human, Mouse, Arabidopsis thaliana, and Drosophila melanogaster cells, and successfully identified differentially expressed genes, clustered them into cohesive functional groups, and constructed novel gene regulatory networks. The RNAMiner web service is available at http://calla.rnet.missouri.edu/rnaminer/index.html.     

Probing the Boundaries of Orthology: The Unanticipated Rapid Evolution of Drosophila centrosomin

The rapid evolution of essential developmental genes and their protein products is both intriguing and problematic. The rapid evolution of gene products with simple protein folds and a lack of well-characterized functional domains typically result in a low discovery rate of orthologous genes. Additionally, in the absence of orthologs it is difficult to study the processes and mechanisms underlying rapid evolution. In this study, we have investigated the rapid evolution of centrosomin (cnn), an essential gene encoding centrosomal protein isoforms required during syncytial development in Drosophila melanogaster. Until recently the rapid divergence of cnn made identification of orthologs difficult and questionable because Cnn violates many of the assumptions underlying models for protein evolution. To overcome these limitations, we have identified a group of insect orthologs and present conserved features likely to be required for the functions attributed to cnn in D. melanogaster. We also show that the rapid divergence of Cnn isoforms is apparently due to frequent coding sequence indels and an accelerated rate of intronic additions and eliminations. These changes appear to be buffered by multi-exon and multi-reading frame maximum potential ORFs, simple protein folds, and the splicing machinery. These buffering features also occur in other genes in Drosophila and may help prevent potentially deleterious mutations due to indels in genes with large coding exons and exon-dense regions separated by small introns. This work promises to be useful for future investigations of cnn and potentially other rapidly evolving genes and proteins.     

Effete,a Drosophila Chromatin-Associated Ubiquitin-Conjugating Enzyme That Affects Telomeric and Heterochromatic Position Effect Variegation

Drosophila telomeres are elongated by the transposition of telomere-specific retrotransposons rather than telomerase activity. Proximal to the terminal transposon array, Drosophila chromosomes contain several kilobases of a complex satellite DNA termed telomere-associated sequences (TASs). Reporter genes inserted into or next to the TAS are silenced through a mechanism called telomere position effect (TPE). TPE is reminiscent of the position effect variegation (PEV) induced by Drosophila constitutive heterochromatin. However, most genes that modulate PEV have no effect on TPE, and systematic searches for TPE modifiers have so far identified only a few dominant suppressors. Surprisingly, only a few of the genes required to prevent telomere fusion have been tested for their effect on TPE. Here, we show that with the exception of the effete (eff; also called UbcD1) mutant alleles, none of the tested mutations at the other telomere fusion genes affects TPE. We also found that mutations in eff, which encodes a class I ubiquitin-conjugating enzyme, act as suppressors of PEV. Thus, eff is one of the rare genes that can modulate both TPE and PEV. Immunolocalization experiments showed that Eff is a major constituent of polytene chromosomes. Eff is enriched at several euchromatic bands and interbands, the TAS regions, and the chromocenter. Our results suggest that Eff associates with different types of chromatin affecting their abilities to regulate gene expression.     

A Novel Nondevelopmental Role of the SAX-7/L1CAM Cell Adhesion Molecule in Synaptic Regulation in Caenorhabditis elegans

The L1CAM family of cell adhesion molecules is a conserved set of single-pass transmembrane proteins that play diverse roles required for proper nervous system development and function. Mutations in L1CAMs can cause the neurological L1 syndrome and are associated with autism and neuropsychiatric disorders. L1CAM expression in the mature nervous system suggests additional functions besides the well-characterized developmental roles. In this study, we demonstrate that the gene encoding the Caenorhabditis elegans L1CAM, sax-7, genetically interacts with gtl-2, as well as with unc-13 and rab-3, genes that function in neurotransmission. These sax-7 genetic interactions result in synthetic phenotypes that are consistent with abnormal synaptic function. Using an inducible sax-7 expression system and pharmacological reagents that interfere with cholinergic transmission, we uncovered a previously uncharacterized nondevelopmental role for sax-7 that impinges on synaptic function.     

TargetCompare: A web interface to compare simultaneous miRNAs targets

MicroRNAs (miRNAs) are small non-coding nucleotide sequences between 17 and 25 nucleotides in length that primarily function in the regulation of gene expression. A since miRNA has thousand of predict targets in a complex, regulatory cell signaling network. Therefore, it is of interest to study multiple target genes simultaneously. Hence, we describe a web tool (developed using Java programming language and MySQL database server) to analyse multiple targets of pre-selected miRNAs. We cross validated the tool in eight most highly expressed miRNAs in the antrum region of stomach. This helped to identify 43 potential genes that are target of at least six of the referred miRNAs. The developed tool aims to reduce the randomness and increase the chance of selecting strong candidate target genes and miRNAs responsible for playing important roles in the studied tissue.

Availability

http://lghm.ufpa.br/targetcompare     

Complete genome sequence of Granulicella mallensis type strain MP5ACTX8T,an acidobacterium from tundra soil

Granulicella mallensis MP5ACTX8^T is a novel species of the genus Granulicella in subdivision 1of Acidobacteria. G. mallensis is of ecological interest being a member of the dominant soil bacterial community active at low temperatures and nutrient limiting conditions in Arctic alpine tundra. G. mallensis is a cold-adapted acidophile and a versatile heterotroph that hydrolyzes a suite of sugars and complex polysaccharides. Genome analysis revealed metabolic versatility with genes involved in metabolism and transport of carbohydrates. These include gene modules encoding the carbohydrate-active enzyme (CAZyme) family involved in breakdown, utilization and biosynthesis of diverse structural and storage polysaccharides including plant based carbon polymers. The genome of Granulicella mallensis MP5ACTX8^T consists of a single replicon of 6,237,577 base pairs (bp) with 4,907 protein-coding genes and 53 RNA genes.