Estimating the Parameters of Selection on Nonsynonymous Mutations in Drosophila pseudoobscura and D. miranda

We present the results of surveys of diversity in sets of >40 X-linked and autosomal loci in samples from natural populations of Drosophila miranda and D. pseudoobscura, together with their sequence divergence from D. affinis. Mean silent site diversity in D. miranda is approximately one-quarter of that in D. pseudoobscura; mean X-linked silent diversity is about three-quarters of that for the autosomes in both species. Estimates of the distribution of selection coefficients against heterozygous, deleterious nonsynonymous mutations from two different methods suggest a wide distribution, with coefficients of variation greater than one, and with the average segregating amino acid mutation being subject to only very weak selection. Only a small fraction of new amino acid mutations behave as effectively neutral, however. A large fraction of amino acid differences between D. pseudoobscura and D. affinis appear to have been fixed by positive natural selection, using three different methods of estimation; estimates between D. miranda and D. affinis are more equivocal. Sources of bias in the estimates, especially those arising from selection on synonymous mutations and from the choice of genes, are discussed and corrections for these applied. Overall, the results show that both purifying selection and positive selection on nonsynonymous mutations are pervasive.SURVEYS of DNA sequence diversity and divergence are shedding light on a number of questions in evolutionary genetics (for recent reviews, see Akey 2009; Sella et al. 2009). Two of the most important questions of this kind concern the distribution of selection coefficients against deleterious mutations affecting protein sequences and the proportion of amino acid sequence differences between related species that have been fixed by positive selection. Several different methods have been proposed for studying each of these questions, using different features of data on polymorphism and divergence at nonsynonymous and silent sites.For example, the parameters of the distribution of selection coefficients against deleterious amino acid mutations have been estimated by contrasting the numbers of nonsynonymous and silent within-species polymorphisms and fixed differences between species (Sawyer and Hartl 1992; Bustamante et al. 2002; Piganeau and Eyre-Walker 2003; Sawyer et al. 2007); by fitting the frequency spectra of nonsynonymous and silent variants to models of selection, mutation, and drift (Akashi 1999; Eyre-Walker et al. 2006; Keightley and Eyre-Walker 2007; Kryukov et al. 2007; Boyko et al. 2008; Eyre-Walker and Keightley 2009); or by comparing levels of nonsynonymous and silent diversities between species with different population sizes (Loewe and Charlesworth 2006; Loewe et al. 2006). The results of these different approaches generally agree in suggesting that there is a wide distribution of selection coefficients against nonsynonymous mutations and that the mean selection coefficient against heterozygous carriers of such mutations is very small. The results imply that a typical individual from a human population carries several hundred weakly deleterious mutations (Eyre-Walker et al. 2006; Kryukov et al. 2007; Boyko et al. 2008); for a typical Drosophila population, with its much higher level of variability, the number is probably an order of magnitude greater (Loewe et al. 2006; Keightley and Eyre-Walker 2007).The presence of this large load of slightly deleterious mutations in human and natural populations, most of which are held at low frequencies by natural selection, has many implications. From the point of view of understanding human genetic disease, it means that we have to face the likelihood that susceptibility to a disease can be influenced by variants at many loci, each with small effects (Kryukov et al. 2007). The pervasive presence of deleterious mutations throughout the genome contributes to inbreeding depression (Charlesworth and Willis 2009) and may mean that the effective population size is reduced by background selection effects, even in regions of the genome with normal levels of genetic recombination (Loewe and Charlesworth 2007). Their presence may contribute so strongly to Hill–Robertson effects (Hill and Robertson 1966; Felsenstein 1974) that they cause severely reduced levels of diversity and adaptation in low-recombination regions of the genome (Charlesworth et al. 2010) and create a selective advantage to maintaining nonzero levels of recombination (Keightley and Otto 2006; Charlesworth et al. 2010). In addition, having an estimate of the distribution of selection coefficients against deleterious nonsynonymous mutations allows their contribution to between-species divergence to be predicted, providing a way of estimating the fraction of fixed nonsynonymous differences caused by positive selection (Loewe et al. 2006; Boyko et al. 2008; Eyre-Walker and Keightley 2009).It is thus important to collect data that shed light on the properties of selection against nonsynonymous mutations in a wide range of systems and also to compare the results from different methods of estimation, since they are subject to different sources of difficulty and biases. In a previous study, we proposed the use of a comparison between two related species with different effective population sizes for this purpose (Loewe and Charlesworth 2006; Loewe et al. 2006), using Drosophila miranda and D. pseudoobscura as material. These are well suited for this type of study, as they are closely related, live together in similar habitats, and yet have very different levels of silent nucleotide diversity, indicating different effective population sizes (N_e). This study was hampered by our inability to compare the same set of loci across the two species and by the small number of loci that could be used. We here present the results of a much larger study of DNA variation at X-linked and autosomal loci for these two species, using D. affinis as a basis for estimating divergence. We compare the results, applying the method of Loewe et al. (2006) with that of Eyre-Walker and Keightley (2009) for estimating the distribution of deleterious selection coefficients and with McDonald–Kreitman test-based methods for estimating the proportion of nonsynonymous differences fixed by positive selection. While broadly confirming the conclusions from earlier studies, we note some possible sources of bias and describe methods for minimizing their effects.     

QTLs and candidate genes for desiccation and abscisic acid content in maize kernels

Background

Kernel moisture at harvest is an important trait since a low value is required to prevent unexpected early germination and ensure seed preservation. It is also well known that early germination occurs in viviparous mutants, which are impaired in abscisic acid (ABA) biosynthesis. To provide some insight into the genetic determinism of kernel desiccation in maize, quantitative trait loci (QTLs) were detected for traits related to kernel moisture and ABA content in both embryo and endosperm during kernel desiccation. In parallel, the expression and mapping of genes involved in kernel desiccation and ABA biosynthesis, were examined to detect candidate genes.

Results

The use of an intermated recombinant inbred line population allowed for precise QTL mapping. For 29 traits examined in an unreplicated time course trial of days after pollination, a total of 78 QTLs were detected, 43 being related to kernel desiccation, 15 to kernel weight and 20 to ABA content. Multi QTL models explained 35 to 50% of the phenotypic variation for traits related to water status, indicating a large genetic control amenable to breeding. Ten of the 20 loci controlling ABA content colocated with previously detected QTLs controlling water status and ABA content in water stressed leaves. Mapping of candidate genes associated with kernel desiccation and ABA biosynthesis revealed several colocations between genes with putative functions and QTLs. Parallel investigation via RT-PCR experiments showed that the expression patterns of the ABA-responsive Rab17 and Rab28 genes as well as the late embryogenesis abundant Emb5 and aquaporin genes were related to desiccation rate and parental allele effect. Database searches led to the identification and mapping of two zeaxanthin epoxidase (ZEP) and five novel 9-cis-epoxycarotenoid dioxygenase (NCED) related genes, both gene families being involved in ABA biosynthesis. The expression of these genes appeared independent in the embryo and endosperm and not correlated with ABA content in either tissue.

Conclusions

A high resolution QTL map for kernel desiccation and ABA content in embryo and endosperm showed several precise colocations between desiccation and ABA traits. Five new members of the maize NCED gene family and another maize ZEP gene were identified and mapped. Among all the identified candidates, aquaporins and members of the Responsive to ABA gene family appeared better candidates than NCEDs and ZEPs.     

MeCP2 Dependent Heterochromatin Reorganization during Neural Differentiation of a Novel Mecp2-Deficient Embryonic Stem Cell Reporter Line

The X-linked Mecp2 is a known interpreter of epigenetic information and mutated in Rett syndrome, a complex neurological disease. MeCP2 recruits HDAC complexes to chromatin thereby modulating gene expression and, importantly regulates higher order heterochromatin structure. To address the effects of MeCP2 deficiency on heterochromatin organization during neural differentiation, we developed a versatile model for stem cell in vitro differentiation. Therefore, we modified murine Mecp2 deficient (Mecp2 ^−/y) embryonic stem cells to generate cells exhibiting green fluorescent protein expression upon neural differentiation. Subsequently, we quantitatively analyzed heterochromatin organization during neural differentiation in wild type and in Mecp2 deficient cells. We found that MeCP2 protein levels increase significantly during neural differentiation and accumulate at constitutive heterochromatin. Statistical analysis of Mecp2 wild type neurons revealed a significant clustering of heterochromatin per nuclei with progressing differentiation. In contrast we found Mecp2 deficient neurons and astroglia cells to be significantly impaired in heterochromatin reorganization. Our results (i) introduce a new and manageable cellular model to study the molecular effects of Mecp2 deficiency, and (ii) support the view of MeCP2 as a central protein in heterochromatin architecture in maturating cells, possibly involved in stabilizing their differentiated state.     

Influence of the APOA5 locus on plasma triglyceride, lipoprotein subclasses, and CVD risk in the Framingham Heart Study

Several polymorphisms in the APOA5 gene have been associated with increased plasma triglyceride (TG) concentrations. However, associations between APOA5 and lipoprotein subclasses, remnant-like particles (RLPs), and cardiovascular disease (CVD) risk have been less explored. We investigated associations of five APOA5 single-nucleotide polymorphisms (SNPs; -1131T>C, -3A>G, 56C>G IVS3+ 476G>A, and 1259T>C) with lipoprotein subfractions and CVD risk in 1,129 men and 1,262 women participating in the Framingham Heart Study. Except for the 56C>G SNP, the other SNPs were in significant linkage disequilibria, resulting in three haplotypes (11111, 22122, and 11211) representing 98% of the population. SNP analyses revealed that the -1131T>C and 56C>G SNPs were significantly associated with higher plasma TG concentrations in both men and women. For RLP and lipoprotein subclasses, we observed gender-specific association for the -1131T>C and 56C>G SNPs. Female carriers of the -1131C allele had higher RLP concentrations, whereas in males, significant associations for RLPs were observed for the 56G allele. Moreover, haplotype analyses confirmed these findings and revealed that the 22122 and 11211 haplotypes exhibited different associations with HDL cholesterol concentrations. In women, the -1131C allele was associated with a higher hazard ratio for CVD (1.85; 95% confidence interval, 1.03-3.34; P = 0.04), in agreement with the association of this SNP with higher RLPs.     

Anthrax toxin: the long and winding road that leads to the kill

The past five years have led to a tremendous increase in our molecular understanding of the mode of action of the anthrax toxin, one of the two main virulence factors produced by Bacillus anthracis. The structures of each of the three components of the toxin--lethal factor (LF), edema factor (EF) and protective antigen (PA)--have been solved not only in their monomeric forms but, depending on the subunit, in a heptameric form, bound to their substrate, co-factor or receptor. The endocytic route followed by the toxin has also been unraveled and the enzymatic mechanisms of EF and LF elucidated.     

Atomic Layer Deposition of Functional Layers for on Chip 3D Li‐Ion All Solid State Microbattery

Nowadays, millimeter scale power sources are key devices for providing autonomy to smart, connected, and miniaturized sensors. However, until now, planar solid state microbatteries do not yet exhibit a sufficient surface energy density. In that context, architectured 3D microbatteries appear therefore to be a good solution to improve the material mass loading while keeping small the footprint area. Beside the design itself of the 3D microbaterry, one important technological barrier to address is the conformal deposition of thin films (lithiated or not) on 3D structures. For that purpose, atomic layer deposition (ALD) technology is a powerful technique that enables conformal coatings of thin film on complex substrate. An original, robust, and highly efficient 3D scaffold is proposed to significantly improve the geometrical surface of miniaturized 3D microbattery. Four functional layers composing the 3D lithium ion microbattery stacking has been successfully deposited on simple and double microtubes 3D templates. In depth synchrotron X‐ray nanotomography and high angle annular dark field transmission electron microscope analyses are used to study the interface between each layer. For the first time, using ALD, anatase TiO₂ negative electrode is coated on 3D tubes with Li₃PO₄ lithium phosphate as electrolyte, opening the way to all solid‐state 3D microbatteries. The surface capacity is significantly increased by the proposed topology (high area enlargement factor – “thick” 3D layer), from 3.5 μA h cm^?2 for a planar layer up to 0.37 mA h cm^?2 for a 3D thin film (105 times higher).