Genotype Imputation Reference Panel Selection Using Maximal Phylogenetic Diversity

The recent dramatic cost reduction of next-generation sequencing technology enables investigators to assess most variants in the human genome to identify risk variants for complex diseases. However, sequencing large samples remains very expensive. For a study sample with existing genotype data, such as array data from genome-wide association studies, a cost-effective approach is to sequence a subset of the study sample and then to impute the rest of the study sample, using the sequenced subset as a reference panel. The use of such an internal reference panel identifies population-specific variants and avoids the problem of a substantial mismatch in ancestry background between the study population and the reference population. To efficiently select an internal panel, we introduce an idea of phylogenetic diversity from mathematical phylogenetics and comparative genomics. We propose the “most diverse reference panel”, defined as the subset with the maximal “phylogenetic diversity”, thereby incorporating individuals that span a diverse range of genotypes within the sample. Using data both from simulations and from the 1000 Genomes Project, we show that the most diverse reference panel can substantially improve the imputation accuracy compared to randomly selected reference panels, especially for the imputation of rare variants. The improvement in imputation accuracy holds across different marker densities, reference panel sizes, and lengths for the imputed segments. We thus propose a novel strategy for planning sequencing studies on samples with existing genotype data.     

Monthly stem increment in relation to climatic variables during 7 years in an East African rainforest

Monthly stem increment of 766 trees was assessed for 7 years in Kakamega Forest, Kenya and related to monthly climatic variables. Mean stem increment of all tree individuals correlated negatively with maximum temperature but not with mean and minimum temperatures. For the precipitation variables sum of precipitation and number of rainy days we found positive correlations. Stem increment of the trees in the early-, mid-, and late-successional groups correlated positively with the number of rainy days. For late-successional trees increment correlated negatively with mean and maximum temperature and positively with all other precipitation variables. For mid-successional trees we found a negative correlation with mean temperature. In addition, the stem increment of most species related positively to precipitation variables and negatively to mean and maximum temperature. In view of the expected increasing temperatures and fewer but heavier rain events, our results suggest that climate change will lead to a reduction in stem increment. The results appertaining to the successional groups imply that early and mid-successional species are better equipped to perform well under the expected future climatic conditions than the late-successional species. This could reduce the role of this East African forest as a carbon store. As the responses to climatic variables were highly group- and species-specific it is likely that climate change will result in a species composition shift, presumably in favour of drought-resistant and heat-tolerating species.     

Intraspecific differentiation of Paramecium novaurelia strains (Ciliophora,Protozoa) inferred from phylogenetic analysis of ribosomal and mitochondrial DNA variation

Paramecium novaurelia Beale and Schneller, 1954, was first found in Scotland and is known to occur mainly in Europe, where it is the most common species of the P. aurelia complex. In recent years, two non-European localities have been described: Turkey and the United States of America. This article presents the analysis of intraspecific variability among 25 strains of P. novaurelia with the application of ribosomal and mitochondrial loci (ITS1-5.8S-ITS2, 5′ large subunit rDNA (5′LSU rDNA) and cytochrome c oxidase subunit 1 (COI) mtDNA). The mean distance observed for all of the studied P. novaurelia sequence pairs was p = 0.008/0.016/0.092 (ITS1-5.8S-ITS2/5′LSU rDNA/COI). Phylogenetic trees (NJ/MP/BI) based on a comparison of all of the analysed sequences show that the studied strains of P. novaurelia form a distinct clade, separate from the P. caudatum outgroup, and are divided into two clusters (A and B) and two branches (C and D). The occurrence of substantial genetic differentiation within P. novaurelia, confirmed by the analysed DNA fragments, indicates a rapid evolution of particular species within the Paramecium genus.     

Multiple nuclear and mitochondrial DNA sequences provide new insights into the phylogeny of South African Lacertids (Lacertidae,Eremiadinae)

Eremiadinae, one of three subfamilies of Lacertidae, are distributed throughout Asia and Africa. Previous phylogenetic studies suggested that one of the main groups of Eremiadinae (the Ethiopian clade) consist of two clades with predominately East‐African and South‐African distribution. Yet, especially the latter one, which includes the genera Pedioplanis, Meroles, Ichnotropis, Tropidosaura and Australolacerta, was not well supported in the molecular phylogenetic analysis. In this study, we analysed the phylogenetic relationships among the genera of the ‘South African clade’ to assess whether this group actually forms a highly supported clade and to address questions concerning the monophyly of the genera. We sequenced sections of the widely used mitochondrial genes coding for 16S rRNA, 12S rRNA and cytochrome b (altogether 2045 bp) as well as the nuclear genes c‐mos, RAG‐1, PRLR, KIF24, EXPH5 and RAG‐2 (altogether 4473 bp). The combined data set increased the support values for several nodes considerably. Yet, the relationships among five major lineages within the ‘South African clade’ are not clearly resolved even with this large data set. We interpret this as a ‘hard polytomy’ due to fast radiation within the South African lacertids. The combined tree based on nine marker genes provides strong support for the ‘South African Clade’ and its sister group relationship with the ‘East African Clade’. Our results confirm the genus Tropidosaura as a monophylum, while Ichnotropis is paraphyletic in our trees: Ichnotropis squamulosa appears more closely related to Meroles than to Ichnotropis capensis. Furthermore, the monophyly of Meroles is questionable as well. Based on our results, I. squamulosa should be transferred from Ichnotropis into the genus Meroles. Also, the two species of Australolacerta (A. australis and A. rupicola) are very distantly related and the genus is perhaps paraphyletic, too. Finally we propose a phylogeographical scenario in the context of palaeoclimatic data and compare it with a previously postulated hypothesis.     

Diversity and selective pressures of anticoagulants in three medicinal leeches (Hirudinida: Hirudinidae,Macrobdellidae)

Although medicinal leeches have long been used as treatment for various ailments because of their potent anticoagulation factors, neither the full diversity of salivary components that inhibit coagulation, nor the evolutionary selection acting on them has been thoroughly investigated. Here, we constructed expressed sequence tag libraries from salivary glands of two species of medicinal hirudinoid leeches, Hirudo verbana and Aliolimnatis fenestrata, and identified anticoagulant‐orthologs through BLASTx searches. The data set then was augmented by the addition of a previously constructed EST library from the macrobdelloid leech Macrobdella decora. The identified orthologs then were compared and contrasted with well‐characterized anticoagulants from a variety of leeches with different feeding habits, including non‐sanguivorous species. Moreover, four different statistical methods for predicting signatures of positive and negative evolutionary pressures were used for 10 rounds each to assess the level and type of selection acting on the molecules as a whole and on specific sites. In total, sequences showing putative BLASTx‐orthology with five and three anticoagulant‐families were recovered in the A. fenestrata and H. verbana EST libraries respectively. Selection pressure analyses predicted high levels of purifying selection across the anticoagulant diversity, although a few isolated sites showed signatures of positive selection. This study represents a first attempt at mapping the anticoagulant repertoires in a comparative fashion across several leech families.     

Importance of rare taxa for bacterial diversity in the rhizosphere of Bt- and conventional maize varieties

Ribosomal 16S rRNA gene pyrosequencing was used to explore whether the genetically modified (GM) Bt-maize hybrid MON 89034 × MON 88017, expressing three insecticidal recombinant Cry proteins of Bacillus thuringiensis, would alter the rhizosphere bacterial community. Fine roots of field cultivated Bt-maize and three conventional maize varieties were analyzed together with coarse roots of the Bt-maize. A total of 547 000 sequences were obtained. Library coverage was 100% at the phylum and 99.8% at the genus rank. Although cluster analyses based on relative abundances indicated no differences at higher taxonomic ranks, genera abundances pointed to variety specific differences. Genera-based clustering depended solely on the 49 most dominant genera while the remaining 461 rare genera followed a different selection. A total of 91 genera responded significantly to the different root environments. As a benefit of pyrosequencing, 79 responsive genera were identified that might have been overlooked with conventional cloning sequencing approaches owing to their rareness. There was no indication of bacterial alterations in the rhizosphere of the Bt-maize beyond differences found between conventional varieties. B. thuringiensis-like phylotypes were present at low abundance (0.1% of Bacteria) suggesting possible occurrence of natural Cry proteins in the rhizospheres. Although some genera indicated potential phytopathogenic bacteria in the rhizosphere, their abundances were not significantly different between conventional varieties and Bt-maize. With an unprecedented sensitivity this study indicates that the rhizosphere bacterial community of a GM maize did not respond abnormally to the presence of three insecticidal proteins in the root tissue.