Accurate Inference of Subtle Population Structure (and Other Genetic Discontinuities) Using Principal Coordinates

Accurate inference of genetic discontinuities between populations is an essential component of intraspecific biodiversity and evolution studies, as well as associative genetics. The most widely-used methods to infer population structure are model-based, Bayesian MCMC procedures that minimize Hardy-Weinberg and linkage disequilibrium within subpopulations. These methods are useful, but suffer from large computational requirements and a dependence on modeling assumptions that may not be met in real data sets. Here we describe the development of a new approach, PCO-MC, which couples principal coordinate analysis to a clustering procedure for the inference of population structure from multilocus genotype data.PCO-MC uses data from all principal coordinate axes simultaneously to calculate a multidimensional "density landscape", from which the number of subpopulations, and the membership within subpopulations, is determined using a valley-seeking algorithm. Using extensive simulations, we show that this approach outperforms a Bayesian MCMC procedure when many loci (e.g. 100) are sampled, but that the Bayesian procedure is marginally superior with few loci (e.g. 10). When presented with sufficient data, PCO-MC accurately delineated subpopulations with population F(st) values as low as 0.03 (G'(st)>0.2), whereas the limit of resolution of the Bayesian approach was F(st) = 0.05 (G'(st)>0.35).We draw a distinction between population structure inference for describing biodiversity as opposed to Type I error control in associative genetics. We suggest that discrete assignments, like those produced by PCO-MC, are appropriate for circumscribing units of biodiversity whereas expression of population structure as a continuous variable is more useful for case-control correction in structured association studies

Impact of Adaptation Currents on Synchronization of Coupled Exponential Integrate-and-Fire Neurons

The ability of spiking neurons to synchronize their activity in a network depends on the response behavior of these neurons as quantified by the phase response curve (PRC) and on coupling properties. The PRC characterizes the effects of transient inputs on spike timing and can be measured experimentally. Here we use the adaptive exponential integrate-and-fire (aEIF) neuron model to determine how subthreshold and spike-triggered slow adaptation currents shape the PRC. Based on that, we predict how synchrony and phase locked states of coupled neurons change in presence of synaptic delays and unequal coupling strengths. We find that increased subthreshold adaptation currents cause a transition of the PRC from only phase advances to phase advances and delays in response to excitatory perturbations. Increased spike-triggered adaptation currents on the other hand predominantly skew the PRC to the right. Both adaptation induced changes of the PRC are modulated by spike frequency, being more prominent at lower frequencies. Applying phase reduction theory, we show that subthreshold adaptation stabilizes synchrony for pairs of coupled excitatory neurons, while spike-triggered adaptation causes locking with a small phase difference, as long as synaptic heterogeneities are negligible. For inhibitory pairs synchrony is stable and robust against conduction delays, and adaptation can mediate bistability of in-phase and anti-phase locking. We further demonstrate that stable synchrony and bistable in/anti-phase locking of pairs carry over to synchronization and clustering of larger networks. The effects of adaptation in aEIF neurons on PRCs and network dynamics qualitatively reflect those of biophysical adaptation currents in detailed Hodgkin-Huxley-based neurons, which underscores the utility of the aEIF model for investigating the dynamical behavior of networks. Our results suggest neuronal spike frequency adaptation as a mechanism synchronizing low frequency oscillations in local excitatory networks, but indicate that inhibition rather than excitation generates coherent rhythms at higher frequencies

Transcranial direct current stimulation generates a transient increase of small-world in brain connectivity: an EEG graph theoretical analysis

Transcranial direct current stimulation (tDCS) is a non-invasive technique able to modulate cortical excitability in a polarity-dependent way. At present, only few studies investigated the effects of tDCS on the modulation of functional connectivity between remote cortical areas. The aim of this study was to investigateâthrough graph theory analysisâhow bipolar tDCS modulate cortical networks high-density EEG recordings were acquired before and after bipolar cathodal, anodal and sham tDCS involving the primary motor and pre-motor cortices of the dominant hemispherein 14 healthy subjects. Results showed that, after bipolar anodal tDCS stimulation, brain networks presented a less evident âsmall worldâ organization with a global tendency to be more random in its functional connections with respect to prestimulus condition in both hemispheres. Results suggest that tDCS globally modulates the cortical connectivity of the brain, modifying the underlying functional organization of the stimulated networks, which might be related to changes in synaptic efficiency of the motor network and related brain areas. This study demonstrated that graph analysis approach to EEG recordings is able to intercept changes in cortical functions mediated by bipolar anodal tDCS mainly involving the dominant M1 and related motor areas. Concluding, tDCS could be an useful technique to help understanding brain rhythms and their topographic functional organization and specificity