ARTIFICIAL SELECTION ON EGG SIZE PERTURBS EARLY PATTERN FORMATION IN DROSOPHILA MELANOGASTER

Pattern formation in Drosophila embryogenesis has been widely investigated as a developmental and evolutionary model of robustness. To ask whether genetic variation for pattern formation is suppressed in this system, artificial selection for divergent egg size was used to challenge the scaling of even‐skipped (eve) pattern formation in mitotic cycle 14 (stage 5) embryos of Drosophila melanogaster. Three‐dimensional confocal imaging revealed shifts in the allometry of eve pair‐rule stripes along both anterior–posterior (A–P) and dorsoventral (D–V) axes as a correlated response to egg size selection, indicating the availability of genetic variation for this buffered trait. Environmental perturbation was not required for the manifestation of this variation. The number of nuclei at the cellular blastoderm stage also changed in response to selection, with large‐egg selected lines having more than 1000 additional nuclei relative to small‐egg lines. This increase in nuclear number in larger eggs does not scale with egg size, however, as nuclear density is inversely correlated with egg length. Nuclear density varies along the A–P axis but does not correlate with the shift in eve stripe allometry between the selection treatments. Despite its macroevolutionary conservation, both eve stripe patterning and blastoderm cell number vary genetically both within and between closely related species.     

The Drosophila eve Insulator Homie Promotes eve Expression and Protects the Adjacent Gene from Repression by Polycomb Spreading

Insulators can block the action of enhancers on promoters and the spreading of repressive chromatin, as well as facilitating specific enhancer-promoter interactions. However, recent studies have called into question whether the activities ascribed to insulators in model transgene assays actually reflect their functions in the genome. The Drosophila even skipped (eve) gene is a Polycomb (Pc) domain with a Pc-group response element (PRE) at one end, flanked by an insulator, an arrangement also seen in other genes. Here, we show that this insulator has three major functions. It blocks the spreading of the eve Pc domain, preventing repression of the adjacent gene, TER94. It prevents activation of TER94 by eve regulatory DNA. It also facilitates normal eve expression. When Homie is deleted in the context of a large transgene that mimics both eve and TER94 regulation, TER94 is repressed. This repression depends on the eve PRE. Ubiquitous TER94 expression is “replaced” by expression in an eve pattern when Homie is deleted, and this effect is reversed when the PRE is also removed. Repression of TER94 is attributable to spreading of the eve Pc domain into the TER94 locus, accompanied by an increase in histone H3 trimethylation at lysine 27. Other PREs can functionally replace the eve PRE, and other insulators can block PRE-dependent repression in this context. The full activity of the eve promoter is also dependent on Homie, and other insulators can promote normal eve enhancer-promoter communication. Our data suggest that this is not due to preventing promoter competition, but is likely the result of the insulator organizing a chromosomal conformation favorable to normal enhancer-promoter interactions. Thus, insulator activities in a native context include enhancer blocking and enhancer-promoter facilitation, as well as preventing the spread of repressive chromatin.     

Conserved boundary elements from the Hox complex of mosquito,Anopheles gambiae

The conservation of hox genes as well as their genomic organization across the phyla suggests that this system of anterior–posterior axis formation arose early during evolution and has come under strong selection pressure. Studies in the split Hox cluster of Drosophila have shown that proper expression of hox genes is dependent on chromatin domain boundaries that prevent inappropriate interactions among different types of cis-regulatory elements. To investigate whether boundary function and their role in regulation of hox genes is conserved in insects with intact Hox clusters, we used an algorithm to locate potential boundary elements in the Hox complex of mosquito, Anopheles gambiae. Several potential boundary elements were identified that could be tested for their functional conservation. Comparative analysis revealed that like Drosophila, the bithorax region in A. gambiae contains an extensive array of boundaries and enhancers organized into domains. We analysed a subset of candidate boundary elements and show that they function as enhancer blockers in Drosophila. The functional conservation of boundary elements from mosquito in fly suggests that regulation of hox genes involving chromatin domain boundaries is an evolutionary conserved mechanism and points to an important role of such elements in key developmentally regulated loci.     

Organization of developmental enhancers in the Drosophila embryo

Most cell-specific enhancers are thought to lack an inherent organization, with critical binding sites distributed in a more or less random fashion. However, there are examples of fixed arrangements of binding sites, such as helical phasing, that promote the formation of higher-order protein complexes on the enhancer DNA template. Here, we investigate the regulatory ‘grammar’ of nearly 100 characterized enhancers for developmental control genes active in the early Drosophila embryo. The conservation of grammar is examined in seven divergent Drosophila genomes. Linked binding sites are observed for particular combinations of binding motifs, including Bicoid–Bicoid, Hunchback–Hunchback, Bicoid–Dorsal, Bicoid–Caudal and Dorsal–Twist. Direct evidence is presented for the importance of Bicoid–Dorsal linkage in the integration of the anterior–posterior and dorsal–ventral patterning systems. Hunchback–Hunchback interactions help explain unresolved aspects of segmentation, including the differential regulation of the eve stripe 3 + 7 and stripe 4 + 6 enhancers. We also present evidence that there is an under-representation of nucleosome positioning sequences in many enhancers, raising the possibility for a subtle higher-order structure extending across certain enhancers. We conclude that grammar of gene control regions is pervasively used in the patterning of the Drosophila embryo.     

Molecular Phylogeny and Functional Genomics of β-Galactoside α2,6-Sialyltransferases That Explain Ubiquitous Expression of st6gal1 Gene in Amniotes

Sialyltransferases are key enzymes in the biosynthesis of sialoglycoconjugates that catalyze the transfer of sialic residue from its activated form to an oligosaccharidic acceptor. β-Galactoside α2,6-sialyltransferases ST6Gal I and ST6Gal II are the two unique members of the ST6Gal family described in higher vertebrates. The availability of genome sequences enabled the identification of more distantly related invertebrates'' st6gal gene sequences and allowed us to propose a scenario of their evolution. Using a phylogenomic approach, we present further evidence of an accelerated evolution of the st6gal1 genes both in their genomic regulatory sequences and in their coding sequence in reptiles, birds, and mammals known as amniotes, whereas st6gal2 genes conserve an ancestral profile of expression throughout vertebrate evolution.     

CRISPR/Cas9 knockout of female-biased genes AeAct-4 or myo-fem in Ae. aegypti results in a flightless phenotype in female,but not male mosquitoes

Aedes aegypti is a vector of dengue, chikungunya, and Zika viruses. Current vector control strategies such as community engagement, source reduction, and insecticides have not been sufficient to prevent viral outbreaks. Thus, interest in novel strategies involving genetic engineering is growing. Female mosquitoes rely on flight to mate with males and obtain a bloodmeal from a host. We hypothesized that knockout of genes specifically expressed in female mosquitoes associated with the indirect flight muscles would result in a flightless female mosquito. Using CRISPR-Cas9 we generated loss-of-function mutations in several genes hypothesized to control flight in mosquitoes, including actin (AeAct-4) and myosin (myo-fem) genes expressed specifically in the female flight muscle. Genetic knockout of these genes resulted in 100% flightless females, with homozygous males able to fly, mate, and produce offspring, albeit at a reduced rate when compared to wild type males. Interestingly, we found that while AeAct-4 was haplosufficient, with most heterozygous individuals capable of flight, this was not the case for myo-fem, where about half of individuals carrying only one intact copy could not fly. These findings lay the groundwork for developing novel mechanisms of controlling Ae. aegypti populations, and our results suggest that this mechanism could be applicable to other vector species of mosquito.     

ZBED Evolution: Repeated Utilization of DNA Transposons as Regulators of Diverse Host Functions

ZBED genes originate from domesticated hAT DNA transposons and encode regulatory proteins of diverse function in vertebrates. Here we reveal the evolutionary relationship between ZBED genes and demonstrate that they are derived from at least two independent domestication events in jawed vertebrate ancestors. We show that ZBEDs form two monophyletic clades, one of which has expanded through several independent duplications in host lineages. Subsequent diversification of ZBED genes has facilitated regulation of multiple diverse fundamental functions. In contrast to known examples of transposable element exaptation, our results demonstrate a novel unprecedented capacity for the repeated utilization of a family of transposable element-derived protein domains sequestered as regulators during the evolution of diverse host gene functions in vertebrates. Specifically, ZBEDs have contributed to vertebrate regulatory innovation through the donation of modular DNA and protein interacting domains. We identify that C7ORF29, ZBED2, 3, 4, and ZBEDX form a monophyletic group together with ZBED6, that is distinct from ZBED1 genes. Furthermore, we show that ZBED5 is related to Buster DNA transposons and is phylogenetically separate from other ZBEDs. Our results offer new insights into the evolution of regulatory pathways, and suggest that DNA transposons have contributed to regulatory complexity during genome evolution in vertebrates.     

Evolution of Stress-Regulated Gene Expression in Duplicate Genes of Arabidopsis thaliana

Due to the selection pressure imposed by highly variable environmental conditions, stress sensing and regulatory response mechanisms in plants are expected to evolve rapidly. One potential source of innovation in plant stress response mechanisms is gene duplication. In this study, we examined the evolution of stress-regulated gene expression among duplicated genes in the model plant Arabidopsis thaliana. Key to this analysis was reconstructing the putative ancestral stress regulation pattern. By comparing the expression patterns of duplicated genes with the patterns of their ancestors, duplicated genes likely lost and gained stress responses at a rapid rate initially, but the rate is close to zero when the synonymous substitution rate (a proxy for time) is >~0.8. When considering duplicated gene pairs, we found that partitioning of putative ancestral stress responses occurred more frequently compared to cases of parallel retention and loss. Furthermore, the pattern of stress response partitioning was extremely asymmetric. An analysis of putative cis-acting DNA regulatory elements in the promoters of the duplicated stress-regulated genes indicated that the asymmetric partitioning of ancestral stress responses are likely due, at least in part, to differential loss of DNA regulatory elements; the duplicated genes losing most of their stress responses were those that had lost more of the putative cis-acting elements. Finally, duplicate genes that lost most or all of the ancestral responses are more likely to have gained responses to other stresses. Therefore, the retention of duplicates that inherit few or no functions seems to be coupled to neofunctionalization. Taken together, our findings provide new insight into the patterns of evolutionary changes in gene stress responses after duplication and lay the foundation for testing the adaptive significance of stress regulatory changes under highly variable biotic and abiotic environments.     

Evolutionary divergence of the paralogs Methoprene tolerant (Met) and germ cell expressed (gce) within the genus Drosophila

Juvenile hormone (JH) signaling underpins both regulatory and developmental pathways in insects. However, the JH receptor is poorly understood. Methoprene tolerant (Met) and germ cell expressed (gce) have been implicated in JH signaling in Drosophila. We investigated the evolution of Met and gce across 12 Drosophila species and found that these paralogs are conserved across at least 63 million years of dipteran evolution. Distinct patterns of selection found using estimates of dN/dS ratios across Drosophila Met and gce coding sequences, along with their incongruent temporal expression profiles in embryonic Drosophila melanogaster, illustrate avenues through which these genes have diverged within the Diptera. Additionally, we demonstrate that the annotated gene CG15032 is the 5′ terminus of gce.In mosquitoes and beetles, a single Met-like homolog displays structural similarity to both Met and gce, and the intron locations are conserved with those of gce. We found that Tribolium and mosquito Met orthologs are assembled from Met- and gce-specific domains in a modular fashion. Our results suggest that Drosophila Met and gce experienced divergent evolutionary pressures following the duplication of an ancestral gce-like gene found in less derived holometabolous insects.     

Genetic regulation of engrailed and wingless in Tribolium segmentation and the evolution of pair-rule segmentation

In Drosophila, primary pair-rule genes establish the parasegmental boundaries and indirectly control the periodic expression of the segment polarity genes engrailed (en) and wingless (wg) via regulation of secondary pair-rule genes. Although orthologs of some Drosophila pair-rule genes are not required for proper segmentation in Tribolium, segmental expression of Tc-en and Tc-wg is conserved. To understand how these segment polarity genes are regulated, we examined the results of expressing one or two pair-rule genes in the absence of the other known pair-rule genes. Expression of one or both of the secondary pair-rule genes, Tc-sloppy-paired (Tc-slp) and Tc-paired (Tc-prd), activated Tc-wg in the absence of the primary pair-rule genes, Tc-even-skipped (Tc-eve), Tc-runt (Tc-run) and Tc-odd-skipped (Tc-odd). Tc-eve alone failed to activate Tc-wg or Tc-en, but in combination with Tc-run or Tc-prd activated Tc-en. These results, interpreted within the pair-rule gene expression patterns, suggest separate models for the genetic regulation of the juxtaposed expression of Tc-wg and Tc-en at odd- and even-numbered parasegmental boundaries, respectively. Conserved interactions between eve and prd at the anterior boundary of odd-numbered parasegments may reflect an ancestral segmentation mechanism that functioned in every segment prior to the evolution of pair-rule segmentation.     

CpG site degeneration triggered by the loss of functional constraint created a highly polymorphic macaque drug-metabolizing gene, CYP1A2

Background

Duplicated genes frequently experience asymmetric rates of sequence evolution. Relaxed selective constraints and positive selection have both been invoked to explain the observation that one paralog within a gene-duplicate pair exhibits an accelerated rate of sequence evolution. In the majority of studies where asymmetric divergence has been established, there is no indication as to which gene copy, ancestral or derived, is evolving more rapidly. In this study we investigated the effect of local synteny (gene-neighborhood conservation) and codon usage on the sequence evolution of gene duplicates in the S. cerevisiae genome. We further distinguish the gene duplicates into those that originated from a whole-genome duplication (WGD) event (ohnologs) versus small-scale duplications (SSD) to determine if there exist any differences in their patterns of sequence evolution.

Results

For SSD pairs, the derived copy evolves faster than the ancestral copy. However, there is no relationship between rate asymmetry and synteny conservation (ancestral-like versus derived-like) in ohnologs. mRNA abundance and optimal codon usage as measured by the CAI is lower in the derived SSD copies relative to ancestral paralogs. Moreover, in the case of ohnologs, the faster-evolving copy has lower CAI and lowered expression.

Conclusions

Together, these results suggest that relaxation of selection for codon usage and gene expression contribute to rate asymmetry in the evolution of duplicated genes and that in SSD pairs, the relaxation of selection stems from the loss of ancestral regulatory information in the derived copy.     

Embryonic expression of the single Tribolium engrailed homolog

We have cloned and sequenced the single Tribolium homolog of the Drosophila engrailed gene. The predicted protein contains a homeobox and several domains conserved among all engrailed genes identified to date. In addition it contains several features specific to the invected homologs of Bombyx and Drosophila, indicating that these features most likely were present in the ancestral gene in the common ancestor of holometabolous insects. We used the cross-reacting monoclonal antibody, 4D9, to follow the expression of the Engrailed protein during segmentation in Tribolium embryos. As in other insects, Engrailed accumulates in the nuclei of cells along the posterior margin of each segment. The first Engrailed stripe appears as the embryonic rudiment condenses. Then as the rudiment elongates into a germ band, Engrailed stripes appear in an anterior to posterior progression, just prior to morphological evidence of the formation of each segment. As in Drosophila (a long germ insect), expression of engrailed in Tribolium (classified as a short germ insect) is preceeded by the expression of several homologous segmentation genes, suggesting that similar genetic regulatory mechanisms are shared by diverse developmental types. © 1994 Wiley-Liss, Inc.     

Mutations in somePolycomb group genes ofDrosophila interfere with regulation of segmentation genes

Mutations in severalPolycomb (Pc) group genes cause maternal-effect or zygotic segmentation defects, suggesting thatPc group genes may regulate the segmentation genes ofDrosophila. We show that individuals doubly heterozygous for mutations inpolyhomeotic and six otherPc group genes show gap, pair rule, and segment polarity segmentation defects. We examined double heterozygous combinations ofPc group and segmentation mutations for enhancement of adult and embryonic segmentation defects.Posterior sex combs andpolyhomeotic interact withKrüppel ² and enhance embryonic phenotypes ofhunchback andknirps, andpolyhomeotic enhanceseven-skipped. Surprisingly, flies carrying duplications ofextra sex combs (esc), that were heterozygous for mutations ofeven-skipped (eve), were extremely subvital. Embryos and surviving adults of this genotype showed strong segmentation defects in even-numbered segments. Antibody studies confirm that expression ofeve is suppressed by duplications ofesc. However,esc duplications have no effect on other gap or pair rule genes tested. To our knowledge, this is only the second triplo-abnormal phenotype associated withPc group genes. Duplications of nine otherPc group genes have no detectable effect oneve. Expression ofengrailed (en) was abnormal in the central nervous systems of mostPc group mutants. These results support a role forPc genes in regulation of some segmentation genes, and suggest thatesc may act differently from otherPc group genes.