On the transition of genetic differentiation from isolation to panmixia: What we can learn from GST and D.

Population genetic differentiation characterizes the repartition of alleles among populations. It is commonly thought that genetic differentiation measures, such as GST and D, should be near zero when allele frequencies are close to their expected value in panmictic populations, and close to one when they are close to their expected value in isolated populations. To analyse those properties, we first derive analytically a reference function f of known parameters that describes how important features of genetic differentiation (e.g. gene diversity, proportion of private alleles, frequency of the most common allele) are close to their expected panmictic and isolation value. We find that the behaviour of function f differs according to three distinct mutation regimes defined by the scaled mutation rate and the number of populations. Then, we compare GST and D to f, and demonstrate that their signal of differentiation strongly depends on the mutation regime. In particular, we show that D captures well the variations of genetic diversity when mutation is weak, otherwise it overestimates it when panmixia is not met. GST detects population differentiation when mutation is intermediate but has a low sensitivity to the variations of genetic diversity when mutation is weak. When mutation is strong the domain of sensitivity of both measures are altered. Finally, we also point out the importance of the number of populations on genetic differentiation measures, and provide recommendations for the use of GST and D.

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.

Population genetic structure and demography analysis in palinurus elephas populations from atlantic and mediterranean waters

Report for the scientific sojourn at the University of Reading, United Kingdom, from January until May 2008. The main objectives have been firstly to infer population structure and parameters in demographic models using a total of 13 microsatellite loci for genotyping approximately 30 individuals per population in 10 Palinurus elephas populations both from Mediterranean and Atlantic waters. Secondly, developing statistical methods to identify discrepant loci, possibly under selection and implement those methods using the R software environment. It is important to consider that the calculation of the probability distribution of the demographic and mutational parameters for a full genetic data set is numerically difficult for complex demographic history (Stephens 2003). The Approximate Bayesian Computation (ABC), based on summary statistics to infer posterior distributions of variable parameters without explicit likelihood calculations, can surmount this difficulty. This would allow to gather information on different demographic prior values (i.e. effective population sizes, migration rate, microsatellite mutation rate, mutational processes) and assay the sensitivity of inferences to demographic priors by assuming different priors.

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.

Trypanosoma cruzi: genetic group with peculiar biochemical and biological behavior

Thirty-two Trypanosoma cruzi strains, isolated from chronic chagasic patients in the northwest of the state of Paraná (Brazil), were analyzed using molecular, biochemical and biological characteristics. Genotypic analysis using randomly amplified polymorphic DNA and simple sequence repeat-anchored polymerase chain reaction amplified profiles showed a large, genetically well-correlated group that contained the majority of the strains and a divergent group that included the PR-150 strain. For glycoconjugate composition, the PR-150 strain was different from the other strains considering the absence or presence of specific bands in aqueous or detergent phases. This strain was also totally different from the others in one out of the six parameters related to in vitro and in vivo biological behavior. We highlight the fact that the PR-150 was totally resistant to benznidazole. For the other biological parameters this strain was not totally distinct from the others, but it showed a peculiar behavior.

Quantitative genetic modeling and inference in the presence of nonignorable missing data.

Natural selection is typically exerted at some specific life stages. If natural selection takes place before a trait can be measured, using conventional models can cause wrong inference about population parameters. When the missing data process relates to the trait of interest, a valid inference requires explicit modeling of the missing process. We propose a joint modeling approach, a shared parameter model, to account for nonrandom missing data. It consists of an animal model for the phenotypic data and a logistic model for the missing process, linked by the additive genetic effects. A Bayesian approach is taken and inference is made using integrated nested Laplace approximations. From a simulation study we find that wrongly assuming that missing data are missing at random can result in severely biased estimates of additive genetic variance. Using real data from a wild population of Swiss barn owls Tyto alba, our model indicates that the missing individuals would display large black spots; and we conclude that genes affecting this trait are already under selection before it is expressed. Our model is a tool to correctly estimate the magnitude of both natural selection and additive genetic variance.

Genetic variability in the 5' UTR and NS5A regions of hepatitis C virus RNA isolated from non-responding and responding patients with chronic HCV genotype 1 infection

Sequence variation among different hepatitis C virus (HCV) isolates has adaptive significance and reflects the modes and intensities of selection mechanisms operating on the virus. In this work, we sought to investigate using classical population genetics parameters, the genetic variability of HCV genotype 1 using the 5' UTR and NS5A regions from treatment non-responding and responding groups of patients. Both regions showed low genetic varia-bility and the 5' UTR showed neutral deviation. No differences were observed in the nonsynonymous/synonymous nucleotide substitution ratio among groups for NS5A. The analysis of molecular variance test of the 5' UTR region showed an 11.94% variation among groups. Phylogenetic analysis showed no correlation between sequence variations and therapeutic responses.

Association analysis of urotensin II gene (UTS2) and flanking regions with biochemical parameters related to insulin resistance.

Background. Urotensin II (UII) is a potent vasoconstrictor peptide, which signals through a G-protein coupled receptor (GPCR) known as GPR14 or urotensin receptor (UTR). UII exerts a broad spectrum of actions in several systems such as vascular cell, heart muscle or pancreas, where it inhibits insulin release. Objective. Given the reported role of UII in insulin secretion, we have performed a genetic association analysis of the UTS2 gene and flanking regions with biochemical parameters related to insulin resistance (fasting glucose, glucose 2 hours after a glucose overload, fasting insulin and insulin resistance estimated as HOMA). Results and Conclusions. We have identified several polymorphisms associated with the analysed clinical traits, not only at the UTS2 gene, but also in thePER3 gene, located upstream from UTS2. Our results are compatible with a role for UII in glucose homeostasis and diabetes although we cannot rule out the possibility that PER3 gene may underlie the reported associations.

Influence of anti-filarial chemotherapy strategies on the genetic structure of Wuchereria bancrofti populations

Lymphatic filarial (LF) parasites have been under anti-filarial drug pressure for more than half a century. Currently, annual mass drug administration (MDA) of diethylcarbamazine (DEC) or ivermectin in combination with albendazole (ALB) have been used globally to eliminate LF. Long-term chemotherapies exert significant pressure on the genetic structure of parasitic populations. We investigated the genetic variation among 210 Wuchereria bancrofti populations that were under three different chemotherapy strategies, namely MDA with DEC alone (group I, n = 74), MDA with DEC and ALB (group II, n = 60) and selective therapy (ST) with DEC (group III, n = 34) to understand the impact of these three drug regimens on the parasite genetic structure. Randomly amplified polymorphic DNA profiles were generated for the three groups of parasite populations; the gene diversity, gene flow and genetic distance values were determined and phylogenetic trees were constructed. Analysis of these parameters indicated that parasite populations under ST with a standard dose of DEC (group III) were genetically more diverse (0.2660) than parasite populations under MDA with DEC alone (group I, H = 0.2197) or with DEC + ALB (group II, H = 0.2317). These results indicate that the MDA may reduce the genetic diversity of W. bancrofti populations when compared to the genetic diversity of parasite populations under ST.

Condition-dependent genetic benefits of extrapair fertilization in female blue tits Cyanistes caeruleus.

In many socially monogamous birds, both partners perform extrapair copulations (EPC). As this behaviour potentially inflicts direct costs on females, they are currently hypothesized to search for genetic benefits for descendants, either as 'good' or 'complementary' genes. Although these hypotheses have found some support, several studies failed to find any beneficial consequence of EPC, and whether this behaviour is adaptive to females is subject to discussion. Here, we test these two hypotheses in a natural population of blue tits by accounting for the effect of most parameters known to potentially affect extrapair fertilization. Results suggest that female body mass affected the type of extrapair genetic benefits obtained. Heavy females obtained extrapair fertilizations when their social male was of low quality (as reflected by sexual display) and produced larger extrapair than within-pair chicks. Lean females obtained extrapair fertilizations when their social mate was genetically similar, thereby producing more heterozygous extrapair chicks. Our results suggest that mating patterns may be condition-dependent.

Sex-biased dispersal in a migratory bat: a characterization using sex-specific demographic parameters.

We studied the noctule bat (Nyctalus noctula), in which the mitochondrial F(ST) is about 10 times that revealed by nuclear markers, to address two questions. We first verified whether random dispersal of one sex is compatible with highly contrasted mitochondrial and nuclear population structures. Using computer simulations, we then assessed the power of multilocus population differentiation tests when the expected population structure departs only slightly from panmixia. Using an island model with sex-specific demographic parameters, we found that random male dispersal is consistent with the population structure observed in the noctule. However, other parameter combinations are also compatible with the data. We computed the minimum sex bias in dispersal (at least 69% of the dispersing individuals are males), a result that would not be available if we had used more classical population genetic models. The power of multilocus population differentiation tests was unexpectedly high, the tests being significant in almost 100% of the replicates, although the observed population structure infered from nuclear markers was extremely low (F(ST) = 0.6%).

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.

Impaired glutathione synthesis in schizophrenia: convergent genetic and functional evidence

Schizophrenia is a complex multifactorial brain disorder with a genetic component. Convergent evidence has implicated oxidative stress and glutathione (GSH) deficits in the pathogenesis of this disease. The aim of the present study was to test whether schizophrenia is associated with a deficit of GSH synthesis. Cultured skin fibroblasts from schizophrenia patients and control subjects were challenged with oxidative stress, and parameters of the rate-limiting enzyme for the GSH synthesis, the glutamate cysteine ligase (GCL), were measured. Stressed cells of patients had a 26% (P = 0.002) decreased GCL activity as compared with controls. This reduction correlated with a 29% (P < 0.001) decreased protein expression of the catalytic GCL subunit (GCLC). Genetic analysis of a trinucleotide repeat (TNR) polymorphism in the GCLC gene showed a significant association with schizophrenia in two independent case-control studies. The most common TNR genotype 7/7 was more frequent in controls [odds ratio (OR) = 0.6, P = 0.003], whereas the rarest TNR genotype 8/8 was three times more frequent in patients (OR = 3.0, P = 0.007). Moreover, subjects with disease-associated genotypes had lower GCLC protein expression (P = 0.017), GCL activity (P = 0.037), and GSH contents (P = 0.004) than subjects with genotypes that were more frequent in controls. Taken together, the study provides genetic and functional evidence that an impaired capacity to synthesize GSH under conditions of oxidative stress is a vulnerability factor for schizophrenia.

Genetic association tests: a method for the joint analysis of family and case-control data.

With the trend in molecular epidemiology towards both genome-wide association studies and complex modelling, the need for large sample sizes to detect small effects and to allow for the estimation of many parameters within a model continues to increase. Unfortunately, most methods of association analysis have been restricted to either a family-based or a case-control design, resulting in the lack of synthesis of data from multiple studies. Transmission disequilibrium-type methods for detecting linkage disequilibrium from family data were developed as an effective way of preventing the detection of association due to population stratification. Because these methods condition on parental genotype, however, they have precluded the joint analysis of family and case-control data, although methods for case-control data may not protect against population stratification and do not allow for familial correlations. We present here an extension of a family-based association analysis method for continuous traits that will simultaneously test for, and if necessary control for, population stratification. We further extend this method to analyse binary traits (and therefore family and case-control data together) and accurately to estimate genetic effects in the population, even when using an ascertained family sample. Finally, we present the power of this binary extension for both family-only and joint family and case-control data, and demonstrate the accuracy of the association parameter and variance components in an ascertained family sample.

Kinetic parameters of silicon uptake by rice cultivars

Silicon is considered an important chemical element for rice, because it can improve tolerance to biotic and abiotic stress. However, in many situations no positive effect of silicon was observed, probably due to genetic factors. The objective of this research was to monitor Si uptake kinetics and identify responses of rice cultivars in terms of Si uptake capacity and use. The experiment was carried out in a greenhouse of the São Paulo State University (UNESP), Brazil. The experiment was arranged in a completely randomized, factorial design with three replications. that consisted of two rice cultivars and two Si levels. Kinetic parameters (Vmax, Km, and Cmin), root morphology variables, dry matter yield, Si accumulation and levels in shoots and roots, uptake efficiency, utilization efficiency, and root/shoot ratio were evaluated. Higher Si concentrations in the nutrient solution did not increase rice dry matter. The development of the low-affinity silicon uptake system of the rice cultivar 'Caiapó' was better than of 'Maravilha'.