Genome-wide prediction of discrete traits using bayesian regressions and machine learning

Background

Genomic selection has gained much attention and the main goal is to increase the predictive accuracy and the genetic gain in livestock using dense marker information. Most methods dealing with the large p (number of covariates) small n (number of observations) problem have dealt only with continuous traits, but there are many important traits in livestock that are recorded in a discrete fashion (e.g. pregnancy outcome, disease resistance). It is necessary to evaluate alternatives to analyze discrete traits in a genome-wide prediction context.

Methods

This study shows two threshold versions of Bayesian regressions (Bayes A and Bayesian LASSO) and two machine learning algorithms (boosting and random forest) to analyze discrete traits in a genome-wide prediction context. These methods were evaluated using simulated and field data to predict yet-to-be observed records. Performances were compared based on the models'' predictive ability.

Results

The simulation showed that machine learning had some advantages over Bayesian regressions when a small number of QTL regulated the trait under pure additivity. However, differences were small and disappeared with a large number of QTL. Bayesian threshold LASSO and boosting achieved the highest accuracies, whereas Random Forest presented the highest classification performance. Random Forest was the most consistent method in detecting resistant and susceptible animals, phi correlation was up to 81% greater than Bayesian regressions. Random Forest outperformed other methods in correctly classifying resistant and susceptible animals in the two pure swine lines evaluated. Boosting and Bayes A were more accurate with crossbred data.

Conclusions

The results of this study suggest that the best method for genome-wide prediction may depend on the genetic basis of the population analyzed. All methods were less accurate at correctly classifying intermediate animals than extreme animals. Among the different alternatives proposed to analyze discrete traits, machine-learning showed some advantages over Bayesian regressions. Boosting with a pseudo Huber loss function showed high accuracy, whereas Random Forest produced more consistent results and an interesting predictive ability. Nonetheless, the best method may be case-dependent and a initial evaluation of different methods is recommended to deal with a particular problem.     

Utility of multilocus genotypes for taxon assignment in stands of closely related European white oaks from Switzerland

Background and Aims

European white oaks (Quercus petraea, Q. pubescens, Q. robur) have long puzzled plant biologists owing to disputed species differentiation. Extensive hybridization or shared ancestry have been proposed as alternative hypotheses to explain why genetic differentiation between these oak species is low. Species delimitation is usually weak and often shows gradual transitions in leaf morphology. Hence, individual identification may be difficult, but remains a critical step for both scientific work and practical management.

Methods

Multilocus genotype data (five nuclear microsatellites) were used from ten Swiss oak stands for taxon identification without a priori grouping of individuals or populations, using model-based Bayesian assignment tests.

Key Results

Three groups best structured the data, indicating that the taxonomical signal was stronger than the spatial signal. Most individuals showed high posterior probabilities for either of three genetic groups that were best circumscribed as taxonomical units. The assignment of a subset of trees, whose taxonomic status had been previously characterized in detail, supported this classification scheme.

Conclusions

Molecular-genetic assignment tests are useful in the identification of species status in critical taxon complexes such as the European white oaks. Such an approach is of practical importance for forest management, e.g. for stand certification or in seed trade to trace the origin of forest products.Key words: Assignment test, Bayesian inference, multilocus genotype, nuclear microsatellites, Quercus sp., species complex     

Phylogeography sheds light on the central–marginal hypothesis in a Mediterranean narrow endemic plant

Background and Aims

Understanding the factors that shape variation in genetic diversity across the geographic ranges of species is an important challenge in the effort to conserve evolutionary processes sustaining biodiversity. The historical influences leading to a central–marginal organization of genetic diversity have been explored for species whose range is known to have expanded from refugia after glacial events. However, this question has rarely been addressed for Mediterranean endemic plants of azonal habitats such as rocky slopes or screes. In this context, this comprehensive study examined molecular and field data from Arenaria provincialis (Caryophyllaceae), a narrow endemic plant of south-eastern France.

Methods

Across the whole geographic range, an investigation was made of whether high levels of abundance and genetic diversity (estimated from amplified fragment length polymorphism markers) are centrally distributed, to evaluate the relevance of the central–marginal hypothesis. Phylogeographic patterns inferred from chloroplast DNA (cpDNA) were used, applying Bayesian methods to test the influence of past biogeographic events. Multivariate analysis combining phylogeographic and ecological data was used to reveal the historical and ecological distinctiveness of populations.

Key Results

Despite the narrow distribution of A. provincialis, a high level of nucleotide variation is found within cpDNA loci, supporting its persistence throughout the Pleistocene period. The area characterized by the highest genetic diversity is centrally located. Structured phylogeography and Bayesian factor analysis supported the hypothesis that the central area of the distribution was the source of both westward and eastward migrations, probably during arid periods of the Pleistocene, and more recently was a crossroads of backward migrations. By contrast, the two areas located today at the range limits are younger, have reduced genetic diversity and are marginal in the ecological gradients.

Conclusions

This study highlights a case of strong population distinctiveness within a narrow range. Phylogeography sheds light on the historical role of the areas centrally situated in the distribution. The current range size and abundance patterns are not sufficient to predict the organization of genetic diversity.     

Nuclear DNA microsatellites reveal genetic variation but a lack of phylogeographical structure in an endangered species, Fraxinus mandshurica, across North-east China

Background and Aims

The widely accepted paradigm that the modern genetic structure of plant species in the northern hemisphere has been largely determined by recolonization from refugia after the last glacial maximum fails to explain the presence of cold-tolerant species at intermediate latitudes. Another generally accepted paradigm is that mountain ridges act as important barriers causing genetic isolation of species, but this too has been challenged in recent studies. The aims of the work reported here were to determine the genetic diversity and distribution patterns of extant natural populations of an endangered cool temperate species, Faxinus mandshurica, and to examine whether these two paradigms are appropriate when applied to this species over a wide geographical scale.

Methods

1435 adult individuals were sampled from 30 natural populations across the main and central range of the species, covering major mountain ranges across North-east China (NEC). Genetic variation was estimated based on nine polymorphic nuclear microsatellite loci. Phylogeographical analyses were employed using various approaches, including Bayesian clustering, spatial analysis of molecular variance, Monmonier''s algorithm, neighbor-joining trees, principal co-ordinate analysis and isolation by distance.

Key Results

Genetic diversity within populations was relatively high, and no significant recent bottlenecks were detected in any of the populations. A significant negative correlation between intra-population genetic diversity and latitude was identified. In contrast, genetic differentiation among all the populations examined was extremely low and no clear geographic genetic structure was identified, with the exception of one distinct population.

Conclusions

The modern genetic structure in this species can be explained by extensive gene flow, an absence of mountains acting as barriers, and the presence of a wide refuge across NEC rather than multiple small refugia. Intra-population genetic variation along latitudes is probably associated with the systematically northward shifts of forest biomes in eastern China during the mid-Holocene. To determine important genetic patterns and identify resources for conservation, however, it will be necessary to examine differentially inherited genetic markers exposed to selection pressures (e.g. chloroplast DNA) and to investigate different generations.Key words: Fraxinus mandshurica, nuclear microsatellites, latitude variation, historical migration, fossil pollen, spatial genetic structure, genetic barriers     

Bayesian neural networks for detecting epistasis in genetic association studies

Background

Discovering causal genetic variants from large genetic association studies poses many difficult challenges. Assessing which genetic markers are involved in determining trait status is a computationally demanding task, especially in the presence of gene-gene interactions.

Results

A non-parametric Bayesian approach in the form of a Bayesian neural network is proposed for use in analyzing genetic association studies. Demonstrations on synthetic and real data reveal they are able to efficiently and accurately determine which variants are involved in determining case-control status. By using graphics processing units (GPUs) the time needed to build these models is decreased by several orders of magnitude. In comparison with commonly used approaches for detecting interactions, Bayesian neural networks perform very well across a broad spectrum of possible genetic relationships.

Conclusions

The proposed framework is shown to be a powerful method for detecting causal SNPs while being computationally efficient enough to handle large datasets.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0368-0) contains supplementary material, which is available to authorized users.     

Early PTSD Symptom Trajectories: Persistence,Recovery, and Response to Treatment: Results from the Jerusalem Trauma Outreach and Prevention Study (J-TOPS)

Context

Uncovering heterogeneities in the progression of early PTSD symptoms can improve our understanding of the disorder''s pathogenesis and prophylaxis.

Objectives

To describe discrete symptom trajectories and examine their relevance for preventive interventions.

Design

Latent Growth Mixture Modeling (LGMM) of data from a randomized controlled study of early treatment. LGMM identifies latent longitudinal trajectories by exploring discrete mixture distributions underlying observable data.

Setting

Hadassah Hospital unselectively receives trauma survivors from Jerusalem and vicinity.

Participants

Adult survivors of potentially traumatic events consecutively admitted to the hospital''s emergency department (ED) were assessed ten days and one-, five-, nine- and fifteen months after ED admission. Participants with data at ten days and at least two additional assessments (n = 957) were included; 125 received cognitive behavioral therapy (CBT) between one and nine months.

Approach

We used LGMM to identify latent parameters of symptom progression and tested the effect of CBT on these parameters. CBT consisted of 12 weekly sessions of either cognitive therapy (n = 41) or prolonged exposure (PE, n = 49), starting 29.8±5.7 days after ED admission, or delayed PE (n = 35) starting at 151.8±42.4 days. CBT effectively reduced PTSD symptoms in the entire sample.

Main Outcome Measure

Latent trajectories of PTSD symptoms; effects of CBT on these trajectories.

Results

Three trajectories were identified: Rapid Remitting (rapid decrease in symptoms from 1- to 5-months; 56% of the sample), Slow Remitting (progressive decrease in symptoms over 15 months; 27%) and Non-Remitting (persistently elevated symptoms; 17%). CBT accelerated the recovery of the Slow Remitting class but did not affect the other classes.

Conclusions

The early course of PTSD symptoms is characterized by distinct and diverging response patterns that are centrally relevant to understanding the disorder and preventing its occurrence. Studies of the pathogenesis of PTSD may benefit from using clustered symptom trajectories as their dependent variables.     

Association between Age at Diagnosis of Graves' Disease and Variants in Genes Involved in Immune Response

Background

Graves'' disease (GD) is a complex disease in which genetic predisposition is modified by environmental factors. The aim of the study was to examine the association between genetic variants in genes encoding proteins involved in immune response and the age at diagnosis of GD.

Methods

735 GD patients and 1216 healthy controls from Poland were included into the study. Eight genetic variants in the HLA-DRB1, TNF, CTLA4, CD40, NFKb, PTPN22, IL4 and IL10 genes were genotyped. Patients were stratified by the age at diagnosis of GD and the association with genotype was analysed.

Results

Polymorphism in the HLA-DRB1, TNF and CTLA4 genes were associated with GD. The carriers of the HLA DRB1*03 allele were more frequent in patients with age at GD diagnosis ≤30 years than in patients with older age at GD diagnosis.

Conclusions

HLADRB1*03 allele is associated with young age at diagnosis of Graves'' disease in polish population.     

Ecological Niche Models and Coalescent Analysis of Gene Flow Support Recent Allopatric Isolation of Parasitoid Wasp Populations in the Mediterranean

Background

The integration of multiple complementary approaches is a powerful way to understand the processes of diversification and speciation. The parasitoid wasp Aphidius transcaspicus Telenga (Hymenoptera: Braconidae) is a parasitoid of Hyalopterus aphids across a wide geographic range. This species shows a remarkable degree of genetic structure among western, central, and eastern Mediterranean population clusters. In this paper we attempt to better characterize this genetic structure.

Methodology/Principal Findings

We use a Bayesian coalescent analysis of gene flow under the Isolation with Migration model using mitochondrial and microsatellite markers together with climate-based ecological niche models to better understand the genetic structure of A. transcaspicus in the Mediterranean. The coalescent analysis revealed low levels of migration among western and eastern Mediterranean populations (Nm<1) that were not statistically distinguishable from zero. Niche models showed that localities within population clusters each occupy areas of continuously high environmental suitability, but are separated from each other by large regions of completely unsuitable habitat that could limit dispersal. Overall, environmental characteristics were similar among the population clusters, though significant differences did emerge.

Conclusions/Significance

These results support contemporary allopatric isolation of Mediterranean populations of A. transcaspicus, which together with previous analyses indicating partial behaviorally mediated reproductive isolation, suggest that the early stages of cryptic speciation may be in progress.     

Genetic consequences of anagenetic speciation in Acer okamotoanum (Sapindaceae) on Ullung Island, Korea

Background and Aim

Anagenesis (also known as phyletic speciation) is an important process of speciation in endemic species of oceanic islands. We investigated genetic variation in Acer okamotoanum, an anagenetically derived species endemic to Ullung Island, South Korea, to infer genetic consequences of anagenesis in comparison with other groups that have undergone cladogenesis (and adaptive radiation).

Methods

We examined genetic variation based on eight polymorphic microsatellite markers from 145 individuals of A. okamotoanum and 134 individuals of its putative progenitor A. mono. We employed standard population genetic analyses, clustering analyses, Bayesian clustering analyses in STRUCTURE and bottleneck analyses.

Key Results

Based on both the Neighbor–Joining tree and Bayesian clustering analyses, clear genetic distinctions were found between the two species. Genetic diversity in terms of allelic richness and heterozygosity shows slightly lower levels in A. okamotoanum in comparison with A. mono. Bayesian clustering analyses showed a relatively high F-value in the cluster of A. okamotoanum, suggesting a strong episode of genetic drift during colonization and speciation. There was no clear evidence of a bottleneck based on allelic frequency distribution and excess of observed heterozygotes, but the M-ratio indicated a historical bottleneck in several populations of A. okamotoanum. No geographical genetic structure within the island was found, and the genetic variation among populations of A. okamotoanum was quite low.

Conclusions

We hypothesized that genetic consequences of oceanic-endemic plants derived via anagenesis would be quite different from those derived via cladogenesis. Populations of A. okamotoanum form a cluster and are clearly differentiated from A. mono, which suggests a single origin for the anagenetically derived island endemic. No pattern of geographical differentiation of populations occurs in A. okamotoanum, which supports the concept of initial founder populations diverging through time by accumulation of mutations in a relatively uniform environment without further specific differentiation.     

Present,past and future of the European rock fern Asplenium fontanum: combining distribution modelling and population genetics to study the effect of climate change on geographic range and genetic diversity

Background and Aims

Climate change is expected to alter the geographic range of many plant species dramatically. Predicting this response will be critical to managing the conservation of plant resources and the effects of invasive species. The aim of this study was to predict the response of temperate homosporous ferns to climate change.

Methods

Genetic diversity and changes in distribution range were inferred for the diploid rock fern Asplenium fontanum along a South–North transect, extending from its putative last glacial maximum (LGM) refugia in southern France towards southern Germany and eastern-central France. This study reconciles observations from distribution models and phylogeographic analyses derived from plastid and nuclear diversity.

Key Results

Genetic diversity distribution and niche modelling propose that genetic diversity accumulates in the LGM climate refugium in southern France with the formation of a diversity gradient reflecting a slow, post-LGM range expansion towards the current distribution range. Evidence supports the fern''s preference for outcrossing, contradicting the expectation that homosporous ferns would populate new sites by single-spore colonization. Prediction of climate and distribution range change suggests that a dramatic loss of range and genetic diversity in this fern is possible. The observed migration is best described by the phalanx expansion model.

Conclusions

The results suggest that homosporous ferns reproducing preferentially by outcrossing accumulate genetic diversity primarily in LGM climate refugia and may be threatened if these areas disappear due to global climate change.     

General Intelligence in Another Primate: Individual Differences across Cognitive Task Performance in a New World Monkey (Saguinus oedipus)

Background

Individual differences in human cognitive abilities show consistently positive correlations across diverse domains, providing the basis for the trait of “general intelligence” (g). At present, little is known about the evolution of g, in part because most comparative studies focus on rodents or on differences across higher-level taxa. What is needed, therefore, are experiments targeting nonhuman primates, focusing on individual differences within a single species, using a broad battery of tasks. To this end, we administered a large battery of tasks, representing a broad range of cognitive domains, to a population of captive cotton-top tamarin monkeys (Saguinus oedipus).

Methodology and Results

Using a Bayesian latent variable model, we show that the pattern of correlations among tasks is consistent with the existence of a general factor accounting for a small but significant proportion of the variance in each task (the lower bounds of 95% Bayesian credibility intervals for correlations between g and task performance all exceed 0.12).

Conclusion

Individual differences in cognitive abilities within at least one other primate species can be characterized by a general intelligence factor, supporting the hypothesis that important aspects of human cognitive function most likely evolved from ancient neural substrates.     

Genetic structure of the tree peony (Paeonia rockii) and the Qinling Mountains as a geographic barrier driving the fragmentation of a large population

Background

Tree peonies are great ornamental plants associated with a rich ethnobotanical history in Chinese culture and have recently been used as an evolutionary model. The Qinling Mountains represent a significant geographic barrier in Asia, dividing mainland China into northern (temperate) and southern (semi–tropical) regions; however, their flora has not been well analyzed. In this study, the genetic differentiation and genetic structure of Paeonia rockii and the role of the Qinling Mountains as a barrier that has driven intraspecific fragmentation were evaluated using 14 microsatellite markers.

Methodology/Principal Findings

Twenty wild populations were sampled from the distributional range of P. rockii. Significant population differentiation was suggested (F_ST value of 0.302). Moderate genetic diversity at the population level (H_S of 0.516) and high population diversity at the species level (H_T of 0.749) were detected. Significant excess homozygosity (F_IS of 0.076) and recent population bottlenecks were detected in three populations. Bayesian clusters, population genetic trees and principal coordinate analysis all classified the P. rockii populations into three genetic groups and one admixed Wenxian population. An isolation-by-distance model for P. rockii was suggested by Mantel tests (r = 0.6074, P<0.001) and supported by AMOVA (P<0.001), revealing a significant molecular variance among the groups (11.32%) and their populations (21.22%). These data support the five geographic boundaries surrounding the Qinling Mountains and adjacent areas that were detected with Monmonier''s maximum-difference algorithm.

Conclusions/Significance

Our data suggest that the current genetic structure of P. rockii has resulted from the fragmentation of a formerly continuously distributed large population following the restriction of gene flow between populations of this species by the Qinling Mountains. This study provides a fundamental genetic profile for the conservation and responsible exploitation of the extant germplasm of this species and for improving the genetic basis for breeding its cultivars.     

Isolated populations of a rare alpine plant show high genetic diversity and considerable population differentiation

Background and Aims

Gene flow and genetic variability within and among alpine plant populations can be greatly influenced by the steep environmental gradients and heterogeneous topography of alpine landscapes. In this study, the effects are examined of natural isolation of alpine habitats on genetic diversity and geographic structure in populations of C. thyrsoides, a rare and isolated European Alpine monocarpic perennial with limited seed dispersal capacity.

Methods

Molecular diversity was analysed for 736 individuals from 32 populations in the Swiss Alps and adjacent Jura mountains using five polymorphic microsatellite loci. Pollen flow was estimated using pollen grain-sized fluorescent powder. In addition, individual-based Bayesian approaches were applied to examine population structure.

Key Results

High within-population genetic diversity (H_E = 0·76) and a relatively low inbreeding coefficient (F_IS = 0·022) were found. Genetic differentiation among populations measured with a standardized measure was considerable (G′_ST = 0·53). A significant isolation-by-distance relationship was found (r = 0·62, P < 0·001) and a significant geographic sub-structure, coinciding with proposed postglacial migration patterns. Altitudinal location and size of populations did not influence molecular variation. Direct measures of pollen flow revealed that insect-mediated pollen dispersal was restricted to short distances within a population.

Conclusions

The natural isolation of suitable habitats for C. thyrsoides restricts gene flow among the populations as expected for a monocarpic species with very limited seed dispersal capacities. The observed high within-population genetic diversity in this rare monocarpic perennial is best explained by its outcrossing behaviour, long-lived individuals and overlapping generations. Despite the high within-population genetic diversity, the considerable genetic differentiation and the clear western–eastern differentiation in this species merits consideration in future conservation efforts.Key words: Alpine plant, Campanula thyrsoides, genetic diversity, gene flow, genetic differentiation, glacial history, G′_ST, habitat isolation, microsatellites, monocarpy, SSR     

Relevance of genetics for conservation policies: the case of Minorcan cork oaks

Background and Aims

Marginal populations of widely distributed species can be of high conservation interest when they hold a significant or unique portion of the genetic diversity of the species. However, such genetic information is frequently lacking. Here the relevance of genetic surveys to develop efficient conservation strategies for such populations is illustrated using cork oak (Quercus suber) from Minorca (Balearic Islands, Spain) as a case study. Cork oak is highly endangered on the island, where no more than 67 individuals live in small, isolated stands in siliceous sites. As a consequence, it was recently granted protected status.

Methods

Two Bayesian clustering approaches were used to analyse the genetic structure of the Minorcan population, on the basis of nuclear microsatellite data. The different groups within the island were also compared with additional island and continental populations surrounding Minorca.

Key Results

Very high genetic diversity was found, with values comparable with those observed in continental parts of the species'' range. Furthermore, the Minorcan oak stands were highly differentiated from one another and were genetically related to different continental populations of France and Spain.

Conclusions

The high levels of genetic diversity and inter-stands differentiation make Minorcan cork oak eligible for specific conservation efforts. The relationship of Minorcan stands to different continental populations in France and Spain probably reflects multiple colonization events. However, discrepancy between chloroplast DNA- and nuclear DNA-based groups does not support a simple scenario of recent introduction. Gene exchanges between neighbouring cork oak stands and with holm oak have created specific and exceptional genetic combinations. They also constitute a wide range of potential genetic resources for research on adaptation to new environmental conditions. Conservation guidelines that take into account these findings are provided.     

A multicenter study confirms CD226 gene association with systemic sclerosis-related pulmonary fibrosis

Introduction

CD226 genetic variants have been associated with a number of autoimmune diseases and recently with systemic sclerosis (SSc). The aim of this study was to test the influence of CD226 loci in SSc susceptibility, clinical phenotypes and autoantibody status in a large multicenter European population.

Methods

A total of seven European populations of Caucasian ancestry were included, comprising 2,131 patients with SSc and 3,966 healthy controls. Three CD226 single nucleotide polymorphisms (SNPs), rs763361, rs3479968 and rs727088, were genotyped using Taqman 5''allelic discrimination assays.

Results

Pooled analyses showed no evidence of association of the three SNPs, neither with the global disease nor with the analyzed subphenotypes. However, haplotype block analysis revealed a significant association for the TCG haplotype (SNP order: rs763361, rs34794968, rs727088) with lung fibrosis positive patients (P_Bonf = 3.18E-02 OR 1.27 (1.05 to 1.54)).

Conclusion

Our data suggest that the tested genetic variants do not individually influence SSc susceptibility but a CD226 three-variant haplotype is related with genetic predisposition to SSc-related pulmonary fibrosis.     

Characterization of Vascular Disease Risk in Postmenopausal Women and Its Association with Cognitive Performance

Objectives

While global measures of cardiovascular (CV) risk are used to guide prevention and treatment decisions, these estimates fail to account for the considerable interindividual variability in pre-clinical risk status. This study investigated heterogeneity in CV risk factor profiles and its association with demographic, genetic, and cognitive variables.

Methods

A latent profile analysis was applied to data from 727 recently postmenopausal women enrolled in the Kronos Early Estrogen Prevention Study (KEEPS). Women were cognitively healthy, within three years of their last menstrual period, and free of current or past CV disease. Education level, apolipoprotein E ε4 allele (APOE4), ethnicity, and age were modeled as predictors of latent class membership. The association between class membership, characterizing CV risk profiles, and performance on five cognitive factors was examined. A supervised random forest algorithm with a 10-fold cross-validation estimator was used to test accuracy of CV risk classification.

Results

The best-fitting model generated two distinct phenotypic classes of CV risk 62% of women were “low-risk” and 38% “high-risk”. Women classified as low-risk outperformed high-risk women on language and mental flexibility tasks (p = 0.008) and a global measure of cognition (p = 0.029). Women with a college degree or above were more likely to be in the low-risk class (OR = 1.595, p = 0.044). Older age and a Hispanic ethnicity increased the probability of being at high-risk (OR = 1.140, p = 0.002; OR = 2.622, p = 0.012; respectively). The prevalence rate of APOE-ε4 was higher in the high-risk class compared with rates in the low-risk class.

Conclusion

Among recently menopausal women, significant heterogeneity in CV risk is associated with education level, age, ethnicity, and genetic indicators. The model-based latent classes were also associated with cognitive function. These differences may point to phenotypes for CV disease risk. Evaluating the evolution of phenotypes could in turn clarify preclinical disease, and screening and preventive strategies.ClinicalTrials.gov {"type":"clinical-trial","attrs":{"text":"NCT00154180","term_id":"NCT00154180"}}NCT00154180     

Language Structure Is Partly Determined by Social Structure

Background

Languages differ greatly both in their syntactic and morphological systems and in the social environments in which they exist. We challenge the view that language grammars are unrelated to social environments in which they are learned and used.

Methodology/Principal Findings

We conducted a statistical analysis of >2,000 languages using a combination of demographic sources and the World Atlas of Language Structures— a database of structural language properties. We found strong relationships between linguistic factors related to morphological complexity, and demographic/socio-historical factors such as the number of language users, geographic spread, and degree of language contact. The analyses suggest that languages spoken by large groups have simpler inflectional morphology than languages spoken by smaller groups as measured on a variety of factors such as case systems and complexity of conjugations. Additionally, languages spoken by large groups are much more likely to use lexical strategies in place of inflectional morphology to encode evidentiality, negation, aspect, and possession. Our findings indicate that just as biological organisms are shaped by ecological niches, language structures appear to adapt to the environment (niche) in which they are being learned and used. As adults learn a language, features that are difficult for them to acquire, are less likely to be passed on to subsequent learners. Languages used for communication in large groups that include adult learners appear to have been subjected to such selection. Conversely, the morphological complexity common to languages used in small groups increases redundancy which may facilitate language learning by infants.

Conclusions/Significance

We hypothesize that language structures are subjected to different evolutionary pressures in different social environments. Just as biological organisms are shaped by ecological niches, language structures appear to adapt to the environment (niche) in which they are being learned and used. The proposed Linguistic Niche Hypothesis has implications for answering the broad question of why languages differ in the way they do and makes empirical predictions regarding language acquisition capacities of children versus adults.     

Investigating the Role of Mitochondrial Haplogroups in Genetic Predisposition to Meningococcal Disease

Background and Aims

Meningococcal disease remains one of the most important infectious causes of death in industrialized countries. The highly diverse clinical presentation and prognosis of Neisseria meningitidis infections are the result of complex host genetics and environmental interactions. We investigated whether mitochondrial genetic background contributes to meningococcal disease (MD) susceptibility.

Methodology/Principal Findings

Prospective controlled study was performed through a national research network on MD that includes 41 Spanish hospitals. Cases were 307 paediatric patients with confirmed MD, representing the largest series of MD patients analysed to date. Two independent sets of ethnicity-matched control samples (CG1 [N = 917]), and CG2 [N = 616]) were used for comparison. Cases and controls underwent mtDNA haplotyping of a selected set of 25 mtDNA SNPs (mtSNPs), some of them defining major European branches of the mtDNA phylogeny. In addition, 34 ancestry informative markers (AIMs) were genotyped in cases and CG2 in order to monitor potential hidden population stratification. Samples of known African, Native American and European ancestry (N = 711) were used as classification sets for the determination of ancestral membership of our MD patients. A total of 39 individuals were eliminated from the main statistical analyses (including fourteen gypsies) on the basis of either non-Spanish self-reported ancestry or the results of AIMs indicating a European membership lower than 95%. Association analysis of the remaining 268 cases against CG1 suggested an overrepresentation of the synonym mtSNP G11719A variant (Pearson''s chi-square test; adjusted P-value = 0.0188; OR [95% CI] = 1.63 [1.22–2.18]). When cases were compared with CG2, the positive association could not be replicated. No positive association has been observed between haplogroup (hg) status of cases and CG1/CG2 and hg status of cases and several clinical variants.

Conclusions

We did not find evidence of association between mtSNPs and mtDNA hgs with MD after carefully monitoring the confounding effect of population sub-structure. MtDNA variability is particularly stratified in human populations owing to its low effective population size in comparison with autosomal markers and therefore, special care should be taken in the interpretation of seeming signals of positive associations in mtDNA case-control association studies.     

Evolutionary diversification and geographical isolation in Dubautia laxa (Asteraceae), a widespread member of the Hawaiian silversword alliance

Background and Aims

The Hawaiian silversword alliance (Asteraceae) is one the best examples of a plant adaptive radiation, exhibiting extensive morphological and ecological diversity. No research within this group has addressed the role of geographical isolation, independent of ecological adaptation, in contributing to taxonomic diversity. The aims of this study were to examine genetic differentiation among subspecies of Dubautia laxa (Asteraceae) to determine if allopatric or sympatric populations and subspecies form distinct genetic clusters to understand better the role of geography in diversification within the alliance.

Methods

Dubautia laxa is a widespread member of the Hawaiian silversword alliance, occurring on four of the five major islands of the Hawaiian archipelago, with four subspecies recognized on the basis of morphological, ecological and geographical variation. Nuclear microsatellites and plastid DNA sequence data were examined. Data were analysed using maximum-likelihood and Bayesian phylogenetic methodologies to identify unique evolutionary lineages.

Key Results

Plastid DNA sequence data resolved two highly divergent lineages, recognized as the Laxa and Hirsuta groups, that are more similar to other members of the Hawaiian silversword alliance than they are to each other. The Laxa group is basal to the young island species of Dubautia, whereas the Hirsuta group forms a clade with the old island lineages of Dubautia and with Argyroxiphium. The divergence between the plastid groups is supported by Bayesian microsatellite clustering analyses, but the degree of nuclear differentiation is not as great. Clear genetic differentiation is only observed between allopatric populations, both within and among islands.

Conclusions

These results indicate that geographical separation has aided diversification in D. laxa, whereas ecologically associated morphological differences are not associated with neutral genetic differentiation. This suggests that, despite the stunning ecological adaptation observed, geography has also played an important role in the Hawaiian silversword alliance plant adaptive radiation.     

Updated three-stage model for the peopling of the Americas

Background

We re-assess support for our three stage model for the peopling of the Americas in light of a recent report that identified nine non-Native American mitochondrial genome sequences that should not have been included in our initial analysis. Removal of these sequences results in the elimination of an early (i.e. ∼40,000 years ago) expansion signal we had proposed for the proto-Amerind population.

Methodology/Findings

Bayesian skyline plot analysis of a new dataset of Native American mitochondrial coding genomes confirms the absence of an early expansion signal for the proto-Amerind population and allows us to reduce the variation around our estimate of the New World founder population size. In addition, genetic variants that define New World founder haplogroups are used to estimate the amount of time required between divergence of proto-Amerinds from the Asian gene pool and expansion into the New World.

Conclusions/Significance

The period of population isolation required for the generation of New World mitochondrial founder haplogroup-defining genetic variants makes the existence of three stages of colonization a logical conclusion. Thus, our three stage model remains an important and useful working hypothesis for researchers interested in the peopling of the Americas and the processes of colonization.