THE USE OF SIMULATIONS IN EVOLUTIONARY POPULATION GENETICS: APPLICATIONS ON HUMANS, OWLS AND VIRTUAL ORGANISMS

Summary (in English) Computer simulations provide a practical way to address scientific questions that would be otherwise intractable. In evolutionary biology, and in population genetics in particular, the investigation of evolutionary processes frequently involves the implementation of complex models, making simulations a particularly valuable tool in the area. In this thesis work, I explored three questions involving the geographical range expansion of populations, taking advantage of spatially explicit simulations coupled with approximate Bayesian computation. First, the neutral evolutionary history of the human spread around the world was investigated, leading to a surprisingly simple model: A straightforward diffusion process of migrations from east Africa throughout a world map with homogeneous landmasses replicated to very large extent the complex patterns observed in real human populations, suggesting a more continuous (as opposed to structured) view of the distribution of modern human genetic diversity, which may play a better role as a base model for further studies. Second, the postglacial evolution of the European barn owl, with the formation of a remarkable coat-color cline, was inspected with two rounds of simulations: (i) determine the demographic background history and (ii) test the probability of a phenotypic cline, like the one observed in the natural populations, to appear without natural selection. We verified that the modern barn owl population originated from a single Iberian refugium and that they formed their color cline, not due to neutral evolution, but with the necessary participation of selection. The third and last part of this thesis refers to a simulation-only study inspired by the barn owl case above. In this chapter, we showed that selection is, indeed, effective during range expansions and that it leaves a distinguished signature, which can then be used to detect and measure natural selection in range-expanding populations. Résumé (en français) Les simulations fournissent un moyen pratique pour répondre à des questions scientifiques qui seraient inabordable autrement. En génétique des populations, l'étude des processus évolutifs implique souvent la mise en oeuvre de modèles complexes, et les simulations sont un outil particulièrement précieux dans ce domaine. Dans cette thèse, j'ai exploré trois questions en utilisant des simulations spatialement explicites dans un cadre de calculs Bayésiens approximés (approximate Bayesian computation : ABC). Tout d'abord, l'histoire de la colonisation humaine mondiale et de l'évolution de parties neutres du génome a été étudiée grâce à un modèle étonnement simple. Un processus de diffusion des migrants de l'Afrique orientale à travers un monde avec des masses terrestres homogènes a reproduit, dans une très large mesure, les signatures génétiques complexes observées dans les populations humaines réelles. Un tel modèle continu (opposé à un modèle structuré en populations) pourrait être très utile comme modèle de base dans l'étude de génétique humaine à l'avenir. Deuxièmement, l'évolution postglaciaire d'un gradient de couleur chez l'Effraie des clocher (Tyto alba) Européenne, a été examiné avec deux séries de simulations pour : (i) déterminer l'histoire démographique de base et (ii) tester la probabilité qu'un gradient phénotypique, tel qu'observé dans les populations naturelles puisse apparaître sans sélection naturelle. Nous avons montré que la population actuelle des chouettes est sortie d'un unique refuge ibérique et que le gradient de couleur ne peux pas s'être formé de manière neutre (sans l'action de la sélection naturelle). La troisième partie de cette thèse se réfère à une étude par simulations inspirée par l'étude de l'Effraie. Dans ce dernier chapitre, nous avons montré que la sélection est, en effet, aussi efficace dans les cas d'expansion d'aire de distribution et qu'elle laisse une signature unique, qui peut être utilisée pour la détecter et estimer sa force.

Visual stimulus-dependent changes in interhemispheric EEG coherence in humans.

We analyzed the coherence of electroencephalographic (EEG) signals recorded symmetrically from the two hemispheres, while subjects (n = 9) were viewing visual stimuli. Considering the many common features of the callosal connectivity in mammals, we expected that, as in our animal studies, interhemispheric coherence (ICoh) would increase only with bilateral iso-oriented gratings located close to the vertical meridian of the visual field, or extending across it. Indeed, a single grating that extended across the vertical meridian significantly increased the EEG ICoh in normal adult subjects. These ICoh responses were obtained from occipital and parietal derivations and were restricted to the gamma frequency band. They were detectable with different EEG references and were robust across and within subjects. Other unilateral and bilateral stimuli, including identical gratings that were effective in anesthetized animals, did not affect ICoh in humans. This fact suggests the existence of regulatory influences, possibly of a top-down kind, on the pattern of callosal activation in conscious human subjects. In addition to establishing the validity of EEG coherence analysis for assaying cortico-cortical connectivity, this study extends to the human brain the finding that visual stimuli cause interhemispheric synchronization, particularly in frequencies of the gamma band. It also indicates that the synchronization is carried out by cortico-cortical connection and suggests similarities in the organization of visual callosal connections in animals and in man.

Sub-cortical and brainstem sites associated with chemo-stimulated increases in ventilation in humans.

We investigated the neural basis for spontaneous chemo-stimulated increases in ventilation in awake, healthy humans. Blood oxygen level dependent (BOLD) functional MRI was performed in nine healthy subjects using T2 weighted echo planar imaging. Brain volumes (52 transverse slices, cortex to high spinal cord) were acquired every 3.9 s. The 30 min paradigm consisted of six, 5-min cycles, each cycle comprising 45 s of hypoxic-isocapnia, 45 s of isooxic-hypercapnia and 45 s of hypoxic-hypercapnia, with 55 s of non-stimulatory hyperoxic-isocapnia (control) separating each stimulus period. Ventilation was significantly (p<0.001) increased during hypoxic-isocapnia, isooxic-hypercapnia and hypoxic-hypercapnia (17.0, 13.8, 24.9 L/min respectively) vs. control (8.4 L/min) and was associated with significant (p<0.05, corrected for multiple comparisons) signal increases within a bilateral network that included the basal ganglia, thalamus, red nucleus, cerebellum, parietal cortex, cingulate and superior mid pons. The neuroanatomical structures identified provide evidence for the spontaneous control of breathing to be mediated by higher brain centres, as well as respiratory nuclei in the brainstem.

Thymic selection generates a large T cell pool recognizing a self-peptide in humans.

The low frequency of self-peptide-specific T cells in the human preimmune repertoire has so far precluded their direct evaluation. Here, we report an unexpected high frequency of T cells specific for the self-antigen Melan-A/MART-1 in CD8 single-positive thymocytes from human histocompatibility leukocyte antigen-A2 healthy individuals, which is maintained in the peripheral blood of newborns and adults. Postthymic replicative history of Melan-A/MART-1-specific CD8 T cells was independently assessed by quantifying T cell receptor excision circles and telomere length ex vivo. We provide direct evidence that the large T cell pool specific for the self-antigen Melan-A/MART-1 is mostly generated by thymic output of a high number of precursors. This represents the only known naive self-peptide-specific T cell repertoire directly accessible in humans.

Deciphering the unusual HLA-A2/Melan-A/MART-1-specific TCR repertoire in humans.

The Melan-A/MART-1(26-35) antigenic peptide is one of the best studied human tumor-associated antigens. It is expressed in healthy melanocytes and malignant melanoma and is recognized by CD8(+) T cells in the context of the MHC class I molecule HLA-A*0201. While an unusually large repertoire of CD8(+) T cells specific for this antigen has been documented, the reasons for its generation have remained elusive. In this issue of the European Journal of Immunology, Pinto et al. [Eur. J. Immunol. 2014. 44: 2811-2821] uncover one important mechanism by comparing the thymic expression of the Melan-A gene to that in the melanocyte lineage. This study shows that medullary thymic epithelial cells (mTECs) dominantly express a truncated Melan-A transcript, the product of misinitiation of transcription. Consequently, the protein product in mTECs lacks the immunodominant epitope spanning residues 26-35, thus precluding central tolerance to this antigen. In contrast, melanocytes and melanoma tumor cells express almost exclusively the full-length Melan-A transcript, thus providing the target antigen for efficient recognition by HLA-A2-restricted CD8(+) T cells. The frequency of these alternative gene transcription modes may be more common than previously appreciated and may represent an important factor modulating the efficiency of central tolerance induction in the thymus.

Plasma ionized magnesium during acute hyperventilation in humans.

1. Respiratory alkalosis accompanies the clinical syndrome of tetany, precipitates cardiac arrhythmias and predisposes to coronary vasoconstriction. Magnesium plays a critical role in the maintenance of membrane function, and magnesium depletion is often associated with cardiac arrhythmias or vasoconstriction. 2. As technology for detecting circulating ionized magnesium (the most interesting form with respect to physiological and biological properties) is now available in the form of new magnesium-selective electrodes, the effect of respiratory alkalosis induced by voluntary overbreathing for 30 min on circulating ionized magnesium was studied in eight healthy subjects. 3. The total plasma magnesium concentration was not modified by hyperventilation. On the contrary, hyperventilation was associated with a significant reduction in the ionized magnesium concentration of 0.05 (0.02-0.15) mmol/l (median and range) and in the free magnesium fraction of 0.06 (0.01-0.19). During hyperventilation the relative intravascular magnesium mass, calculated from changes in total plasma magnesium concentration and haematocrit, decreased significantly. 4. It is concluded that acute overbreathing reduces the circulating ionized magnesium concentration and the intravascular magnesium mass. It is therefore conceivable that extracellular magnesium deficiency is at least a subsidiary cause of the syndrome of tetany and the cardiac complications that are precipitated by hyperventilation.

Regulation of gastric and pancreatic lipase secretion by CCK and cholinergic mechanisms in humans.

Gastric lipase (HGL) contributes significantly to fat digestion. However, little is known about its neurohormonal regulation in humans. We studied the role of CCK and cholinergic mechanisms in the postprandial regulation of HGL and pancreatic lipase (HPL) secretion in six healthy subjects. Gastric emptying of a mixed meal and outputs of HGL, pepsin, acid, and HPL were determined with a double-indicator technique. Three experiments were performed in random order: intravenous infusion of 1) placebo, 2) low-dose atropine (5 micrograms.kg-.h-1), and 3) the CCK-A receptor antagonist loxiglumide (22 mumol.kg-.h-1). Atropine decreased postprandial outputs of HGL, pepsin, gastric acid, and HPL (P < 0.03) while slowing gastric emptying (P < 0.05). Loxiglumide markedly increased the secretion of HGL, pepsin, and acid while distinctly reducing HPL outputs and accelerating gastric emptying (P < 0.03). Plasma CCK and gastrin levels increased during loxiglumide infusion (P < 0.03). Atropine enhanced gastrin but not CCK release. Postprandial HGL, pepsin, and acid secretion are under positive cholinergic but negative CCK control, whereas HPL is stimulated by cholinergic and CCK mechanisms. We conclude that CCK and cholinergic mechanisms have an important role in the coordination of HGL and HPL secretion to optimize digestion of dietary lipids in humans.

A new electrophysiological index for working memory load in humans.

Working memory, the ability to store and simultaneously manipulate information, is affected in several neuropsychiatric disorders which lead to severe cognitive and functional deficits. An electrophysiological marker for this process could help identify early cerebral function abnormalities. In subjects performing working memory-specific n-back tasks, event-related potential analysis revealed a positive-negative waveform (PNwm) component modulated in amplitude by working memory load. It occurs in the expected time range for this process, 140-280 ms after stimulus onset, superimposed on the classical P200 and N200 components. Independent Component Analysis extracted two functional components with latencies and topographical scalp distributions similar to the PNwm. Our results imply that the PNwm represents a new electrophysiological index for working memory load in humans.

Simulation of amoxicillin pharmacokinetics in humans for the prevention of streptococcal endocarditis in rats.

The pharmacokinetic determinants of successful antibiotic prophylaxis of endocarditis are not precisely known. Differences in half-lives of antibiotics between animals and humans preclude extrapolation of animal results to human situations. To overcome this limitation, we have mimicked in rats the amoxicillin kinetics in humans following a 3-g oral dose (as often used for prophylaxis of endocarditis) by delivering the drug through a computerized pump. Rats with catheter-induced vegetations were challenged with either of two strains of antibiotic-tolerant viridans group streptococci. Antibiotics were given either through the pump (to simulate the whole kinetic profile during prophylaxis in humans) or as an intravenous bolus which imitated only the peak level of amoxicillin (18 mg/liter) in human serum. Prophylaxis by intravenous bolus was inoculum dependent and afforded a limited protection only in rats challenged with the minimum inoculum size infecting > or = 90% of untreated controls. In contrast, simulation of kinetics in humans significantly protected animals challenged with 10 to 100 times the inoculum of either of the test organisms infecting > or = 90% of untreated controls. Thus, simulation of the profiles of amoxicillin prophylaxis in human serum was more efficacious than mere imitation of the transient peak level in rats. This confirms previous studies suggesting that the duration for which the serum amoxicillin level remained detectable (not only the magnitude of the peak) was an important parameter in successful prophylaxis of endocarditis. The results also suggest that single-dose prophylaxis with 3 g of amoxicillin in humans might be more effective than predicted by conventional animal models in which only peak levels of antibiotic in human serum were stimulated.

Carbon isotopic ratio analysis by gas chromatography/combustion/isotope ratio mass spectrometry for the detection of gamma-hydroxybutyric acid (GHB) administration to humans.

Since GHB (gamma-hydroxybutyric acid) is naturally produced in the human body, clinical and forensic toxicologists must be able to discriminate between endogenous levels and a concentration resulting from exposure. To suggest an alternative to the use of interpretative concentration cut-offs, the detection of exogenous GHB in urine specimens was investigated by means of gas chromatography/combustion/isotope ratio mass spectrometry (GC/C/IRMS). GHB was isolated from urinary matrix by successive purification on Oasis MCX and Bond Elute SAX solid-phase extraction (SPE) cartridges prior to high-performance liquid chromatography (HPLC) fractioning using an Atlantis dC18 column eluted with a mixture of formic acid and methanol. Subsequent intramolecular esterification of GHB leading to the formation of gamma-butyrolactone (GBL) was carried out to avoid introduction of additional carbon atoms for carbon isotopic ratio analysis. A precision of 0.3 per thousand was determined using this IRMS method for samples at GHB concentrations of 10 mg/L. The (13)C/(12)C ratios of GHB in samples of subjects exposed to the drug ranged from -32.1 to -42.1 per thousand, whereas the results obtained for samples containing GHB of endogenous origin at concentration levels less than 10 mg/L were in the range -23.5 to -27.0 per thousand. Therefore, these preliminary results show that a possible discrimination between endogenous and exogenous GHB can be made using carbon isotopic ratio analyses.

Stimulation-induced increases of astrocytic oxidative metabolism in rats and humans investigated with 1-11C-acetate

The purpose of this study was to investigate astrocytic oxidative metabolism using 1-(11)C-acetate. 1-(11)C-acetate kinetics were evaluated in the rat somatosensory cortex using a beta-scintillator during different manipulations (test-retest, infraorbital nerve stimulation, and administration of acetazolamide or dichloroacetate). In humans a visual activation paradigm was used and kinetics were measured with positron emission tomography. Data were analyzed using a one-tissue compartment model. The following features supported the hypothesis that washout of radiolabel (k(2)) is because of (11)C-CO(2) and therefore related to oxygen consumption (CMRO(2)): (1) the onset of (11)C washout was delayed; (2)k(2) was not affected by acetazolamide-induced blood flow increase; (3)k(2) demonstrated a significant increase during stimulation in rats (from 0.014+/-0.007 to 0.027+/-0.006 per minute) and humans (from 0.016+/-0.010 to 0.026+/-0.006 per minute); and (4) dichloroacetate led to a substantial decrease of k(2). In the test-retest experiments K(1) and k(2) were very stable. In summary, 1-(11)C-acetate seems a promising tracer to investigate astrocytic oxidative metabolism in vivo. If the washout rate indeed represents the production of (11)C-CO(2), then its increase during stimulation would point to a substantially higher astrocytic oxidative metabolism during brain activation. However, the quantitative relationship between k(2) and CMRO(2) needs to be determined in future experiments.

Renal and hormonal responses to direct renin inhibition with aliskiren in healthy humans.

BACKGROUND: Pharmacological interruption of the renin-angiotensin system focuses on optimization of blockade. As a measure of intrarenal renin activity, we have examined renal plasma flow (RPF) responses in a standardized protocol. Compared with responses with angiotensin-converting enzyme inhibition (rise in RPF approximately 95 mL x min(-1) x 1.73 m(-2)), greater renal vasodilation with angiotensin receptor blockers (approximately 145 mL x min(-1) x 1.73 m(-2)) suggested more effective blockade. We predicted that blockade with the direct oral renin inhibitor aliskiren would produce renal vascular responses exceeding those induced by angiotensin-converting enzyme inhibitors and angiotensin receptor blockers. METHODS AND RESULTS: Twenty healthy normotensive subjects were studied on a low-sodium (10 mmol/d) diet, receiving separate escalating doses of aliskiren. Six additional subjects received captopril 25 mg as a low-sodium comparison and also received aliskiren on a high-sodium (200 mmol/d) diet. RPF was measured by clearance of para-aminohippurate. Aliskiren induced a remarkable dose-related renal vasodilation in low-sodium balance. The RPF response was maximal at the 600-mg dose (197+/-27 mL x min(-1) x 1.73 m(-2)) and exceeded responses to captopril (92+/-20 mL x min(-1) x 1.73 m(-2); P<0.01). Furthermore, significant residual vasodilation was observed 48 hours after each dose (P<0.01). The RPF response on a high-sodium diet was also higher than expected (47+/-17 mL x min(-1) x 1.73 m(-2)). Plasma renin activity and angiotensin levels were reduced in a dose-related manner. As another functional index of the effect of aliskiren, we found significant natriuresis on both diets. CONCLUSIONS: Renal vasodilation in healthy people with the potent renin inhibitor aliskiren exceeded responses seen previously with angiotensin-converting enzyme inhibitors and angiotensin receptor blockers. The effects were longer lasting and were associated with significant natriuresis. These results indicate that aliskiren may provide more complete and thus more effective blockade of the renin-angiotensin system.