Probabilistic base calling of Solexa sequencing data.

BACKGROUND: Solexa/Illumina short-read ultra-high throughput DNA sequencing technology produces millions of short tags (up to 36 bases) by parallel sequencing-by-synthesis of DNA colonies. The processing and statistical analysis of such high-throughput data poses new challenges; currently a fair proportion of the tags are routinely discarded due to an inability to match them to a reference sequence, thereby reducing the effective throughput of the technology. RESULTS: We propose a novel base calling algorithm using model-based clustering and probability theory to identify ambiguous bases and code them with IUPAC symbols. We also select optimal sub-tags using a score based on information content to remove uncertain bases towards the ends of the reads. CONCLUSION: We show that the method improves genome coverage and number of usable tags as compared with Solexa's data processing pipeline by an average of 15%. An R package is provided which allows fast and accurate base calling of Solexa's fluorescence intensity files and the production of informative diagnostic plots.

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone.

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long-standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well-studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best-supported scenario was a long period of allopatric isolation during the first three-quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.

Impact of social structure variation on population genetics, social conflicts and life histories in ants

Summary The evolution of social structures and breeding systems in animals is a complex process that combines ecological, genetical and social factors. This thesis sheds light on important changes in population genetics, life-history and social behavior that are associated with variation in social structure in ants. The socially polymorphic ant Formica selysi was chosen as the model organism because single- and multiple-queen colonies occur in close proximity within a single large population. The shift from single- to multiple-queen colonies is generally associated with profound changes in dispersal behavior and mode of colony founding. In chapter 1, we examine the genetic consequences of variation in social structure at both the colony and population levels. A detailed microsatellite analysis reveals that both colony types have similar mating systems, with few or no queen turnover. Furthermore, the complete lack of genetic differentiation observed between single- and multiple-queen colonies provides no support to the hypothesis that change in queen number leads to restricted gene flow between social forms. Besides changes in the genetic composition of the colony, the variation in the number of queens per colony is associated with changes in a network of behavioral and life-history traits that have been described as forming a "polygyny syndrome". In chapter 2, we demonstrate that multiple-queen colonies profoundly differ from single-queen ones in terms of size, nest density and lifespan of colonies, in weight of queens produced, as well as in allocation to reproductive individuals relative to workers. These multifaceted changes in life-history traits can provide various fitness benefits to members of multiple-queen colonies. Increasing the number of queens in a colony usually results in a decreased level of aggression towards non-nestmates. The phenotype matching hypothesis predicts that, compared to single-queen colonies, multiple-queen colonies have more diverse genetically-derived cues used for recognition, resulting in a lower ability to discriminate non-nestmates. In sharp contrast to this hypothesis, we show in chapter 3 that single- and multiple-queen colonies exhibit on average similar levels of aggression. Moreover, stronger aggression is recorded between colonies of different social structure than between colonies of the same social structure. Several hypotheses propose that the evolution of multiple-queen colonies is at least partly due to benefits resulting from an increase in colony genetic diversity. The task-efficiency hypothesis holds that genetic variation improves task performance due to a more complete or more sensitive expression of the genetically-based division of labor. In .chapter 4, we evaluate if higher colony genetic diversity increases worker size polymorphism and thus may improve division of labor. We show that despite the fact that worker size has a heritable component, higher levels of genetic diversity do not result in more polymorphic workers. The smaller size and lower polymorphism levels of workers of multiple-queen colonies compared to single-queen ones further indicate that an increase in colony genetic diversity does not increase worker size polymorphism but might improve colony homeostasis. In chapter 5, we provide clear evidence for an ongoing conflict between queens and workers on sex allocation, as predicted by kin selection theory. Our data show that queens of F. selysi strongly influence colony sex allocation by biasing the sex ratio of their eggs. However, there is also evidence that workers eliminated some male brood, resulting in a population sex-investment ratio that is between the queens' and workers' equilibria. Résumé L'évolution des structures sociales et systèmes d'accouplement chez les animaux est un processus complexe combinant à la fois des facteurs écologiques, génétiques et sociaux. Cette thèse met en lumière des changements importants dans la génétique des populations, les traits d'histoire de vie et les comportements sociaux qui sont associés à des variations de structure sociale chez les fourmis. Durant ce travail, nous avons étudié une population de Formica selysi composée à la fois de colonies à une reine et de colonies à plusieurs reines. La transition de colonie à une reine à colonie à plusieurs reines est généralement associée à des changements profonds dans le comportement de dispersion ainsi que le mode de fondation des sociétés. Dans le chapitre 1, nous examinons les conséquences génétiques de la variation de structure sociale tant au niveau de la colonie qu'au niveau de la population. Une analyse détaillée à l'aide de marqueurs microsatellites nous révèle que les deux types de colonies ont des systèmes d'accouplements similaires avec peu ou pas de renouvellement de reines. L'absence totale de différenciation génétique entre les colonies à une et à plusieurs reines n'apporte aucun support à l'hypothèse selon laquelle un changement dans le nombre de reines conduit à un flux de gènes restreint entre les deux formes sociales. A côté de changements dans la composition génétique de la colonie, la variation du nombre de reines dans une colonie est associée à une multitude de changements comportementaux et de traits d'histoire de vie qui ont été décrits comme formant un "syndrome polygyne". Dans le chapitre 2, nous démontrons que les colonies à plusieurs reines diffèrent profondément des colonies à une reine en terme de taille, densité de nids, longévité des colonies, poids des nouvelles reines produites ainsi que dans l'allocation entre les individus reproducteurs et les ouvrières. Ces changements multiples dans les traits d'histoire de vie peuvent apporter des bénéfices variés en terme de fitness aux colonies à plusieurs reines. L'augmentation du nombre de reines dans une colonie est généralement associée à une baisse du degré d'agressivité envers les fourmis étrangères au nid. L'hypothèse "phénotype matching" prédit que les colonies à plusieurs reines ont une plus grande diversité dans les facteurs d'origine génétique utilisés pour la reconnaissance, résultant en une capacité diminuée à discriminer une fourmi étrangère au nid. Contrairement à cette hypothèse, nous montrons dans le chapitre 3 que les colonies à une et à plusieurs reines ont des niveaux d'agressivité similaires. De plus, une agressivité accrue est observée entre colonies de structures sociales différentes comparée à des colonies de même structure sociale. Plusieurs hypothèses ont proposé que l'évolution de colonies ä plusieurs reines soit en partie due aux bénéfices résultant d'une augmentation de la diversité génétique dans la colonie. L'hypothèse "task efficiency" prédit que la diversité génétique améliore l'efficacité à effectuer certaines tâches grâce à une expression plus complète et plus souple d'une division du travail génétiquement déterminée. Nous évaluons dans le chapitre 4 si un accroissement de la diversité génétique augmente le polymorphisme de taille des ouvrières, d'où peut ainsi découler une meilleure division du travail. Nous montrons qu'en dépit du fait que la taille des ouvrières soit un caractère héritable, une forte diversité génétique ne se traduit pas par un plus fort polymorphisme chez les ouvrières. Les ouvrières de colonies à plusieurs reines sont plus petites et moins polymorphes que celles des colonies à une seule reine. Dans le chapitre 5, nous démontrons l'existence d'un conflit ouvert entre reines et ouvrières à propos de l'allocation dans les sexes, comme le prédit la théorie de la sélection de parentèle. Nos données révèlent que les reines de F. selysi influencent fortement l'allocation dans les sexes en biaisant la sexe ratio des oeufs. Cependant, certains indices indiquent que les ouvrières éliminent une partie du couvain mâle, ce qui a pour effet d'avoir un investissement dans les sexes au niveau de la population intermédiaire entre les intérêts des reines et des ouvrières.

Constant p53 pathway inactivation in a large series of soft tissue sarcomas with complex genetics.

Alterations of the p53 pathway are among the most frequent aberrations observed in human cancers. We have performed an exhaustive analysis of TP53, p14, p15, and p16 status in a large series of 143 soft tissue sarcomas, rare tumors accounting for around 1% of all adult cancers, with complex genetics. For this purpose, we performed genomic studies, combining sequencing, copy number assessment, and expression analyses. TP53 mutations and deletions are more frequent in leiomyosarcomas than in undifferentiated pleomorphic sarcomas. Moreover, 50% of leiomyosarcomas present TP53 biallelic inactivation, whereas most undifferentiated pleomorphic sarcomas retain one wild-type TP53 allele (87.2%). The spectrum of mutations between these two groups of sarcomas is different, particularly with a higher rate of complex mutations in undifferentiated pleomorphic sarcomas. Most tumors without TP53 alteration exhibit a deletion of p14 and/or lack of mRNA expression, suggesting that p14 loss could be an alternative genotype for direct TP53 inactivation. Nevertheless, the fact that even in tumors altered for TP53, we could not detect p14 protein suggests that other p14 functions, independent of p53, could be implicated in sarcoma oncogenesis. In addition, both p15 and p16 are frequently codeleted or transcriptionally co-inhibited with p14, essentially in tumors with two wild-type TP53 alleles. Conversely, in TP53-altered tumors, p15 and p16 are well expressed, a feature not incompatible with an oncogenic process.

Effects of long-term averaging of quantitative blood pressure traits on the detection of genetic associations.

Blood pressure (BP) is a heritable, quantitative trait with intraindividual variability and susceptibility to measurement error. Genetic studies of BP generally use single-visit measurements and thus cannot remove variability occurring over months or years. We leveraged the idea that averaging BP measured across time would improve phenotypic accuracy and thereby increase statistical power to detect genetic associations. We studied systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP), and pulse pressure (PP) averaged over multiple years in 46,629 individuals of European ancestry. We identified 39 trait-variant associations across 19 independent loci (p < 5 × 10(-8)); five associations (in four loci) uniquely identified by our LTA analyses included those of SBP and MAP at 2p23 (rs1275988, near KCNK3), DBP at 2q11.2 (rs7599598, in FER1L5), and PP at 6p21 (rs10948071, near CRIP3) and 7p13 (rs2949837, near IGFBP3). Replication analyses conducted in cohorts with single-visit BP data showed positive replication of associations and a nominal association (p < 0.05). We estimated a 20% gain in statistical power with long-term average (LTA) as compared to single-visit BP association studies. Using LTA analysis, we identified genetic loci influencing BP. LTA might be one way of increasing the power of genetic associations for continuous traits in extant samples for other phenotypes that are measured serially over time.

Population pharmacokinetic analysis and pharmacogenetics of raltegravir in HIV-positive and healthy individuals.

The objectives of this study were to characterize raltegravir (RAL) population pharmacokinetics in HIV-positive (HIV(+)) and healthy individuals, identify influential factors, and search for new candidate genes involved in UDP glucuronosyltransferase (UGT)-mediated glucuronidation. The pharmacokinetic analysis was performed with NONMEM. Genetic association analysis was performed with PLINK using the relative bioavailability as the phenotype. Simulations were performed to compare once- and twice-daily regimens. A 2-compartment model with first-order absorption adequately described the data. Atazanavir, gender, and bilirubin levels influenced RAL relative bioavailability, which was 30% lower in HIV(+) than in healthy individuals. UGT1A9*3 was the only genetic variant possibly influencing RAL pharmacokinetics. The majority of RAL pharmacokinetic variability remains unexplained by genetic and nongenetic factors. Owing to the very large variability, trough drug levels might be very low under the standard dosing regimen, raising the question of a potential relevance of therapeutic drug monitoring of RAL in some situations.

Swain: Stata module to correct the SEM chi-square overidentification test in small sample sizes or complex models. Statistical Software Components S457617, Boston College Department of Economics.

Swain corrects the chi-square overidentification test (i.e., likelihood ratio test of fit) for structural equation models whethr with or without latent variables. The chi-square statistic is asymptotically correct; however, it does not behave as expected in small samples and/or when the model is complex (cf. Herzog, Boomsma, & Reinecke, 2007). Thus, particularly in situations where the ratio of sample size (n) to the number of parameters estimated (p) is relatively small (i.e., the p to n ratio is large), the chi-square test will tend to overreject correctly specified models. To obtain a closer approximation to the distribution of the chi-square statistic, Swain (1975) developed a correction; this scaling factor, which converges to 1 asymptotically, is multiplied with the chi-square statistic. The correction better approximates the chi-square distribution resulting in more appropriate Type 1 reject error rates (see Herzog & Boomsma, 2009; Herzog, et al., 2007).

Adaptive strategy for the statistical analysis of connectomes.

We study an adaptive statistical approach to analyze brain networks represented by brain connection matrices of interregional connectivity (connectomes). Our approach is at a middle level between a global analysis and single connections analysis by considering subnetworks of the global brain network. These subnetworks represent either the inter-connectivity between two brain anatomical regions or by the intra-connectivity within the same brain anatomical region. An appropriate summary statistic, that characterizes a meaningful feature of the subnetwork, is evaluated. Based on this summary statistic, a statistical test is performed to derive the corresponding p-value. The reformulation of the problem in this way reduces the number of statistical tests in an orderly fashion based on our understanding of the problem. Considering the global testing problem, the p-values are corrected to control the rate of false discoveries. Finally, the procedure is followed by a local investigation within the significant subnetworks. We contrast this strategy with the one based on the individual measures in terms of power. We show that this strategy has a great potential, in particular in cases where the subnetworks are well defined and the summary statistics are properly chosen. As an application example, we compare structural brain connection matrices of two groups of subjects with a 22q11.2 deletion syndrome, distinguished by their IQ scores.

Estimating the basic reproductive number from viral sequence data.

Epidemiological processes leave a fingerprint in the pattern of genetic structure of virus populations. Here, we provide a new method to infer epidemiological parameters directly from viral sequence data. The method is based on phylogenetic analysis using a birth-death model (BDM) rather than the commonly used coalescent as the model for the epidemiological transmission of the pathogen. Using the BDM has the advantage that transmission and death rates are estimated independently and therefore enables for the first time the estimation of the basic reproductive number of the pathogen using only sequence data, without further assumptions like the average duration of infection. We apply the method to genetic data of the HIV-1 epidemic in Switzerland.

Adaptation and the genetics of social behaviour.

In recent years much progress has been made towards understanding the selective forces involved in the evolution of social behaviour including conflicts over reproduction among group members. Here, I argue that an important additional step necessary for advancing our understanding of the resolution of potential conflicts within insect societies is to consider the genetics of the behaviours involved. First, I discuss how epigenetic modifications of behaviour may affect conflict resolution within groups. Second, I review known natural polymorphisms of social organization to demonstrate that a lack of consideration of the genetic mechanisms involved may lead to erroneous explanations of the adaptive significance of behaviour. Third, I suggest that, on the basis of recent genetic studies of sexual conflict in Drosophila, it is necessary to reconsider the possibility of within-group manipulation by means of chemical substances (i.e. pheromones). Fourth, I address the issue of direct versus indirect genetic effects, which is of particular importance for the study of behaviour in social groups. Fifth, I discuss the issue of how a genetic influence on dominance hierarchies and reproductive division of labour can have secondary effects, for example in the evolution of promiscuity. Finally, because the same sets of genes (e.g. those implicated in chemical signalling and the responses that are triggered) may be used even in species as divergent as ants, cooperative breeding birds and primates, an integration of genetic mechanisms into the field of social evolution may also provide unifying ideas.

An investigation of the potential of DIP-STR markers for DNA mixture analyses

The genetic characterization of unbalanced mixed stains remains an important area where improvement is imperative. In fact, using the standard tools of forensic DNA profiling (i.e., STR markers), the profile of the minor contributor in mixed DNA stains cannot be successfully detected if its quantitative share of DNA is less than 10% of the mixed trace. This is due to the fact that the major contributor's profile "masks" that of the minor contributor. Besides known remedies to this problem, such as Y-STR analysis, a new compound genetic marker that consists of a Deletion/Insertion Polymorphism (DIP) linked to a Short Tandem Repeat (STR) polymorphism, has recently been developed and proposed [1]. These novel markers are called DIP-STR markers. This paper compares, from a statistical and forensic perspective, the potential usefulness of these novel DIP-STR markers (i) with traditional STR markers in cases of moderately unbalanced mixtures, and (ii) with Y-STR markers in cases of female-male mixtures. This is done through a comparison of the distribution of 100,000 likelihood ratio values obtained using each method on simulated mixtures. This procedure is performed assuming, in turn, the prosecution's and the defence's point of view.

Genetics of essential tremor.

Essential tremor (ET) is a prevalent condition manifesting with progressive action tremor. Although ET was traditionally viewed as a sporadic disease, a significant proportion of cases report a positive family history of tremor. Autosomal dominant inheritance can be demonstrated in many families. Previously, genome-wide linkage studies in families mapped three loci for ET, hereditary essential tremor-1 (ETM1), ETM2 and ETM3. However, no causal mutation has been replicated in candidate genes within these loci, including dopamine D3 receptor (DRD3) and HS1-binding protein 3 (HS1BP3). Recently, the first genome-wide association study in ET followed by replication studies conducted in diverse populations identified a significant association between the leucine-rich repeat and Ig domain containing 1 gene (LINGO1) SNP rs9652490 and risk for ET Although further novel variants were indentified in LINGO1 and its paralog LINGO2 that may be associated with risk for ET, the pathogenic mechanisms involved remain elusive. Given the possibility that ET as a complex trait may be influenced by the combined effects of rare variants, novel high-throughput technologies sequencing all exons across the genome (exome sequencing) or the whole genome (genome sequencing) may become crucial in understanding/deciphering the genetic background of ET.

Association of a common vitamin D-binding protein polymorphism with inflammatory bowel disease.

OBJECTIVE: Inflammatory bowel diseases (IBDs), Crohn's disease, and ulcerative colitis (UC), are multifactorial disorders, characterized by chronic inflammation of the intestine. A number of genetic components have been proposed to contribute to IBD pathogenesis. In this case-control study, we investigated the association between two common vitamin D-binding protein (DBP) genetic variants and IBD susceptibility. These two single nucleotide polymorphisms (SNPs) in exon 11 of the DBP gene, at codons 416 (GAT>GAG; Asp>Glu) and 420 (ACG>AAG; Thr>Lys), have been previously suggested to play roles in the etiology of other autoimmune diseases. METHODS: Using TaqMan SNP technology, we have genotyped 884 individuals (636 IBD cases and 248 non-IBD controls) for the two DBP variants. RESULTS: On statistical analysis, we observed that the DBP 420 variant Lys is less frequent in IBD cases than in non-IBD controls (allele frequencies, P=0.034; homozygous carrier genotype frequencies, P=0.006). This inverse association between the DBP 420 Lys and the disease remained significant, when non-IBD participants were compared with UC (homozygous carrier genotype frequencies, P=0.022) or Crohn's disease (homozygous carrier genotype frequencies, P=0.016) patients separately. Although the DBP position 416 alone was not found to be significantly associated with IBD, the haplotype DBP_2, consisting of 416 Asp and 420 Lys, was more frequent in the non-IBD population, particularly notably when compared with the UC group (Odds ratio, 4.390). CONCLUSION: Our study adds DBP to the list of potential genes that contribute to the complex genetic etiology of IBD, and further emphasizes the association between vitamin D homeostasis and intestinal inflammation.

Evolutionary genetics of arbuscular mycorrhizal fungi - causes and consequences of multiple genomes in Glomus intraradices

Résumé Les champignons endomycorhiziens arbusculaires (CEA) forment des symbioses avec la plupart des plantes terrestres. Les CEA influencent la croissance des plantes et la biodiversité. Ils sont supposés avoir évolué de manière asexuée pendant au moins 400 millions d'années et aucune diversification morphologique majeure n'a été constatée. Pour ces raisons, les CEA sont considérés comme d'anciens asexués. Très peu d'espèces sont connues actuellement. Les individus de ces champignons contiennent des noyaux génétiquement différents dans un cytoplasme continu. La signification évolutive, la variabilité et la maintenance des génomes multiples au sein des individus sont inconnues. Ce travail a démontré qu'une population du CEA Glomus intraradices est génétiquement très variable. Nous avons conclu que les plantes hôtes plutôt que la différenciation géographique devraient être responsables de cette grande diversité. Puis nous avons cherché l'existence de recombinaison entre génotypes dans une population. Nous avons détecté un groupe recombinant au sein de la population, ce qui met en doute l'état d'anciens asexués des CEA. Nous avons également détecté l'occurrence de fusions d'hyphes et l'échange de noyaux entre isolats génétiquement différents. La descendance hybride issue de cet échange était viable et distincte phénotypiquement des isolats parentaux. En résumé, ce travail identifie des événements cruciaux dans le cycle de vie des CEA qui ont le potentiel d'influencer l'évolution de génomes multiples. L'étude des conséquences de ces événements sur les interactions avec les plantes hôtes pourrait éclaircir significativement la compréhension de la symbiose entre plantes et CEA. Abstract Arbuscular mycorrhizal fungi (AMF) are important symbionts of most land plants. AMF influence plant growth and biodiversity. Very few extant species are described. AMF are thought to have evolved asexually for at least 400 million years and no major morphological diversification has occurred. Due to these reasons, they were termed `ancient asexuals'. Fungal individuals harbour genetically different nuclei in a continuous cytoplasm. The variability, maintenance and evolutionary significance of multiple genomes within individuals are unknown. This work showed that a population of the AMF Glomus intraradices harbours very high genetic diversity. We concluded that host plants rather than geographic differentiation were responsible for this diversity. Furthermore, we investigated whether recombination occurred among genotypes of a G. intraradices population. The identification of a core group of recombining genotypes in the population refutes the assumption of ancient asexuality in AMF. We found that genetically different isolates can form hyphal fusions and exchange nuclei. The hybrid progeny produced by the exchange was viable and phenotypically distinct from the parental isolates. Taken together, this work provided evidence for key events in the AMF life cycle, that influence the evolution of multiple genomes. Studying the consequences of these events on the interaction with host plants may significantly further the understanding of the AMF-plant symbiosis.