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

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

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

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

Specific In Vivo Staining of Astrocytes in the Whole Brain after Intravenous Injection of Sulforhodamine Dyes

Fluorescent staining of astrocytes without damaging or interfering with normal brain functions is essential for intravital microscopy studies. Current methods involved either transgenic mice or local intracerebral injection of sulforhodamine 101. Transgenic rat models rarely exist, and in mice, a backcross with GFAP transgenic mice may be difficult. Local injections of fluorescent dyes are invasive. Here, we propose a non-invasive, specific and ubiquitous method to stain astrocytes in vivo. This method is based on iv injection of sulforhodamine dyes and is applicable on rats and mice from postnatal age to adulthood. The astrocytes staining obtained after iv injection was maintained for nearly half a day and showed no adverse reaction on astrocytic calcium signals or electroencephalographic recordings in vivo. The high contrast of the staining facilitates the image processing and allows to quantify 3D morphological parameters of the astrocytes and to characterize their network. Our method may become a reference for in vivo staining of the whole astrocytes population in animal models of neurological disorders

Epilepsy in Dcx Knockout Mice Associated with Discrete Lamination Defects and Enhanced Excitability in the Hippocampus

Patients with Doublecortin (DCX) mutations have severe cortical malformations associated with mental retardation and epilepsy. Dcx knockout (KO) mice show no major isocortical abnormalities, but have discrete hippocampal defects. We questioned the functional consequences of these defects and report here that Dcx KO mice are hyperactive and exhibit spontaneous convulsive seizures. Changes in neuropeptide Y and calbindin expression, consistent with seizure occurrence, were detected in a large proportion of KO animals, and convulsants, including kainate and pentylenetetrazole, also induced seizures more readily in KO mice. We show that the dysplastic CA3 region in KO hippocampal slices generates sharp wave-like activities and possesses a lower threshold for epileptiform events. Video-EEG monitoring also demonstrated that spontaneous seizures were initiated in the hippocampus. Similarly, seizures in human patients mutated for DCX can show a primary involvement of the temporal lobe. In conclusion, seizures in Dcx KO mice are likely to be due to abnormal synaptic transmission involving heterotopic cells in the hippocampus and these mice may therefore provide a useful model to further study how lamination defects underlie the genesis of epileptiform activities

Human and Non-Human Primate Genomes Share Hotspots of Positive Selection

Among primates, genome-wide analysis of recent positive selection is currently limited to the human species because it requires extensive sampling of genotypic data from many individuals. The extent to which genes positively selected in human also present adaptive changes in other primates therefore remains unknown. This question is important because a gene that has been positively selected independently in the human and in other primate lineages may be less likely to be involved in human specific phenotypic changes such as dietary habits or cognitive abilities. To answer this question, we analysed heterozygous Single Nucleotide Polymorphisms (SNPs) in the genomes of single human, chimpanzee, orangutan, and macaque individuals using a new method aiming to identify selective sweeps genome-wide. We found an unexpectedly high number of orthologous genes exhibiting signatures of a selective sweep simultaneously in several primate species, suggesting the presence of hotspots of positive selection. A similar significant excess is evident when comparing genes positively selected during recent human evolution with genes subjected to positive selection in their coding sequence in other primate lineages and identified using a different test. These findings are further supported by comparing several published human genome scans for positive selection with our findings in non-human primate genomes. We thus provide extensive evidence that the co-occurrence of positive selection in humans and in other primates at the same genetic loci can be measured with only four species, an indication that it may be a widespread phenomenon. The identification of positive selection in humans alongside other primates is a powerful tool to outline those genes that were selected uniquely during recent human evolution