A higher mutational burden in females supports a "female protective model" in neurodevelopmental disorders.

Increased male prevalence has been repeatedly reported in several neurodevelopmental disorders (NDs), leading to the concept of a "female protective model." We investigated the molecular basis of this sex-based difference in liability and demonstrated an excess of deleterious autosomal copy-number variants (CNVs) in females compared to males (odds ratio [OR] = 1.46, p = 8 × 10(-10)) in a cohort of 15,585 probands ascertained for NDs. In an independent autism spectrum disorder (ASD) cohort of 762 families, we found a 3-fold increase in deleterious autosomal CNVs (p = 7 × 10(-4)) and an excess of private deleterious single-nucleotide variants (SNVs) in female compared to male probands (OR = 1.34, p = 0.03). We also showed that the deleteriousness of autosomal SNVs was significantly higher in female probands (p = 0.0006). A similar bias was observed in parents of probands ascertained for NDs. Deleterious CNVs (>400 kb) were maternally inherited more often (up to 64%, p = 10(-15)) than small CNVs < 400 kb (OR = 1.45, p = 0.0003). In the ASD cohort, increased maternal transmission was also observed for deleterious CNVs and SNVs. Although ASD females showed higher mutational burden and lower cognition, the excess mutational burden remained, even after adjustment for those cognitive differences. These results strongly suggest that females have an increased etiological burden unlinked to rare deleterious variants on the X chromosome. Carefully phenotyped and genotyped cohorts will be required for identifying the symptoms, which show gender-specific liability to mutational burden.

A corneal dystrophy associated with transforming growth factor beta-induced Gly623Asp mutation an amyloidogenic phenotype.

PURPOSE: To present the light and electron microscopic findings of a unique corneal dystrophy never before described in a German family carrying the Gly623Asp Mutation of the TGFBI gene with late clinical onset. DESIGN: Experimental study. PARTICIPANTS: Four affected and 6 nonaffected family members. METHODS: Slit-lamp examination, photographic documentation, and isolation of genomic DNA from peripheral blood leucocytes obtained from each family member examined. Exons 3, 4, 5, and 11 to 14 of the TGFBI gene were amplified and sequenced in these family members. Five corneal buttons of 3 affected siblings were excised at the time of penetrating keratoplasty. Light and electron microscopic examination were performed including immunohistochemistry with antibodies against keratoepithelin (KE) 2 and 15. MAIN OUTCOME MEASURES: Clinical and histologic characteristics of corneal opacification in affected patients and presence of coding region changes in the TGFBI gene. RESULTS: The specimens showed destructive changes in Bowman's layer and the adjacent stroma. Patchy Congo red-positive amyloid deposits were found within the epithelium in 1 cornea, in Bowman's layer and in the anterior stroma of all specimens also showing KE2, but not KE15, immunostaining. Electron microscopy revealed deposits mainly located in the anterior stroma and Bowman's layer and in small amounts in the basal area of some epithelial cells. The destroyed areas were strongly Alcian blue-positive, the Masson Trichrome stain proved mainly negative for the deposits. All affected but none of the unaffected family members had a heterozygous missense mutation in exon 14 of the TGFBI gene (G-->A transition at nucleotide 1915) replacing glycin by aspartic acid amino acid (Gly623Asp) at position 623 of the KE protein. CONCLUSIONS: In contrast with the patient carrying the Gly623Asp mutation of the TGFBI gene described by Afshari et al, our cases presented with Salzmann's nodular degeneration-like clinical features and their specimens contained KE2-positive amyloid. The reason for this now "meeting the expectation histologic phenotype" is unclear. The histologic findings emphasize that this is a unique corneal dystrophy, which shares no clinical characteristics with Reis-Bücklers' dystrophy and should be treated as a distinct entity. FINANCIAL DISCLOSURE(S): The authors have no proprietary or commercial interest in any materials discussed in this article.

Stable transmission of targeted gene modification using single-stranded oligonucleotides with flanking LNAs.

Targeted mutagenesis directed by oligonucleotides (ONs) is a promising method for manipulating the genome in higher eukaryotes. In this study, we have compared gene editing by different ONs on two new target sequences, the eBFP and the rd1 mutant photoreceptor betaPDE cDNAs, which were integrated as single copy transgenes at the same genomic site in 293T cells. Interestingly, antisense ONs were superior to sense ONs for one target only, showing that target sequence can by itself impart strand-bias in gene editing. The most efficient ONs were short 25 nt ONs with flanking locked nucleic acids (LNAs), a chemistry that had only been tested for targeted nucleotide mutagenesis in yeast, and 25 nt ONs with phosphorothioate linkages. We showed that LNA-modified ONs mediate dose-dependent target modification and analyzed the importance of LNA position and content. Importantly, when using ONs with flanking LNAs, targeted gene modification was stably transmitted during cell division, which allowed reliable cloning of modified cells, a feature essential for further applications in functional genomics and gene therapy. Finally, we showed that ONs with flanking LNAs aimed at correcting the rd1 stop mutation could promote survival of photoreceptors in retinas of rd1 mutant mice, suggesting that they are also active in vivo.

Le mélanome et les polymorphismes associés au dysfonctionnement du système immunitaire

La recherche biomédicale profite de plus en plus au développement des techniques de séquençage et d'analyse de l'ADN. Les coûts du séquençage ont drastiquement baissés au cours de ces dernières années et les genomes-wides associations studies (GWAS) ont révolutionné l'approche de la recherche génétique en mettant en évidence associations et single-nucleotide-polymorphisms (SNPs) qui pourraient être importantes pour la susceptibilité à développer des maladies dites communes. La majorité des cancers appartiennent à cette définition de maladie commune, ils sont généralement causés par une accumulation de lésions/mutations de l'ADN aboutissant à une perte de contrôle de la prolifération et du cycle cellulaire. Ces mutations peuvent être héréditaires, acquises ou une combinaison des deux. Dans la plupart des cancers communs (cancers qui n'ont pas une hérédité familiale importante) les mutations de l'ADN sont souvent amenées par des facteurs tels que inflammation chronique, tabac, virus, exposition aux radiations, aux agents chimiques. Ceci est le cas pour le mélanome également, un cancer de la peau qui est corrélé à l'exposition des rayons UV solaires ou artificiels. Une hypothèse largement acceptée aujourd'hui est que les tumeurs, à travers leur accumulation progressive de mutations somatiques et d'anomalies chromosomiques, finissent par échapper au contrôle exercé par le système immunitaire. Il est par conséquence imaginable que des polymorphismes naturels puissent renforcer ou affaiblir la capacité du système immunitaire à freiner voir arrêter la progression tumorale.

Design and analysis of ChIP-Seq experiments for nuclear factor I DNA-binding proteins

SUMMARY : Eukaryotic DNA interacts with the nuclear proteins using non-covalent ionic interactions. Proteins can recognize specific nucleotide sequences based on the sterical interactions with the DNA and these specific protein-DNA interactions are the basis for many nuclear processes, e.g. gene transcription, chromosomal replication, and recombination. New technology termed ChIP-Seq has been recently developed for the analysis of protein-DNA interactions on a whole genome scale and it is based on immunoprecipitation of chromatin and high-throughput DNA sequencing procedure. ChIP-Seq is a novel technique with a great potential to replace older techniques for mapping of protein-DNA interactions. In this thesis, we bring some new insights into the ChIP-Seq data analysis. First, we point out to some common and so far unknown artifacts of the method. Sequence tag distribution in the genome does not follow uniform distribution and we have found extreme hot-spots of tag accumulation over specific loci in the human and mouse genomes. These artifactual sequence tags accumulations will create false peaks in every ChIP-Seq dataset and we propose different filtering methods to reduce the number of false positives. Next, we propose random sampling as a powerful analytical tool in the ChIP-Seq data analysis that could be used to infer biological knowledge from the massive ChIP-Seq datasets. We created unbiased random sampling algorithm and we used this methodology to reveal some of the important biological properties of Nuclear Factor I DNA binding proteins. Finally, by analyzing the ChIP-Seq data in detail, we revealed that Nuclear Factor I transcription factors mainly act as activators of transcription, and that they are associated with specific chromatin modifications that are markers of open chromatin. We speculate that NFI factors only interact with the DNA wrapped around the nucleosome. We also found multiple loci that indicate possible chromatin barrier activity of NFI proteins, which could suggest the use of NFI binding sequences as chromatin insulators in biotechnology applications. RESUME : L'ADN des eucaryotes interagit avec les protéines nucléaires par des interactions noncovalentes ioniques. Les protéines peuvent reconnaître les séquences nucléotidiques spécifiques basées sur l'interaction stérique avec l'ADN, et des interactions spécifiques contrôlent de nombreux processus nucléaire, p.ex. transcription du gène, la réplication chromosomique, et la recombinaison. Une nouvelle technologie appelée ChIP-Seq a été récemment développée pour l'analyse des interactions protéine-ADN à l'échelle du génome entier et cette approche est basée sur l'immuno-précipitation de la chromatine et sur la procédure de séquençage de l'ADN à haut débit. La nouvelle approche ChIP-Seq a donc un fort potentiel pour remplacer les anciennes techniques de cartographie des interactions protéine-ADN. Dans cette thèse, nous apportons de nouvelles perspectives dans l'analyse des données ChIP-Seq. Tout d'abord, nous avons identifié des artefacts très communs associés à cette méthode qui étaient jusqu'à présent insoupçonnés. La distribution des séquences dans le génome ne suit pas une distribution uniforme et nous avons constaté des positions extrêmes d'accumulation de séquence à des régions spécifiques, des génomes humains et de la souris. Ces accumulations des séquences artéfactuelles créera de faux pics dans toutes les données ChIP-Seq, et nous proposons différentes méthodes de filtrage pour réduire le nombre de faux positifs. Ensuite, nous proposons un nouvel échantillonnage aléatoire comme un outil puissant d'analyse des données ChIP-Seq, ce qui pourraient augmenter l'acquisition de connaissances biologiques à partir des données ChIP-Seq. Nous avons créé un algorithme d'échantillonnage aléatoire et nous avons utilisé cette méthode pour révéler certaines des propriétés biologiques importantes de protéines liant à l'ADN nommés Facteur Nucléaire I (NFI). Enfin, en analysant en détail les données de ChIP-Seq pour la famille de facteurs de transcription nommés Facteur Nucléaire I, nous avons révélé que ces protéines agissent principalement comme des activateurs de transcription, et qu'elles sont associées à des modifications de la chromatine spécifiques qui sont des marqueurs de la chromatine ouverte. Nous pensons que lés facteurs NFI interagir uniquement avec l'ADN enroulé autour du nucléosome. Nous avons également constaté plusieurs régions génomiques qui indiquent une éventuelle activité de barrière chromatinienne des protéines NFI, ce qui pourrait suggérer l'utilisation de séquences de liaison NFI comme séquences isolatrices dans des applications de la biotechnologie.

Decreasing population selection rates of resistance mutation K65R over time in HIV-1 patients receiving combination therapy including tenofovir.

OBJECTIVES: The use of tenofovir is highly associated with the emergence of mutation K65R, which confers broad resistance to nucleoside/nucleotide analogue reverse transcriptase inhibitors (NRTIs), especially when tenofovir is combined with other NRTIs also selecting for K65R. Although recent HIV-1 treatment guidelines discouraging these combinations resulted in reduced K65R selection with tenofovir, updated information on the impact of currently recommended regimens on the population selection rate of K65R is presently lacking. METHODS: In this study, we evaluated changes over time in the selection rate of resistance mutation K65R in a large population of 2736 HIV-1-infected patients failing combination antiretroviral treatment between 2002 and 2010. RESULTS: The K65R resistance mutation was detected in 144 patients, a prevalence of 5.3%. A large majority of observed K65R cases were explained by the use of tenofovir, reflecting its wide use in clinical practice. However, changing patterns over time in NRTIs accompanying tenofovir resulted in a persistent decreasing probability of K65R selection by tenofovir-based therapy. The currently recommended NRTI combination tenofovir/emtricitabine was associated with a low probability of K65R emergence. For any given dual NRTI combination including tenofovir, higher selection rates of K65R were consistently observed with a non-nucleoside reverse transcriptase inhibitor than with a protease inhibitor as the third agent. DISCUSSION: Our finding of a stable time trend of K65R despite elevated use of tenofovir illustrates increased potency of current HIV-1 therapy including tenofovir.

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Sequence conservation in Plasmodium falciparum alpha-helical coiled coil domains proposed for vaccine development.

BACKGROUND: The availability of the P. falciparum genome has led to novel ways to identify potential vaccine candidates. A new approach for antigen discovery based on the bioinformatic selection of heptad repeat motifs corresponding to alpha-helical coiled coil structures yielded promising results. To elucidate the question about the relationship between the coiled coil motifs and their sequence conservation, we have assessed the extent of polymorphism in putative alpha-helical coiled coil domains in culture strains, in natural populations and in the single nucleotide polymorphism data available at PlasmoDB. METHODOLOGY/PRINCIPAL FINDINGS: 14 alpha-helical coiled coil domains were selected based on preclinical experimental evaluation. They were tested by PCR amplification and sequencing of different P. falciparum culture strains and field isolates. We found that only 3 out of 14 alpha-helical coiled coils showed point mutations and/or length polymorphisms. Based on promising immunological results 5 of these peptides were selected for further analysis. Direct sequencing of field samples from Papua New Guinea and Tanzania showed that 3 out of these 5 peptides were completely conserved. An in silico analysis of polymorphism was performed for all 166 putative alpha-helical coiled coil domains originally identified in the P. falciparum genome. We found that 82% (137/166) of these peptides were conserved, and for one peptide only the detected SNPs decreased substantially the probability score for alpha-helical coiled coil formation. More SNPs were found in arrays of almost perfect tandem repeats. In summary, the coiled coil structure prediction was rarely modified by SNPs. The analysis revealed a number of peptides with strictly conserved alpha-helical coiled coil motifs. CONCLUSION/SIGNIFICANCE: We conclude that the selection of alpha-helical coiled coil structural motifs is a valuable approach to identify potential vaccine targets showing a high degree of conservation.

Repeated unidirectional introgression towards Populus balsamifera in contact zones of exotic and native poplars.

As the evolutionary significance of hybridization is largely dictated by its extent beyond the first generation, we broadly surveyed patterns of introgression across a sympatric zone of two native poplars (Populus balsamifera, Populus deltoides) in Quebec, Canada within which European exotic Populus nigra and its hybrids have been extensively planted since the 1800s. Single nucleotide polymorphisms (SNPs) that appeared fixed within each species were characterized by DNA-sequencing pools of pure individuals. Thirty-five of these diagnostic SNPs were employed in a high-throughput assay that genotyped 635 trees of different age classes, sampled from 15 sites with various degrees of anthropogenic disturbance. The degree of admixture within sampled trees was then assessed through Bayesian clustering of genotypes. Hybrids were present in seven of the populations, with 2.4% of all sampled trees showing spontaneous admixture. Sites with hybrids were significantly more disturbed than pure stands, while hybrids comprised both immature juveniles and trees of reproductive age. All three possible F1s were detected. Advanced-generation hybrids were consistently biased towards P. balsamifera regardless of whether hybridization had occurred with P. deltoides or P. nigra. Gene exchange between P. deltoides and P. nigra was not detected beyond the F1 generation; however, detection of a trihybrid demonstrates that even this apparent reproductive isolation does not necessarily result in an evolutionary dead end. Collectively, results demonstrate the natural fertility of hybrid poplars and suggest that introduced genes could potentially affect the genetic integrity of native trees, similar to that arising from introgression between natives.

IL28B polymorphism is significantly correlated with ifn anti-viral effectiveness already on first day of pegylated interferon-alpha and ribavirin therapy of chronic HCV infection

Background and Aims: Recently, single nucleotide polymorphisms (SNPs) in IL28B were shown to correlate with response to pegylated interferon-a (IFN) and ribavirin therapy of chronic HCV infection. However, the cause for the SNPs effect on therapy response and its application for direct anti-viral (DAV) treatment are not clear. Here, we analyze early HCV kinetics as function of IL28B SNPs to determine its specific effect on viral dynamics. Methods: IL28B SNPs rs8099917, rs12979860 and rs12980275 were genotyped in 252 chronically HCV infected Caucasian naïve patients (67% HCV genotype 1, 28% genotype 2-3) receiving peginterferonalfa- 2a (180 mg/qw) plus ribavirin (1000-1200 mg/qd) in the DITTO study. HCV-RNA was measured (LD = 50 IU/ml) frequently during first 28 days. Results: RVR was achieved in 33% of genotype 1 patients with genotype CC at rs12979860 versus 12-16% for genotypes TT and CT (P < 0.03). Significant (P < 0.001) difference in viral decline was observed already at day 1 (see Figure). First phase decline was significantly (P < 0.001) larger in patients with genotype CC (2.0 log) than for TT and CT genotypes (0.6 and 0.8), indicating IFN anti-viral effectiveness in blocking virion production of 99% versus 75-84%. There was no significant association between second phase slope and rs12979860 genotype in patients with a first phase decline larger than 1 log. HCV kinetics as function of IL28b SNP. The same trend (not shown) was observed for HCV genotype 2-3 patients with different SNP genotype distribution that may indicate differential selection pressure as function of HCV genotype. Similar results were observed for SNPs rs8099917 and rs12980275, with a strong linkage disequilibrium among the 3 loci allowing to define the composite haplotype best associated with IFN effectiveness. Conclusions: IFN effectiveness in blocking virion production/ release is strongly affected by IL28B SNPs, but not other viral dynamic properties such as infected cell loss rate. Thus, IFN based therapy, as standard-of-care or in combination with DAV, should consider IL28B SNPs for prediction and personalized treatment, while response to pure DAV treatment may be less affected by IL28B SNPs. Additional analyses are undergoing to pinpoint the SNP effect on IFN anti-viral effectiveness.

A mega-analysis of genome-wide association studies for major depressive disorder.

Prior genome-wide association studies (GWAS) of major depressive disorder (MDD) have met with limited success. We sought to increase statistical power to detect disease loci by conducting a GWAS mega-analysis for MDD. In the MDD discovery phase, we analyzed more than 1.2 million autosomal and X chromosome single-nucleotide polymorphisms (SNPs) in 18 759 independent and unrelated subjects of recent European ancestry (9240 MDD cases and 9519 controls). In the MDD replication phase, we evaluated 554 SNPs in independent samples (6783 MDD cases and 50 695 controls). We also conducted a cross-disorder meta-analysis using 819 autosomal SNPs with P<0.0001 for either MDD or the Psychiatric GWAS Consortium bipolar disorder (BIP) mega-analysis (9238 MDD cases/8039 controls and 6998 BIP cases/7775 controls). No SNPs achieved genome-wide significance in the MDD discovery phase, the MDD replication phase or in pre-planned secondary analyses (by sex, recurrent MDD, recurrent early-onset MDD, age of onset, pre-pubertal onset MDD or typical-like MDD from a latent class analyses of the MDD criteria). In the MDD-bipolar cross-disorder analysis, 15 SNPs exceeded genome-wide significance (P<5 × 10(-8)), and all were in a 248 kb interval of high LD on 3p21.1 (chr3:52 425 083-53 822 102, minimum P=5.9 × 10(-9) at rs2535629). Although this is the largest genome-wide analysis of MDD yet conducted, its high prevalence means that the sample is still underpowered to detect genetic effects typical for complex traits. Therefore, we were unable to identify robust and replicable findings. We discuss what this means for genetic research for MDD. The 3p21.1 MDD-BIP finding should be interpreted with caution as the most significant SNP did not replicate in MDD samples, and genotyping in independent samples will be needed to resolve its status.

Host innate immune responses to microbial pathogens.

Sepsis is among the leading causes of death worldwide and its incidence is increasing. Defined as the host response to infection, sepsis is a clinical syndrome considered to be the expression of a dysregulated immune reaction induced by danger signals that may lead to organ failure and death. Remarkable progresses have been made in our understanding of the molecular basis of host defenses in recent years. The host defense response is initiated by innate immune sensors of danger signals designated under the collective name of pattern-recognition receptors. Members of the family of microbial sensors include the complement system, the Toll-like receptors, the nucleotide-binding oligomerization domainlike receptors, the RIG-I-like helicases and the C-type lectin receptors. Ligand-activated pattern-recognition receptors kick off a cascade of intracellular events resulting in the expression of co-stimulatory molecules and release of effector molecules playing a fundamental role in the initiation of the innate and adaptive immune responses. Fine tuning of proinflammatory and anti-inflammatory reactions is critical for keeping the innate immune response in check. Overwhelming or dysregulated responses induced by infectious stimuli may have dramatic consequences for the host as shown by the profound derangements observed in sepsis. Unfortunately, translational research approaches aimed at the development of therapies targeting newly identified innate immune pathways have not held their promises. Indeed, all recent clinical investigations of adjunctive anti-sepsis treatments had little, if any, impact on morbidity and all-cause mortality of sepsis. Dissecting the mechanisms underlying the transition from infection to sepsis is essential for solving the sepsis enigma. Important components of the puzzle have already been identified, but the hunt must go on in the laboratory and at the bedside.

Disease-causing mutations improving the branch site and polypyrimidine tract: pseudoexon activation of LINE-2 and antisense Alu lacking the poly(T)-tail.

Cryptic exons or pseudoexons are typically activated by point mutations that create GT or AG dinucleotides of new 5' or 3' splice sites in introns, often in repetitive elements. Here we describe two cases of tetrahydrobiopterin deficiency caused by mutations improving the branch point sequence and polypyrimidine tracts of repeat-containing pseudoexons in the PTS gene. In the first case, we demonstrate a novel pathway of antisense Alu exonization, resulting from an intronic deletion that removed the poly(T)-tail of antisense AluSq. The deletion brought a favorable branch point sequence within proximity of the pseudoexon 3' splice site and removed an upstream AG dinucleotide required for the 3' splice site repression on normal alleles. New Alu exons can thus arise in the absence of poly(T)-tails that facilitated inclusion of most transposed elements in mRNAs by serving as polypyrimidine tracts, highlighting extraordinary flexibility of Alu repeats in shaping intron-exon structure. In the other case, a PTS pseudoexon was activated by an A>T substitution 9 nt upstream of its 3' splice site in a LINE-2 sequence, providing the first example of a disease-causing exonization of the most ancient interspersed repeat. These observations expand the spectrum of mutational mechanisms that introduce repetitive sequences in mature transcripts and illustrate the importance of intronic mutations in alternative splicing and phenotypic variability of hereditary disorders.

Immunogold labeling and cerium cytochemistry of the enzyme ecto-5'-nucleotidase in promastigote forms of Leishmania species

We have applied both enzyme cytochemistry and immunological labeling techniques to characterize the enzyme 5'-nucleotidase (5'-Nase), at the ultrastructural level, in promastigote forms of four Leishmania species: Leishmania amazonensis, Leishmania mexicana, Leishmania donovani and Leishmania chagasi. The cerium phosphate staining was localized at the surface of the cell body, the flagellum and the flagellar pocket membranes of all the parasites studied. The immunogold labelling technique confirmed these results. In this report we localized 5'-Nase in L. chagasi and L. amazonensis which have been implicated respectively in visceral and cutaneous forms of leishmaniasis. In addition, we confirmed the localization of this phosphomonoesterase in the other two species studied. The superior quality of the images, obtained with both methodologies, confirms that these parasites possess mechanisms capable of hydrolyzing nucleotide monophosphates, and that the expression of 5'-Nase is associated with the outer surface of the plasma membrane.

Detecting epistasis with restricted response patterns in pairs of biallelic loci.

Well-established examples of genetic epistasis between a pair of loci typically show characteristic patterns of phenotypic distributions in joint genotype tables. However, inferring epistasis given such data is difficult due to the lack of power in commonly used approaches, which decompose the epistatic patterns into main plus interaction effects followed by testing the interaction term. Testing additive-only or all terms may have more power, but they are sensitive to nonepistatic patterns. Alternatively, the epistatic patterns of interest can be enumerated and the best matching one is found by searching through the possibilities. Although this approach requires multiple testing correction over possible patterns, each pattern can be fitted with a regression model with just one degree of freedom and thus the overall power can still be high, if the number of possible patterns is limited. Here we compare the power of the linear decomposition and pattern search methods, by applying them to simulated data generated under several patterns of joint genotype effects with simple biological interpretations. Interaction-only tests are the least powerful; while pattern search approach is the most powerful if the range of possibilities is restricted, but still includes the true pattern.