Gene organization of the mating type regions in the ectomycorrhizal fungus Laccaria bicolor reveals distinct evolution between the two mating type loci

In natural conditions, basidiomycete ectomycorrhizal fungi such as Laccaria bicolor are typically in the dikaryotic state when forming symbioses with trees, meaning that two genetically different individuals have to fuse or 'mate'. Nevertheless, nothing is known about the molecular mechanisms of mating in these ecologically important fungi. Here, advantage was taken of the first sequenced genome of the ectomycorrhizal fungus, Laccaria bicolor, to determine the genes that govern the establishment of cell-type identity and orchestrate mating. The L. bicolor mating type loci were identified through genomic screening. The evolutionary history of the genomic regions that contained them was determined by genome-wide comparison of L. bicolor sequences with those of known tetrapolar and bipolar basidiomycete species, and by phylogenetic reconstruction of gene family history. It is shown that the genes of the two mating type loci, A and B, are conserved across the Agaricales, but they are contained in regions of the genome with different evolutionary histories. The A locus is in a region where the gene order is under strong selection across the Agaricales. By contrast, the B locus is in a region where the gene order is likely under a low selection pressure but where gene duplication, translocation and transposon insertion are frequent.

The murine interleukin-2 receptor gamma chain gene: organization, chromosomal localization and expression in the adult thymus.

Defects in the interleukin-2 receptor gamma (IL-2R gamma) chain in the man result in an X-linked severe combined immunodeficiency, SCIDX1, characterized by an absence of T-cell differentiation. This phenotype may result from pertubations in IL-2, IL-4-, IL-7- or IL-15-mediated signaling, as the IL-2R gamma chain forms an integral component of these receptor systems. We have isolated and characterized cDNA and genomic clones for the murine IL-2R gamma. The gene (Il2rg) is well conserved between mouse and man with respect to overall structure and size, and contains regions of high conservation in the promoter region as well. Il2rg maps to mouse X chromosome region 40, in a region of synteny with human Xq12-13.1. We have also explored the expression of the IL-2R gamma during thymocyte development. IL-2R gamma transcripts are detected in the earliest thymocyte precursor cells and persist throughout intrathymic development into the mature peripheral compartment. Genomic clones for the murine IL-2R gamma will allow for further studies on the regulation and function of this gene in vivo.

Organization of lin genes and IS6100 among different strains of hexachlorocyclohexane-degrading Sphingomonas paucimobilis: evidence for horizontal gene transfer.

The organization of lin genes and IS6100 was studied in three strains of Sphingomonas paucimobilis (B90A, Sp+, and UT26) which degraded hexachlorocyclohexane (HCH) isomers but which had been isolated at different geographical locations. DNA-DNA hybridization data revealed that most of the lin genes in these strains were associated with IS6100, an insertion sequence classified in the IS6 family and initially found in Mycobacterium fortuitum. Eleven, six, and five copies of IS6100 were detected in B90A, Sp+, and UT26, respectively. IS6100 elements in B90A were sequenced from five, one, and one regions of the genomes of B90A, Sp+, and UT26, respectively, and were found to be identical. DNA-DNA hybridization and DNA sequencing of cosmid clones also revealed that S. paucimobilis B90A contains three and two copies of linX and linA, respectively, compared to only one copy of these genes in strains Sp+ and UT26. Although the copy number and the sequence of the remaining genes of the HCH degradative pathway (linB, linC, linD, and linE) were nearly the same in all strains, there were striking differences in the organization of the linA genes as a result of replacement of portions of DNA sequences by IS6100, which gave them a strange mosaic configuration. Spontaneous deletion of linD and linE from B90A and of linA from Sp+ occurred and was associated either with deletion of a copy of IS6100 or changes in IS6100 profiles. The evidence gathered in this study, coupled with the observation that the G+C contents of the linA genes are lower than that of the remaining DNA sequence of S. paucimobilis, strongly suggests that all these strains acquired the linA gene through horizontal gene transfer mediated by IS6100. The association of IS6100 with the rest of the lin genes further suggests that IS6100 played a role in shaping the current lin gene organization.

A dual functional origin of transfer in the ICEclc genomic island of Pseudomonas knackmussii B13.

Genomic islands (GEIs) are large DNA segments, present in most bacterial genomes, that are most likely acquired via horizontal gene transfer. Here, we study the self-transfer system of the integrative and conjugative element ICEclc of Pseudomonas knackmussii B13, which stands model for a larger group of ICE/GEI with syntenic core gene organization. Functional screening revealed that unlike conjugative plasmids and other ICEs ICEclc carries two separate origins of transfer, with different sequence context but containing a similar repeat motif. Conjugation experiments with GFP-labelled ICEclc variants showed that both oriTs are used for transfer and with indistinguishable efficiencies, but that having two oriTs results in an estimated fourfold increase of ICEclc transfer rates in a population compared with having a single oriT. A gene for a relaxase essential for ICEclc transfer was also identified, but in vivo strand exchange assays suggested that the relaxase processes both oriTs in a different manner. This unique dual origin of transfer system might have provided an evolutionary advantage for distribution of ICE, a hypothesis that is supported by the fact that both oriT regions are conserved in several GEIs related to ICEclc.

Regulation of bistable ICEclc transfer in pseudomonas knackmussii B13

L'élément génétique intégratif et conjugatif auto-transférable de 103 kb qui se trouve dans le génome de Pseudomonas knackmussii B13 (ICEc/c) confère la capacité de dégrader le 3-chlorobenzoate et le 2-aminophénol. L'élément ICE c/c peut être transféré par conjugaison de la souche B13 à diverses bêta- et gamma- protéobactéries. Seule une sous-population de 3 à 5% des cellules transfère l'élément, les cellules dites "compétentes pour le transfert". L'acquisition de la compétence pour le transfert est vraisemblablement la conséquence d'une régulation bistable, conduisant une partie des cellules au transfert de l'élément ICE c/c tandis que, dans les autres, l'élément reste quiescent et ne se transfère pas. À ce jour, les mécanismes et les acteurs moléculaires qui régulent l'activation bistable de l'élément sont restés inconnus. Mon travail de doctorat visait à identifier les éléments bistables du régulon de la compétence pour le transfert et d'analyser les fondements moléculaires de la bistabilité de l'élément ICE c/c chez P. knackmussii. Le premier chapitre introduit le thème du transfert génétique horizontal avec un accent particulier sur les éléments intégratifs et conjugatifs (ICE) et ICEcIc. L'état actuel des connaissances sur l'organisation génétique, la régulation, l'intégration et le transfert de différents modèles de ICEs est exposé en détail. En outre, je m'étends sur les phénomènes d'hétérogénéité et de bistabilité phénotyplques, qu'on peut distinguer dans une population isogénique dans des conditions de culture homogènes, et qui sont susceptibles de jouer un rôle dans le transfert de l'élément ICE c/c, dans la mesure où il ne s'active et n'est transférable que dans une très petite sous-population de cellules. Dans le chapitre 2, je présente une analyse globale des régions promotrices minimales des gènes appartenant au régulon de la compétence pour le transfert de l'élément ICE c/c. Nous avons étudié les caractéristiques d'expression des promoteurs et, s'ils s'avéraient bistables, leur activation dans le temps par comparaison avec le mutant lntB13. Pour ce faire, nous avons utilisé des fusions de promoteurs avec des gènes rapporteurs et testé l'expression bistable chez P. knackmussii par microscopie à épifluorescence. Pour six promoteurs présentant une expression bistable, nous avons employé de la microscopie temporelle pour déterminer la chronologie de leur expression par rapport à Pint et PinR. Parmi eux, nous avons identifié deux gènes exprimés précocement et trois gènes exprimés tardivement dans le processus d'acquisition de la compétence de transfert. Dans le chapitre 3, j'expose une analyse d'expression génétique pour l'un des groupes de gènes dont la transcription est la plus élevée dans la région conservée de ICE c/c, les gènes orf81655-orf68241 contenus dans une région de 14 kb. Nous montrons d'abord que cet opéron fait partie du même régulon bistable que intB13 et inrR et analysons les caractéristiques génétiques qui conduisent à une transcription élevée. Nous étudions les fonctions biologiques de ce groupe de gènes par des délétlons ciblées et montrons que certaines d'entre elles empêchent le transfert de l'élément. Nous approfondissons la caractérlsatlon de I'orf8l655 en construisant une fusion transcrlptionnelle avec le gène codant pour la protéine fluorescente verte (egfp) (en utilisant le système minl-Tn5). L'expression de Vorf81655 dans des cellules individuelles est comparée au signal mesuré par hybridation in situ en fluorescence (FISH) sur le ARN messager du gène. En utilisant FISH, des délétlons du promoteur et de l'analyse directe de transcription, nous avons localisé la région promotrice du groupe de gènes. En outre, nous avons utilisé des mutations dirigées pour comprendre la bistabilité de cette région promotrice, caractérisée par une transcription très élevée et une traduction lente de l'ARN messager.  Dans le chapitre 4, nous nous efforçons de comprendre comment la bistabilité est générée au sein du régulon te de l'élément ICE c/c. Pour ce faire, nous avons tenté de reconstituer une expression bistable, dans un hôte qui ne présente pas de bistabilité naturellement, à partir d'éléments génétiques individuels. L'hôte choisi est Pseudomonas putida dans lequel nous avons introduit une copie unique de Pint, PinR ou PaipA fusionnés à la egfp, construits qui permettent d'observer l'apparition de bistabilité. Nous avons ensuite construit différents assemblages de composants génétiques de l'élément ICE c/c, en nous concentrant sur la région parA-inrR. En effet, nous avons pu démontrer qu'une expression bistable apparaît dans P. putida grâce à ces éléments en l'absence de l'élément ICE c/c complet. À noter que la plupart des construits génétiques activent PaipA ou P|,,R, mais qu'un seul recrée la bistabilité de Pint, ce qui suggère que la région parA-inrR permet à la fois d'engendrer la bistabilité et d'opérer la transition entre les promoteurs précoces et les promoteurs tardifs du régulon de la bistabilité. Dans le chapitre 5, nous concluons sur une discussion de la pertinence de nos résultats et sur de futures perspectives de recherche. -- The 103-kb self-transmissible integrative and conjugative element (ICE) of Pseudomonas knackmussii B13 (ICEc/c) confers the capacity to degrade 3- chlorobenzoate and 2-aminophenol. ICEc/c can be conjugated from strain B13 to a variety of Beta- and Gammaproteobacteria. Interestingly, ICE c/c transfer is observed in a subpopulatlon of cells (3-5%) only, the so-called 'transfer competent' cells. The formation of transfer competence (tc) is thought to be the consequence of a 'bistable' decision, which forces those cells to follow the developmental path which leads to ICEc/c transfer, whereas in others ICE c/c remains silent and does not transfer. So far, the mechanisms and molecular partners generating this bistable transfer activation in cells of P. knackmussii B13 remain mostly unidentified. This thesis aimed at understanding the extent of the tc bistability regulon and to dissect the molecular basis of bistabillty formation of ICEc/c in P. knackmussii. The first chapter is a general Introduction on horizontal gene transfer (HGT) with particular emphasis on ICEs and ICE c/c. The emphasis is made on the current knowledge about the HGT gene organization, regulation and specific integration and transfer aspects of the different ICEs models. Furthermore, I focus on the phenomena of phenotypic heterogeneity and bistability (the property of two distinguishable phenotypes existing within an isogenic population under homogeneous conditions), which may play a particular role in ICEc/c behaviour, since ICE activation and transfer only occurs in a very small subpopulation of cells. In Chapter Two, I focus on a global analysis of the different core promoters that might belong to the ICEc/c tc pathway regulon. We studied both expression patterns of ICEc/c promoters and, once being identified as "bistable", their temporal activation compared to that of intB13. In order to do this, we used promoter reporter fusions and tested blstability expression in P. knackmussii using epifluorescence microscopy. For the 6 promoters that showed bistable expression, we used time-lapse microscopy to study the timing of promoter expression in comparison to that of P,,,t or PlnR. We could establish two "early" and 3 "late" phase promoters in the process of transfer competence. In Chapter Three, I focused my attention on analysis of gene expression of one of the most highly transcribed gene clusters in the conserved core region of ICEc/c, a 14-kb gene cluster formed by the genes orf81655-orf68241. First we showed that this operon is part of the same bistability 'regulon' as intB13 and inrR, and analysed the genetic features that lead to high transcription. We studied the potential biological function of this cluster for ICE c/c by making specific gene deletions, showing that some interrupt ICEc/c transfer. We further analysed the orfdl655 promoter by constructing transcriptional egfp fusion reporter strains using the miniTn5 delivery system. Expression of the orf81655 promoter in single cells was compared to signals measured by Fluorescence In Situ Hybridization (FISH) on orfSl655 mRNA. We localized the promoter region of the gene cluster using FISH, promoter deletions, and by direct transcript analysis. We further used site-directed mutagenesis to understand the bistability character of the promoter region and the extremely high transcription but low translation from this mRNA. In Chapter Four, we set out to understand how bistability is generated in the tc pathway of ICEc/c. For this we tried rebuilding bistable expression from ICEc/c individual gene components in a host, which normally does not display bistability. As host we used P. putida without ICEc/c but with a single copy Pint-, PlnR- or PalpA- egfp fusion that enabled us to verify bistability formation. Subsequently, we built different assemblages of ICEc/c gene components, focusing on the parA-inrR region. Indeed, we found that bistable expression can be build from those components in P. putida without ICEc/c. Interestingly, most genetic constructs activated PaipA or PlnR, but only one resulted in bistable activation of PinT. This suggests that the parA-inrR region acts as a bistability "generator", but also as a bistability "relay" from early to late promoters in the tc pathway hierarchy. In the final fifth chapter, we conclude with a discussion of the relevance of the present thesis and the resulting perspectives for future studies.

Genomic organization, fine-mapping, and expression of the human islet-brain 1 (IB1)/c-Jun-amino-terminal kinase interacting protein-1 (JIP-1) gene.

Islet-brain 1 (IB1), a regulator of the pancreatic beta-cell function in the rat, is homologous to JIP-1, a murine inhibitor of c-Jun amino-terminal kinase (JNK). Whether IB1 and JIP-1 are present in humans was not known. We report the sequence of the 2133-bp human IB1 cDNA, the expression, structure, and fine-mapping of the human IB1 gene, and the characterization of an IB1 pseudogene. Human IB1 is 94% identical to rat IB1. The tissue-specific expression of IB1 in human is similar to that observed in rodent. The IB1 gene contains 12 exons and maps to chromosome 11 (11p11.2-p12), a region that is deleted in DEFECT-11 syndrome. Apart from an IB1 pseudogene on chromosome 17 (17q21), no additional IB1-related gene was found in the human genome. Our data indicate that the sequence and expression pattern of IB1 are highly conserved between rodent and human and provide the necessary tools to investigate whether IB1 is involved in human diseases.

Complex organization of CTF/NF-I, C/EBP, and HNF3 binding sites within the promoter of the liver-specific vitellogenin gene.

Vitellogenin genes are expressed specifically in the liver of female oviparous vertebrates under the strict control of estrogen. To explain this tissue-specific expression, we performed a detailed analysis of the Xenopus laevis vitellogenin gene B1 promoter by DNase I footprinting and gel mobility-shift assays. We characterized five binding sites for the ubiquitous factor CTF/NF-I. Two of these sites are close to the TATA-box, whereas the others are located on both sides of the estrogen responsive unit formed by two imperfect estrogen response elements. Moreover two liver-enriched factors, C/EBP and HNF3, were found to interact with multiple closely spaced proximal promoter elements in the first 100 base pairs upstream of the TATA-box. To confirm the physiological significance of this in vitro analysis, in vivo DNase I footprinting experiments were carried out using the ligation-mediated polymerase chain reaction technique. The various cis-elements characterized in vitro as binding sites for known transcription factors and more particularly for liver-enriched transcription factors are efficiently recognized in vivo as well, suggesting that they play an important role in the control of the liver-specific vitellogenin gene B1 expression.

Genomic organization and gene expression in a chromosomal region of Leishmania major.

Little is known about the relation between the genome organization and gene expression in Leishmania. Bioinformatic analysis can be used to predict genes and find homologies with known proteins. A model was proposed, in which genes are organized into large clusters and transcribed from only one strand, in the form of large polycistronic primary transcripts. To verify the validity of this model, we studied gene expression at the transcriptional, post-transcriptional and translational levels in a unique locus of 34kb located on chr27 and represented by cosmid L979. Sequence analysis revealed 115 ORFs on either DNA strand. Using computer programs developed for Leishmania genes, only nine of these ORFs, localized on the same strand, were predicted to code for proteins, some of which show homologies with known proteins. Additionally, one pseudogene, was identified. We verified the biological relevance of these predictions. mRNAs from nine predicted genes and proteins from seven were detected. Nuclear run-on analyses confirmed that the top strand is transcribed by RNA polymerase II and suggested that there is no polymerase entry site. Low levels of transcription were detected in regions of the bottom strand and stable transcripts were identified for four ORFs on this strand not predicted to be protein-coding. In conclusion, the transcriptional organization of the Leishmania genome is complex, raising the possibility that computer predictions may not be comprehensive.

Genome organization and gene expression shape the transposable element distribution in the Drosophila melanogaster euchromatin.

The distribution of transposable elements (TEs) in a genome reflects a balance between insertion rate and selection against new insertions. Understanding the distribution of TEs therefore provides insights into the forces shaping the organization of genomes. Past research has shown that TEs tend to accumulate in genomic regions with low gene density and low recombination rate. However, little is known about the factors modulating insertion rates across the genome and their evolutionary significance. One candidate factor is gene expression, which has been suggested to increase local insertion rate by rendering DNA more accessible. We test this hypothesis by comparing the TE density around germline- and soma-expressed genes in the euchromatin of Drosophila melanogaster. Because only insertions that occur in the germline are transmitted to the next generation, we predicted a higher density of TEs around germline-expressed genes than soma-expressed genes. We show that the rate of TE insertions is greater near germline- than soma-expressed genes. However, this effect is partly offset by stronger selection for genome compactness (against excess noncoding DNA) on germline-expressed genes. We also demonstrate that the local genome organization in clusters of coexpressed genes plays a fundamental role in the genomic distribution of TEs. Our analysis shows that-in addition to recombination rate-the distribution of TEs is shaped by the interaction of gene expression and genome organization. The important role of selection for compactness sheds a new light on the role of TEs in genome evolution. Instead of making genomes grow passively, TEs are controlled by the forces shaping genome compactness, most likely linked to the efficiency of gene expression or its complexity and possibly their interaction with mechanisms of TE silencing.

Current status of a model system: the gene Gp-9 and its association with social organization in fire ants.

The Gp-9 gene in fire ants represents an important model system for studying the evolution of social organization in insects as well as a rich source of information relevant to other major evolutionary topics. An important feature of this system is that polymorphism in social organization is completely associated with allelic variation at Gp-9, such that single-queen colonies (monogyne form) include only inhabitants bearing B-like alleles while multiple-queen colonies (polygyne form) additionally include inhabitants bearing b-like alleles. A recent study of this system by Leal and Ishida (2008) made two major claims, the validity and significance of which we examine here. After reviewing existing literature, analyzing the methods and results of Leal and Ishida (2008), and generating new data from one of their study sites, we conclude that their claim that polygyny can occur in Solenopsis invicta in the U.S.A. in the absence of expression of the b-like allele Gp-9(b) is unfounded. Moreover, we argue that available information on insect OBPs (the family of proteins to which GP-9 belongs), on the evolutionary/population genetics of Gp-9, and on pheromonal/behavioral control of fire ant colony queen number fails to support their view that GP-9 plays no role in the chemosensory-mediated communication that underpins regulation of social organization. Our analyses lead us to conclude that there are no new reasons to question the existing consensus view of the Gp-9 system outlined in Gotzek and Ross (2007).

The knock-out of ARP3a gene affects F-actin cytoskeleton organization altering cellular tip growth, morphology and development in moss Physcomitrella patens.

The seven subunit Arp2/3 complex is a highly conserved nucleation factor of actin microfilaments. We have isolated the genomic sequence encoding a putative Arp3a protein of the moss Physcomitrella patens. The disruption of this ARP3A gene by allele replacement has generated loss-of-function mutants displaying a complex developmental phenotype. The loss-of function of ARP3A gene results in shortened, almost cubic chloronemal cells displaying affected tip growth and lacking differentiation to caulonemal cells. In moss arp3a mutants, buds differentiate directly from chloronemata to form stunted leafy shoots having differentiated leaves similar to wild type. Yet, rhizoids never differentiate from stem epidermal cells. To characterize the F-actin organization in the arp3a-mutated cells, we disrupted ARP3A gene in the previously described HGT1 strain expressing conditionally the GFP-talin marker. In vivo observation of the F-actin cytoskeleton during P. patens development demonstrated that loss-of-function of Arp3a is associated with the disappearance of specific F-actin cortical structures associated with the establishment of localized cellular growth domains. Finally, we show that constitutive expression of the P. patens Arp3a and its Arabidopsis thaliana orthologs efficiently complement the mutated phenotype indicating a high degree of evolutionary conservation of the Arp3 function in land plants.

Comparison of the organization and fine structure of a chicken and a Xenopus laevis vitellogenin gene.

The structural organization and the coding nucleotide sequence of the Xenopus laevis A2 and the chicken major vitellogenin genes have been compared. Both genes show the same exon-intron organization. However, the degree of homology between the nucleotide and derived amino acid sequences varies extensively along the genes. Several of the 35 exons are quite similar, and a unique cysteine motif in the lipovitellin II domain is conserved between the two genes. In contrast, one internal region is quite divergent. Part of this region encodes phosvitin, which appears to have evolved rapidly by both point mutations and duplications of serines or short other amino acid stretches. On the basis of these observations, we discuss the possible mechanism of evolution of phosvitin in vertebrates.

A simple genetic basis for complex social behaviour mediates widespread gene expression differences.

A remarkable social polymorphism is controlled by a single Mendelian factor in the fire ant Solenopsis invicta. A genomic element marked by the gene Gp-9 determines whether workers tolerate one or many fertile queens in their colony. Gp-9 was recently shown to be part of a supergene with two nonrecombining variants, SB and Sb. SB/SB and SB/Sb queens differ in how they initiate new colonies, and in many physiological traits, for example odour and maturation rate. To understand how a single genetic element can affect all these traits, we used a microarray to compare gene expression patterns between SB/SB and SB/Sb queens of three different age classes: 1-day-old unmated queens, 11-day-old unmated queens and mated, fully reproductive queens collected from mature field colonies. The number of genes that were differentially expressed between SB/SB and SB/Sb queens of the same age class was smallest in 1-day-old queens, maximal in 11-day-old queens and intermediate in reproductive queens. Gene ontology analysis showed that SB/SB queens upregulate reproductive genes faster than SB/Sb queens. For all age classes, genes inside the supergene were overrepresented among the differentially expressed genes. Consistent with the hypothesized greater number of transposons in the Sb supergene, 13 transposon genes were upregulated in SB/Sb queens. Viral genes were also upregulated in SB/Sb mature queens, consistent with the known greater parasite load in colonies headed by SB/Sb queens compared with colonies headed by SB/SB queens. Eighteen differentially expressed genes between reproductive queens were involved in chemical signalling. Our results suggest that many genes in the supergene are involved in regulating social organization and queen phenotypes in fire ants.

Linking melanism to brain development: expression of a melanism-related gene in barn owl feather follicles covaries with sleep ontogeny.

BACKGROUND: Intra-specific variation in melanocyte pigmentation, common in the animal kingdom, has caught the eye of naturalists and biologists for centuries. In vertebrates, dark, eumelanin pigmentation is often genetically determined and associated with various behavioral and physiological traits, suggesting that the genes involved in melanism have far reaching pleiotropic effects. The mechanisms linking these traits remain poorly understood, and the potential involvement of developmental processes occurring in the brain early in life has not been investigated. We examined the ontogeny of rapid eye movement (REM) sleep, a state involved in brain development, in a wild population of barn owls (Tyto alba) exhibiting inter-individual variation in melanism and covarying traits. In addition to sleep, we measured melanistic feather spots and the expression of a gene in the feather follicles implicated in melanism (PCSK2). RESULTS: As in mammals, REM sleep declined with age across a period of brain development in owlets. In addition, inter-individual variation in REM sleep around this developmental trajectory was predicted by variation in PCSK2 expression in the feather follicles, with individuals expressing higher levels exhibiting a more precocial pattern characterized by less REM sleep. Finally, PCSK2 expression was positively correlated with feather spotting. CONCLUSIONS: We demonstrate that the pace of brain development, as reflected in age-related changes in REM sleep, covaries with the peripheral activation of the melanocortin system. Given its role in brain development, variation in nestling REM sleep may lead to variation in adult brain organization, and thereby contribute to the behavioral and physiological differences observed between adults expressing different degrees of melanism.

Molecular variation at a candidate gene implicated in the regulation of fire ant social behavior

The fire ant Solenopsis invicta and its close relatives display an important social polymorphism involving differences in colony queen number. Colonies are headed by either a single reproductive queen (monogyne form) or multiple queens (polygyne form). This variation in social organization is associated with variation at the gene Gp-9, with monogyne colonies harboring only B-like allelic variants and polygyne colonies always containing b-like variants as well. We describe naturally occurring variation at Gp-9 in fire ants based on 185 full-length sequences, 136 of which were obtained from S. invicta collected over much of its native range. While there is little overall differentiation between most of the numerous alleles observed, a surprising amount is found in the coding regions of the gene, with such substitutions usually causing amino acid replacements. This elevated coding-region variation may result from a lack of negative selection acting to constrain amino acid replacements over much of the protein, different mutation rates or biases in coding and non-coding sequences, negative selection acting with greater strength on non-coding than coding regions, and/or positive selection acting on the protein. Formal selection analyses provide evidence that the latter force played an important role in the basal b-like lineages coincident with the emergence of polygyny. While our data set reveals considerable paraphyly and polyphyly of S. invicta sequences with respect to those of other fire ant species, the b-like alleles of the socially polymorphic species are monophyletic. An expanded analysis of colonies containing alleles of this clade confirmed the invariant link between their presence and expression of polygyny. Finally, our discovery of several unique alleles bearing various combinations of b-like and B-like codons allows us to conclude that no single b-like residue is completely predictive of polygyne behavior and, thus, potentially causally involved in its expression. Rather, all three typical b-like residues appear to be necessary.