Interaction between bisphenol A and dietary sugar affects global gene transcription in Drosophila melanogaster

Human exposure to environmental toxins is a public health issue. The microarray data available in the Gene Expression Omnibus database under accession number GSE55655 and GSE55670GSE55655GSE55670 show the isolated and combined effects of dietary sugar and two organic compounds present in a variety of plastics [bisphenol A (BPA) and Bis(2-ethylhexyl) phthalate (DEHP)] on global gene expression in Drosophila melanogaster. The study was carried out with samples collected from flies exposed to these compounds for a limited period of time (48 h) in the adult stage, or throughout the entire development of the insect. The arrays were normalized using the limma/Bioconductor package. Differential expression was inferred using linear models in limma and BAGEL. The data show that each compound had its unique consequences to gene expression, and that the individual effect of each organic compound is maximized with the joint ingestion of dietary sugar.

Mito-nuclear interactions as drivers of gene movement on and off the X-chromosome

Background

Mito-nuclear gene interactions regulate energy conversion, and are fundamental to eukaryotes. Generally, mito-nuclear coadaptation would be most efficient if the interacting nuclear genes were X-linked, because this maximizes the probability of favorable mito-nuclear allelic combinations co-transmitting across generations. Thus, under a coadaptation (CA) hypothesis, nuclear genes essential for mitochondrial function might be under selection to relocate to the X-chromosome. However, maternal inheritance predisposes the mitochondrial DNA (mtDNA) to accumulate variation that, while male-harming, is benign to females. Numerous nuclear genes were recently reported in Drosophila melanogaster, which exhibit male-specific patterns of differential expression when placed alongside different mtDNA haplotypes, suggesting that nuclear genes are sensitive to an underlying male-specific mitochondrial mutation load. These genes are thus candidates for involvement in mito-nuclear interactions driven by sexual conflict (SC), and selection might have moved them off the X-chromosome to facilitate an optimal evolutionary counter-response, through males, to the presence of male-harming mtDNA mutations. Furthermore, the presence of male-harming mtDNA mutations could exert selection for modifiers on the Y-chromosome, thus placing these mito-sensitive nuclear genes at the center of an evolutionary tug-of-war between mitochondrion and Y-chromosome.

We test these hypotheses by examining the chromosomal distributions of three distinct sets of mitochondrial-interacting nuclear genes in D. melanogaster; the first is a list of genes with mitochondrial annotations by Gene Ontologies, the second is a list comprising the core evolutionary-conserved mitochondrial proteome, and the third is a list of genes involved in male-specific responses to maternally-inherited mitochondrial variation and which might be putative targets of Y-chromosomal regulation.

Results

Genes with mitochondrial annotations and genes representing the mitochondrial proteome do not exhibit statistically-significant biases in chromosomal representation. However, genes exhibiting sex-specific sensitivity to mtDNA are under-represented on the X-chromosome, over-represented among genes known to be sensitive to Y-chromosomal variation, and among genes previously associated with male fitness, but under-represented among genes associated with direct sexual antagonism.

Conclusions

Our results are consistent with the SC hypothesis, suggesting that mitochondrial mutational pressure selects for gene movement off-the-X, hence enabling mito-nuclear coadaptation to proceed along trajectories that result in optimized fitness in both sexes.

How do y-chromosomes modulate genome-wide epigenetic States: genome folding, chromatin sinks, and gene expression

The Y chromosomes of Drosophila melanogaster and D. simulans contain only a handful of protein-coding genes, which are related to sperm mobility and reproductive fitness. Despite low or absent protein coding polymorphism, the Drosophila Y chromosome has been associated with natural phenotypic variation, including variation in the expression of hundreds to thousands of genes located on autosomes and on the X chromosome. Polymorphisms present in the large blocks of heterochromatin and consisting of differences in the amounts and kinds of sequences for satellite DNA and transposable elements may be the source of this modulation. Here we review the evidence and discuss mechanisms for global epigenetic regulation by repetitious elements in the Y chromosome. We also discuss how the discovery of this new function impacts the current knowledge about Y chromosome origin, its current dynamics, and future fate.

Characterizing the Retinoblastoma 1 locus: putative elements for Rb1 regulation by in silico analysis

Limited understanding of the Rb1 locus hinders genetic and epigenetic analyses of Retinoblastoma, a childhood cancer of the nervous systems. In this study, we used in silico tools to investigate and review putative genetic and epigenetic elements of the Rb1 gene. We report transcription start sites, CpG islands, and regulatory moieties that are likely to influence transcriptional states of this gene. These might contribute genetic and epigenetic information modulating tissue-specific transcripts and expression levels of Rb1. The elements we identified include tandem repeats that reside within or next to CpG islands near Rb1's transcriptional start site, and that are likely to be polymorphic among individuals. Our analyses highlight the complexity of this gene and suggest opportunities and limitations for future studies of retinoblastoma, genetic counseling, and the accurate identification of patients at greater risk of developing the malignancy.

Natural variation of the Y chromosome suppresses sex ratio distortion and modulates testis-specific gene expression in Drosophila simulans

X-linked sex-ratio distorters that disrupt spermatogenesis can cause a deficiency in functional Y-bearing sperm and a female-biased sex ratio. Y-linked modifiers that restore a normal sex ratio might be abundant and favored when a X-linked distorter is present. Here we investigated natural variation of Y-linked suppressors of sex-ratio in the Winters systems and the ability of these chromosomes to modulate gene expression in Drosophila simulans. Seventy-eight Y chromosomes of worldwide origin were assayed for their resistance to the X-linked sex-ratio distorter gene Dox. Y chromosome diversity caused males to sire ∼63% to ∼98% female progeny. Genome-wide gene expression analysis revealed hundreds of genes differentially expressed between isogenic males with sensitive (high sex ratio) and resistant (low sex ratio) Y chromosomes from the same population. Although the expression of about 75% of all testis-specific genes remained unchanged across Y chromosomes, a subset of post-meiotic genes was upregulated by resistant Y chromosomes. Conversely, a set of accessory gland-specific genes and mitochondrial genes were downregulated in males with resistant Y chromosomes. The D. simulans Y chromosome also modulated gene expression in XXY females in which the Y-linked protein-coding genes are not transcribed. The data suggest that the Y chromosome might exert its regulatory functions through epigenetic mechanisms that do not require the expression of protein-coding genes. The gene network that modulates sex ratio distortion by the Y chromosome is poorly understood, other than that it might include interactions with mitochondria and enriched for genes expressed in post-meiotic stages of spermatogenesis.

More epigenetic hits than meets the eye: microRNAs and genes associated with the tumorigenesis of retinoblastoma

Retinoblastoma (RB), a childhood neoplasia of the retinoblasts, can occur unilaterally or bilaterally, with one or multiple foci per eye. RB is associated with somatic loss of function of both alleles of the tumor suppressor gene RB1. Hereditary forms emerge due to germline loss of function mutations in RB1 alleles. RB has long been the prototypic “model” cancer ever since Knudson's “two-hit” hypothesis. However, a simple two-hit model for RB is challenged by an increasing number of studies documenting additional hits that contribute to RB development. Here we review the genetics and epigenetics of RB with a focus on the role of small non-coding RNAs (microRNAs) and on novel findings indicating the relevance of DNA methylation in the development and prognosis of this neoplasia. Studies point to an elaborated landscape of genetic and epigenetic complexity, in which a number of events and pahtways play crucial roles in the origin and prognosis of RB. These include roles for microRNAs, inprinted loci, and parent-of-origin contributions to RB1 regulation and RB progression. This complexity is also manifested in the structure of the RB1 locus itself: it includes numerous repetitive DNA segments and retrotransposon insertion elements, some of which are actively transcribed from the RB1 locus. Altogether, we conclude that RB1 loss of function represents the tip of an iceberg of events that determine RB development, progression, severity, and disease risk. Comprehensive assessment of personalized RB risk will require genetic and epigenetic evaluations beyond RB1 protein coding sequences.

Copy-number variation: the balance between gene dosage and expression in Drosophila melanogaster

Copy-number variants (CNVs) reshape gene structure, modulate gene expression, and contribute to significant phenotypic variation. Previous studies have revealed CNV patterns in natural populations of Drosophila melanogaster and suggested that selection and mutational bias shape genomic patterns of CNV. Although previous CNV studies focused on heterogeneous strains, here, we established a number of second-chromosome substitution lines to uncover CNV characteristics when homozygous. The percentage of genes harboring CNVs is higher than found in previous studies. More CNVs are detected in homozygous than heterozygous substitution strains, suggesting the comparative genomic hybridization arrays underestimate CNV owing to heterozygous masking. We incorporated previous gene expression data collected from some of the same substitution lines to investigate relationships between CNV gene dosage and expression. Most genes present in CNVs show no evidence of increased or diminished transcription, and the fraction of such dosage-insensitive CNVs is greater in heterozygotes. More than 70% of the dosage-sensitive CNVs are recessive with undetectable effects on transcription in heterozygotes. A deficiency of singletons in recessive dosage-sensitive CNVs supports the hypothesis that most CNVs are subject to negative selection. On the other hand, relaxed purifying selection might account for the higher number of protein–protein interactions in dosage-insensitive CNVs than in dosage-sensitive CNVs. Dosage-sensitive CNVs that are upregulated and downregulated coincide with copy-number increases and decreases. Our results help clarify the relation between CNV dosage and gene expression in the D. melanogaster genome.

Ribosomal DNA deletions modulate genome-wide gene expression: "rDNA-sensitive" genes and natural variation

The ribosomal rDNA gene array is an epigenetically-regulated repeated gene locus. While rDNA copy number varies widely between and within species, the functional consequences of subtle copy number polymorphisms have been largely unknown. Deletions in the Drosophila Y-linked rDNA modifies heterochromatin-induced position effect variegation (PEV), but it has been unknown if the euchromatic component of the genome is affected by rDNA copy number. Polymorphisms of naturally occurring Y chromosomes affect both euchromatin and heterochromatin, although the elements responsible for these effects are unknown. Here we show that copy number of the Y-linked rDNA array is a source of genome-wide variation in gene expression. Induced deletions in the rDNA affect the expression of hundreds to thousands of euchromatic genes throughout the genome of males and females. Although the affected genes are not physically clustered, we observed functional enrichments for genes whose protein products are located in the mitochondria and are involved in electron transport. The affected genes significantly overlap with genes affected by natural polymorphisms on Y chromosomes, suggesting that polymorphic rDNA copy number is an important determinant of gene expression diversity in natural populations. Altogether, our results indicate that subtle changes to rDNA copy number between individuals may contribute to biologically relevant phenotypic variation.

Fine-scale genetic mapping of a hybrid sterility factor between Drosophila simulans and D. mauritiana: the varied and elusive functions of "speciation genes"

Background

Hybrid male sterility (HMS) is a usual outcome of hybridization between closely related animal species. It arises because interactions between alleles that are functional within one species may be disrupted in hybrids. The identification of genes leading to hybrid sterility is of great interest for understanding the evolutionary process of speciation. In the current work we used marked P-element insertions as dominant markers to efficiently locate one genetic factor causing a severe reduction in fertility in hybrid males of Drosophila simulans and D. mauritiana.

Results

Our mapping effort identified a region of 9 kb on chromosome 3, containing three complete and one partial coding sequences. Within this region, two annotated genes are suggested as candidates for the HMS factor, based on the comparative molecular characterization and public-source information. Gene Taf1 is partially contained in the region, but yet shows high polymorphism with four fixed non-synonymous substitutions between the two species. Its molecular functions involve sequence-specific DNA binding and transcription factor activity. Gene agt is a small, intronless gene, whose molecular function is annotated as methylated-DNA-protein-cysteine S-methyltransferase activity. High polymorphism and one fixed non-synonymous substitution suggest this is a fast evolving gene. The gene trees of both genes perfectly separate D. simulans and D. mauritiana into monophyletic groups. Analysis of gene expression using microarray revealed trends that were similar to those previously found in comparisons between whole-genome hybrids and parental species.

Conclusions

The identification following confirmation of the HMS candidate gene will add another case study leading to understanding the evolutionary process of hybrid incompatibility.

The Opossum genome reveals further evidence for regulatory evolution in mammalian diversification

The first marsupial genome to be sequenced emphasizes the importance of noncoding sequences in mammalian evolution.

The sequencing of the euchromatic genome of a marsupial, the opossum Monodelphis domestica, identifies shared and unique features of marsupial and placental genomes and reveals a prominent role for the evolution of non-protein-coding elements.

Polymorphic Y chromosomes harbor cryptic variation with manifold functional consequences

The paucity of polymorphisms in single-copy genes on the Y chromosome of Drosophila contrasts with data indicating that this chromosome has polymorphic phenotypic effects on sex ratio, temperature sensitivity, behavior, and fitness. We show that the Y chromosome of D. melanogaster harbors substantial genetic diversity in the form of polymorphisms for genetic elements that differentially affect the expression of hundreds of X-linked and autosomal genes. The affected genes are more highly expressed in males, more meagerly expressed in females, and more highly divergent between species. Functionally, they affect microtubule stability, lipid and mitochondrial metabolism, and the thermal sensitivity of spermatogenesis. Our findings provide a mechanism for adaptive phenotypic variation associated with the Y chromosome.

Genetic Properties Influencing the Evolvability of Gene Expression

Identifying the properties of gene networks that influence their evolution is a fundamental research goal. However, modes of evolution cannot be inferred solely from the distribution of natural variation, because selection interacts with demography and mutation rates to shape polymorphism and divergence. We estimated the effects of naturally occurring mutations on gene expression while minimizing the effect of natural selection. We demonstrate that sensitivity of gene expression to mutations increases with both increasing trans-mutational target size and the presence of a TATA box. Genes with greater sensitivity to mutations are also more sensitive to systematic environmental perturbations and stochastic noise. These results provide a mechanistic basis for gene expression evolvability that can serve as a foundation for realistic models of regulatory evolution.

Small mammals of Chapada dos Veadeiros National Park (Cerrado of Central Brazil): ecologic, karyologic, and taxonomic considerations

This work is based on a survey of small mammals carried out in the Chapada dos Veadeiros National Park, a natural reserve located in the mountains of the Planalto Central Goiano in the Cerrado of Central Brazil. The 227 specimens collected represented six marsupial and 13 rodent species. Taxonomic, karyologic, and ecologic considerations are present and discussed in the present work. Our data reflected the faunal heterogeneity with respect to both elevation and vegetation because only eight of the 19 species were collected at both high and low elevations. The composition of the small mammal fauna of the park is influenced by predominance of forest formations at low elevations and cerrado with rupestrian areas at high elevations. Presence of endemic species and one undescribed demonstrated that the cerrado has an endemic fauna and a little known diversity of small mammals.

Male sterility at extreme temperatures: a significant but neglected phenomenon for understanding Drosophila climatic adaptations

The thermal range for viability is quite variable among Drosophila species and it has long been known that these variations are correlated with geographic distribution: temperate species are on average more cold tolerant but more heat sensitive than tropical species. At both ends of their viability range, sterile males have been observed in all species investigated so far. This symmetrical phenomenon restricts the temperature limits within which permanent cultures can be kept in the laboratory. Thermal heat sterility thresholds are very variable across species from 23 degrees C in heat sensitive species up to 31 degrees C in heat tolerant species. In Drosophila melanogaster, genetic variations are observed among geographic populations. Tropical populations are more tolerant to heat induced sterility and recover more rapidly than temperate ones. A genetic analysis revealed that about 50% of the difference observed between natural populations was due to the Y chromosome. Natural populations have not reached a selection limit, however: thermal tolerance was still increased by keeping strains at a high temperature, close to the sterility threshold. On the low temperature side, a symmetrical reverse phenomenon seems to exist: temperate populations are more tolerant to cold than tropical ones. Compared to Mammals, drosophilids exhibit two major differences: first, male sterility occurs not only at high temperature, but also at a low temperature; second, sterility thresholds are not evolutionarily constrained, but highly variable. Altogether, significant and sometimes major genetic variations have been observed between species, between geographic races of the same species, and even between strains kept in the laboratory under different thermal regimes. In each case, it is easily argued that the observed variations correspond to adaptations to climatic conditions, and that male sterility is a significant component of fitness and a target of natural selection.

Evolution of proteins and gene expression levels are coupled in Drosophila and are independently associated with mRNA abundance, protein length, and number of protein-protein interactions

Organismic evolution requires that variation at distinct hierarchical levels and attributes be coherently integrated, often in the face of disparate environmental and genetic pressures. A central part of the evolutionary analysis of biological systems remains to decipher the causal connections between organism-wide (or genome-wide) attributes (e.g., mRNA abundance, protein length, codon bias, recombination rate, genomic position, mutation rate, etc) as well as their role-together with mutation, selection, and genetic drift-in shaping patterns of evolutionary variation in any of the attributes themselves. Here we combine genome-wide evolutionary analysis of protein and gene expression data to highlight fundamental relationships among genomic attributes and their associations with the evolution of both protein sequences and gene expression levels. Our results show that protein divergence is positively coupled with both gene expression polymorphism and divergence. We show moreover that although the number of protein-protein interactions in Drosophila is negatively associated with protein divergence as well as gene expression polymorphism and divergence, protein-protein interactions cannot account for the observed coupling between regulatory and structural evolution. Furthermore, we show that proteins with higher rates of amino acid substitutions tend to have larger sizes and tend to be expressed at lower mRNA abundances, whereas genes with higher levels of gene expression divergence and polymorphism tend to have shorter sizes and tend to be expressed at higher mRNA abundances. Finally, we show that protein length is negatively associated with both number of protein-protein interactions and mRNA abundance and that interacting proteins in Drosophila show similar amounts of divergence. We suggest that protein sequences and gene expression are subjected to similar evolutionary dynamics, possibly because of similarity in the fitness effect (i.e., strength of stabilizing selection) of disruptions in a gene's protein sequence or its mRNA expression. We conclude that, as more and better data accumulate, understanding the causal connections among biological traits and how they are integrated over time to constrain or promote structural and regulatory evolution may finally become possible.

Phylogenetic footprinting reveals extensive conservation of Sonic Hedgehog (SHH) regulatory elements

Sonic Hedgehog (SHH) plays a fundamental role in numerous developmental processes including morphogenesis of limbs, nervous system, and teeth. Using a Bayesian alignment algorithm for phylogenetic footprinting we analyzed approximately 28 kb of noncoding DNA in the SHH locus of human and mouse. This showed that the length of conserved noncoding sequences (4196 nt) shared by these species was approximately 3 times larger than the SHH coding sequence (1386 nt). Most segments were located in introns (53%) or within 2-kb regions upstream (16%) or downstream (20%) of the first and last SHH codon. Even though regions more than 2 kb upstream or downstream of the first and last SHH codon represented 57% (16 kb) of the sequence compared, they accounted for only 11% (494 nt) of the total length of conserved noncoding segments. One region of 650 nt downstream of SHH was identified as a putative scaffold/matrix attachment region (SMAR). Human-mouse analysis was complemented by sequencing in apes, monkeys, rodents, and bats, thus further confirming the evolutionary conservation of some segments. Gel-shift assays indicated that conserved segments are targeted by nuclear proteins and showed differences between two cell types that expressed different levels of SHH, namely human endothelial cells and breast cancer cells. The relevance of these findings with respect to regulation of SHH expression during normal and pathologic development is discussed.

Rates of divergence in gene expression profiles of primates, mice, and flies: stabilizing selection and variability among functional categories

The extent to which natural selection shapes phenotypic variation has long been a matter of debate among those studying organic evolution. We studied the patterns of gene expression polymorphism and divergence in several datasets that ranged from comparisons between two very closely related laboratory strains of mice to comparisons across a considerably longer time scale, such as between humans and chimpanzees, two species of mice, and two species of Drosophila. The results were analyzed and interpreted in view of neutral models of phenotypic evolution. Our analyses used a number of metrics to show that most mRNA levels are evolutionary stable, changing little across the range of taxonomic distances compared. This implies that, overall, widespread stabilizing selection on transcription levels has prevented greater evolutionary changes in mRNA levels. Nevertheless, the range of rates of divergence is large with highly significant differences in the rate and patterns of transcription divergence across functional classes defined on the basis of the gene ontology annotation (primates and mice datasets) or on the basis of the pattern of sex-biased gene expression (Drosophila). Moreover, rates of divergence of sex-biased genes in the contrast between Drosophila species show a distinct pattern from that observed in the contrast between populations of D. melanogaster. Hence, we discuss the time scale of the changes observed and its consequences for the relationship between variation in gene expression within and between species. Finally, we argue that differences in mRNA levels of the magnitudes observed herein could be explained by a remarkably small number of generations of directional selection.

The lowest diploid number (2n = 16) yet found in any primate: Callicebus lugens (Humboldt, 1811)

Morphologic, molecular and karyologic analyses of Callicebus lugens (Humboldt, 1811) of known geographic origin supported the proposition that this is a valid species. Morphologic and morphometric analyses showed evident differences between C. lugens and two other related taxa of the same group (Callicebus purinus and Callicebus torquatus). Cytochrome b DNA analyses (maximum parsimony, neighbour joining and maximum likelihood) were congruent in showing a strong association between C. lugens and Callicebus sp. of the torquatus group in one branch and a sister branch further divided into two clades: one with species of the personatus group and another, with species of the moloch group. Karyotypic analysis showed that C. lugens has the lowest diploid chromosome number of the primate order (2n = 16). Comparisons with other congeneric species clearly supported the proposition that C. lugens is karyotypically similar to others of the torquatus group.

Gene assignment in Ateles paniscus chamek (Platyrrhini, Primates). Allocation of 18 markers of human syntenic groups 1, 2, 7, 14, 15, 17 and 22

Eighteen markers allocated to human syntenic groups 1, 2, 7, 14, 15, 17 and 22 were assigned to the chromosome complement of the neotropical primate Ateles paniscus chamek. These new allocations and existing gene charts in this species were compared with chromosome painting patterns produced by human chromosome probes in the congeneric species A teles geoffroyi and with available data on the human genome and gene mapping. These comparisons showed congruent findings in Ateles and provided good evidence of how several human syntenic groups were evolutionarily rearranged.