Genomic Variation and Its Impact on Gene Expression in Drosophila melanogaster

Understanding the relationship between genetic and phenotypic variation is one of the great outstanding challenges in biology. To meet this challenge, comprehensive genomic variation maps of human as well as of model organism populations are required. Here, we present a nucleotide resolution catalog of single-nucleotide, multi-nucleotide, and structural variants in 39 Drosophila melanogaster Genetic Reference Panel inbred lines. Using an integrative, local assembly-based approach for variant discovery, we identify more than 3.6 million distinct variants, among which were more than 800,000 unique insertions, deletions (indels), and complex variants (1 to 6,000 bp). While the SNP density is higher near other variants, we find that variants themselves are not mutagenic, nor are regions with high variant density particularly mutation-prone. Rather, our data suggest that the elevated SNP density around variants is mainly due to population-level processes. We also provide insights into the regulatory architecture of gene expression variation in adult flies by mapping cis-expression quantitative trait loci (cis-eQTLs) for more than 2,000 genes. Indels comprise around 10% of all cis-eQTLs and show larger effects than SNP cis-eQTLs. In addition, we identified two-fold more gene associations in males as compared to females and found that most cis-eQTLs are sex-specific, revealing a partial decoupling of the genomic architecture between the sexes as well as the importance of genetic factors in mediating sex-biased gene expression. Finally, we performed RNA-seq-based allelic expression imbalance analyses in the offspring of crosses between sequenced lines, which revealed that the majority of strong cis-eQTLs can be validated in heterozygous individuals.     

Hfq Is a Global Regulator That Controls the Pathogenicity of Staphylococcus aureus

The Hfq protein is reported to be an RNA chaperone, which is involved in the stress response and the virulence of several pathogens. In E. coli, Hfq can mediate the interaction between some sRNAs and their target mRNAs. But it is controversial whether Hfq plays an important role in S. aureus. In this study, we found that the deletion of hfq gene in S. aureus 8325-4 can increase the surface carotenoid pigments. The hfq mutant was more resistant to oxidative stress but the pathogenicity of the mutant was reduced. We reveal that the Hfq protein can be detected only in some S. aureus strains. Using microarray and qRT-PCR, we identified 116 genes in the hfq mutant which had differential expression from the wild type, most of which are related to the phenotype and virulence of S. aureus. Among the 116 genes, 49 mRNAs can specifically bind Hfq protein, which indicates that Hfq also acts as an RNA binding protein in S. aureus. Our data suggest that Hfq protein of S. aureus is a multifunctional regulator involved in stress and virulence.     

Impact of Elevated Plasma Serotonin on Global Gene Expression of Murine Megakaryocytes

Background

Serotonin (5-HT) is a biogenic amine that also acts as a mitogen and a developmental signal early in rodent embryogenesis. Genetic and pharmacological disruption of 5-HT signaling causes various diseases and disorders via mediating central nervous system, cardiovascular system, and serious abnormalities on a growing embryo. Today, neither the effective modulators on 5-HT signaling pathways nor the genes affected by 5-HT signal are well known yet.

Methodology/Principal Findings

In an attempt to identify the genes altered by 5-HT signaling pathways, we analyzed the global gene expression via the Illumina array platform using the mouse WG-6 v2.0 Expression BeadChip containing 45,281 probe sets representing 30,854 genes in megakaryocytes isolated from mice infused with 5-HT or saline. We identified 723 differentially expressed genes of which 706 were induced and 17 were repressed by elevated plasma 5-HT.

Conclusions/Significance

Hierarchical gene clustering analysis was utilized to represent relations between groups and clusters. Using gene ontology mining tools and canonical pathway analyses, we identified multiple biological pathways that are regulated by 5-HT: (i) cytoskeletal remodeling, (ii) G-protein signaling, (iii) vesicular transport, and (iv) apoptosis and survival. Our data encompass the first extensive genome-wide based profiling in the progenitors of platelets in response to 5-HT elevation in vivo.     

NAD+-Glycohydrolase Promotes Intracellular Survival of Group A Streptococcus

A global increase in invasive infections due to group A Streptococcus (S. pyogenes or GAS) has been observed since the 1980s, associated with emergence of a clonal group of strains of the M1T1 serotype. Among other virulence attributes, the M1T1 clone secretes NAD⁺-glycohydrolase (NADase). When GAS binds to epithelial cells in vitro, NADase is translocated into the cytosol in a process mediated by streptolysin O (SLO), and expression of these two toxins is associated with enhanced GAS intracellular survival. Because SLO is required for NADase translocation, it has been difficult to distinguish pathogenic effects of NADase from those of SLO. To resolve the effects of the two proteins, we made use of anthrax toxin as an alternative means to deliver NADase to host cells, independently of SLO. We developed a novel method for purification of enzymatically active NADase fused to an amino-terminal fragment of anthrax toxin lethal factor (LFn-NADase) that exploits the avid, reversible binding of NADase to its endogenous inhibitor. LFn-NADase was translocated across a synthetic lipid bilayer in vitro in the presence of anthrax toxin protective antigen in a pH-dependent manner. Exposure of human oropharyngeal keratinocytes to LFn-NADase in the presence of protective antigen resulted in cytosolic delivery of NADase activity, inhibition of protein synthesis, and cell death, whereas a similar construct of an enzymatically inactive point mutant had no effect. Anthrax toxin-mediated delivery of NADase in an amount comparable to that observed during in vitro infection with live GAS rescued the defective intracellular survival of NADase-deficient GAS and increased the survival of SLO-deficient GAS. Confocal microscopy demonstrated that delivery of LFn-NADase prevented intracellular trafficking of NADase-deficient GAS to lysosomes. We conclude that NADase mediates cytotoxicity and promotes intracellular survival of GAS in host cells.     

Impact of Gene Dosage on Gene Expression,Biological Processes and Survival in Cervical Cancer: A Genome-Wide Follow-Up Study

We investigated the role of tumor copy number (CN)–altered genome (CN-AG) in the carcinogenesis of cervical cancer (CC), especially its effect on gene expression, biological processes, and patient survival. Fifty-nine human papillomavirus 16 (HPV16)-positive CCs were investigated with microarrays–31 for mapping CN-AG and 55 for global gene expression, with 27 CCs in common. Five-year survival was investigated in 55 patients. Deletions and amplifications >2.5 Mb were defined as CN alterations. The %CN-AG varied from 0 to 32.2% (mean = 8.1±8.9). Tumors were classified as low (mean = 0.5±0.6, n = 11), medium (mean = 5.4±2.4, n = 10), or high (mean = 19.2±6.6, n = 10) CN. The highest %CN-AG was found in 3q, which contributed an average of 55% of all CN alterations. Genome-wide, only 5.3% of CN-altered genes were deregulated directly by gene dosage. In contrast, the rate in fully duplicated 3q was twice as high. Amplification of 3q explained 23.2% of deregulated genes in whole tumors (r² = 0.232, p = 0.006; analysis of variance), including genes located in 3q and other chromosomes. A total of 862 genes were deregulated exclusively in high-CN tumors, but only 22.9% were CN altered. This suggests that the remaining genes are not deregulated directly by gene dosage, but by mechanisms induced in trans by CN-altered genes. Anaphase-promoting complex/cyclosome (APC/C)-dependent proteasome proteolysis, glycolysis, and apoptosis were upregulated, whereas cell adhesion and angiogenesis were downregulated exclusively in high-CN tumors. The high %CN-AG and upregulated gene expression profile of APC/C-dependent proteasome proteolysis were associated with poor patient survival (p<0.05, log-rank test). Along with glycolysis, they were linearly associated with FIGO stage (r>0.38, p<0.01, Spearman test). Therefore, inhibition of APC/C-dependent proteasome proteolysis and glycolysis could be useful for CC treatment. However, whether they are indispensable for tumor growth remains to be demonstrated.     

Unveiling the Intracellular Survival Gene Kit of Trypanosomatid Parasites

Trypanosomatids are unicellular protozoans of medical and economical relevance since they are the etiologic agents of infectious diseases in humans as well as livestock. Whereas Trypanosoma cruzi and different species of Leishmania are obligate intracellular parasites, Trypanosoma brucei and other trypanosomatids develop extracellularly throughout their entire life cycle. After their genomes have been sequenced, various comparative genomic studies aimed at identifying sequences involved with host cell invasion and intracellular survival have been described. However, for only a handful of genes, most of them present exclusively in the T. cruzi or Leishmania genomes, has there been any experimental evidence associating them with intracellular parasitism. With the increasing number of published complete genome sequences of members of the trypanosomatid family, including not only different Trypanosoma and Leishmania strains and subspecies but also trypanosomatids that do not infect humans or other mammals, we may now be able to contemplate a slightly better picture regarding the specific set of parasite factors that defines each organism''s mode of living and the associated disease phenotypes. Here, we review the studies concerning T. cruzi and Leishmania genes that have been implicated with cell invasion and intracellular parasitism and also summarize the wealth of new information regarding the mode of living of intracellular parasites that is resulting from comparative genome studies that are based on increasingly larger trypanosomatid genome datasets.     

Faster-X Evolution of Gene Expression in Drosophila

DNA sequences on X chromosomes often have a faster rate of evolution when compared to similar loci on the autosomes, and well articulated models provide reasons why the X-linked mode of inheritance may be responsible for the faster evolution of X-linked genes. We analyzed microarray and RNA–seq data collected from females and males of six Drosophila species and found that the expression levels of X-linked genes also diverge faster than autosomal gene expression, similar to the “faster-X” effect often observed in DNA sequence evolution. Faster-X evolution of gene expression was recently described in mammals, but it was limited to the evolutionary lineages shortly following the creation of the therian X chromosome. In contrast, we detect a faster-X effect along both deep lineages and those on the tips of the Drosophila phylogeny. In Drosophila males, the dosage compensation complex (DCC) binds the X chromosome, creating a unique chromatin environment that promotes the hyper-expression of X-linked genes. We find that DCC binding, chromatin environment, and breadth of expression are all predictive of the rate of gene expression evolution. In addition, estimates of the intraspecific genetic polymorphism underlying gene expression variation suggest that X-linked expression levels are not under relaxed selective constraints. We therefore hypothesize that the faster-X evolution of gene expression is the result of the adaptive fixation of beneficial mutations at X-linked loci that change expression level in cis. This adaptive faster-X evolution of gene expression is limited to genes that are narrowly expressed in a single tissue, suggesting that relaxed pleiotropic constraints permit a faster response to selection. Finally, we present a conceptional framework to explain faster-X expression evolution, and we use this framework to examine differences in the faster-X effect between Drosophila and mammals.