Inference of chromosomal inversion dynamics from Pool-Seq data in natural and laboratory populations of Drosophila melanogaster.

Sequencing of pools of individuals (Pool-Seq) represents a reliable and cost-effective approach for estimating genome-wide SNP and transposable element insertion frequencies. However, Pool-Seq does not provide direct information on haplotypes so that, for example, obtaining inversion frequencies has not been possible until now. Here, we have developed a new set of diagnostic marker SNPs for seven cosmopolitan inversions in Drosophila melanogaster that can be used to infer inversion frequencies from Pool-Seq data. We applied our novel marker set to Pool-Seq data from an experimental evolution study and from North American and Australian latitudinal clines. In the experimental evolution data, we find evidence that positive selection has driven the frequencies of In(3R)C and In(3R)Mo to increase over time. In the clinal data, we confirm the existence of frequency clines for In(2L)t, In(3L)P and In(3R)Payne in both North America and Australia and detect a previously unknown latitudinal cline for In(3R)Mo in North America. The inversion markers developed here provide a versatile and robust tool for characterizing inversion frequencies and their dynamics in Pool-Seq data from diverse D. melanogaster populations.

Molecular cophylogenetic relationships between European bats and their ectoparasitic mites (Acari, Spinturnicidae).

Cospeciation between host-parasite species is generally thought to result in mirror-image congruent phylogenies. Incongruence can be explained by mechanisms such as host switching, duplication, failure to speciate and sorting events. To investigate the level of association in the host-parasite relationship between Spinturnicid mites and their bat hosts, we constructed the phylogenetic tree of the genus Spinturnix (Acari, Mesostigmata) and compared it to the host phylogeny. We sequenced 938bp of the mitochondrial 16S rDNA and Cytochrome Oxydase subunit I (COI) genes among eleven morphospecies of Spinturnix collected on 20 European Vespertilionid and Rhinolophid bat species. Phylogenetic reconstruction of hosts and parasites showed statistical evidence for cospeciation and suggested that their evolutionary history involved also failure to speciate events and host switches. The latter seem to be mainly promoted by similar roosting habits of the host. As currently understood, host associations of Spinturnicid mites likely results from a complex interaction between the phylogenetic history of the host and the behaviour and the ecology of both parasite and host.

Spatio-temporal population genetic structure of the parasitic mite Spinturnix bechsteini is shaped by its own demography and the social system of its bat host.

Information about the population genetic structures of parasites is important for an understanding of parasite transmission pathways and ultimately the co-evolution with their hosts. If parasites cannot disperse independently of their hosts, a parasite's population structure will depend upon the host's spatial distribution. Geographical barriers affecting host dispersal can therefore lead to structured parasite populations. However, how the host's social system affects the genetic structure of parasite populations is largely unknown. We used mitochondrial DNA (mtDNA) to describe the spatio-temporal population structure of a contact-transmitted parasitic wing mite (Spinturnix bechsteini) and compared it to that of its social host, the Bechstein's bat (Myotis bechsteinii). We observed no genetic differentiation between mites living on different bats within a colony. This suggests that mites can move freely among bats of the same colony. As expected in case of restricted inter-colony dispersal, we observed a strong genetic differentiation of mites among demographically isolated bat colonies. In contrast, we found a strong genetic turnover between years when we investigated the temporal variation of mite haplotypes within colonies. This can be explained with mite dispersal occuring between colonies and bottlenecks of mite populations within colonies. The observed absence of isolation by distance could be the result from genetic drift and/or from mites dispersing even between remote bat colonies, whose members may meet at mating sites in autumn or in hibernacula in winter. Our data show that the population structure of this parasitic wing mite is influenced by its own demography and the peculiar social system of its bat host.

Sex-chromosome differentiation parallels postglacial range expansion in European tree frogs (Hyla arborea).

Occasional XY recombination is a proposed explanation for the sex-chromosome homomorphy in European tree frogs. Numerous laboratory crosses, however, failed to detect any event of male recombination, and a detailed survey of NW-European Hyla arborea populations identified male-specific alleles at sex-linked loci, pointing to the absence of XY recombination in their recent history. Here, we address this paradox in a phylogeographic framework by genotyping sex-linked microsatellite markers in populations and sibships from the entire species range. Contrasting with postglacial populations of NW Europe, which display complete absence of XY recombination and strong sex-chromosome differentiation, refugial populations of the southern Balkans and Adriatic coast show limited XY recombination and large overlaps in allele frequencies. Geographically and historically intermediate populations of the Pannonian Basin show intermediate patterns of XY differentiation. Even in populations where X and Y occasionally recombine, the genetic diversity of Y haplotypes is reduced below the levels expected from the fourfold drop in copy numbers. This study is the first in which X and Y haplotypes could be phased over the distribution range in a species with homomorphic sex chromosomes; it shows that XY-recombination patterns may differ strikingly between conspecific populations, and that recombination arrest may evolve rapidly (<5000 generations).

MHC-genotype of progeny influenced by parental infection.

In a previous series of in vitro fertilization experiments with mice we found non-random combination of major histocompatibility complex (MHC) haplotypes in the very early embryos. Our results suggested that two selection mechanisms were operating: (i) the eggs selected specific sperm; and (ii) the second meiotic division in the eggs was influenced by the type of sperm that entered the egg. Furthermore, the proportion of MHC-heterozygous embryos varied over time, suggesting that non-random fertilization was dependent on an external factor that changed over time. As a higher frequency of heterozygous individuals correlated with an uncontrolled epidemic by MHV (mouse hepatitis virus), we suggested that MHV-infection might have influenced the outcome of fertilization. Here, we present an experiment that tests this hypothesis. We infected randomly chosen mice with MHV and sham-infected control mice five days before pairing. We recovered the two-cell embryos from the oviduct, cultured them until the blastocyst stage, and determined the genotype of each resulting blastocyst by polymerase chain reaction. We found the pattern that we expected from our previous experiments: virus-infected mice produced more MHC-heterozygous embryos than sham-infected ones. This suggests that parents are able to promote specific combinations of MHC-haplotypes during fertilization according to the presence or absence of a viral infection.

Identification and potential origin of invasive clawed frogs Xenopus (Anura: Pipidae) in Sicily based on mitochondrial and nuclear DNA

African clawed frogs of the widespread polytypic species Xenopus laevis Daudin, 1802 (ranging large parts of sub-Saharan Africa) have been spreading since the 1940s, and have established reproductive populations in Europe, Asia and the Americas, where they can have negative impact as competitors of native amphibians and as disease vectors for chytridomycosis or ranaviruses. Here we use two mitochondrial (cytochrome b, 16S rDNA) and one nuclear (RAG 1: Recombination Associated Gene 1) DNA markers to infer the potential origin of invasive clawed frogs from Sicily that represent the largest invasive population in Europe. Identical mtDNA haplotypes match with those of Xenopus laevis, and Sicilian clawed frogs very probably belong to a lineage from the Cape Region of South Africa, most likely originating from a laboratory stock. Nuclear data support this conclusion. Identical mtDNA sequences (cyt b, 16S) of frogs sampled across their range in Sicily suggest the occurrence of a single source population and a potential bottleneck at their release, but faster evolving multilocus nuclear data (microsatellites, SNPs) on the population genetics would be important in the future to better support this hypothesis

Ultra high throughput sequencing excludes MDH1 as candidate gene for RP28-linked retinitis pigmentosa.

PURPOSE: Mutations in IDH3B, an enzyme participating in the Krebs cycle, have recently been found to cause autosomal recessive retinitis pigmentosa (arRP). The MDH1 gene maps within the RP28 arRP linkage interval and encodes cytoplasmic malate dehydrogenase, an enzyme functionally related to IDH3B. As a proof of concept for candidate gene screening to be routinely performed by ultra high throughput sequencing (UHTs), we analyzed MDH1 in a patient from each of the two families described so far to show linkage between arRP and RP28. METHODS: With genomic long-range PCR, we amplified all introns and exons of the MDH1 gene (23.4 kb). PCR products were then sequenced by short-read UHTs with no further processing. Computer-based mapping of the reads and mutation detection were performed by three independent software packages. RESULTS: Despite the intrinsic complexity of human genome sequences, reads were easily mapped and analyzed, and all algorithms used provided the same results. The two patients were homozygous for all DNA variants identified in the region, which confirms previous linkage and homozygosity mapping results, but had different haplotypes, indicating genetic or allelic heterogeneity. None of the DNA changes detected could be associated with the disease. CONCLUSIONS: The MDH1 gene is not the cause of RP28-linked arRP. Our experimental strategy shows that long-range genomic PCR followed by UHTs provides an excellent system to perform a thorough screening of candidate genes for hereditary retinal degeneration.

Plastid genome characterisation in Brassica and Brassicaceae using a new set of nine SSRs.

We report a new set of nine primer pairs specifically developed for amplification of Brassica plastid SSR markers. The wide utility of these markers is demonstrated for haplotype identification and detection of polymorphism in B. napus, B. nigra, B. oleracea, B. rapa and in related genera Arabidopsis, Camelina, Raphanus and Sinapis. Eleven gene regions (ndhB-rps7 spacer, rbcL-accD spacer, rpl16 intron, rps16 intron, atpB-rbcL spacer, trnE-trnT spacer, trnL intron, trnL-trnF spacer, trnM-atpE spacer, trnR-rpoC2 spacer, ycf3-psaA spacer) were sequenced from a range of Brassica and related genera for SSR detection and primer design. Other sequences were obtained from GenBank/EMBL. Eight out of nine selected SSR loci showed polymorphism when amplified using the new primers and a combined analysis detected variation within and between Brassica species, with the number of alleles detected per locus ranging from 5 (loci MF-6, MF-1) to 11 (locus MF-7). The combined SSR data were used in a neighbour-joining analysis (SMM, D (DM) distances) to group the samples based on the presence and absence of alleles. The analysis was generally able to separate plastid types into taxon-specific groups. Multi-allelic haplotypes were plotted onto the neighbour joining tree. A total number of 28 haplotypes were detected and these differentiated 22 of the 41 accessions screened from all other accessions. None of these haplotypes was shared by more than one species and some were not characteristic of their predicted type. We interpret our results with respect to taxon differentiation, hybridisation and introgression patterns relating to the 'Triangle of U'.

Polymorphisms in toll-like receptor 4 and toll-like receptor 9 influence viral load in a seroincident cohort of HIV-1-infected individuals.

OBJECTIVES: Toll-like receptors (TLRs) are innate immune sensors that are integral to resisting chronic and opportunistic infections. Mounting evidence implicates TLR polymorphisms in susceptibilities to various infectious diseases, including HIV-1. We investigated the impact of TLR single nucleotide polymorphisms (SNPs) on clinical outcome in a seroincident cohort of HIV-1-infected volunteers. DESIGN: We analyzed TLR SNPs in 201 antiretroviral treatment-naive HIV-1-infected volunteers from a longitudinal seroincident cohort with regular follow-up intervals (median follow-up 4.2 years, interquartile range 4.4). Participants were stratified into two groups according to either disease progression, defined as peripheral blood CD4(+) T-cell decline over time, or peak and setpoint viral load. METHODS: Haplotype tagging SNPs from TLR2, TLR3, TLR4, and TLR9 were detected by mass array genotyping, and CD4(+) T-cell counts and viral load measurements were determined prior to antiretroviral therapy initiation. The association of TLR haplotypes with viral load and rapid progression was assessed by multivariate regression models using age and sex as covariates. RESULTS: Two TLR4 SNPs in strong linkage disequilibrium [1063 A/G (D299G) and 1363 C/T (T399I)] were more frequent among individuals with high peak viral load compared with low/moderate peak viral load (odds ratio 6.65, 95% confidence interval 2.19-20.46, P < 0.001; adjusted P = 0.002 for 1063 A/G). In addition, a TLR9 SNP previously associated with slow progression was found less frequently among individuals with high viral setpoint compared with low/moderate setpoint (odds ratio 0.29, 95% confidence interval 0.13-0.65, P = 0.003, adjusted P = 0.04). CONCLUSION: This study suggests a potentially new role for TLR4 polymorphisms in HIV-1 peak viral load and confirms a role for TLR9 polymorphisms in disease progression.

A Novel Homozygous R764H Mutation in CRB1 Causes Autosomal Recessive Retinitis Pigmentosa in a Consanguineous Family

Purpose: Retinitis pigmentosa (RP; MIM 268000) is a hereditary disease characterized by poor night vision and progressive loss of photoreceptors, eventually leading to blindness. This degenerative process primarily affects peripheral vision due to the loss of rods. Autosomal recessive RP (arRP) is clinically and genetically heterogeneous. It has been associated with mutations in different genes, including CRB1 (Crumbs homolog 1). The aim of this study was to determine the causative gene in a Tunisian patient with arRP born to non consanguineous parents.Methods: Four accessible family members were included. They underwent full ophthalmic examination with best corrected Snellen visual acuity, fundus photography and fluoroangiography. Haplotype analyses were used to test linkage in the family to 20 arRP loci, including ABCA4, LRAT, USH2A, RP29, CERKL, CNGA1, CNGB1, CRB1, EYS, RP28, MERTK, NR2E3, PDE6A, PDE6B, RGR, RHO, RLBP1, TULP1. All exons and intron-exon junctions of candidate genes not excluded by haplotype analysis were PCR amplified and directly sequenced.Results: A 39 aged affected member was individualized. Best corrected visual acuity was OR: 20/63, OS: 20/80. Visual loss began at the third decade. Funduscopic examination and FA revealed typical advanced RP changes with bone spicule-shaped pigment deposits in the posterior pole and the mild periphery along with retinal atrophy, narrowing of the vessels and waxy optic discs. Haplotypes analysis revealed homozygosity with microsatellites markers D1S412 and D1S413 on chromosome 1q31.3. These markers flanked the CRB1 gene. Our results excluded linkage of all the other arRP loci/ genes tested. Sequencing of the 12 coding exons and splice sites of CRB1 gene disclosed a homozygous missense mutation in exon 7 at nucleotide c.(2291 G>A), resulting in an Arg to Hist substitution (p.R764H).Conclusions: R764H is a novel mutation associated with CRB1-related arRP. Previously, an R764C mutation was observed. Extending the mutation spectrum of CRB1 with additional families is important for genotype-phenotype correlations.

Chromosomal versus mitochondrial DNA evolution: tracking the evolutionary history of the southwestern european populations of the Sorex araneus group (Mammalia Insectivora)

The shrews of the Sorer araneus group have undergone a spectacular chromosome evolution. The karyotype of Sorer granarius is generally considered ancestral to those of Sorer coronatus and S. araneus. However, a sequence of 777 base pairs of the cytochrome b gene of the mitochondrial DNA (mtDNA) produces a quite different picture: S. granarius is closely related to the populations of S. araneus from the Pyrenees and from the northwestern Alps, whereas S. coronatus and S. araneus from Italy and the southern Alps represent two well-separated lineages. It is suggested that mtDNA and chromosomal evolution are in this case largely independant processes. Whereas mtDNA haplotypes are closely linked to the geographical history of the populations, chromosomal mutations were probably transmitted from one population to another. Available data suggest that the impressive chromosome polymorphism of this group is quite a recent phenomenon.

Parallel evolution of facial stripe patterns in the Neolamprologus brichardi/pulcher species complex endemic to Lake Tanganyika.

Colour pattern diversity can be due to random processes or to natural or sexual selection. Consequently, similarities in colour patterns are not always correlated with common ancestry, but may result from convergent evolution under shared selection pressures or drift. Neolamprologus brichardi and Neolamprologus pulcher have been described as two distinct species based on differences in the arrangement of two dark bars on the operculum. Our study uses DNA sequences of the mitochondrial control region to show that relatedness of haplotypes disagrees with species assignment based on head colour pattern. This suggests repeated parallel evolution of particular stripe patterns. The complete lack of shared haplotypes between populations of the same or different phenotypes reflects strong philopatric behaviour, possibly induced by the cooperative breeding mode in which offspring remain in their natal territory and serve as helpers until they disperse to nearby territories or take over a breeding position. Concordant phylogeographic patterns between N. brichardi/N. pulcher populations and other rock-dwelling cichlids suggest that the same colonization routes have been taken by sympatric species and that these routes were affected by lake level fluctuations in the past.

A polymorphic microsatellite in the tumor necrosis factor alpha promoter identifies an allele unique to the NZW mouse strain.

We have amplified a (CA)n:(GT)n microsatellite from the TNF promoters of a panel of mouse strains using the polymerase chain reaction. The length of the microsatellites was polymorphic, with eight alleles observed among 15 inbred strains bearing seven distinct H-2 haplotypes, and four outbred strains. In B10 congenic strains, the TNF allele detected by microsatellite polymorphism segregated with the MHC, and in recombinant haplotypes (NOD, NZW), it segregated with H-2D. The TNF allele found in the NZW strain (H-2z) was distinct from those of all other haplotypes, consistent with the hypothesis that this strain may carry a genetic defect in TNF production.

The vitamin D receptor gene bAt (CCA) haplotype impairs the response to pegylated-interferon/ribavirin-based therapy in chronic hepatitis C patients.

BACKGROUND: Chronic hepatitis C infection is a major cause of end-stage liver disease. Therapy outcome is influenced by 25-OH vitamin D deficiency. To further address this observation, our study investigates the impact of the vitamin D receptor (NR1I1) haplotype and combined effects of plasma vitamin D levels in a well-described cohort of hepatitis C patients. METHODS: A total of 155 chronic hepatitis C patients were recruited from the Swiss Hepatitis C Cohort Study for NR1I1 genotyping and plasma 25-OH vitamin D level measurement. NR1I1 genotype data and combined effects of plasma 25-OH vitamin D level were analysed regarding therapy response (sustained virological response). RESULTS: A strong association was observed between therapy non-response and the NR1I1 CCA (bAt) haplotype consisting of rs1544410 (BsmI) C, rs7975232 (ApaI) C and rs731236 (TaqI) A alleles. Of the HCV patients carrying the CCA haplotype, 50.3% were non-responders (odds ratio [OR] 1.69, 95% CI 1.07, 2.67; P=0.028). A similar association was observed for the combinational CCCCAA genotype (OR 2.94, 95% CI 1.36, 6.37; P=0.007). The combinational CCCCAA genotype was confirmed as an independent risk factor for non-response in multivariate analysis (OR 2.50, 95% CI 1.07, 5.87; P=0.034). Analysing combined effects, a significant impact of low 25-OH vitamin D levels on sustained virological response were only seen in patients with the unfavourable NR1I1 CCA (bAt) haplotype (OR for non-SVR 3.55; 95% CI 1.005, 12.57; P=0.049). CONCLUSIONS: NR1I1 vitamin D receptor polymorphisms influence response to pegylated-interferon/ribavirin-based therapy in chronic hepatitis C and exert an additive genetic predisposition to previously described low 25-OH vitamin D serum levels.

Crohn's disease-associated polymorphism within the PTPN2 gene affects muramyl-dipeptide-induced cytokine secretion and autophagy.

BACKGROUND: The single nucleotide polymorphism (SNP) rs2542151 within the gene locus region encoding protein tyrosine phosphatase non-receptor type 2 (PTPN2) has been associated with Crohn's disease (CD), ulcerative colitis (UC), type-I diabetes, and rheumatoid arthritis. We have previously shown that PTPN2 regulates mitogen-activated protein kinase (MAPK) signaling and cytokine secretion in human THP-1 monocytes and intestinal epithelial cells (IEC). Here, we studied whether intronic PTPN2 SNP rs1893217 regulates immune responses to the nucleotide-oligomerization domain 2 (NOD2) ligand, muramyl-dipeptide (MDP). MATERIALS AND METHODS: Genomic DNA samples from 343 CD and 663 non-IBD control patients (male and female) from a combined German, Swiss, and Polish cohort were genotyped for the presence of the PTPN2 SNPs, rs2542151, and rs1893217. PTPN2-variant rs1893217 was introduced into T(84) IEC or THP-1 cells using a lentiviral vector. RESULTS: We identified a novel association between the genetic variant, rs1893217, located in intron 7 of the PTPN2 gene and CD. Human THP-1 monocytes carrying this variant revealed increased MAPK activation as well as elevated mRNA expression of T-bet transcription factor and secretion of interferon-γ in response to the bacterial wall component, MDP. In contrast, secretion of interleukin-8 and tumor necrosis factor were reduced. In both, T(84) IEC and THP-1 monocytes, autophagosome formation was impaired. CONCLUSIONS: We identified a novel CD-associated PTPN2 variant that modulates innate immune responses to bacterial antigens. These findings not only provide key insights into the effects of a functional mutation on a clinically relevant gene, but also reveal how such a mutation could contribute to the onset of disease.