Using Classical Population Genetics Tools with Heterochroneous Data: Time Matters!

BACKGROUND:New polymorphism datasets from heterochroneous data have arisen thanks to recent advances in experimental and microbial molecular evolution, and the sequencing of ancient DNA (aDNA). However, classical tools for population genetics analyses do not take into account heterochrony between subsets, despite potential bias on neutrality and population structure tests. Here, we characterize the extent of such possible biases using serial coalescent simulations. METHODOLOGY/PRINCIPAL FINDINGS:We first use a coalescent framework to generate datasets assuming no or different levels of heterochrony and contrast most classical population genetic statistics. We show that even weak levels of heterochrony ( approximately 10% of the average depth of a standard population tree) affect the distribution of polymorphism substantially, leading to overestimate the level of polymorphism theta, to star like trees, with an excess of rare mutations and a deficit of linkage disequilibrium, which are the hallmark of e.g. population expansion (possibly after a drastic bottleneck). Substantial departures of the tests are detected in the opposite direction for more heterochroneous and equilibrated datasets, with balanced trees mimicking in particular population contraction, balancing selection, and population differentiation. We therefore introduce simple corrections to classical estimators of polymorphism and of the genetic distance between populations, in order to remove heterochrony-driven bias. Finally, we show that these effects do occur on real aDNA datasets, taking advantage of the currently available sequence data for Cave Bears (Ursus spelaeus), for which large mtDNA haplotypes have been reported over a substantial time period (22-130 thousand years ago (KYA)). CONCLUSIONS/SIGNIFICANCE:Considering serial sampling changed the conclusion of several tests, indicating that neglecting heterochrony could provide significant support for false past history of populations and inappropriate conservation decisions. We therefore argue for systematically considering heterochroneous models when analyzing heterochroneous samples covering a large time scale

Molecular Evolution of Regulatory Genes in Spruces from Different Species and Continents: Heterogeneous Patterns of Linkage Disequilibrium and Selection but Correlated Recent Demographic Changes

Genes involved in transcription regulation may represent valuable targets in association genetics studies because of their key roles in plant development and potential selection at the molecular level. Selection and demographic signatures at the sequence level were investigated for five regulatory genes belonging to the knox-I family (KN1, KN2, KN3, KN4) and the HD-Zip III family (HB-3) in three Picea species affected by post-glacial recolonization in North America and Europe. To disentangle neutral and selective forces and estimate linkage disequilibrium (LD) on a gene basis, complete or nearly complete gene sequences were analysed. Nucleotide variation within species, haplotype structure, LD, and neutrality tests, in addition to coalescent simulations based on Tajima’s D and Fay and Wu’s H, were estimated. Nucleotide diversity was generally low in all species (average π = 0.002–0.003) and much heterogeneity was seen in LD and selection signatures among genes and species. Most of the genes harboured an excess of both rare and frequent alleles in the three species. Simulations showed that this excess was significantly higher than that expected under neutrality and a bottleneck during the Last Glacial Maximum followed by population expansion at the Pleistocene/Holocene boundary or shortly after best explains the correlated sequence patterns. These results indicate that despite recent large demographic changes in the three boreal species from two continents, species-specific selection signatures could still be detected from the analysis of nearly complete regulatory gene sequences. Such different signatures indicate differential subfunctionalization of gene family members in the three congeneric species

GABAergic Neuron Deficit As An Idiopathic Generalized Epilepsy Mechanism: The Role Of BRD2 Haploinsufficiency In Juvenile Myoclonic Epilepsy

Idiopathic generalized epilepsy (IGE) syndromes represent about 30% of all epilepsies. They have strong, but elusive, genetic components and sex-specific seizure expression. Multiple linkage and population association studies have connected the bromodomain-containing gene BRD2 to forms of IGE. In mice, a null mutation at the homologous Brd2 locus results in embryonic lethality while heterozygous Brd2+/− mice are viable and overtly normal. However, using the flurothyl model, we now show, that compared to the Brd2+/+ littermates, Brd2+/− males have a decreased clonic, and females a decreased tonic-clonic, seizure threshold. Additionally, long-term EEG/video recordings captured spontaneous seizures in three out of five recorded Brd2+/− female mice. Anatomical analysis of specific regions of the brain further revealed significant differences in Brd2+/− vs +/+ mice. Specifically, there were decreases in the numbers of GABAergic (parvalbumin- or GAD67-immunopositive) neurons along the basal ganglia pathway, i.e., in the neocortex and striatum of Brd2+/− mice, compared to Brd2+/+ mice. There were also fewer GABAergic neurons in the substantia nigra reticulata (SNR), yet there was a minor, possibly compensatory increase in the GABA producing enzyme GAD67 in these SNR cells. Further, GAD67 expression in the superior colliculus and ventral medial thalamic nucleus, the main SNR outputs, was significantly decreased in Brd2+/− mice, further supporting GABA downregulation. Our data show that the non-channel-encoding, developmentally critical Brd2 gene is associated with i) sex-specific increases in seizure susceptibility, ii) the development of spontaneous seizures, and iii) seizure-related anatomical changes in the GABA system, supporting BRD2's involvement in human IGE

Recent and Recurrent Selective Sweeps of the Antiviral RNAi Gene Argonaute-2 in Three Species of Drosophila

Antagonistic host–parasite interactions can drive rapid adaptive evolution in genes of the immune system, and such arms races may be an important force shaping polymorphism in the genome. The RNA interference pathway gene Argonaute-2 (AGO2) is a key component of antiviral defense in Drosophila, and we have previously shown that genes in this pathway experience unusually high rates of adaptive substitution. Here we study patterns of genetic variation in a 100-kbp region around AGO2 in three different species of Drosophila. Our data suggest that recent independent selective sweeps in AGO2 have reduced genetic variation across a region of more than 50 kbp in Drosophila melanogaster, D. simulans, and D. yakuba, and we estimate that selection has fixed adaptive substitutions in this gene every 30–100 thousand years. The strongest signal of recent selection is evident in D. simulans, where we estimate that the most recent selective sweep involved an allele with a selective advantage of the order of 0.5–1% and occurred roughly 13–60 Kya. To evaluate the potential consequences of the recent substitutions on the structure and function of AGO2, we used fold-recognition and homology-based modeling to derive a structural model for the Drosophila protein, and this suggests that recent substitutions in D. simulans are overrepresented at the protein surface. In summary, our results show that selection by parasites can consistently target the same genes in multiple species, resulting in areas of the genome that have markedly reduced genetic diversity

Awareness of cognitive decline trajectories in asymptomatic individuals at risk for AD

Background: Lack of awareness of cognitive decline (ACD) is common in late-stage Alzheimer’s disease (AD). Recent studies showed that ACD can also be reduced in the early stages. Methods: We described different trends of evolution of ACD over 3 years in a cohort of memory-complainers and their association to amyloid burden and brain metabolism. We studied the impact of ACD at baseline on cognitive scores’ evolution and the association between longitudinal changes in ACD and in cognitive score. Results: 76.8% of subjects constantly had an accurate ACD (reference class). 18.95% showed a steadily heightened ACD and were comparable to those with accurate ACD in terms of demographic characteristics and AD biomarkers. 4.25% constantly showed low ACD, had significantly higher amyloid burden than the reference class, and were mostly men. We found no overall effect of baseline ACD on cognitive scores’ evolution and no association between longitudinal changes in ACD and in cognitive scores. Conclusions: ACD begins to decrease during the preclinical phase in a group of individuals, who are of great interest and need to be further characterized. Trial registration: The present study was conducted as part of the INSIGHT-PreAD study. The identification number of INSIGHT-PreAD study (ID-RCB) is 2012-A01731-42

Genetic Impact of a Severe El Niño Event on Galápagos Marine Iguanas (Amblyrhynchus cristatus)

The El Niño-Southern Oscillation (ENSO) is a major source of climatic disturbance, impacting the dynamics of ecosystems worldwide. Recent models predict that human-generated rises in green-house gas levels will cause an increase in the strength and frequency of El Niño warming events in the next several decades, highlighting the need to understand the potential biological consequences of increased ENSO activity. Studies have focused on the ecological and demographic implications of El Niño in a range of organisms, but there have been few systematic attempts to measure the impact of these processes on genetic diversity in populations. Here, we evaluate whether the 1997–1998 El Niño altered the genetic composition of Galápagos marine iguana populations from eleven islands, some of which experienced mortality rates of up to 90% as a result of El Niño warming. Specifically, we measured the temporal variation in microsatellite allele frequencies and mitochondrial DNA diversity (mtDNA) in samples collected before (1991/1993) and after (2004) the El Niño event. Based on microsatellite data, only one island (Marchena) showed signatures of a genetic bottleneck, where the harmonic mean of the effective population size (Ne) was estimated to be less than 50 individuals during the period between samplings. Substantial decreases in mtDNA variation between time points were observed in populations from just two islands (Marchena and Genovesa). Our results suggests that, for the majority of islands, a single, intense El Niño event did not reduce marine iguana populations to the point where substantial neutral genetic diversity was lost. In the case of Marchena, simultaneous changes to both nuclear and mitochondrial DNA variation may also be the result of a volcanic eruption on the island in 1991. Therefore, studies that seek to evaluate the genetic impact of El Niño must also consider the confounding or potentially synergistic effect of other environmental and biological forces shaping populations

Low Cognitive Awareness, but Not Complaint, is a Good Marker of Preclinical Alzheimer's Disease

Subjective cognitive decline (SCD) may result from many conditions, including Alzheimer's disease (AD)

Large-Scale Selective Sweep among Segregation Distorter Chromosomes in African Populations of Drosophila melanogaster

Segregation Distorter (SD) is a selfish, coadapted gene complex on chromosome 2 of Drosophila melanogaster that strongly distorts Mendelian transmission; heterozygous SD/SD+ males sire almost exclusively SD-bearing progeny. Fifty years of genetic, molecular, and theory work have made SD one of the best-characterized meiotic drive systems, but surprisingly the details of its evolutionary origins and population dynamics remain unclear. Earlier analyses suggested that the SD system arose recently in the Mediterranean basin and then spread to a low, stable equilibrium frequency (1–5%) in most natural populations worldwide. In this report, we show, first, that SD chromosomes occur in populations in sub-Saharan Africa, the ancestral range of D. melanogaster, at a similarly low frequency (∼2%), providing evidence for the robustness of its equilibrium frequency but raising doubts about the Mediterranean-origins hypothesis. Second, our genetic analyses reveal two kinds of SD chromosomes in Africa: inversion-free SD chromosomes with little or no transmission advantage; and an African-endemic inversion-bearing SD chromosome, SD-Mal, with a perfect transmission advantage. Third, our population genetic analyses show that SD-Mal chromosomes swept across the African continent very recently, causing linkage disequilibrium and an absence of variability over 39% of the length of the second chromosome. Thus, despite a seemingly stable equilibrium frequency, SD chromosomes continue to evolve, to compete with one another, or evade suppressors in the genome

Chromosomal Rearrangements between Serotype A and D Strains in Cryptococcus neoformans

Cryptococcus neoformans is a major human pathogenic fungus that can cause meningoencephalitis in immunocompromised hosts. It contains two divergent varieties, var. grubii (serotype A) and var. neoformans (serotype D), as well as hybrids (serotype AD) between these two varieties. In this study, we investigated the extent of chromosomal rearrangements between the two varieties, estimated the effects of chromosomal rearrangements on recombination frequencies, and surveyed the potential polymorphisms of the rearrangements among natural strains of the three serotypes. Through the analyses of two sequenced genomes from strains H99 (representing var. grubii) and JEC21 (representing var. neoformans), we revealed a total of 32 unambiguous chromosome rearrangements, including five translocations, nine simple inversions, and 18 complex rearrangements. Our analyses identified that overall, rearranged regions had recombination frequencies about half of those around syntenic regions. Using a direct PCR screening strategy, we examined the potential polymorphisms of 11 rearrangements among 64 natural C. neoformans strains from five countries. We found no polymorphism within var. neoformans and very limited polymorphism within var. grubii. However, strains of serotype AD showed significant polymorphism, consistent with their hybrid origins coupled with differential loss of heterozygosity. We discuss the implications of these results on the genome structure, ecology, and evolution of C. neoformans

Convergent Processing of Both Positive and Negative Motivational Signals by the VTA Dopamine Neuronal Populations

Dopamine neurons in the ventral tegmental area (VTA) have been traditionally studied for their roles in reward-related motivation or drug addiction. Here we study how the VTA dopamine neuron population may process fearful and negative experiences as well as reward information in freely behaving mice. Using multi-tetrode recording, we find that up to 89% of the putative dopamine neurons in the VTA exhibit significant activation in response to the conditioned tone that predict food reward, while the same dopamine neuron population also respond to the fearful experiences such as free fall and shake events. The majority of these VTA putative dopamine neurons exhibit suppression and offset-rebound excitation, whereas ∼25% of the recorded putative dopamine neurons show excitation by the fearful events. Importantly, VTA putative dopamine neurons exhibit parametric encoding properties: their firing change durations are proportional to the fearful event durations. In addition, we demonstrate that the contextual information is crucial for these neurons to respectively elicit positive or negative motivational responses by the same conditioned tone. Taken together, our findings suggest that VTA dopamine neurons may employ the convergent encoding strategy for processing both positive and negative experiences, intimately integrating with cues and environmental context