The complete linkage disequilibrium test: a test that points to causative mutations underlying quantitative traits

<p>Abstract</p> <p>Background</p> <p>Genetically, SNP that are in complete linkage disequilibrium with the causative SNP cannot be distinguished from the causative SNP. The Complete Linkage Disequilibrium (CLD) test presented here tests whether a SNP is in complete LD with the causative mutation or not. The performance of the CLD test is evaluated in 1000 simulated datasets.</p> <p>Methods</p> <p>The CLD test consists of two steps i.e. analysis I and analysis II. Analysis I consists of an association analysis of the investigated region. The log-likelihood values from analysis I are next ranked in descending order and in analysis II the CLD test evaluates differences in log-likelihood ratios between the best and second best markers. Under the null-hypothesis distribution, the best SNP is in greater LD with the QTL than the second best, while under the alternative-CLD-hypothesis, the best SNP is alike-in-state with the QTL. To find a significance threshold, the test was also performed on data excluding the causative SNP. The 5^th, 10^thand 50^thhighest T_CLDvalue from 1000 replicated analyses were used to control the type-I-error rate of the test at p = 0.005, p = 0.01 and p = 0.05, respectively.</p> <p>Results</p> <p>In a situation where the QTL explained 48% of the phenotypic variance analysis I detected a QTL in 994 replicates (p = 0.001), where 972 were positioned in the correct QTL position. When the causative SNP was excluded from the analysis, 714 replicates detected evidence of a QTL (p = 0.001). In analysis II, the CLD test confirmed 280 causative SNP from 1000 simulations (p = 0.05), i.e. power was 28%. When the effect of the QTL was reduced by doubling the error variance, the power of the test reduced relatively little to 23%. When sequence data were used, the power of the test reduced to 16%. All SNP that were confirmed by the CLD test were positioned in the correct QTL position.</p> <p>Conclusions</p> <p>The CLD test can provide evidence for a causative SNP, but its power may be low in situations with closely linked markers. In such situations, also functional evidence will be needed to definitely conclude whether the SNP is causative or not.</p

Accuracy of breeding values of 'unrelated' individuals predicted by dense SNP genotyping

<p>Abstract</p> <p>Background</p> <p>Recent developments in SNP discovery and high throughput genotyping technology have made the use of high-density SNP markers to predict breeding values feasible. This involves estimation of the SNP effects in a training data set, and use of these estimates to evaluate the breeding values of other 'evaluation' individuals. Simulation studies have shown that these predictions of breeding values can be accurate, when training and evaluation individuals are (closely) related. However, many general applications of genomic selection require the prediction of breeding values of 'unrelated' individuals, i.e. individuals from the same population, but not particularly closely related to the training individuals.</p> <p>Methods</p> <p>Accuracy of selection was investigated by computer simulation of small populations. Using scaling arguments, the results were extended to different populations, training data sets and genome sizes, and different trait heritabilities.</p> <p>Results</p> <p>Prediction of breeding values of unrelated individuals required a substantially higher marker density and number of training records than when prediction individuals were offspring of training individuals. However, when the number of records was 2*N_e*L and the number of markers was 10*N_e*L, the breeding values of unrelated individuals could be predicted with accuracies of 0.88 – 0.93, where N_eis the effective population size and L the genome size in Morgan. Reducing this requirement to 1*N_e*L individuals, reduced prediction accuracies to 0.73–0.83.</p> <p>Conclusion</p> <p>For livestock populations, 1N_eL requires about ~30,000 training records, but this may be reduced if training and evaluation animals are related. A prediction equation is presented, that predicts accuracy when training and evaluation individuals are related. For humans, 1N_eL requires ~350,000 individuals, which means that human disease risk prediction is possible only for diseases that are determined by a limited number of genes. Otherwise, genotyping and phenotypic recording need to become very common in the future.</p

Adaptation of gastrointestinal nematode parasites to host genotype: single locus simulation models

Background: Breeding livestock for improved resistance to disease is an increasingly important selection goal. However, the risk of pathogens adapting to livestock bred for improved disease resistance is difficult to quantify. Here, we explore the possibility of gastrointestinal worms adapting to sheep bred for low faecal worm egg count using computer simulation. Our model assumes sheep and worm genotypes interact at a single locus, such that the effect of an A allele in sheep is dependent on worm genotype, and the B allele in worms is favourable for parasitizing the A allele sheep but may increase mortality on pasture. We describe the requirements for adaptation and test if worm adaptation (1) is slowed by non-genetic features of worm infections and (2) can occur with little observable change in faecal worm egg count. Results: Adaptation in worms was found to be primarily influenced by overall worm fitness, viz. the balance between the advantage of the B allele during the parasitic stage in sheep and its disadvantage on pasture. Genetic variation at the interacting locus in worms could be from de novo or segregating mutations, but de novo mutations are rare and segregating mutations are likely constrained to have (near) neutral effects on worm fitness. Most other aspects of the worm infection we modelled did not affect the outcomes. However, the host-controlled mechanism to reduce faecal worm egg count by lowering worm fecundity reduced the selection pressure on worms to adapt compared to other mechanisms, such as increasing worm mortality. Temporal changes in worm egg count were unreliable for detecting adaptation, despite the steady environment assumed in the simulations. Conclusions: Adaptation of worms to sheep selected for low faecal worm egg count requires an allele segregating in worms that is favourable in animals with improved resistance but less favourable in other animals. Obtaining alleles with this specific property seems unlikely. With support from experimental data, we conclude that selection for low faecal worm egg count should be stable over a short time frame (e.g. 20 years). We are further exploring model outcomes with multiple loci and comparing outcomes to other control strategies

A two-domain elevator mechanism for sodium/proton antiport

Sodium/proton (Na+/H+) antiporters, located at the plasma membrane in every cell, are vital for cell homeostasis1. In humans, their dysfunction has been linked to diseases, such as hypertension, heart failure and epilepsy, and they are well-established drug targets2. The best understood model system for Na+/H+ antiport is NhaA from Escherichia coli1, 3, for which both electron microscopy and crystal structures are available4, 5, 6. NhaA is made up of two distinct domains: a core domain and a dimerization domain. In the NhaA crystal structure a cavity is located between the two domains, providing access to the ion-binding site from the inward-facing surface of the protein1, 4. Like many Na+/H+ antiporters, the activity of NhaA is regulated by pH, only becoming active above pH 6.5, at which point a conformational change is thought to occur7. The only reported NhaA crystal structure so far is of the low pH inactivated form4. Here we describe the active-state structure of a Na+/H+ antiporter, NapA from Thermus thermophilus, at 3 Å resolution, solved from crystals grown at pH 7.8. In the NapA structure, the core and dimerization domains are in different positions to those seen in NhaA, and a negatively charged cavity has now opened to the outside. The extracellular cavity allows access to a strictly conserved aspartate residue thought to coordinate ion binding1, 8, 9 directly, a role supported here by molecular dynamics simulations. To alternate access to this ion-binding site, however, requires a surprisingly large rotation of the core domain, some 20° against the dimerization interface. We conclude that despite their fast transport rates of up to 1,500 ions per second3, Na+/H+ antiporters operate by a two-domain rocking bundle model, revealing themes relevant to secondary-active transporters in general

Two-loop Yang-Mills diagrams from superstring amplitudes

Starting from the superstring amplitude describing interactions among D-branes with a constant world-volume field strength, we present a detailed analysis of how the open string degeneration limits reproduce the corresponding field theory Feynman diagrams. A key ingredient in the string construction is represented by the twisted (Prym) super differentials, as their periods encode the information about the background field. We provide an efficient method to calculate perturbatively the determinant of the twisted period matrix in terms of sets of super-moduli appropriate to the degeneration limits. Using this result we show that there is a precise one-to-one correspondence between the degeneration of different factors in the superstring amplitudes and one-particle irreducible Feynman diagrams capturing the gauge theory effective action at the two-loop level.Comment: 42 pages plus appendices, 10 figure

SDSS-IV MaNGA: Spatially resolved star formation histories in galaxies as a function of galaxy mass and type

© 2016 The Authors. We study the internal gradients of stellar population propertieswithin 1.5 Re for a representative sample of 721 galaxies, with stellar masses ranging between 109M⊙ and 1011.5M⊙ from the SDSS-IV MaNGA Integral-Field-Unit survey. Through the use of our full spectral fitting code FIREFLY, we derive light- and mass-weighted stellar population properties and their radial gradients, as well as full star formation and metal enrichment histories. We also quantify the impact that different stellar population models and full spectral fitting routines have on the derived stellar population properties and the radial gradient measurements. In our analysis, we find that age gradients tend to be shallow for both early-type and late-type galaxies. Mass-weighted age gradients of early-types arepositive (~0.09 dex/Re) pointing to 'outside- in' progression of star formation, while late-type galaxies have negative light-weighted age gradients (~-0.11 dex/Re), suggesting an 'inside-out' formation of discs. We detect negative metallicity gradients in both early- and late-type galaxies, but these are significantly steeper in late-types, suggesting that the radial dependence of chemical enrichment processes and the effect of gas inflow and metal transport are far more pronounced in discs. Metallicity gradients of both morphological classes correlate with galaxy mass, with negative metallicity gradients becoming steeper with increasing galaxy mass. The correlation with mass is stronger for late-type galaxies, with a slope of d(∇[Z/H])/d(logM) ~ -0.2 ± 0.05, compared to d(∇[Z/H])/d(logM) ~ -0.05 ± 0.05 for early-types. This result suggests that the merger history plays a relatively small role in shaping metallicity gradients of galaxies

The predictive value of early behavioural assessments in pet dogs: a longitudinal study from neonates to adults

Studies on behavioural development in domestic dogs are of relevance for matching puppies with the right families, identifying predispositions for behavioural problems at an early stage, and predicting suitability for service dog work, police or military service. The literature is, however, inconsistent regarding the predictive value of tests performed during the socialisation period. Additionally, some practitioners use tests with neonates to complement later assessments for selecting puppies as working dogs, but these have not been validated. We here present longitudinal data on a cohort of Border collies, followed up from neonate age until adulthood. A neonate test was conducted with 99 Border collie puppies aged 2–10 days to assess activity, vocalisations when isolated and sucking force. At the age of 40–50 days, 134 puppies (including 93 tested as neonates) were tested in a puppy test at their breeders' homes. All dogs were adopted as pet dogs and 50 of them participated in a behavioural test at the age of 1.5 to 2 years with their owners. Linear mixed models found little correspondence between individuals' behaviour in the neonate, puppy and adult test. Exploratory activity was the only behaviour that was significantly correlated between the puppy and the adult test. We conclude that the predictive validity of early tests for predicting specific behavioural traits in adult pet dogs is limited

Quantifying the relative roles of selective and neutral processes in defining eukaryotic microbial communities

We have a limited understanding of the relative contributions of different processes that regulate microbial communities, which are crucial components of both natural and agricultural ecosystems. The contributions of selective and neutral processes in defining community composition are often confounded in field studies because as one moves through space, environments also change. Managed ecosystems provide an excellent opportunity to control for this and evaluate the relative strength of these processes by minimising differences between comparable niches separated at different geographic scales. We use next-generation sequencing to characterize the variance in fungal communities inhabiting adjacent fruit, soil and bark in comparable vineyards across 1000 kms in New Zealand. By compartmentalizing community variation, we reveal that niche explains at least four times more community variance than geographic location. We go beyond merely demonstrating that different communities are found in both different niches and locations by quantifying the forces that define these patterns. Overall, selection unsurprisingly predominantly shapes these microbial communities, but we show the balance of neutral processes also have a significant role in defining community assemblage in eukaryotic microbes

The GABA transporter 1 (SLC6A1): a novel candidate gene for anxiety disorders

Recent evidence suggests that the GABA transporter 1 (GAT-1; SLC6A1) plays a role in the pathophysiology and treatment of anxiety disorders. In order to understand the impact of genetic variation within SLC6A1 on pathological anxiety, we performed a case–control association study with anxiety disorder patients with and without syndromal panic attacks. Using the method of sequential addition of cases, we found that polymorphisms in the 5′ flanking region of SLC6A1 are highly associated with anxiety disorders when considering the severity of syndromal panic attacks as phenotype covariate. Analysing the effect size of the association, we observed a constant increase in the odds ratio for disease susceptibility with an increase in panic severity (OR ~ 2.5 in severely affected patients). Nominally significant association effects were observed considering the entire patient sample. These data indicate a high load of genetic variance within SLC6A1 on pathological anxiety and highlight GAT-1 as a promising target for treatment of anxiety disorders with panic symptoms

Unusually Long Palindromes Are Abundant in Mitochondrial Control Regions of Insects and Nematodes

BACKGROUND: Palindromes are known to be involved in a variety of biological processes. In the present investigation we carried out a comprehensive analysis of palindromes in the mitochondrial control regions (CRs) of several animal groups to study their frequency, distribution and architecture to gain insights into the origin of replication of mtDNA. METHODOLOGY/PRINCIPAL FINDINGS: Many species of Arthropoda, Nematoda, Mollusca and Annelida harbor palindromes and inverted repeats (IRs) in their CRs. Lower animals like cnidarians and higher animal groups like chordates are almost devoid of palindromes and IRs. The study revealed that palindrome occurrence is positively correlated with the AT content of CRs, and that IRs are likely to give rise to longer palindromes. CONCLUSIONS/SIGNIFICANCE: The present study attempts to explain possible reasons and gives in silico evidence for absence of palindromes and IRs from CR of vertebrate mtDNA and acquisition and retention of the same in insects. Study of CRs of different animal phyla uncovered unique architecture of this locus, be it high abundance of long palindromes and IRs in CRs of Insecta and Nematoda, or short IRs of 10–20 nucleotides with a spacer region of 12–14 bases in subphylum Chelicerata, or nearly complete of absence of any long palindromes and IRs in Vertebrata, Cnidaria and Echinodermata