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.

Genomic evidence for role of inversion 3RP of Drosophila melanogaster in facilitating climate change adaptation.

Chromosomal inversion polymorphisms are common in animals and plants, and recent models suggest that alternative arrangements spread by capturing different combinations of alleles acting additively or epistatically to favour local adaptation. It is also thought that inversions typically maintain favoured combinations for a long time by suppressing recombination between alternative chromosomal arrangements. Here, we consider patterns of linkage disequilibrium and genetic divergence in an old inversion polymorphism in Drosophila melanogaster (In(3R)Payne) known to be associated with climate change adaptation and a recent invasion event into Australia. We extracted, karyotyped and sequenced whole chromosomes from two Australian populations, so that changes in the arrangement of the alleles between geographically separated tropical and temperate areas could be compared. Chromosome-wide linkage disequilibrium (LD) analysis revealed strong LD within the region spanned by In(3R)Payne. This genomic region also showed strong differentiation between the tropical and the temperate populations, but no differentiation between different karyotypes from the same population, after controlling for chromosomal arrangement. Patterns of differentiation across the chromosome arm and in gene ontologies were enhanced by the presence of the inversion. These data support the notion that inversions are strongly selected by bringing together combinations of genes, but it is still not clear if such combinations act additively or epistatically. Our data suggest that climatic adaptation through inversions can be dynamic, reflecting changes in the relative abundance of different forms of an inversion and ongoing evolution of allelic content within an inversion.

Genomic evidence for adaptive inversion clines in Drosophila melanogaster.

Clines in chromosomal inversion polymorphisms-presumably driven by climatic gradients-are common but there is surprisingly little evidence for selection acting on them. Here we address this long-standing issue in Drosophila melanogaster by using diagnostic single nucleotide polymorphism (SNP) markers to estimate inversion frequencies from 28 whole-genome Pool-seq samples collected from 10 populations along the North American east coast. Inversions In(3L)P, In(3R)Mo, and In(3R)Payne showed clear latitudinal clines, and for In(2L)t, In(2R)NS, and In(3R)Payne the steepness of the clinal slopes changed between summer and fall. Consistent with an effect of seasonality on inversion frequencies, we detected small but stable seasonal fluctuations of In(2R)NS and In(3R)Payne in a temperate Pennsylvanian population over 4 years. In support of spatially varying selection, we observed that the cline in In(3R)Payne has remained stable for >40 years and that the frequencies of In(2L)t and In(3R)Payne are strongly correlated with climatic factors that vary latitudinally, independent of population structure. To test whether these patterns are adaptive, we compared the amount of genetic differentiation of inversions versus neutral SNPs and found that the clines in In(2L)t and In(3R)Payne are maintained nonneutrally and independent of admixture. We also identified numerous clinal inversion-associated SNPs, many of which exhibit parallel differentiation along the Australian cline and reside in genes known to affect fitness-related traits. Together, our results provide strong evidence that inversion clines are maintained by spatially-and perhaps also temporally-varying selection. We interpret our data in light of current hypotheses about how inversions are established and maintained.

Restricted gene flow at specific parts of the shrew genome in chromosomal hybrid zones.

The species and races of the shrews of the Sorex araneus group exhibit a broad range of chromosomal polymorphisms. European taxa of this group are parapatric and form contact or hybrid zones that span an extraordinary variety of situations, ranging from absolute genetic isolation to almost free gene flow. This variety seems to depend for a large part on the chromosome composition of populations, which are primarily differentiated by various Robertsonian fusions of a subset of acrocentric chromosomes. Previous studies suggested that chromosomal rearrangements play a causative role in the speciation process. In such models, gene flow should be more restricted for markers on chromosomes involved in rearrangements than on chromosomes common in both parent species. In the present study, we address the possibility of such differential gene flow in the context of two genetically very similar but karyotypically different hybrid zones between species of the S. araneus group using microsatellite loci mapped to the chromosome arm level. Interspecific genetic structure across rearranged chromosomes was in general larger than across common chromosomes. However, the difference between the two classes of chromosomes was only significant in the hybrid zone where the complexity of hybrids is expected to be larger. These differences did not distinguish populations within species. Therefore, the rearranged chromosomes appear to affect the reproductive barrier between karyotypic species, although the strength of this effect depends on the complexity of the hybrids produced.

Effect of single nucleotide polymorphisms in cytochrome P450 isoenzyme and N-acetyltransferase 2 genes on the metabolism of artemisinin-based combination therapies in malaria patients from Cambodia and Tanzania.

The pharmacogenetics of antimalarial agents are poorly known, although the application of pharmacogenetics might be critical in optimizing treatment. This population pharmacokinetic-pharmacogenetic study aimed at assessing the effects of single nucleotide polymorphisms (SNPs) in cytochrome P450 isoenzyme genes (CYP, namely, CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, and CYP3A5) and the N-acetyltransferase 2 gene (NAT2) on the pharmacokinetics of artemisinin-based combination therapies in 150 Tanzanian patients treated with artemether-lumefantrine, 64 Cambodian patients treated with artesunate-mefloquine, and 61 Cambodian patients treated with dihydroartemisinin-piperaquine. The frequency of SNPs varied with the enzyme and the population. Higher frequencies of mutant alleles were found in Cambodians than Tanzanians for CYP2C9*3, CYP2D6*10 (100C → T), CYP3A5*3, NAT2*6, and NAT2*7. In contrast, higher frequencies of mutant alleles were found in Tanzanians for CYP2D6*17 (1023C → T and 2850C → T), CYP3A4*1B, NAT2*5, and NAT2*14. For 8 SNPs, no significant differences in frequencies were observed. In the genetic-based population pharmacokinetic analyses, none of the SNPs improved model fit. This suggests that pharmacogenetic data need not be included in appropriate first-line treatments with the current artemisinin derivatives and quinolines for uncomplicated malaria in specific populations. However, it cannot be ruled out that our results represent isolated findings, and therefore more studies in different populations, ideally with the same artemisinin-based combination therapies, are needed to evaluate the influence of pharmacogenetic factors on the clearance of antimalarials.

Evaluation of the reconstruction limits of a frequency-independent crosshole georadar waveform inversion scheme in the presence of dispersion

Waveform tomographic imaging of crosshole georadar data is a powerful method to investigate the shallow subsurface because of its ability to provide images of pertinent petrophysical parameters with extremely high spatial resolution. All current crosshole georadar waveform inversion strategies are based on the assumption of frequency-independent electromagnetic constitutive parameters. However, in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behavior. In this paper, we evaluate synthetically the reconstruction limits of a recently published crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. Our results indicate that, when combined with a source wavelet estimation procedure that provides a means of partially accounting for the frequency-dependent effects through an "effective" wavelet, the inversion algorithm performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Genome-wide patterns of latitudinal differentiation among populations of Drosophila melanogaster from North America.

Understanding the genetic underpinnings of adaptive change is a fundamental but largely unresolved problem in evolutionary biology. Drosophila melanogaster, an ancestrally tropical insect that has spread to temperate regions and become cosmopolitan, offers a powerful opportunity for identifying the molecular polymorphisms underlying clinal adaptation. Here, we use genome-wide next-generation sequencing of DNA pools ('pool-seq') from three populations collected along the North American east coast to examine patterns of latitudinal differentiation. Comparing the genomes of these populations is particularly interesting since they exhibit clinal variation in a number of important life history traits. We find extensive latitudinal differentiation, with many of the most strongly differentiated genes involved in major functional pathways such as the insulin/TOR, ecdysone, torso, EGFR, TGFβ/BMP, JAK/STAT, immunity and circadian rhythm pathways. We observe particularly strong differentiation on chromosome 3R, especially within the cosmopolitan inversion In(3R)Payne, which contains a large number of clinally varying genes. While much of the differentiation might be driven by clinal differences in the frequency of In(3R)P, we also identify genes that are likely independent of this inversion. Our results provide genome-wide evidence consistent with pervasive spatially variable selection acting on numerous loci and pathways along the well-known North American cline, with many candidates implicated in life history regulation and exhibiting parallel differentiation along the previously investigated Australian cline.

Impact of genetic polymorphisms in cytomegalovirus glycoprotein B on outcomes in solid-organ transplant recipients with cytomegalovirus disease.

BACKGROUND:It is unknown whether specific viral polymorphisms affect in vivo therapeutic response in patients with cytomegalovirus (CMV) disease. Polymorphisms in the CMV glycoprotein B (gB) gene allow discrimination of 4 distinct genotypes (gB1-gB4). We assessed the influence of gB genotypes on the clinical and virologic outcome of CMV disease. METHODS:Solid-organ transplant recipients enrolled in a multicenter trial of CMV disease treatment (VICTOR study) were included in this study. CMV gB genotyping was performed using quantitative real-time polymerase chain reaction at day 0 (start of antiviral therapy). RESULTS:Among 239 patients with CMV disease, the prevalence of gB strain types was 26% for gB1, 10% for gB2, 10% for gB3, and 5% for gB4, whereas mixed infections were present in 49%. Donor-seropositive/recipient-seropositive patients were more likely to have mixed gB infection than donor-seropositive/recipient-seronegative patients (40% vs. 12%; P = .001). Median baseline viral loads were higher and time to viral eradication was longer ( P = .006 and P = .026 , respectively) for mixed infection versus infection with a single genotype. In a multivariate model, mixed gB infection was a significant predictor of failure to eradicate virus by day 21 (mixed vs single genotype; odds ratio, 2.66; 95% confidence interval, 1.31-5.38; P = .007 ) after controlling for baseline viral load, CMV serostatus at baseline, ganciclovir resistance, and antiviral treatment. No effect of gB genotype was seen on virologic or clinical CMV recurrence. CONCLUSIONS:No specific gB genotype appears to confer a specific CMV virulence advantage. However, mixed gB genotype infections are associated with higher viral loads and delayed viral clearance.

Recurrent loss of heterozygosity in 1p36 associated with TNFRSF14 mutations in IRF4 translocation negative pediatric follicular lymphomas.

Pediatric follicular lymphoma is a rare disease that differs genetically and clinically from its adult counterpart. With the exception of pediatric follicular lymphoma with IRF4-translocation, the genetic events associated with these lymphomas have not yet been defined. We applied array-comparative genomic hybridization and molecular inversion probe assay analyses to formalin-fixed paraffin-embedded tissues from 18 patients aged 18 years and under with IRF4 translocation negative follicular lymphoma. All evaluable cases lacked t(14;18). Only 6 of 16 evaluable cases displayed chromosomal imbalances with gains or amplifications of 6pter-p24.3 (including IRF4) and deletion and copy number neutral-loss of heterozygosity in 1p36 (including TNFRSF14) being most frequent. Sequencing of TNFRSF14 located in the minimal region of loss in 1p36.32 showed nine mutations in 7 cases from our series. Two subsets of pediatric follicular lymphoma were delineated according to the presence of molecular alterations, one with genomic aberrations associated with higher grade and/or diffuse large B-cell lymphoma component and more widespread disease, and another one lacking genetic alterations associated with more limited disease.

ABCB1 and cytochrome P450 polymorphisms: clinical pharmacogenetics of clozapine.

To examine the genetic factors influencing clozapine kinetics in vivo, 75 patients treated with clozapine were genotyped for CYPs and ABCB1 polymorphisms and phenotyped for CYP1A2 and CYP3A activity. CYP1A2 activity and dose-corrected trough steady-state plasma concentrations of clozapine correlated significantly (r = -0.61; P = 1 x 10), with no influence of the CYP1A2*1F genotype (P = 0.38). CYP2C19 poor metabolizers (*2/*2 genotype) had 2.3-fold higher (P = 0.036) clozapine concentrations than the extensive metabolizers (non-*2/*2). In patients comedicated with fluvoxamine, a strong CYP1A2 inhibitor, clozapine and norclozapine concentrations correlate with CYP3A activity (r = 0.44, P = 0.075; r = 0.63, P = 0.007, respectively). Carriers of the ABCB1 3435TT genotype had a 1.6-fold higher clozapine plasma concentrations than noncarriers (P = 0.046). In conclusion, this study has shown for the first time a significant in vivo role of CYP2C19 and the P-gp transporter in the pharmacokinetics of clozapine. CYP1A2 is the main CYP isoform involved in clozapine metabolism, with CYP2C19 contributing moderately, and CYP3A4 contributing only in patients with reduced CYP1A2 activity. In addition, ABCB1, but not CYP2B6, CYP2C9, CYP2D6, CYP3A5, nor CYP3A7 polymorphisms, influence clozapine pharmacokinetics.

Interleukin-28B polymorphisms, IP-10, viral load and age predict the outcome of therapy in genotype 1 hepatitis C virus under real-life conditions

Background and Aims: IL28B polymorphisms, interferon (IFN)-gamma inducible protein-10 (IP-10) levels and the homeostasis model assessment of insulin resistance (HOMA-IR) score have been reported to predict rapid (RVR) and sustained (SVR) virological response in chronic hepatitis C (CHC), but it is not known whether these factors represent independent, clinically useful predictors. The aim of the study was to assess factors (including IL28B polymorphisms, IP-10 levels and HOMA-IR score) independently predicting response to therapy in CHC under real life conditions.Methods: Multivariate analysis of factors predicting RVR and SVR in 280 consecutive, treatment-naive CHC patients treated with pegylated IFN alpha and ribavirin in a prospective multicenter study.Results: Independent predictors of RVR were HCV RNA < 400,000 IU/ml (OR11.37; 95% CI 3.03-42.6), rs12980275 AA (vs. AG/GG) (OR 7.09; 1.97-25.56) and IP-10 (OR 0.04; 0.003-0.56) in HCV genotype 1 patients and lower baseline γ-glutamyl-transferase levels (OR = 0.02; 0.0009-0.31) in HCV genotype 3 patients. Independent predictors of SVR were rs12980275 AA (OR 9.68; 3.44-27.18), age < 40 yrs (OR = 4.79; 1.50-15.34) and HCV RNA < 400,000 IU/ml (OR 2.74; 1.03-7.27) in HCV genotype 1 patients and rs12980275 AA (OR = 6.26; 1.98-19.74) and age < 40 yrs (OR 5.37; 1.54-18.75) in the 88 HCV genotype 1 patients without a RVR. RVR was by itself predictive of SVR in HCV genotype 1 patients (32 of 33, 97%; OR 33.0; 4.06-268.32) and the only independent predictor of SVR in HCV genotype 2 (OR 9.0, 1.72-46.99; p=0.009) or 3 patients (OR 7.8, 1.43-42.67; p=0.01).Conclusions: In HCV genotype 1 patients, IL28B polymorphisms, HCV RNA load and IP-10 independently predict RVR. The combination of IL28B polymorphisms, HCV RNA level and age may yield more accurate pretreatment prediction of SVR. HOMA-IR score is not associated with viral response.

Selected polymorphisms in sex hormone-related genes, circulating sex hormones and risk of endometrial cancer.

BACKGROUND: The role of estrogen and progesterone in the development of endometrial cancer is well documented. Few studies have examined the association of genetic variants in sex hormone-related genes with endometrial cancer risk. METHODS: We conducted a case-control study nested within three cohorts to examine the association of endometrial cancer risk with polymorphisms in hormone-related genes among 391 cases (92% postmenopausal at diagnosis) and 712 individually-matched controls. We also examined the association of these polymorphisms with circulating levels of sex hormones and SHBG in a cross-sectional analysis including 596 healthy postmenopausal women at blood donation (controls from this nested case-control study and from a nested case-control study of breast cancer in one of the three cohorts). RESULTS: Adjusting for endometrial cancer risk factors, the A allele of rs4775936 in CYP19 was significantly associated (OR(per allele)=1.22, 95% CI=1.01-1.47, p(trend)=0.04), while the T allele of rs10046 was marginally associated with increased risk of endometrial cancer (OR(per allele)=1.20, 95% CI=0.99-1.45, p(trend)=0.06). PGR rs1042838 was also marginally associated with risk (OR(per allele)=1.25, 95% CI=0.96-1.61, p(trend)=0.09). No significant association was found for the other polymorphisms, i.e. CYP1B1 rs1800440 and rs1056836, UGT1A1 rs8175347, SHBG rs6259 and ESR1 rs2234693. Rs8175347 was significantly associated with postmenopausal levels of estradiol, free estradiol and estrone and rs6259 with SHBG and estradiol. CONCLUSION: Our findings support an association between genetic variants in CYP19, and possibly PGR, and risk of endometrial cancer.