Dynamics and function of intron sequences of the wingless gene during the evolution of the Drosophila genus

Conserved Critical Evolutionary Gene Structures in Orthologs

Unravelling gene structure requires the identification and understanding of the constraints that are often associated with the evolutionary history and functional domains of genes. We speculated in this manuscript with the possibility of the existence in orthologs of an emergent highly conserved gene structure that might explain their coordinated evolution during speciation events and their parental function. Here, we will address the following issues: (1) is there any conserved hypothetical structure along ortholog gene sequences? (2) If any, are such conserved structures maintained and conserved during speciation events? The data presented show evidences supporting this hypothesis. We have found that, (1) most orthologs studied share highly conserved compositional structures not observed previously. (2) While the percent identity of nucleotide sequences of orthologs correlates with the percent identity of composon sequences, the number of emergent compositional structures conserved during speciation does not correlate with the percent identity. (3) A broad range of species conserves the emergent compositional stretches. We will also discuss the concept of critical gene structure.

Patterns of exon-intron architecture variation of genes in eukaryotic genomes

Background

The origin and importance of exon-intron architecture comprises one of the remaining mysteries of gene evolution. Several studies have investigated the variations of intron length, GC content, ordinal position in a gene and divergence. However, there is little study about the structural variation of exons and introns.

Results

We investigated the length, GC content, ordinal position and divergence in both exons and introns of 13 eukaryotic genomes, representing plant and animal. Our analyses revealed that three basic patterns of exon-intron variation were present in nearly all analyzed genomes (P < 0.001 in most cases): an ordinal reduction of length and divergence in both exon and intron, a co-variation between exon and its flanking introns in their length, GC content and divergence, and a decrease of average exon (or intron) length, GC content and divergence as the total exon numbers of a gene increased. In addition, we observed that the shorter introns had either low or high GC content, and the GC content of long introns was intermediate.

Conclusion

Although the factors contributing to these patterns have not been identified, our results provide three important clues: common factor(s) exist and may shape both exons and introns; the ordinal reduction patterns may reflect a time-orderly evolution; and the larger first and last exons may be splicing-required. These clues provide a framework for elucidating mechanisms involved in the organization of eukaryotic genomes and particularly in building exon-intron structures.

Regulation of the feedback antagonist naked cuticle by Wingless signaling

Signaling pathways usually activate transcriptional targets in a cell type-specific manner. Notable exceptions are pathway-specific feedback antagonists, which serve to restrict the range or duration of the signal. These factors are often activated by their respective pathways in a broad array of cell types. For example, the Wnt ligand Wingless (Wg) activates the naked cuticle (nkd) gene in all tissues examined throughout Drosophila development. How does the nkd gene respond in such an unrestricted manner to Wg signaling? Analysis in cell culture revealed regions of the nkd locus that contain Wg response elements (WREs) that are directly activated by the pathway via the transcription factor TCF. In flies, Wg signaling activates these WREs in multiple tissues, in distinct but overlapping patterns. These WREs are necessary and largely sufficient for nkd expression in late stage larval tissues, but only contribute to part of the embryonic expression pattern of nkd. These results demonstrate that nkd responsiveness to Wg signaling is achieved by several WREs which are broadly (but not universally) activated by the pathway. The existence of several WREs in the nkd locus may have been necessary to allow the Wg signaling-Nkd feedback circuit to remain intact as Wg expression diversified during animal evolution.

Detecting the limits of regulatory element conservation and divergence estimation using pairwise and multiple alignments

Background

Molecular evolutionary studies of noncoding sequences rely on multiple alignments. Yet how multiple alignment accuracy varies across sequence types, tree topologies, divergences and tools, and further how this variation impacts specific inferences, remains unclear.

Results

Here we develop a molecular evolution simulation platform, CisEvolver, with models of background noncoding and transcription factor binding site evolution, and use simulated alignments to systematically examine multiple alignment accuracy and its impact on two key molecular evolutionary inferences: transcription factor binding site conservation and divergence estimation. We find that the accuracy of multiple alignments is determined almost exclusively by the pairwise divergence distance of the two most diverged species and that additional species have a negligible influence on alignment accuracy. Conserved transcription factor binding sites align better than surrounding noncoding DNA yet are often found to be misaligned at relatively short divergence distances, such that studies of binding site gain and loss could easily be confounded by alignment error. Divergence estimates from multiple alignments tend to be overestimated at short divergence distances but reach a tool specific divergence at which they cease to increase, leading to underestimation at long divergences. Our most striking finding was that overall alignment accuracy, binding site alignment accuracy and divergence estimation accuracy vary greatly across branches in a tree and are most accurate for terminal branches connecting sister taxa and least accurate for internal branches connecting sub-alignments.

Conclusion

Our results suggest that variation in alignment accuracy can lead to errors in molecular evolutionary inferences that could be construed as biological variation. These findings have implications for which species to choose for analyses, what kind of errors would be expected for a given set of species and how multiple alignment tools and phylogenetic inference methods might be improved to minimize or control for alignment errors.

Evolution of cis-regulatory sequence and function in Diptera

Cis-regulatory sequences direct patterns of gene expression essential for development and physiology. Evolutionary changes in these sequences contribute to phenotypic divergence. Despite their importance, cis-regulatory regions remain one of the most enigmatic features of the genome. Patterns of sequence evolution can be used to identify cis-regulatory elements, but the power of this approach depends upon the relationship between sequence and function. Comparative studies of gene regulation among Diptera reveal that divergent sequences can underlie conserved expression, and that expression differences can evolve despite largely similar sequences. This complex structure-function relationship is the primary impediment for computational identification and interpretation of cis-regulatory sequences. Biochemical characterization and in vivo assays of cis-regulatory sequences on a genomic-scale will relieve this barrier.

A 3' cis-regulatory region controls wingless expression in the Drosophila eye and leg primordia

The precise regulation of wingless (wg) expression in the Drosophila eye disc is key to control the anteroposterior and dorsoventral patterning of this disc. Here, we identify an eye disc-specific wg cis-regulatory element that functions as a regulatory rheostat. Pannier (Pnr), a transcription factor previously proposed to act as an upstream activator of wg, is sufficient to activate the eye disc enhancer but required for wg expression only in the peripodial epithelium of the disc. We propose that this regulation of wg by Pnr appeared associated to the development of the peripodial epithelium in higher dipterans and was added to an existing mechanism regulating the deployment of wingless in the dorsal region of the eye primordium. In addition, our analysis identifies a separate ventral disc enhancer that lies adjacent to the eye-specific one, and thus altogether, they define a 1-kb genomic region where disc-specific enhancers of the wg gene are located.