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.

A scaled difference chi-square test statistic for moment structure analysis

A family of scaling corrections aimed to improve the chi-square approximation of goodness-of-fit test statistics in small samples, large models, and nonnormal data was proposed in Satorra and Bentler (1994). For structural equations models, Satorra-Bentler's (SB) scaling corrections are available in standard computer software. Often, however, the interest is not on the overall fit of a model, but on a test of the restrictions that a null model say ${\cal M}_0$ implies on a less restricted one ${\cal M}_1$. If $T_0$ and $T_1$ denote the goodness-of-fit test statistics associated to ${\cal M}_0$ and ${\cal M}_1$, respectively, then typically the difference $T_d = T_0 - T_1$ is used as a chi-square test statistic with degrees of freedom equal to the difference on the number of independent parameters estimated under the models ${\cal M}_0$ and ${\cal M}_1$. As in the case of the goodness-of-fit test, it is of interest to scale the statistic $T_d$ in order to improve its chi-square approximation in realistic, i.e., nonasymptotic and nonnormal, applications. In a recent paper, Satorra (1999) shows that the difference between two Satorra-Bentler scaled test statistics for overall model fit does not yield the correct SB scaled difference test statistic. Satorra developed an expression that permits scaling the difference test statistic, but his formula has some practical limitations, since it requires heavy computations that are notavailable in standard computer software. The purpose of the present paper is to provide an easy way to compute the scaled difference chi-square statistic from the scaled goodness-of-fit test statistics of models ${\cal M}_0$ and ${\cal M}_1$. A Monte Carlo study is provided to illustrate the performance of the competing statistics.

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'.

Seventy-five genetic loci influencing the human red blood cell.

Anaemia is a chief determinant of global ill health, contributing to cognitive impairment, growth retardation and impaired physical capacity. To understand further the genetic factors influencing red blood cells, we carried out a genome-wide association study of haemoglobin concentration and related parameters in up to 135,367 individuals. Here we identify 75 independent genetic loci associated with one or more red blood cell phenotypes at P < 10(-8), which together explain 4-9% of the phenotypic variance per trait. Using expression quantitative trait loci and bioinformatic strategies, we identify 121 candidate genes enriched in functions relevant to red blood cell biology. The candidate genes are expressed preferentially in red blood cell precursors, and 43 have haematopoietic phenotypes in Mus musculus or Drosophila melanogaster. Through open-chromatin and coding-variant analyses we identify potential causal genetic variants at 41 loci. Our findings provide extensive new insights into genetic mechanisms and biological pathways controlling red blood cell formation and function.

Genetic evidence for a role of adiponutrin in the metabolism of apolipoprotein B-containing lipoproteins.

Adiponutrin (PNPLA3) is a predominantly liver-expressed transmembrane protein with phospholipase activity that is regulated by fasting and feeding. Recent genome-wide association studies identified PNPLA3 to be associated with hepatic fat content and liver function, thus pointing to a possible involvement in the hepatic lipoprotein metabolism. The aim of this study was to examine the association between two common variants in the adiponutrin gene and parameters of lipoprotein metabolism in 23,274 participants from eight independent West-Eurasian study populations including six population-based studies [Bruneck (n = 800), KORA S3/F3 (n = 1644), KORA S4/F4 (n = 1814), CoLaus (n = 5435), SHIP (n = 4012), Rotterdam (n = 5967)], the SAPHIR Study as a healthy working population (n = 1738) and the Utah Obesity Case-Control Study including a group of 1037 severely obese individuals (average BMI 46 kg/m2) and 827 controls from the same geographical region of Utah. We observed a strong additive association of a common non-synonymous variant within adiponutrin (rs738409) with age-, gender-, and alanine-aminotransferase-adjusted lipoprotein concentrations: each copy of the minor allele decreased levels of total cholesterol on average by 2.43 mg/dl (P = 8.87 x 10(-7)), non-HDL cholesterol levels by 2.35 mg/dl (P = 2.27 x 10(-6)) and LDL cholesterol levels by 1.48 mg/dl (P = 7.99 x 10(-4)). These associations remained significant after correction for multiple testing. We did not observe clear evidence for associations with HDL cholesterol or triglyceride concentrations. In conclusion, our study suggests that adiponutrin is involved in the metabolism of apoB-containing lipoproteins.

Genetic isolation of insular populations of the Maghrebian bat, Myotis punicus, in the Mediterranean Basin

We investigate the population genetic structure of the Maghrebian bat, Myotis punicus, between the mainland and islands to assess the island colonization pattern and current gene flow between nearby islands and within the mainland. Location North Africa and the Mediterranean islands of Corsica and Sardinia. Methods We sequenced part of the control region (HVII) of 79 bats across 11 colonies. The phylogeographical pattern was assessed by analysing molecular diversity indices, examining differentiation among populations and estimating divergence time. In addition, we genotyped 182 bats across 10 colonies at seven microsatellite loci. We used analysis of molecular variance and a Bayesian approach to infer nuclear population structure. Finally, we estimated sex-specific dispersal between Corsica and Sardinia. Results Mitochondrial analyses indicated that colonies between Corsica, Sardinia and North Africa are highly differentiated. Within islands there was no difference between colonies, while at the continental level Moroccan and Tunisian populations were highly differentiated. Analyses with seven microsatellite loci showed a similar pattern. The sole difference was the lack of nuclear differentiation between populations in North Africa, suggesting a male-biased dispersal over the continental area. The divergence time of Sardinian and Corsican populations was estimated to date back to the early and mid-Pleistocene. Main conclusions Island colonization by the Maghrebian bats seems to have occurred in a stepping-stone manner and certainly pre-dated human colonization. Currently, open water seems to prevent exchange of bats between the two islands, despite their ability to fly and the narrowness of the strait of Bonifacio. Corsican and Sardinian populations are thus currently isolated from any continental gene pool and must therefore be considered as different evolutionarily significant units (ESU).

Genetic diversity and relationship in American and African oil palm as revealed by RFLP and AFLP molecular markers

The objective of this work was to evaluate the genetic diversity, its organization and the genetic relationships within oil palm (Elaeis oleifera (Kunth) Cortés, from America, and E. guineensis (Jacq.), from Africa) germplasm using Restriction Fragment Length Polymorphism (RFLP) and Amplified Fragment Length Polymorphism (AFLP). In complement to a previous RFLP study on 241 E. oleifera accessions, 38 E. guineensis accessions were analyzed using the same 37 cDNA probes. These accessions covered a large part of the geographical distribution areas of these species in America and Africa. In addition, AFLP analysis was performed on a sub-set of 40 accessions of E. oleifera and 22 of E. guineensis using three pairs of enzyme/primer combinations. Data were subjected to Factorial Analysis of Correspondence (FAC) and cluster analysis, with parameters of genetic diversity being also studied. Results appeared congruent between RFLP and AFLP. In the E. oleifera, AFLP confirmed the strong structure of genetic diversity revealed by RFLP, according to geographical origin of the studied material, with the identification of the same four distinct genetic groups: Brazil, French Guyana/Surinam, Peru, north of Colombia/Central America. Both markers revealed that genetic divergence between the two species is of the same magnitude as that among provenances of E. oleifera. This finding is in discrepancy with the supposed early tertiary separation of the two species.