Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Effects of genetic alterations in arbuscular mycorrhizal fungi on plant growth and on their response to a change of host

Abstract :The majority of land plants form the symbiosis with arbuscular mycorrhizal fungi (AMF). The AM symbiosis has existed for hundreds of millions of years but little or no specificity seems to have co- evolved between the partners and only about 200 morphospecies of AMF are known. The fungi supply the plants most notably with phosphate in exchange for carbohydrates. The fungi improve plant growth, protect them against pathogens and herbivores and the symbiosis plays a key role in ecosystem productivity and plant diversity. The fungi are coenocytic, grow clonally and no sexual stage in their life cycle is known. For these reasons, they are presumed ancient asexuals. Evidence suggests that AMF contain populations of genetically different nucleotypes coexisting in a common cytoplasm. Consequently, the nucleotype content of new clonal offspring could potentially be altered by segregation of nuclei at spore formation and by genetic exchange between different AMF. Given the importance of AMF, it is surprising that remarkably little is known about the genetics and genomics of the fungi.The main goal of this thesis was to investigate the combined effects of plant species differences and of genetic exchange and segregation in AMF on the symbiosis. This work showed that single spore progeny can receive a different assortment of nucleotypes compared to their parent and compared to other single spore progeny. This is the first direct evidence that segregation occurs in AMF. We then showed that both genetic exchange and segregation can lead to new progeny that differentially alter plant growth compared to their parents. We also found that genetic exchange and segregation can lead to different development of the fungus during the establishment of the symbiosis. Finally, we found that a shift of host species can differentially alter the phenotypes and genotypes of AMF progeny obtained by genetic exchange and segregation compared to their parents.Overall, this study confirms the multigenomic state of the AMF Glomus intraradices because our findings are possible only if the fungus contains genetically different nuclei. We demonstrated the importance of the processes of genetic exchange and segregation to produce, in a very short time span, new progeny with novel symbiotic effects. Moreover, our results suggest that different host species could affect the fate of different nucleotypes following genetic exchange and segregation in AMF, and can potentially contribute to the maintenance of genetic diversity within AMF individuals. This work brings new insights into understanding how plants and fungi have coevolved and how the genetic diversity in AMF can be maintained. We recommend that the intra-ir1dividual AMF diversity and these processes should be considered in future research on this symbiosis.Résumé :La majorité des plantes terrestres forment des symbioses avec les champignons endomycorhiziens arbusculaires (CEA). Cette symbiose existe depuis plusieurs centaines de millions d'années mais peu ou pas de spécificité semble avoir co-évoluée entre les partenaires et seulement 200 morpho-espèces de CEA sont connues. Le champignon fournit surtout aux plantes du phosphate en échange de carbohydrates. Le champignon augmente la croissance des plantes, les protège contre des pathogènes et herbivores et la symbiose joue un rôle clé dans la productivité des écosystèmes et de la diversité des plantes. Les CEA sont coenocytiques, se reproduisent clonalement et aucune étape sexuée n'est connue dans leur cycle de vie. Pour ces raisons, ils sont présumés comme anciens asexués. Des preuves suggèrent que les CEA ont des populations de nucleotypes différents coexistant dans un cytoplasme commun. Par conséquent, le contenu en nucleotype des nouveaux descendants clonaux pourrait être altéré par la ségrégation des noyaux lors de la fonnation des spores et par l'échange génétique entre différents CEA. Etant donné l'importance des CEA, il est surprenant que si peu soit connu sur la génétique et la génomique du champignon.Le principal but de cette thèse a été d'étudier les effets combinés de différentes espèces de plantes et des mécanismes d'échange génétique et de ségrégation chez les CEA sur la symbiose. Ce travail a montré que chaque nouvelle spore produite pouvait recevoir un assortiment différent de noyaux comparé au parent ou comparé à d'autres nouvelles spores. Ceci est la première preuve directe que la ségrégation peut se produire chez les CEA. Nous avons ensuite montré qu'à la fois l'échange génétique et la ségrégation pouvaient mener à de nouveaux descendants qui altèrent différemment la croissance des plantes, comparé à leurs parents. Nous avons également trouvé que l'échange génétique et la ségrégation pouvaient entraîner des développements différents du champignon pendant l'établissement de la symbiose. Pour finir, nous avons trouvé qu'un changement d'espèce de l'hôte pouvait altérer différemment les phénotypes et génotypes des descendants issus d'échange génétique et de ségrégation, comparé à leurs parents.Globalement, cette étude confirme l'état multigénomique du CEA Glumus intraradices car nous résultats sont possibles seulement si le champignon possède des noyaux génétiquement différents. Nous avons démontrés l'importance des mécanismes d'échange génétique et de ségrégation pour produire en très peu de temps de nouveaux descendants ayant des effets symbiotiques nouveaux. De plus, nos résultats suggèrent que différentes espèces de plantes peuvent agir sur le devenir des nucleotypes après l'échange génétique et la ségrégation chez les CEA, et pourraient contribuer à la maintenance de la diversité génétique au sein d'un même CEA. Ce travail apporte des éléments nouveaux pour comprendre comment les plantes et les champignons ont coévolué et comment la diversité génétique chez les CEA peut être maintenue. Nous recommandons de considérer la diversité génétique intra-individuelle des CEA et ces mécanismes lors de futures recherches sur cette symbiose.

Differentiation and genetic analysis of Rhodnius prolixus and Rhodnius colombiensis by rDNA and RAPD amplification

Domiciliated Rhodnius prolixus and sylvatic R. colombiensis were analyzed in order to confirm their genetic divergence and verify the risk that the latter represents in the domiciliation process, and to provide tools for identifying the sources of possible reinfestation by triatomines in human dwellings allowing control programs to be undertaken. Comparison of random amplified polymorphic DNA amplification patterns and cluster analysis suggests reproductive discontinuity between the two species. The calculated statistical F value of 0.24 and effective migration rate of 0.6 individuals per generation are insufficient to maintain genetic homogeneity between them and confirm the absence of present genetic flow. R. colombiensis presents higher intrapopulation variability. Polymerase chain reaction of ribosomal DNA supports these findings. The low genetic flow between the two species implies that R. colombiensis do not represent an epidemiological risk for the domiciliary transmission of Trypanosoma cruzi in the Tolima Department. The lower variability of the domiciliated R. prolixus could result in greater susceptibility to the use of pesticides in control programs.

Genetic and environmental components of phenotypic variation in immune response and body size of a colonial bird, Delichon urbica (the house martin).

Directional selection for parasite resistance is often intense in highly social host species. Using a partial cross-fostering experiment we studied environmental and genetic variation in immune response and morphology in a highly colonial bird species, the house martin (Delichon urbica). We manipulated intensity of infestation of house martin nests by the haematophagous parasitic house martin bug Oeciacus hirundinis either by spraying nests with a weak pesticide or by inoculating them with 50 bugs. Parasitism significantly affected tarsus length, T cell response, immunoglobulin and leucocyte concentrations. We found evidence of strong environmental effects on nestling body mass, body condition, wing length and tarsus length, and evidence of significant additive genetic variance for wing length and haematocrit. We found significant environmental variance, but no significant additive genetic variance in immune response parameters such as T cell response to the antigenic phytohemagglutinin, immunoglobulins, and relative and absolute numbers of leucocytes. Environmental variances were generally greater than additive genetic variances, and the low heritabilities of phenotypic traits were mainly a consequence of large environmental variances and small additive genetic variances. Hence, highly social bird species such as the house martin, which are subject to intense selection by parasites, have a limited scope for immediate microevolutionary response to selection because of low heritabilities, but also a limited scope for long-term response to selection because evolvability as indicated by small additive genetic coefficients of variation is weak.

Antiretroviral resistance and genetic diversity of human immunodeficiency virus type 1 isolates from the Federal District, Central Brazil

In the context of universal access to antiretroviral therapy, the surveillance of human immunodeficiency virus type 1 (HIV-1) genetic diversity and resistance becomes pivotal. In this work our purpose was to describe the genetic variability; prevalence of drug-resistance mutations; and genotypic resistance profiles in HIV-1 infected individuals under antiretroviral treatment, from the Federal District, Brasília, Central Brazil. The entire viral protease and codons 19 to 234 of the reverse transcriptase gene from 45 HIV-1 isolates were amplified and sequenced for subtyping and genotyping. By phylogenetic analysis, 96% of the samples clustered with subtype B and the remaining 4% with HIV-1 subtype F sequences. One major protease inhibitor resistance-associated mutation, I50V, was detected in 38% of the samples. Minor mutations were also found at the protease gene: L10I/V (7%), K20M (2%), M36I (11%), L63P (20%), A71T (2%), and V77I (7%). Many mutations associated with reduced susceptibility to nucleoside or non-nucleoside reverse transcriptase inhibitors were detected: M41L (11%), E44D (4%), D67N (11%), T69D (2%), K70R (11%), L74V (2%), L100I (4%), K103N (18%), V118I (9%), Y181C (11%), M184V (18%), G190A (4%), T215Y (4%), and K219E (4%). This study has shown that 84% of the studied population from the Federal District, showing evidences of therapy failure, presented viral genomic mutations associated with drug resistance. The main antiretrovirals to which this population showed resistance were the PI amprenavir (38%), the NNRTIs delavirdine, nevirapine (31%), and efavirenz (24%), and the NRTIs lamivudine (18%), abacavir, and zidovudine (13%).

Genetic-morphometric variation in Culex quinquefasciatus from Brazil and La Plata, Argentina

Variation among natural populations of Culex (Culex) quinquefasciatus Say is associated with different vectorial capacities. The species Cx. quinquefasciatus is present in the equatorial, tropical and subtropical zones in the Brazilian territory, with intermediate forms between Cx. quinquefasciatus and Culex pipiens occurring in regions of latitudes around 33°-35°S. Herein, we studied geographically distinct populations of Cx. quinquefasciatus by genetic characterization and analysis of intra-specific wing morphometrics. After morphological analysis, molecular characterization of Cx. quinquefasciatus and intermediate forms was performed by polymerase chain reaction of the polymorphic nuclear region of the second intron of the acetylcholinesterase locus. Additionally, the morphology of adult female wings collected from six locations was analyzed. Wing centroid sizes were significantly different between some geographical pairs. Mean values of R2/R2+3 differed significantly after pairwise comparisons. The overall wing shape represented by morphometric characters could be divided into two main groupings. Our data suggest that Brazilian samples are morphologically and genetically distinct from the Argentinean samples and also indicated a morphological distinction between northern and southern populations of Brazilian Cx. quinquefasciatus. We suggest that wing morphology may be used for preliminary assessment of population structure of Cx. quinquefasciatusin Brazil.

Genetic and morphological characterisation of a new species of the genus Hysterothylacium (Nematoda) from Paralichthys isosceles Jordan, 1890 (Pisces: Teleostei) of the Neotropical Region, state of Rio de Janeiro, Brazil

Taking into account the difficulties of taxonomic identification of larval anisakid nematodes based on morphological characters, genetic analyses were performed, together with those usually applied, in order to identify anisakid larvae found in the flounder Paralichthys isosceles from the littoral of the state of Rio de Janeiro, Brazil. The analysis of 1,820 larvae revealed a new species, similar to Hysterothylacium MD, Hysterothylacium 2, Hysterothylacium KB and Hysterothylacium sp regarding the absence of the larval tooth, an excretory pore situated below the nerve ring level, and slender lateral alae. Moreover, the new species differs from Hysterothylacium fortalezae and Hysterothylacium reliquens with regard to the number and size of spines present on the tail end and from Hysterothylacium patagonicus by the absence of interlabia. The maximum parsimony and neighbour joining tree topologies based on the 18S ribosomal DNA gene, complete internal transcribed spacer region and cytochrome oxidase 2 (COII) gene demonstrated that the Brazilian larvae belong to Raphidascarididae and represent a unique genetic entity, confirmed as a new Hysterothylacium species. Furthermore, the new species presents COII genetic signatures and shares polymorphisms with Raphidascarididae members. This is the first description of a new anisakid species from Brazil through the integration of morphological and molecular taxonomy data.

Genetic diversity of NS3 protease from Brazilian HCV isolates and possible implications for therapy with direct-acting antiviral drugs

The hepatitis C virus (HCV) NS3 protease has been one of the molecular targets of new therapeutic approaches. Its genomic sequence variability in Brazilian HCV isolates is poorly documented. To obtain more information on the magnitude of its genetic diversity, 114 Brazilian HCV samples were sequenced and analysed together with global reference sequences. Genetic distance (d) analyses revealed that subtype 1b had a higher degree of heterogeneity (d = 0.098) than subtypes 1a (d = 0.060) and 3a (d = 0.062). Brazilian isolates of subtype 1b were distributed in the phylogenetic tree among sequences from other countries, whereas most subtype 1a and 3a sequences clustered into a single branch. Additional characterisation of subtype 1a in clades 1 and 2 revealed that all but two Brazilian subtype 1a sequences formed a distinct and strongly supported (approximate likelihood-ratio test = 93) group of sequences inside clade 1. Moreover, this subcluster inside clade 1 presented an unusual phenotypic characteristic in relation to the presence of resistance mutations for macrocyclic inhibitors. In particular, the mutation Q80K was found in the majority of clade 1 sequences, but not in the Brazilian isolates. These data demonstrate that Brazilian HCV subtypes display a distinct pattern of genetic diversity and reinforce the importance of sequence information in future therapeutic approaches.

Genetic diversity of EBV-encoded LMP1 in the Swiss HIV Cohort Study and implication for NF-Κb activation.

Epstein-Barr virus (EBV) is associated with several types of cancers including Hodgkin's lymphoma (HL) and nasopharyngeal carcinoma (NPC). EBV-encoded latent membrane protein 1 (LMP1), a multifunctional oncoprotein, is a powerful activator of the transcription factor NF-κB, a property that is essential for EBV-transformed lymphoblastoid cell survival. Previous studies reported LMP1 sequence variations and induction of higher NF-κB activation levels compared to the prototype B95-8 LMP1 by some variants. Here we used biopsies of EBV-associated cancers and blood of individuals included in the Swiss HIV Cohort Study (SHCS) to analyze LMP1 genetic diversity and impact of sequence variations on LMP1-mediated NF-κB activation potential. We found that a number of variants mediate higher NF-κB activation levels when compared to B95-8 LMP1 and mapped three single polymorphisms responsible for this phenotype: F106Y, I124V and F144I. F106Y was present in all LMP1 isolated in this study and its effect was variant dependent, suggesting that it was modulated by other polymorphisms. The two polymorphisms I124V and F144I were present in distinct phylogenetic groups and were linked with other specific polymorphisms nearby, I152L and D150A/L151I, respectively. The two sets of polymorphisms, I124V/I152L and F144I/D150A/L151I, which were markers of increased NF-κB activation in vitro, were not associated with EBV-associated HL in the SHCS. Taken together these results highlighted the importance of single polymorphisms for the modulation of LMP1 signaling activity and demonstrated that several groups of LMP1 variants, through distinct mutational paths, mediated enhanced NF-κB activation levels compared to B95-8 LMP1.

Division of labour and worker size polymorphism in ant colonies: the impact of social and genetic factors

Division of labour among workers is central to the organisation and ecological success of insect societies. If there is a genetic component to worker size, morphology or task preference, an increase in colony genetic diversity arising from the presence of multiple breeders per colony might improve division of labour. We studied the genetic basis of worker size and task preference in Formica selysi, an ant species that shows natural variation in the number of mates per queen and the number of queens per colony. Worker size had a heritable component in colonies headed by a doubly mated queen (h(2)=0.26) and differed significantly among matrilines in multiple-queen colonies. However, higher levels of genetic diversity did not result in more polymorphic workers across single- or multiple-queen colonies. In addition, workers from multiple-queen colonies were consistently smaller and less polymorphic than workers from single-queen colonies. The relationship between task, body size and genetic lineage appeared to be complex. Foragers were significantly larger than brood-tenders, which may provide energetic or ergonomic advantages to the colony. Task specialisation was also often associated with genetic lineage. However, genetic lineage and body size were often correlated with task independently of each other, suggesting that the allocation of workers to tasks is modulated by multiple factors. Overall, these results indicate that an increase in colony genetic diversity does not increase worker size polymorphism but might improve colony homeostasis.

Genetic and phenotypic architecture of metabolic syndrome-associated components in dyslipidemic and normolipidemic subjects: the GEMS Study.

Atherogenic dyslipidemia, manifest by low HDL-cholesterol and high TG levels, is an important component of ATP-III defined metabolic syndrome. Here, we dissected the phenotypic and genetic architecture of these traits by assessing their relationships with other metabolically relevant measures, including plasma adipo-cytokines, highly sensitive C-reactive protein (hsCRP) and LDL particle size, in a large family data set (n=2800) and in an independent set of dyslipidemic cases (n=716) and normolipidemic controls (n=1073). We explored the relationships among these phenotypes using variable clustering and then estimated their genetic heritabilities and cross-trait correlations. In families, four clusters explained 61% of the total variance, with one adiposity-related cluster (including hsCRP), one BP-related cluster, and two lipid-related clusters (HDL-C, TG, adiponectin and LDL particle size; apoB and non-HDL-C). A similar structure was observed in dyslipidemic cases and normolipidemic controls. The genetic correlations in the families largely paralleled the phenotype clustering results, suggesting that common genes having pleiotropic effects contributed to the correlations observed. In summary, our analyses support a model of metabolic syndrome with two major components, body fat and lipids, each with two subcomponents, and quantifies their degree of overlap with each other and with metabolic-syndrome related measures (adipokines, LDL particle size and hsCRP).

Genetic and phenotypic population divergence on a microgeographic scale in brown trout.

Salmonid populations of many rivers are rapidly declining. One possible explanation is that habitat fragmentation increases genetic drift and reduces the populations' potential to adapt to changing environmental conditions. We measured the genetic and eco-morphological diversity of brown trout (Salmo trutta) in a Swiss stream system, using multivariate statistics and Bayesian clustering. We found large genetic and phenotypic variation within only 40 km of stream length. Eighty-eight percent of all pairwise F(ST) comparisons and 50% of the population comparisons in body shape were significant. High success rates of population assignment tests confirmed the distinctiveness of populations in both genotype and phenotype. Spatial analysis revealed that divergence increased with waterway distance, the number of weirs, and stretches of poor habitat between sampling locations, but effects of isolation-by-distance and habitat fragmentation could not be fully disentangled. Stocking intensity varied between streams but did not appear to erode genetic diversity within populations. A lack of association between phenotypic and genetic divergence points to a role of local adaptation or phenotypically plastic responses to habitat heterogeneity. Indeed, body shape could be largely explained by topographic stream slope, and variation in overall phenotype matched the flow regimes of the respective habitats.

Linguistic and maternal genetic diversity are not correlated in native mexicans

Mesoamerica, defined as the broad linguistic and cultural area from middle southern Mexico to Costa Rica, might have played a pivotal role during the colonization of theAmerican continent. It has been suggested that the Mesoamerican isthmus could have played an important role in severely restricting prehistorically gene flow between North and SouthAmerica. Although the Native American component has been already described in admixedMexican populations, few studies have been carried out in native Mexican populations. In thisstudy we present mitochondrial DNA (mtDNA) sequence data for the first hypervariable region (HVR-I) in 477 unrelated individuals belonging to eleven different native populations from Mexico. Almost all the Native Mexican mtDNAs could be classified into the four pan-Amerindian haplogroups (A2, B2, C1 and D1); only three of them could be allocated to the rare Native American lineage D4h3. Their haplogroup phylogenies are clearly star-like, as expected from relatively young populations that have experienced diverse episodes of genetic drift (e.g. extensive isolation, genetic drift and founder effects) and posterior population expansions. In agreement with this observation is the fact that Native Mexican populations show a high degree of heterogeneity in their patterns of haplogroup frequencies. HaplogroupX2a was absent in our samples, supporting previous observations where this clade was only detected in the American northernmost areas. The search for identical sequences in the American continent shows that, although Native Mexican populations seem to show a closer relationship to North American populations, they cannot be related to a single geographical region within the continent. Finally, we did not find significant population structure on the maternal lineages when considering the four main and distinct linguistic groups represented in our Mexican samples (Oto-Manguean, Uto-Aztecan, Tarascan, and Mayan), suggesting that genetic divergence predates linguistic diversification in Mexico.

Origin and genetic differentiation of Three Mexican Native Groups (Purépechas, Triquis and Mayas): contribution of CODIS-STRs to the history of human populations of Mesoamerica

BACKGROUND: CODIS-STRs in Native Mexican groups have rarely been analysed for human identification and anthropological purposes. AIM:To analyse the genetic relationships and population structure among three Native Mexican groups from Mesoamerica.SUBJECTS AND METHODS: 531 unrelated Native individuals from Mexico were PCR-typed for 15 and 9 autosomal STRs (Identifiler™ and Profiler™ kits, respectively), including five population samples: Purépechas (Mountain, Valley and Lake), Triquis and Yucatec Mayas. Previously published STR data were included in the analyses. RESULTS:Allele frequencies and statistical parameters of forensic importance were estimated by population. The majority of Native groups were not differentiated pairwise, excepting Triquis and Purépechas, which was attributable to their relative geographic and cultural isolation. Although Mayas, Triquis and Purépechas-Mountain presented the highest number of private alleles, suggesting recurrent gene flow, the elevated differentiation of Triquis indicates a different origin of this gene flow. Interestingly, Huastecos and Mayas were not differentiated, which is in agreement with the archaeological hypothesis that Huastecos represent an ancestral Maya group. Interpopulation variability was greater in Natives than in Mestizos, both significant.CONCLUSION: Although results suggest that European admixture has increased the similarity between Native Mexican groups, the differentiation and inconsistent clustering by language or geography stresses the importance of serial founder effect and/or genetic drift in showing their present genetic relationships.

Peak and persistent excess of genetic diversity following an abrupt migration increase.

Genetic diversity is essential for population survival and adaptation to changing environments. Demographic processes (e.g., bottleneck and expansion) and spatial structure (e.g., migration, number, and size of populations) are known to shape the patterns of the genetic diversity of populations. However, the impact of temporal changes in migration on genetic diversity has seldom been considered, although such events might be the norm. Indeed, during the millions of years of a species' lifetime, repeated isolation and reconnection of populations occur. Geological and climatic events alternately isolate and reconnect habitats. We analytically document the dynamics of genetic diversity after an abrupt change in migration given the mutation rate and the number and sizes of the populations. We demonstrate that during transient dynamics, genetic diversity can reach unexpectedly high values that can be maintained over thousands of generations. We discuss the consequences of such processes for the evolution of species based on standing genetic variation and how they can affect the reconstruction of a population's demographic and evolutionary history from genetic data. Our results also provide guidelines for the use of genetic data for the conservation of natural populations.