Gene copy number polymorphisms in an arbuscular mycorrhizal fungal population.

Gene copy number polymorphism was studied in a population of the arbuscular mycorrhizal fungus Glomus intraradices by using a quantitative PCR approach on four different genomic regions. Variation in gene copy number was found for a pseudogene and for three ribosomal genes, providing conclusive evidence for a widespread occurrence of macromutational events in the population.

Fine mapping and functional analysis of the multiple sclerosis risk gene CD6.

CD6 has recently been identified and validated as risk gene for multiple sclerosis (MS), based on the association of a single nucleotide polymorphism (SNP), rs17824933, located in intron 1. CD6 is a cell surface scavenger receptor involved in T-cell activation and proliferation, as well as in thymocyte differentiation. In this study, we performed a haptag SNP screen of the CD6 gene locus using a total of thirteen tagging SNPs, of which three were non-synonymous SNPs, and replicated the recently reported GWAS SNP rs650258 in a Spanish-Basque collection of 814 controls and 823 cases. Validation of the six most strongly associated SNPs was performed in an independent collection of 2265 MS patients and 2600 healthy controls. We identified association of haplotypes composed of two non-synonymous SNPs [rs11230563 (R225W) and rs2074225 (A257V)] in the 2(nd) SRCR domain with susceptibility to MS (P max(T) permutation = 1×10(-4)). The effect of these haplotypes on CD6 surface expression and cytokine secretion was also tested. The analysis showed significantly different CD6 expression patterns in the distinct cell subsets, i.e. - CD4(+) naïve cells, P = 0.0001; CD8(+) naïve cells, P<0.0001; CD4(+) and CD8(+) central memory cells, P = 0.01 and 0.05, respectively; and natural killer T (NKT) cells, P = 0.02; with the protective haplotype (RA) showing higher expression of CD6. However, no significant changes were observed in natural killer (NK) cells, effector memory and terminally differentiated effector memory T cells. Our findings reveal that this new MS-associated CD6 risk haplotype significantly modifies expression of CD6 on CD4(+) and CD8(+) T cells.

Molecular identification of Giardia duodenalis in Ecuador by polymerase chain reaction-restriction fragment length polymorphism

The aim of this study was to determine the genetic diversity of Giardia duodenalis present in a human population living in a northern Ecuadorian rain forest. All Giardia positive samples (based on an ELISA assay) were analysed using a semi-nested polymerase chain reaction-restriction fragment length polymorphism assay that targets the glutamate dehydrogenase (gdh) gene; those amplified were subsequently genotyped using NlaIV and RsaI enzymes. The gdh gene was successfully amplified in 74 of 154 ELISA positive samples; 69 of the 74 samples were subsequently genotyped. Of these 69 samples, 42 (61%) were classified as assemblage B (26 as BIII and 16 as BIV), 22 (32%) as assemblage A (3 as AI and 19 as AII) and five (7%) as mixed AII and BIII types. In this study site we observe similar diversity in genotypes to other regions in Latin America, though in contrast to some previous studies, we found similar levels of diarrheal symptoms in those individuals infected with assemblage B compared with those infected with assemblage A.

Evidence of new risk genetic factor to systemic lupus erythematosus: the UBASH3A gene

The ubiquitin associated and Src-homology 3 (SH3) domain containing A (UBASH3a) is a suppressor of T-cell receptor signaling, underscoring antigen presentation to T-cells as a critical shared mechanism of diseases pathogenesis. The aim of the present study was to determine whether the UBASH3a gene influence the susceptibility to systemic lupus erythematosus (SLE) in Caucasian populations. We evaluated five UBASH3a polymorphisms (rs2277798, rs2277800, rs9976767, rs13048049 and rs17114930), using TaqMan® allelic discrimination assays, in a discovery cohort that included 906 SLE patients and 1165 healthy controls from Spain. The SNPs that exhibit statistical significance difference were evaluated in a German replication cohort of 360 SLE patients and 379 healthy controls. The case-control analysis in the Spanish population showed a significant association between the rs9976767 and SLE (Pc = 9.9E-03 OR = 1.21 95%CI = 1.07-1.37) and a trend of association for the rs2277798 analysis (P = 0.09 OR = 0.9 95%CI = 0.79-1.02). The replication in a German cohort and the meta-analysis confirmed that the rs9976767 (Pc = 0.02; Pc = 2.4E-04, for German cohort and meta-analysis, respectively) and rs2277798 (Pc = 0.013; Pc = 4.7E-03, for German cohort and meta-analysis, respectively) UBASH3a variants are susceptibility factors for SLE. Finally, a conditional regression analysis suggested that the most likely genetic variation responsible for the association was the rs9976767 polymorphism. Our results suggest that UBASH3a gene plays a role in the susceptibility to SLE. Moreover, our study indicates that UBASH3a can be considered as a common genetic factor in autoimmune diseases.

Genetic characterisation of Langerin gene in human immunodeficiency virus-1-infected women from Bahia, Brazil

Studies on human genetic variations are a useful source of knowledge about human immunodeficiency virus (HIV)-1 infection. The Langerin protein, found at the surface of Langerhans cells, has an important protective role in HIV-1 infection. Differences in Langerin function due to host genetic factors could influence susceptibility to HIV-1 infection. To verify the frequency of mutations in the Langerin gene, 118 samples from HIV-1-infected women and 99 samples from HIV-1-uninfected individuals were selected for sequencing of the promoter and carbohydrate recognition domain (CRD)-encoding regions of the Langerin gene. Langerin promoter analysis revealed two single nucleotide polymorphisms (SNPs) and one mutation in both studied groups, which created new binding sites for certain transcription factors, such as NFAT5, HOXB9.01 and STAT6.01, according to MatInspector software analysis. Three SNPs were observed in the CRD-encoding region in HIV-1-infected and uninfected individuals: p.K313I, c.941C>T and c.983C>T. This study shows that mutations in the Langerin gene are present in the analysed populations at different genotypic and allelic frequencies. Further studies should be conducted to verify the role of these mutations in HIV-1 susceptibility.

Identification of a novel splice-site mutation in the CYP1A2 gene.

AIMS: To identify the molecular basis for a low CYP1A2 metabolic status, as determined by a caffeine phenotyping test, in a 71-year-old, nonsmoking, Caucasian woman who presented with very high clozapine concentrations despite being administered a standard dose of the drug. METHODS: The nucleotide sequence of the 7 exons, exon-intron boundaries and 5'-flanking region of the CYP1A2 gene was analysed by direct sequencing. RESULTS: Only one heterozygous point mutation was identified in the donor splice site of intron 6 (3534G > A) of CYP1A2. This mutation could cause abnormal RNA splicing and therefore lead to a truncated nonfunctional enzyme. No other carrier of this mutation was identified in a population of 100 unrelated healthy Caucasians. CONCLUSIONS: This is the first report of a splice-site mutation affecting the CYP1A2 gene. This polymorphism is a likely explanation for the low CYP1A2 activity associated with high clozapine concentrations in this patient.

Polymorphism at a sex-linked transcription cofactor in European tree frogs (Hyla arborea): Sex-antagonistic selection or neutral processes?

Nascent sex chromosomes offer a unique opportunity to investigate the evolutionary fate of genesrecently trapped in non-recombining segments. A housekeeping gene (MED15) was recently shown to lie on the nascent sex-chromosomes of the European tree frog (Hyla arborea), with different alleles fixed on the X and the Y chromosomes. Here we document a polymorphism (glutamine deletion) in the X copy of the gene, and use population surveys and experimental crosses to test whether this polymorphism is neutral or maintained by sex-antagonistic selection. Tadpoles from parents of known genotypes revealed significant discrepancies from Mendelian inheritance, suggesting possible sex-antagonistic effects under laboratory conditions. Quantitatively, however, these effects did not meet the conditions for polymorphism maintenance. Furthermore, field estimates of female genotypic frequencies did not differ from Hardy-Weinberg equilibrium and allelic frequencies on the X chromosome did not differ between sexes. In conclusion, although sex antagonistic effects cannot be excluded given the laboratory conditions, the X-linked polymorphism under study appears neutral in the wild. Alternatively, sex-antagonistic selection might still account for the fixation of a male specific allele on the Y chromosome.

Gene-lifestyle interaction and type 2 diabetes: the EPIC interact case-cohort study.

BACKGROUND Understanding of the genetic basis of type 2 diabetes (T2D) has progressed rapidly, but the interactions between common genetic variants and lifestyle risk factors have not been systematically investigated in studies with adequate statistical power. Therefore, we aimed to quantify the combined effects of genetic and lifestyle factors on risk of T2D in order to inform strategies for prevention. METHODS AND FINDINGS The InterAct study includes 12,403 incident T2D cases and a representative sub-cohort of 16,154 individuals from a cohort of 340,234 European participants with 3.99 million person-years of follow-up. We studied the combined effects of an additive genetic T2D risk score and modifiable and non-modifiable risk factors using Prentice-weighted Cox regression and random effects meta-analysis methods. The effect of the genetic score was significantly greater in younger individuals (p for interaction  = 1.20×10-4). Relative genetic risk (per standard deviation [4.4 risk alleles]) was also larger in participants who were leaner, both in terms of body mass index (p for interaction  = 1.50×10-3) and waist circumference (p for interaction  = 7.49×10-9). Examination of absolute risks by strata showed the importance of obesity for T2D risk. The 10-y cumulative incidence of T2D rose from 0.25% to 0.89% across extreme quartiles of the genetic score in normal weight individuals, compared to 4.22% to 7.99% in obese individuals. We detected no significant interactions between the genetic score and sex, diabetes family history, physical activity, or dietary habits assessed by a Mediterranean diet score. CONCLUSIONS The relative effect of a T2D genetic risk score is greater in younger and leaner participants. However, this sub-group is at low absolute risk and would not be a logical target for preventive interventions. The high absolute risk associated with obesity at any level of genetic risk highlights the importance of universal rather than targeted approaches to lifestyle intervention.

Functional and molecular evolutionary analysis of the "Brevix Radix" gene family in mono-and dicotyledonous plants

ABSTRACTIn contrast to animals, plants cannot move from their place of birth and, therefore, need to adapt to their particular habitat in order to survive. Thus, plant development is remarkably plastic, making plants an ideal system for the isolation of genes that account for intraspecific natural variation and possibly environmental adaptation. However, to date, this approach mostly identified null alleles and missed mutations with subtle effects. For instance, BREVIS RADIX (BRX) has been isolated as a key regulator of root growth through a naturally occurring loss-of-function allele in the Arabidopsis thaliana accession Uk-1 and is the founding member of a highly-conserved plant-specific gene family.In this work, we show that a strong selective pressure is acting on the BRX gene family and dates back before the monocot-dicot divergence. However, functional diversification is observed mainly in dicotyledon BRX family genes and is correlated with acceleration in the evolutionary rates in the N-terminal regions. Population genetic data revealed that BRX is highly conserved across Arabidopsis accessions and presents signatures of adaptation. Interestingly, a seven amino acid deletion polymorphism in BRX sequence was found in a few accessions, which seems to be responsible for their enhanced primary root growth. Nevertheless, BRX might not only be active in the root, as suggested by its expression in the shoot. Indeed, leaves and cotyledons of brx mutants are significantly smaller than wild- type. This phenotype is a direct consequence of the absence of BRX function in the shoot rather than an indirect effect of an altered root system growth. Interestingly, cotyledons of brx plants reflect the same physiological defects as the root. Moreover, phenotypes in BRX gain-of-function plants, such as epinastic leaves and increased epidermal cell size, could be associated with an increase in leaf brassinosteroid content.Collectively, these results indicate that BRX contributes to local adaptation by ubiquitously regulating plant growth, probably through the modulation of brassinosteroid biosynthesis.RÉSUMÉContrairement à la plupart des animaux, les plantes ne peuvent se mouvoir et doivent ainsi s'adapter à leur environnement pour survivre. Pour cette raison, elles représentent un système idéal pour l'identification de gènes contribuant à la variation naturelle intra- spécifique, ainsi qu'à l'adaptation. Cependant, cette approche a, jusqu'à présent, surtout permis d'isoler des allèles nuls et non des mutations conférant des effets plus subtiles. C'est le cas du gène Β REVIS RADIX (BRX), un régulateur clé de la croissance racinaire, qui a été identifié grâce à un allèle non-fonctionnel présent dans l'accession naturelle d'Arabidopsis thaliana Uk-1. BRX et ses homologues des plantes mono- et dicotylédones forment une famille très conservée et spécifique aux plantes.Dans ce travail, nous démontrons que la famille de gènes BRX est soumise à une forte pression de sélection qui remonte avant la divergence entre mono- et dicotylédones. Cependant, une diversification fonctionnelle a été observée chez les gènes des dicotylédones et corrèle avec une accélération de la vitesse d'évolution dans leur région N- terminale. Une analyse génétique de différentes accessions naturelles d'Arabidopsis a révélé que BRX est hautement conservé et présente des signatures d'adaptation. Remarquablement, un polymorphisme de délétion de sept acides aminés a été détecté dans quelques accessions et a pour conséquence une plus forte croissance de la racine primaire. Néanmoins, il semble que le rôle de BRX ne se limite pas qu'à la racine, comme indiqué par son expression dans les parties aériennes de la plante. En effet, les mutants brx présentent des cotylédons et des feuilles significativement plus petits que le type sauvage, une conséquence directe de l'absence d'activité de BRX dans ces organes. Nous avons aussi noté que les cotylédons des mutants brx, à l'instar des racines, ont une perception altérée de l'auxine et peuvent être complémentés par l'application exogène de brassinostéroïdes. De plus, dans des plantes présentant un gain de fonction BRX, les feuilles sont épinastiques et les cellules de leur épiderme plus grandes. Ces phénotypes sont accompagnés d'une augmentation de la concentration de brassinostéroïdes dans les feuilles. Conjointement, ces résultats démontrent que BRX contribue à une adaptation locale de la plante par la régulation générale de sa croissance, probablement en modulant la biosynthèse des brassinostéroïdes.

Blood pressure loci identified with a gene-centric array.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

Intronic variants in the NFKB1 gene may influence hearing forecast in patients with unilateral sensorineural hearing loss in Meniere's disease.

Meniere's disease is an episodic vestibular syndrome associated with sensorineural hearing loss (SNHL) and tinnitus. Patients with MD have an elevated prevalence of several autoimmune diseases (rheumatoid arthritis, systemic lupus erythematosus, ankylosing spondylitis and psoriasis), which suggests a shared autoimmune background. Functional variants of several genes involved in the NF-κB pathway, such as REL, TNFAIP3, NFKB1 and TNIP1, have been associated with two or more immune-mediated diseases and allelic variations in the TLR10 gene may influence bilateral affectation and clinical course in MD. We have genotyped 716 cases of MD and 1628 controls by using the ImmunoChip, a high-density genotyping array containing 186 autoimmune loci, to explore the association of immune system related-loci with sporadic MD. Although no single nucleotide polymorphism (SNP) reached a genome-wide significant association (p<10(-8)), we selected allelic variants in the NF-kB pathway for further analyses to evaluate the impact of these SNPs in the clinical outcome of MD in our cohort. None of the selected SNPs increased susceptibility for MD in patients with uni or bilateral SNHL. However, two potential regulatory variants in the NFKB1 gene (rs3774937 and rs4648011) were associated with a faster hearing loss progression in patients with unilateral SNHL. So, individuals with unilateral MD carrying the C allele in rs3774937 or G allele in rs4648011 had a shorter mean time to reach hearing stage 3 (>40 dB HL) (log-rank test, corrected p values were p = 0.009 for rs3774937 and p = 0.003 for rs4648011, respectively). No variants influenced hearing in bilateral MD. Our data support that the allelic variants rs3774937 and rs4648011 can modify hearing outcome in patients with MD and unilateral SNHL.

A molecular platform for the diagnosis of multidrug-resistant and pre-extensively drug-resistant tuberculosis based on single nucleotide polymorphism mutations present in Colombian isolates of Mycobacterium tuberculosis

Developing a fast, inexpensive, and specific test that reflects the mutations present in Mycobacterium tuberculosis isolates according to geographic region is the main challenge for drug-resistant tuberculosis (TB) control. The objective of this study was to develop a molecular platform to make a rapid diagnosis of multidrug-resistant (MDR) and extensively drug-resistant TB based on single nucleotide polymorphism (SNP) mutations present in therpoB, katG, inhA,ahpC, and gyrA genes from Colombian M. tuberculosis isolates. The amplification and sequencing of each target gene was performed. Capture oligonucleotides, which were tested before being used with isolates to assess the performance, were designed for wild type and mutated codons, and the platform was standardised based on the reverse hybridisation principle. This method was tested on DNA samples extracted from clinical isolates from 160 Colombian patients who were previously phenotypically and genotypically characterised as having susceptible or MDR M. tuberculosis. For our method, the kappa index of the sequencing results was 0,966, 0,825, 0,766, 0,740, and 0,625 forrpoB, katG, inhA,ahpC, and gyrA, respectively. Sensitivity and specificity were ranked between 90-100% compared with those of phenotypic drug susceptibility testing. Our assay helps to pave the way for implementation locally and for specifically adapted methods that can simultaneously detect drug resistance mutations to first and second-line drugs within a few hours.

Association of the solute carrier family 11 member 1 gene polymorphisms with susceptibility to leprosy in a Brazilian sample

Natural resistance-associated macrophage protein 1/solute carrier family 11 member 1 gene (Nramp1/Slc11a1) is a gene that controls the susceptibility of inbred mice to intracellular pathogens. Polymorphisms in the human Slc11a1/Nramp1 gene have been associated with host susceptibility to leprosy. This study has evaluated nine polymorphisms of the Slc11a1/Nramp1 gene [(GT)n, 274C/T, 469+14G/C, 577-18G/A, 823C/T, 1029 C/T, 1465-85G/A, 1703G/A, and 1729+55del4] in 86 leprosy patients (67 and 19 patients had the multibacillary and the paucibacillary clinical forms of the disease, respectively), and 239 healthy controls matched by age, gender, and ethnicity. The frequency of allele 2 of the (GT)n polymorphism was higher in leprosy patients [p = 0.04, odds ratio (OR) = 1.49], whereas the frequency of allele 3 was higher in the control group (p = 0.03; OR = 0.66). Patients carrying the 274T allele (p = 0.04; OR = 1.49) and TT homozygosis (p = 0.02; OR = 2.46), such as the 469+14C allele (p = 0.03; OR = 1.53) of the 274C/T and 469+14G/C polymorphisms, respectively, were more frequent in the leprosy group. The leprosy and control groups had similar frequency of the 577-18G/A, 823C/T, 1029C/T, 1465-85G/A, 1703G/A, and 1729+55del4 polymorphisms. The 274C/T polymorphism in exon 3 and the 469+14G/C polymorphism in intron 4 were associated with susceptibility to leprosy, while the allele 2 and 3 of the (GT)n polymorphism in the promoter region were associated with susceptibility and protection to leprosy, respectively.

Dominant PRPF31 Mutations Are Hypostatic to a Recessive CNOT3 Polymorphism in Retinitis Pigmentosa: A Novel Phenomenon of "Linked Trans-Acting Epistasis".

Mutations in PRPF31 are responsible for autosomal dominant retinitis pigmentosa (adRP, RP11 form) and affected families show nonpenetrance. Differential expression of the wildtype PRPF31 allele is responsible for this phenomenon: coinheritance of a mutation and a higher expressing wildtype allele provide protection against development of disease. It has been suggested that a major modulating factor lies in close proximity to the wildtype PRPF31 gene on Chromosome 19, implying that a cis-acting factor directly alters PRPF31 expression. Variable expression of CNOT3 is one determinant of PRPF31 expression. This study explored the relationship between CNOT3 (a trans-acting factor) and its paradoxical cis-acting nature in relation to RP11. Linkage analysis on Chromosome 19 was performed in mutation-carrying families, and the inheritance of the wildtype PRPF31 allele in symptomatic-asymptomatic sibships was assessed-confirming that differential inheritance of wildtype chromosome 19q13 determines the clinical phenotype (P < 2.6 × 10(-7) ). A theoretical model was constructed that explains the apparent conflict between the linkage data and the recent demonstration that a trans-acting factor (CNOT3) is a major nonpenetrance factor: we propose that this apparently cis-acting effect arises due to the intimate linkage of CNOT3 and PRPF31 on Chromosome 19q13-a novel mechanism that we have termed "linked trans-acting epistasis."

Single loci with major effects in fire ants : the effects of the complementary sex-determining locus and the green beard supergene on gene expression

Etant données la complexité et la redondance des réseaux de gènes influençant de nombreux phénotypes, l'étude des rares cas d'un locus unique ayant des effets importants sur de nombreux phénotypes peut fournir des informations cruciales sur l'évolution des traits complexes. Nous avons séquencé le génome de la fourmi de feu Solenopsis invicta pour étudier comment l'expression des gènes détermine les effets majeurs et étendus de deux loci uniques sur le phénotype. Le premier locus concerne la détermination du sexe par le modèle des allèles complémentaires. Ce locus est connu pour déterminer le sexe chez tous les hyménoptères mais n'a été caractérisé que chez les abeilles. Les hétérozygotes pour ce locus se développent en reines diploïdes (ou ouvrières stériles) alors que les homozygotes se développent en mâles diploïdes incapables de produire du sperme et les hémizygotes en mâles haploïdes fertiles. Nous avons comparé l'expression des gènes entre les reines et les deux types de mâles au stade pupe, ainsi que 1 et 11 jours après l'émergence. Nous avons trouvé un changement prononcé de l'expression des gènes chez les mâles diploïdes, passant de très proche de celle des reines au stade pupe à identique aux mâles haploïdes 11 jours après l'émergence. Cela signifie que les mâles diploïdes sont condamnés à être stériles parce que les effets après émergence du locus de détermination du sexe ne per¬mettent pas d'effacer les effets de la ploïdie sur l'expression des gènes pendant le stade pupe, quand la spermatogénèse prend place. Le second locus aux effets majeurs que nous avons étudié est le supergène dit "green beard", qui consiste en 616 gènes couvrant 55% d'un chromosome (13 Mb) et est caractérisé par une absence de recombinaison entre les deux variants du supergène : "Social B" et "Social b" (SB et Sb). Au travers de l'effet "green beard", par lequel les ouvrières avec le supergène Sb discriminent favorablement les reines qui partagent ce supergène de façon perceptible, le génotype des reines fondatrices au niveau de ce supergène détermine l'organisation de la colonie : soit elle contient une seule reine SB/SB, soit plusieurs reines SB/Sb. Nous avons montré que le chromosome Sb a évolué comme le chromosome Y, accumulant probablement des allèles favorables dans des colonies avec plusieurs reines mais défavorables dans des colonies avec une seule reine (cf. gènes sexuellement antagonistes), ainsi que des transposons et des séquences répéti¬tives. Nous avons également montré que le polymorphisme du supergène cause de grandes différences d'expression chez les ouvrières et particulièrement les reines mais pas chez les mâles. Pour comprendre comment le polymorphisme du supergène chez les reines peut affecter l'organisation de la colonie, nous avons comparé l'expression entre les génotypes SB/SB et SB/Sb chez des reines vierges (1 et 11 jours) et des reines matures. Nous avons montré que les reines SB/SB sur-régulent des gènes impliqués dans la reproduction, expli-quant pourquoi elle grandissent plus rapidement et peuvent fonder des colonies de façon indépendante, tandis que les reines SB/Sb (qui ne peuvent fonder une nouvelle colonie) sur-régulent des gènes de signalement chimique qui affectent l'organisation des colonies par l'effet "green beard". - Given the complexity and redundancy of the gene networks that underlie many pheno- types, the study of rare cases of a single locus having major effects on many phenotypes can give powerful insights into the evolution of complex traits. We sequenced the genome of Solenopsis invicta fire ants to study how gene expression mediates the widespread major effects of two single loci on phenotype. The first is the complementary sex-determining locus, which is known to exist in most Hymenoptera despite being characterized only for honeybees. Heterozygotes at this locus become diploid queens (or sterile workers), homozy¬gotes become aspermic diploid males, and hemizygotes become fertile haploid males. We compared gene expression between queens and both types of males in pupae and 1 and 11 days after eclosion. We found a pronounced shift in gene expression in diploid males, from being nearly identical to queens as pupae to identical to haploid males 11 days after eclosion. This means that diploid males are condemned to sterility because the overriding effects of the sex locus after eclosion cannot undo the ploidy effects on expression during the pupal stage, when spermatogenesis must be completed. The second locus with major ef¬fects that we studied was the so-called "green beard" supergene, which consists of 616 genes encompassing 55% of one chromosome (13 Mb), without recombination between the two variants "Social B" and "Social b" (SB and Sb) supergene. Through the green beard effect, i.e. workers with the Sb supergene discriminating in favor of queens who perceptibly share this supergene, the founding queen's genotype at the supergene determines colony organi¬zation: either headed by a single SB/SB queen or many SB/Sb queens. We show that the Sb chromosome evolved like a Y-chromosome, probably accumulating alleles beneficial in multi-queen colonies but disadvantageous in single-queen colonies (cf. sexually antagonistic genes), as well as transposons and repetitive sequences. We also show that the polymor¬phism of the supergene causes widespread expression differences in workers and especially queens but not in males. To understand how the polymorphism at the supergene in queen can transform colony organization, we compared the expression between SB/SB and SB/Sb virgin queens (1 and 11 days) and mother queens. We show that SB/SB queens up-regulate genes involved in reproduction, explaining why they mature faster and can found colonies independently, while SB/Sb queens (which cannot found colonies) up-regulate chemical signaling genes that can transform colonies through the green beard effect.