A theory of oligogenic adaptation of a quantitative trait.

Rapid phenotypic adaptation is widespread in nature, but the underlying genetic dynamics remain controversial. Whereas population genetics envisages sequential beneficial substitutions, quantitative genetics assumes a collective response through subtle shifts in allele frequencies. This dichotomy of a monogenic and a highly polygenic view of adaptation raises the question of a middle ground, as well as the factors controlling the transition. Here, we consider an additive quantitative trait with equal locus effects under Gaussian stabilizing selection that adapts to a new trait optimum after an environmental change. We present an analytical framework based on Yule branching processes to describe how phenotypic adaptation is achieved by collective changes in allele frequencies at the underlying loci. In particular, we derive an approximation for the joint allele-frequency distribution conditioned on the trait mean as a comprehensive descriptor of the adaptive architecture. Depending on the model parameters, this architecture reproduces the well-known patterns of sequential, monogenic sweeps, or of subtle, polygenic frequency shifts. Between these endpoints, we observe oligogenic architecture types that exhibit characteristic patterns of partial sweeps. We find that a single compound parameter, the population-scaled background mutation rate Θbg, is the most important predictor of the type of adaptation, while selection strength, the number of loci in the genetic basis, and linkage only play a minor role.

Polygenic adaptation: From sweeps to subtle frequency shifts.

Evolutionary theory has produced two conflicting paradigms for the adaptation of a polygenic trait. While population genetics views adaptation as a sequence of selective sweeps at single loci underlying the trait, quantitative genetics posits a collective response, where phenotypic adaptation results from subtle allele frequency shifts at many loci. Yet, a synthesis of these views is largely missing and the population genetic factors that favor each scenario are not well understood. Here, we study the architecture of adaptation of a binary polygenic trait (such as resistance) with negative epistasis among the loci of its basis. The genetic structure of this trait allows for a full range of potential architectures of adaptation, ranging from sweeps to small frequency shifts. By combining computer simulations and a newly devised analytical framework based on Yule branching processes, we gain a detailed understanding of the adaptation dynamics for this trait. Our key analytical result is an expression for the joint distribution of mutant alleles at the end of the adaptive phase. This distribution characterizes the polygenic pattern of adaptation at the underlying genotype when phenotypic adaptation has been accomplished. We find that a single compound parameter, the population-scaled background mutation rate Θbg, explains the main differences among these patterns. For a focal locus, Θbg measures the mutation rate at all redundant loci in its genetic background that offer alternative ways for adaptation. For adaptation starting from mutation-selection-drift balance, we observe different patterns in three parameter regions. Adaptation proceeds by sweeps for small Θbg ≲ 0.1, while small polygenic allele frequency shifts require large Θbg ≳ 100. In the large intermediate regime, we observe a heterogeneous pattern of partial sweeps at several interacting loci.

Bounds to parapatric speciation: A Dobzhansky-Muller incompatibility model involving autosomes, X chromosomes, and mitochondria.

We investigate the conditions for the origin and maintenance of postzygotic isolation barriers, so called (Bateson-)Dobzhansky-Muller incompatibilities or DMIs, among populations that are connected by gene flow. Specifically, we compare the relative stability of pairwise DMIs among autosomes, X chromosomes, and mitochondrial genes. In an analytical approach based on a continent-island framework, we determine how the maximum permissible migration rates depend on the genomic architecture of the DMI, on sex bias in migration rates, and on sex-dependence of allelic and epistatic effects, such as dosage compensation. Our results show that X-linkage of DMIs can enlarge the migration bounds relative to autosomal DMIs or autosome-mitochondrial DMIs, in particular in the presence of dosage compensation. The effect is further strengthened with male-biased migration. This mechanism might contribute to a higher density of DMIs on the X chromosome (large X-effect) that has been observed in several species clades. Furthermore, our results agree with empirical findings of higher introgression rates of autosomal compared to X-linked loci.

Between- and within-site variability of fMRI localizations.

This study provides first data about the spatial variability of fMRI sensorimotor localizations when investigating the same subjects at different fMRI sites. Results are comparable to a previous patient study. We found a median between-site variability of about 6 mm independent of task (motor or sensory) and experimental standardization (high or low). An intraclass correlation coefficient analysis using data quality measures indicated a major influence of the fMRI site on variability. In accordance with this, within-site localization variability was considerably lower (about 3 mm). We conclude that the fMRI site is a considerable confound for localization of brain activity. However, when performed by experienced clinical fMRI experts, brain pathology does not seem to have a relevant impact on the reliability of fMRI localizations. Hum Brain Mapp 37:2151-2160, 2016. © 2016 Wiley Periodicals, Inc.

Variability of clinical functional MR imaging results: a multicenter study.

PURPOSE: To investigate intersite variability of clinical functional magnetic resonance (MR) imaging, including influence of task standardization on variability and use of various parameters to inform the clinician whether the reliability of a given functional localization is high or low. MATERIALS AND METHODS: Local ethics committees approved the study; all participants gave written informed consent. Eight women and seven men (mean age, 40 years) were prospectively investigated at three experienced functional MR sites with 1.5- (two sites) or 3-T (one site) MR. Nonstandardized motor and highly standardized somatosensory versions of a frequently requested clinical task (localization of the primary sensorimotor cortex) were used. Perirolandic functional MR variability was assessed (peak activation variability, center of mass [COM] variability, intraclass correlation values, overlap ratio [OR], activation size ratio). Data quality measures for functional MR images included percentage signal change (PSC), contrast-to-noise ratio (CNR), and head motion parameters. Data were analyzed with analysis of variance and a correlation analysis. RESULTS: Localization of perirolandic functional MR activity differed by 8 mm (peak activity) and 6 mm (COM activity) among sites. Peak activation varied up to 16.5 mm (COM range, 0.4-16.5 mm) and 45.5 mm (peak activity range, 1.8-45.5 mm). Signal strength (PSC, CNR) was significantly lower for the somatosensory task (mean PSC, 1.0% ± 0.5 [standard deviation]; mean CNR, 1.2 ± 0.4) than for the motor task (mean PSC, 2.4% ± 0.8; mean CNR, 2.9 ± 0.9) (P < .001, both). Intersite variability was larger with low signal strength (negative correlations between signal strength and peak activation variability) even if the task was highly standardized (mean OR, 22.0% ± 18.9 [somatosensory task] and 50.1% ± 18.8 [motor task]). CONCLUSION: Clinical practice and clinical functional MR biomarker studies should consider that the center of task-specific brain activation may vary up to 16.5 mm, with the investigating site, and should maximize functional MR signal strength and evaluate reliability of local results with PSC and CNR.

New type of cortical neuroplasticity after nerve repair in brachial plexus lesions.

BACKGROUND: In brachial plexus avulsion, a recent technique connects the ending of the disrupted musculocutaneous nerve to the side of the intact phrenic nerve to regain elbow flexion. This requires the phrenic nerve to perform a new double function: independent control of breathing and elbow flexion. Neuroplastic changes associated with acquisition of double nerve functions have not yet been investigated. OBJECTIVE: To evaluate neuroplastic changes associated with acquisition of double nerve functions in a monofunctional nerve (phrenic nerve). DESIGN: Clinical and functional magnetic resonance imaging investigations during arm movements, forced inspiration, and motor control tasks. SETTING: Investigations at the Medical University of Vienna, Vienna, Austria. PARTICIPANTS: Three healthy control subjects, 2 patients with phrenic nerve end-to-side coaptation, and 1 control patient with C7 end-to-end coaptation (same clinical presentation but phrenic nerve unchanged). RESULTS: Clinical documentation showed that both patients with phrenic nerve end-to-side coaptation were able to control the diaphragm and the biceps independently via the same phrenic nerve. In contrast to all controls, both patients with phrenic nerve end-to-side coaptation activated the cortical diaphragm areas with flexion of the diseased arm. CONCLUSION: Our functional magnetic resonance imaging data indicate that the patient's cortical diaphragm areas reorganize in such a way that independent control of breathing and elbow flexion is possible with the same neuronal population.

An fMRI marker for peripheral nerve regeneration.

Final outcome after surgical repair of peripheral nerve transections varies. Here, we present the first longitudinal functional magnetic resonance imaging (fMRI) observation of cortical somatosensory reorganization patterns after surgery. A 43-year-old man presented with isolated complete transection of the right median nerve and underwent immediate epineural end-to-end coaptation. Applying standardized vibrotactile median nerve stimulation, 3 T brain activation maps were evaluated at 1, 7, 15 weeks and 1 year after surgery. Initially, the affected hemisphere showed no primary activation but increased frontoparietal activity. After 1 year, primary activation had recovered, and frontoparietal activity was decreased relative to the nonaffected hemisphere. Based on these longitudinal fMRI patterns, we propose a new marker for restoration of somatosensory function, which may not be provided by electrophysiological methods.

How much are clinical fMRI reports influenced by standard postprocessing methods? An investigation of normalization and region of interest effects in the medial temporal lobe.

Recent evidence has indicated that standard postprocessing methods such as template-based region of interest (ROI) definition and normalization of individual brains to a standard template may influence final outcome of functional magnetic resonance imaging investigations. Here, we provide the first comprehensive investigation into whether ROI definition and normalization may also change the clinical interpretation of patient data. A series of medial temporal lobe epilepsy patients were investigated with a clinical memory paradigm and individually delineated as well as template-based ROIs. Different metrics for activation quantification were applied. Results show that the application of template-based ROIs can significantly change the clinical interpretation of individual patient data. This relates to sensitivity for brain activation and hemispheric dominance. We conclude that individual ROIs should be defined on nontransformed functional data and that use of more than one metric for activation quantification is beneficial.