Fine time scaling of purifying selection on human nonsynonymous mtDNA mutations based on the worldwide population tree and mother-child pairs

A high-resolution mtDNA phylogenetic tree allowed us to look backward in time to investigate purifying selection. Purifying selection was very strong in the last 2,500 years, continuously eliminating pathogenic mutations back until the end of the Younger Dryas (∼11,000 years ago), when a large population expansion likely relaxed selection pressure. This was preceded by a phase of stable selection until another relaxation occurred in the out-of-Africa migration. Demography and selection are closely related: expansions led to relaxation of selection and higher pathogenicity mutations significantly decreased the growth of descendants. The only detectible positive selection was the recurrence of highly pathogenic nonsynonymous mutations (m.3394T>C-m.3397A>G-m.3398T>C) at interior branches of the tree, preventing the formation of a dinucleotide STR (TATATA) in the MT-ND1 gene. At the most recent time scale in 124 mother-children transmissions, purifying selection was detectable through the loss of mtDNA variants with high predicted pathogenicity. A few haplogroup-defining sites were also heteroplasmic, agreeing with a significant propensity in 349 positions in the phylogenetic tree to revert back to the ancestral variant. This nonrandom mutation property explains the observation of heteroplasmic mutations at some haplogroup-defining sites in sequencing datasets, which may not indicate poor quality as has been claimed.

Mitochondrial targeting adaptation of the hominoid-specific glutamate dehydrogenase driven by positive Darwinian selection.

Many new gene copies emerged by gene duplication in hominoids, but little is known with respect to their functional evolution. Glutamate dehydrogenase (GLUD) is an enzyme central to the glutamate and energy metabolism of the cell. In addition to the single, GLUD-encoding gene present in all mammals (GLUD1), humans and apes acquired a second GLUD gene (GLUD2) through retroduplication of GLUD1, which codes for an enzyme with unique, potentially brain-adapted properties. Here we show that whereas the GLUD1 parental protein localizes to mitochondria and the cytoplasm, GLUD2 is specifically targeted to mitochondria. Using evolutionary analysis and resurrected ancestral protein variants, we demonstrate that the enhanced mitochondrial targeting specificity of GLUD2 is due to a single positively selected glutamic acid-to-lysine substitution, which was fixed in the N-terminal mitochondrial targeting sequence (MTS) of GLUD2 soon after the duplication event in the hominoid ancestor approximately 18-25 million years ago. This MTS substitution arose in parallel with two crucial adaptive amino acid changes in the enzyme and likely contributed to the functional adaptation of GLUD2 to the glutamate metabolism of the hominoid brain and other tissues. We suggest that rapid, selectively driven subcellular adaptation, as exemplified by GLUD2, represents a common route underlying the emergence of new gene functions.

The Effect of Inoculum Size on Selection of In Vitro Resistance to Vancomycin, Daptomycin, and Linezolid in Methicillin-Resistant Staphylococcus aureus.

OBJECTIVES: The inoculum effect (IE) is an increase in the minimum inhibitory concentration (MIC) at high bacterial densities. The effect of three inoculum sizes on the selection of resistance to vancomycin, daptomycin, and linezolid was investigated in methicillin-resistant Staphylococcus aureus (MRSA). METHODS: Low (10(4) CFU/ml), medium (10(6) CFU/ml), and high (10(8) CFU/ml) inocula of MRSA were exposed to twofold increasing concentrations of either drug during 15 days of cycling. MICs for low (MICL), medium (MICM), and high (MICH) inocula were determined daily. Conventional MICs were measured at days 1, 5, 10, and 15. Experiments were performed in triplicate. RESULTS: At the beginning of the experiment a small IE was observed for vancomycin (MICL=1 μg/ml, MICM=1-2 μg/ml, and MICH=2 μg/ml) and a significant IE for daptomycin (MICL=0.25 μg/ml, MICM=0.25-0.5 μg/ml, and MICH=2 μg/ml). Linezolid exhibited no IE at low and medium inocula (MICL=1 μg/ml and MICM=1-2 μg/ml), but with the high inoculum, concentrations up to 2,048 μg/ml did not fully inhibit visual growth. During cycling, increase of MIC was observed for all antibiotics. At day 15, MICL, MICM, and MICH of vancomycin were 2-4, 4-8, and 4-16 μg/ml and of daptomycin were 0.5-2, 8-128, and 64-256 μg/ml, respectively. MICL and MICM of linezolid were 1 and 2-4 μg/ml, respectively. Conventional MICs showed vancomycin and daptomycin selection of resistance since day 5 depending on the inocula. No selection of linezolid resistance was observed. CONCLUSIONS: Our results showed the importance of the inoculum size in the development of resistance. Measures aimed at lowering the inoculum at the site of infection should be used whenever possible in parallel to antimicrobial therapy.

Parallel evolution of facial stripe patterns in the Neolamprologus brichardi/pulcher species complex endemic to Lake Tanganyika.

Colour pattern diversity can be due to random processes or to natural or sexual selection. Consequently, similarities in colour patterns are not always correlated with common ancestry, but may result from convergent evolution under shared selection pressures or drift. Neolamprologus brichardi and Neolamprologus pulcher have been described as two distinct species based on differences in the arrangement of two dark bars on the operculum. Our study uses DNA sequences of the mitochondrial control region to show that relatedness of haplotypes disagrees with species assignment based on head colour pattern. This suggests repeated parallel evolution of particular stripe patterns. The complete lack of shared haplotypes between populations of the same or different phenotypes reflects strong philopatric behaviour, possibly induced by the cooperative breeding mode in which offspring remain in their natal territory and serve as helpers until they disperse to nearby territories or take over a breeding position. Concordant phylogeographic patterns between N. brichardi/N. pulcher populations and other rock-dwelling cichlids suggest that the same colonization routes have been taken by sympatric species and that these routes were affected by lake level fluctuations in the past.

Neutral and selection-driven decay of sexual traits in asexual stick insects.

Environmental shifts and lifestyle changes may result in formerly adaptive traits becoming non-functional or maladaptive. The subsequent decay of such traits highlights the importance of natural selection for adaptations, yet its causes have rarely been investigated. To study the fate of formerly adaptive traits after lifestyle changes, we evaluated sexual traits in five independently derived asexual lineages, including traits that are specific to males and therefore not exposed to selection. At least four of the asexual lineages retained the capacity to produce males that display normal courtship behaviours and are able to fertilize eggs of females from related sexual species. The maintenance of male traits may stem from pleiotropy, or from these traits only regressing via drift, which may require millions of years to generate phenotypic effects. By contrast, we found parallel decay of sexual traits in females. Asexual females produced altered airborne and contact signals, had modified sperm storage organs, and lost the ability to fertilize their eggs, impeding reversals to sexual reproduction. Female sexual traits were decayed even in recently derived asexuals, suggesting that trait changes following the evolution of asexuality, when they occur, proceed rapidly and are driven by selective processes rather than drift.

C4 Photosynthesis evolved in grasses via parallel adaptive genetic changes.

Phenotypic convergence is a widespread and well-recognized evolutionary phenomenon. However, the responsible molecular mechanisms remain often unknown mainly because the genes involved are not identified. A well-known example of physiological convergence is the C4 photosynthetic pathway, which evolved independently more than 45 times [1]. Here, we address the question of the molecular bases of the C4 convergent phenotypes in grasses (Poaceae) by reconstructing the evolutionary history of genes encoding a C4 key enzyme, the phosphoenolpyruvate carboxylase (PEPC). PEPC genes belong to a multigene family encoding distinct isoforms of which only one is involved in C4 photosynthesis [2]. By using phylogenetic analyses, we showed that grass C4 PEPCs appeared at least eight times independently from the same non-C4 PEPC. Twenty-one amino acids evolved under positive selection and converged to similar or identical amino acids in most of the grass C4 PEPC lineages. This is the first record of such a high level of molecular convergent evolution, illustrating the repeatability of evolution. These amino acids were responsible for a strong phylogenetic bias grouping all C4 PEPCs together. The C4-specific amino acids detected must be essential for C4 PEPC enzymatic characteristics, and their identification opens new avenues for the engineering of the C4 pathway in crops.

Review. Characterization and selection of hexaploid wheats containing resistance to Heterodera avenae or Mayetiola destructor introgressed from Aegilops

Cereal cyst nematode (CCN, Heterodera avenae) and Hessian fly (HF, Mayetiola destructor) are two major pests affecting wheat crops worldwide including important cereal areas of Spain. Aegilops ventricosa and Ae. triuncialis were used as donors in a strategy to introduce resistance genes (RG) for these pests in hexaploid wheat (Triticum aestivum L.). Two 42 chromosomes introgression lines have been derived from Ae. ventricosa: H-93-8 and H-93-33 carrying genes Cre2 and H27 conferring resistance to CCN and HF, respectively. Line TR-3531 with 42 chromosomes has been derived from Ae. triuncialis and carries RGs conferring resistance for CCN (Cre7) and for HF (H30). Alien material has been incorporated in lines H-93 by chromosomal substitution and recombination, while in line TR-3531 homoeologous recombination affecting small DNA fragments has played a major role. It has been demonstrated that Cre2, Cre7, H27 and H30 are major single dominant genes and not allelic of other previously described RGs. Biochemical and molecular-biology studies of the defense mechanism triggered by Cre2 and Cre7 have revealed specific induction of peroxidase and other antioxidant enzymes. In parallel to these basic studies advanced lines carrying resistance genes for CNN and/or HF have been developed. Selection was done using molecular markers for eventually «pyramiding» resistance genes. Several isozyme and RAPD markers have been described and, currently, new markers based on transposable elements and NBS-LRR sequences are being developed. At present, two advanced lines have already been included at the Spanish Catalogue of Commercial Plant Varieties.

Spatially varying selection shapes life history clines among populations of Drosophila melanogaster from sub-Saharan Africa.

Clines in life history traits, presumably driven by spatially varying selection, are widespread. Major latitudinal clines have been observed, for example, in Drosophila melanogaster, an ancestrally tropical insect from Africa that has colonized temperate habitats on multiple continents. Yet, how geographic factors other than latitude, such as altitude or longitude, affect life history in this species remains poorly understood. Moreover, most previous work has been performed on derived European, American and Australian populations, but whether life history also varies predictably with geography in the ancestral Afro-tropical range has not been investigated systematically. Here, we have examined life history variation among populations of D. melanogaster from sub-Saharan Africa. Viability and reproductive diapause did not vary with geography, but body size increased with altitude, latitude and longitude. Early fecundity covaried positively with altitude and latitude, whereas lifespan showed the opposite trend. Examination of genetic variance-covariance matrices revealed geographic differentiation also in trade-off structure, and QST -FST analysis showed that life history differentiation among populations is likely shaped by selection. Together, our results suggest that geographic and/or climatic factors drive adaptive phenotypic differentiation among ancestral African populations and confirm the widely held notion that latitude and altitude represent parallel gradients.

Genomics of clinal variation in Drosophila: disentangling the interactions of selection and demography.

Clines in phenotypes and genotype frequencies across environmental gradients are commonly taken as evidence for spatially varying selection. Classical examples include the latitudinal clines in various species of Drosophila, which often occur in parallel fashion on multiple continents. Today, genomewide analysis of such clinal systems provides a fantastic opportunity for unravelling the genetics of adaptation, yet major challenges remain. A well-known but often neglected problem is that demographic processes can also generate clinality, independent of or coincident with selection. A closely related issue is how to identify true genic targets of clinal selection. In this issue of Molecular Ecology, three studies illustrate these challenges and how they might be met. Bergland et al. report evidence suggesting that the well-known parallel latitudinal clines in North American and Australian D. melanogaster are confounded by admixture from Africa and Europe, highlighting the importance of distinguishing demographic from adaptive clines. In a companion study, Machado et al. provide the first genomic comparison of latitudinal differentiation in D. melanogaster and its sister species D. simulans. While D. simulans is less clinal than D. melanogaster, a significant fraction of clinal genes is shared between both species, suggesting the existence of convergent adaptation to clinaly varying selection pressures. Finally, by drawing on several independent sources of evidence, Bo?ičević et al. identify a functional network of eight clinal genes that are likely involved in cold adaptation. Together, these studies remind us that clinality does not necessarily imply selection and that separating adaptive signal from demographic noise requires great effort and care.

Efficient Peer Selection in P2P JXTA-Based Platforms

Peer-reviewed

Scalable Feature Selection Applications for Genome-Wide Association Studies of Complex Diseases

Personalized medicine will revolutionize our capabilities to combat disease. Working toward this goal, a fundamental task is the deciphering of geneticvariants that are predictive of complex diseases. Modern studies, in the formof genome-wide association studies (GWAS) have afforded researchers with the opportunity to reveal new genotype-phenotype relationships through the extensive scanning of genetic variants. These studies typically contain over half a million genetic features for thousands of individuals. Examining this with methods other than univariate statistics is a challenging task requiring advanced algorithms that are scalable to the genome-wide level. In the future, next-generation sequencing studies (NGS) will contain an even larger number of common and rare variants. Machine learning-based feature selection algorithms have been shown to have the ability to effectively create predictive models for various genotype-phenotype relationships. This work explores the problem of selecting genetic variant subsets that are the most predictive of complex disease phenotypes through various feature selection methodologies, including filter, wrapper and embedded algorithms. The examined machine learning algorithms were demonstrated to not only be effective at predicting the disease phenotypes, but also doing so efficiently through the use of computational shortcuts. While much of the work was able to be run on high-end desktops, some work was further extended so that it could be implemented on parallel computers helping to assure that they will also scale to the NGS data sets. Further, these studies analyzed the relationships between various feature selection methods and demonstrated the need for careful testing when selecting an algorithm. It was shown that there is no universally optimal algorithm for variant selection in GWAS, but rather methodologies need to be selected based on the desired outcome, such as the number of features to be included in the prediction model. It was also demonstrated that without proper model validation, for example using nested cross-validation, the models can result in overly-optimistic prediction accuracies and decreased generalization ability. It is through the implementation and application of machine learning methods that one can extract predictive genotype–phenotype relationships and biological insights from genetic data sets.

Neurofuzzy mixture of experts network parallel learning and model construction algorithms

A connection between a fuzzy neural network model with the mixture of experts network (MEN) modelling approach is established. Based on this linkage, two new neuro-fuzzy MEN construction algorithms are proposed to overcome the curse of dimensionality that is inherent in the majority of associative memory networks and/or other rule based systems. The first construction algorithm employs a function selection manager module in an MEN system. The second construction algorithm is based on a new parallel learning algorithm in which each model rule is trained independently, for which the parameter convergence property of the new learning method is established. As with the first approach, an expert selection criterion is utilised in this algorithm. These two construction methods are equivalent in their effectiveness in overcoming the curse of dimensionality by reducing the dimensionality of the regression vector, but the latter has the additional computational advantage of parallel processing. The proposed algorithms are analysed for effectiveness followed by numerical examples to illustrate their efficacy for some difficult data based modelling problems.

Optimization Of Groundwater Remediation With New Efficient Parallel Algorithm For Global Optimization

Application of optimization algorithm to PDE modeling groundwater remediation can greatly reduce remediation cost. However, groundwater remediation analysis requires a computational expensive simulation, therefore, effective parallel optimization could potentially greatly reduce computational expense. The optimization algorithm used in this research is Parallel Stochastic radial basis function. This is designed for global optimization of computationally expensive functions with multiple local optima and it does not require derivatives. In each iteration of the algorithm, an RBF is updated based on all the evaluated points in order to approximate expensive function. Then the new RBF surface is used to generate the next set of points, which will be distributed to multiple processors for evaluation. The criteria of selection of next function evaluation points are estimated function value and distance from all the points known. Algorithms created for serial computing are not necessarily efficient in parallel so Parallel Stochastic RBF is different algorithm from its serial ancestor. The application for two Groundwater Superfund Remediation sites, Umatilla Chemical Depot, and Former Blaine Naval Ammunition Depot. In the study, the formulation adopted treats pumping rates as decision variables in order to remove plume of contaminated groundwater. Groundwater flow and contamination transport is simulated with MODFLOW-MT3DMS. For both problems, computation takes a large amount of CPU time, especially for Blaine problem, which requires nearly fifty minutes for a simulation for a single set of decision variables. Thus, efficient algorithm and powerful computing resource are essential in both cases. The results are discussed in terms of parallel computing metrics i.e. speedup and efficiency. We find that with use of up to 24 parallel processors, the results of the parallel Stochastic RBF algorithm are excellent with speed up efficiencies close to or exceeding 100%.

Genetic associations among average annual productivity, growth traits, and stayability: A parallel between Nelore and composite beef cattle

This study was conducted to examine the relationship among average annual productivity of the cow (PRODAM), yearling weight (YW), postweaning BW gain (PWG), scrotal circumference (SC), and stayability in the herd for at least 6 yr (STAY) of Nelore and composite beef cattle. Measurements were taken on animals born between 1980 and 2010 on 70 farms located in 7 Brazilian states. Estimates of heritability and genetic and environmental correlations were obtained by Bayesian approach with 5-trait animal models. Genetic trends were estimated by regressing means of estimated breeding values by year of birth. The heritability estimates were between 0.14 and 0.47. Estimates of genetic correlation among female traits (PRODAM and STAY) and growth traits ranged from-0.02 to 0.30. Estimates of genetic correlations ranged from 0.23 to 0.94 among growth traits indicating that selection for these traits could be successful in tropical breeding programs. Genetic correlations among all traits were favorable and simultaneous selection for growth, productivity, and stayability is therefore possible. Genetic correlation between PRODAM and STAY was 0.99 and 0.85 for Nelore and composite cattle, respectively. Therefore, PRODAM and STAY might be influenced by many of the same genes. The inclusion of PRODAM instead of STAY as a selection criterion seems to be more advantageous for tropical breeding programs because the generation interval required to obtain accurate estimates of genetic merit for PRODAM is shorter. Average annual genetic changes were greater in Nelore than in composite cattle. This was not unexpected because the breeding program of composite cattle included a large number of farms, different production environments, and genetic level of the herds and breeds. Thus, the selection process has become more difficult in this population. © 2013 American Society of Animal Science. All rights reserved.