Structural and Functional Diversity of Peptide Toxins from Tarantula Haplopelma hainanum (Ornithoctonus hainana) Venom Revealed by Transcriptomic, Peptidomic, and Patch Clamp Approaches

Spider venom is a complex mixture of bioactive peptides to subdue their prey. Early estimates suggested that over 400 venom peptides are produced per species. In order to investigate the mechanisms responsible for this impressive diversity, transcriptomics based on second generation high throughput sequencing was combined with peptidomic assays to characterize the venom of the tarantula Haplopelma hainanum. The genes expressed in the venom glands were identified, and the bioactivity of their protein products was analyzed using the patch clamp technique. A total of 1,136 potential toxin precursors were identified that clustered into 90 toxin groups, of which 72 were novel. The toxin peptides clustered into 20 cysteine scaffolds that included between 4 and 12 cysteines, and 14 of these groups were newly identified in this spider. Highly abundant toxin peptide transcripts were present and resulted from hypermutation and/or fragment insertion/deletion. In combination with variable post-translational modifications, this genetic variability explained how a limited set of genes can generate hundreds of toxin peptides in venom glands. Furthermore, the intraspecies venom variability illustrated the dynamic nature of spider venom and revealed how complex components work together to generate diverse bioactivities that facilitate adaptation to changing environments, types of prey, and milking regimes in captivity.

"The usual suspects"- analysis of transcriptome sequences reveals deviating B gene activity in C. vulgaris bud bloomers

BACKGROUND

The production of heather (Calluna vulgaris) in Germany is highly dependent on cultivars with mutated flower morphology, the so-called diplocalyx bud bloomers. So far, this unique flower type of C. vulgaris has not been reported in any other plant species. The flowers are characterised by an extremely extended flower attractiveness, since the flower buds remain closed throughout the complete flowering season. The flowers of C. vulgaris bud bloomers are male sterile, because the stamens are absent. Furthermore, petals are converted into sepals. Therefore the diplocalyx bud bloomer flowers consist of two whorls of sepals directly followed by the gynoecium.

RESULTS

A broad comparison was undertaken to identify genes differentially expressed in the bud flowering phenotype and in the wild type of C. vulgaris. Transcriptome sequence reads were generated using 454 sequencing of two flower type specific cDNA libraries. In total, 360,000 sequence reads were obtained, assembled to 12,200 contigs, functionally mapped, and annotated. Transcript abundances were compared and 365 differentially expressed genes detected. Among these differentially expressed genes, Calluna vulgaris PISTILLATA (CvPI) which is the orthologue of the Arabidopsis B gene PISTILLATA (PI) was considered as the most promising candidate gene. Quantitative Reverse Transcription Polymerase Chain Reaction (qRT PCR) was performed to analyse the gene expression levels of two C. vulgaris B genes CvPI and Calluna vulgaris APETALA 3 (CvAP3) in both flower types. CvAP3 which is the orthologue of the Arabidopsis B gene APETALA 3 (AP3) turned out to be ectopically expressed in sepals of wild type and bud bloomer flowers. CvPI expression was proven to be reduced in the bud blooming flowers.

CONCLUSIONS

Differential expression patterns of the B-class genes CvAP3 and CvPI were identified to cause the characteristic morphology of C. vulgaris flowers leading to the following hypotheses: ectopic expression of CvAP3 is a convincing explanation for the formation of a completely petaloid perianth in both flower types. In C. vulgaris, CvPI is essential for determination of petal and stamen identity. The characteristic transition of petals into sepals potentially depends on the observed deficiency of CvPI and CvAP3 expression in bud blooming flowers.

Bacterial community structure in freshwater springs infested with the invasive plant species Hydrilla verticillata

The phylogenetic composition and physiological profiles of bacterial communities in freshwater springs were evaluated during the blooming and non-blooming stages of the invasive plant species, Hydrilla verticillata. Community-level physiological profiles (CLPPs) and pyrosequencing of 16S rRNA gene amplicons were used to study potential Hydrilla mediated shifts in the physiological potential and phylogenetic composition of the bacterial community in infested systems. The results of CLPP revealed that the microbes in the Hydrilla invaded sites utilized less substrates during blooming periods than during nonblooming periods of the plant. Spearman's rank correlation analysis showed some relationships between the relative abundances of bacterial taxa and the Biolog substrate utilization pattern. The relative abundance of the identified taxa showed some striking differences based on the blooming status of Hydrilla and to a lesser extent on site variation. The relative abundance of Actinobacteria, Bacteriodetes, and Verrucomicrobia was generally higher during Hydrilla blooms, while Deltaproteobacteria was generally higher during non-blooming stages of Hydrilla. The detected genera also varied based on the blooming stages of the plant. Based on the findings, it appears that Hydrilla alters the phylogenetic composition and structure of the bacterial community during the blooming stage.

Transcriptome and proteome profiling of adventitious root development in hybrid larch (Larix kaempferi × Larix olgensis)

BACKGROUND

Hybrids of larch (Larix kaempferi × Larix olgensis) are important afforestation species in northeastern China. They are routinely propagated via rooted stem cuttings. Despite the importance of rooting, little is known about the regulation of adventitious root development in larch hybrids. 454 GS FLX Titanium technology represents a new method for characterizing the transcriptomes of non-model species. This method can be used to identify differentially expressed genes, and then two-dimensional difference gel electrophoresis (2D-DIGE) and matrix-assisted laser desorption-ionization time-of-flight mass spectrometry (MALDI-TOF/TOF MS) analyses can be used to analyze their corresponding proteins. In this study, we analyzed semi-lignified cuttings of two clones of L. kaempferi × L. olgensis with different rooting capacities to study the molecular basis of adventitious root development.

RESULTS

We analyzed two clones; clone 25-5, with strong rooting capacity, and clone 23-12, with weak rooting capacity. We constructed four cDNA libraries from 25-5 and 23-12 at two development stages. Sequencing was conducted using the 454 pyrosequencing platform. A total of 957832 raw reads was produced; 95.07% were high-quality reads, and were assembled into 45137 contigs and 61647 singletons. The functions of the unigenes, as indicated by their Gene Ontology annotation, included diverse roles in the molecular functions, biological processes, and cellular component categories. We analyzed 75 protein spots (-fold change ≥ 2, P ≤ 0.05) by 2D-DIGE, and identified the differentially expressed proteins using MALDI-TOF/TOF MS. A joint analysis of transcriptome and proteome showed genes related to two pathways, polyamine synthesis and stress response, might play an important role on adventitious root development.

CONCLUSIONS

These results provide fundamental and important information for research on the molecular mechanism of adventitious root development. We also demonstrated for the first time the combined use of two important technologies as a powerful approach to advance research on non-model, but otherwise important, larch species.

Isolation and characterization of 16 polymorphic microsatellite loci for Euchiloglanis kishinouyei

A total of 16 novel polymorphic microsatellite marker loci were isolated from a genomic library of Euchiloglanis kishinouyei and further characterized using a sample from a wild population consisting of 40 individuals. The number of alleles among loci ranged from three to 17 and the observed and expected heterozygosities ranged from 0·030 to 0·950 and 0·163 to 0·980, respectively. The average polymorphic information content (PIC) of all loci was 0·450.

The Antibody Germline/Maturation Hypothesis, Elicitation of Broadly Neutralizing Antibodies Against HIV-1 and Cord Blood IgM Repertoires

We have previously observed that all known potent broadly neutralizing antibodies (bnAbs) against HIV-1 are highly divergent from their putative germline predecessors in contrast to bnAbs against viruses causing acute infections such as henipaviruses and SARS CoV, which are much less divergent from their germline counterparts. Consequently, we have hypothesized that germline antibodies may not bind to the HIV-1 envelope glycoprotein (Env) because they are so different compared to the highly somatically mutated HIV-1-specific bnAbs. We have further hypothesized that the immunogenicity of highly conserved epitopes on the HIV-1 envelope glycoproteins (Envs) may be reduced or eliminated by their very weak or absent interactions with germline antibodies and immune responses leading to the elicitation of bnAbs may not be initiated and/or sustained. Even if such responses are initiated, the maturation pathways are so extraordinarily complex that prolonged periods of time may be required for elicitation of bnAbs with defined unique sequences. We provided the initial evidence supporting this antibody germline/maturation hypothesis, which prompted a number of studies to design vaccine immunogens that could bind putative germline predecessors of known bnAbs and to explore complex B cell lineages. However, guiding the immune system through the exceptionally complex antibody maturation pathways to elicit known bnAbs remains a major challenge. Here, we discuss studies exploring the antibody germline/maturation hypothesis as related to elicitation of bnAbs against HIV-1 and present our recent data demonstrating the existence of germline-like precursors of VRC01 antibodies in a human cord blood IgM library.

Bidirectional genomic exchange between Helicobacter pylori strains from a family in Coventry, United Kingdom

The human gastric pathogen Helicobacter pylori is characterised by a high mutation rate and frequent recombination during mixed infection, which result in extensive genetic diversity and rapid allelic diversification. Mixed infections are believed to be much more common in regions with a high H. pylori prevalence than in industrialised countries. To better understand the genomic flexibility of H. pylori in a low prevalence region, we used 454 sequencing technology to investigate whole genome sequences of H. pylori strains isolated from members of three generations of a family living in Coventry, UK. The genomes of four H. pylori strains isolated from a grandfather, two of his sons and one grandson were sequenced. Three of these genomes showed a high overall sequence similarity, suggesting a recent common ancestor. The genomes differed by 316-336 SNPs, and recombination events (imports) resulted in 170-251 clusters of polymorphisms (CNPs). Imports were particularly frequent in genes encoding Helicobacter outer membrane proteins, suggesting an adaptation of the strains to their individual host. The fourth strain differed substantially from these three highly related strains but still shared long fragments of identical sequence, which most likely reflect imports from the highly related family variants. The data show extensive bidirectional exchange of DNA between the strains isolated from the family members, illustrating both the convergence and divergence effect that recombination can lead to. Detailed analysis of the distribution of SNPs and imports permits to draw up a complex scenario of the transmission history involving infection with at least two, and probably more separate strains. This complexity and the resulting high frequency of recombination were unexpected for an industrialised country where the prevalence of H. pylori infection has strongly declined in recent decades.

Species richness of arbuscular mycorrhizal fungi: associations with grassland plant richness and biomass

Although experiments show a positive association between vascular plant and arbuscular mycorrhizal fungal (AMF) species richness, evidence from natural ecosystems is scarce. Furthermore, there is little knowledge about how AMF richness varies with belowground plant richness and biomass. We examined relationships among AMF richness, above- and belowground plant richness, and plant root and shoot biomass in a native North American grassland. Root-colonizing AMF richness and belowground plant richness were detected from the same bulk root samples by 454-sequencing of the AMF SSU rRNA and plant trnL genes. In total we detected 63 AMF taxa. Plant richness was 1.5 times greater belowground than aboveground. AMF richness was significantly positively correlated with plant species richness, and more strongly with below- than aboveground plant richness. Belowground plant richness was positively correlated with belowground plant biomass and total plant biomass, whereas aboveground plant richness was positively correlated only with belowground plant biomass. By contrast, AMF richness was negatively correlated with belowground and total plant biomass. Our results indicate that AMF richness and plant belowground richness are more strongly related with each other and with plant community biomass than with the plant aboveground richness measures that have been almost exclusively considered to date.

The complete mitochondrial genome of the sandbar shark Carcharhinus plumbeus

The sandbar shark, Carcharhinus plumbeus, a major representative species in shark fisheries worldwide is now considered vulnerable to overfishing. A pool of 774,234 Roche 454 shotgun sequences from one individual were assembled into a 16,706 bp mitogenome with 33× average coverage depth. It comprised 13 protein coding genes, 22 transfer RNA's, 2 ribosomal genes and 2 non-coding regions, typical of a vertebrate mitogenome. As expected for sharks, an A-T nucleotide bias was evident. This adds to rapidly growing number of mitogenome assemblies for the economically important Carcharhinidae family. The C. plumbeus mitogenome will assist researchers, fisheries and conservation managers interested in shark molecular systematics, phylogeography, conservation genetics, population and stock structure.

Characterization of the reniform nematode genome by shotgun sequencing

The reniform nematode (RN), a major agricultural pest particularly on cotton in the United States, is among the major plant-parasitic nematodes for which limited genomic information exists. In this study, over 380 Mb of sequence data were generated from pooled DNA of four adult female RNs and assembled into 67,317 contigs, including 25,904 (38.5%) predicted coding contigs and 41,413 (61.5%) noncoding contigs. Most of the characterized repeats were of low complexity (88.9%), and 0.9% of the contigs matched with 53.2% of GenBank ESTs. The most frequent Gene Ontology (GO) terms for molecular function and biological process were protein binding (32%) and embryonic development (20%). Further analysis showed that 741 (1.1%), 94 (0.1%), and 169 (0.25%) RN genomic contigs matched with 1328 (13.9%), 1480 (5.4%), and 1330 (7.4%) supercontigs of Meloidogyne incognita, Brugia malayi, and Pristionchus pacificus, respectively. Chromosome 5 of Caenorhabditis elegans had the highest number of hits to the RN contigs. Seven putative detoxification genes and three carbohydrate-active enzymes (CAZymes) involved in cell wall degradation were studied in more detail. Additionally, kinases, G protein-coupled receptors, and neuropeptides functioning in physiological, developmental, and regulatory processes were identified in the RN genome.

Development of microsatellite markers for the coastal shrub Scaevola taccada (Goodeniaceae)

PREMISE OF THE STUDY

Microsatellite markers were developed for the coastal shrub species Scaevola taccada to estimate the population genetic structure, which may reflect different seed dispersal patterns. •

METHODS AND RESULTS

Thirteen microsatellite primer sets were developed for S. taccada using 454 pyrosequencing. The primer sets were tested on 64 individuals sampled from two populations in Japan. Fragments were amplified using the primers, with one to 10 alleles per locus, and the expected heterozygosity ranged from 0.00 to 0.85. •

CONCLUSIONS

These results indicate the utility of markers in S. taccada for broad estimations of the population genetic structure of this species.

Massively parallel pyrosequencing of the mitochondrial genome with the 454 methodology in forensic genetics

RESULTS

of sequencing of whole mitochondrial genome, HV1 and HV2 DNA with the second generation system (SGS) Roche 454 GS Junior were compared with results of Sanger sequencing and SNP typing with SNaPshot single base extension detected with MALDI-TOF and capillary electrophoresis. We investigated the performance of the software analysis of the data, reproducibility, ability to sequence homopolymeric regions, detection of mixtures and heteroplasmy as well as the implications of the depth of coverage. We found full reproducibility between samples sequenced twice with SGS. We found close to full concordance between the mtDNA sequences of 26 samples obtained with (1) the 454 SGS method using a depth of coverage above 100 and (2) Sanger sequencing and SNP typing. The discrepancies were primarily observed in homopolymeric regions. The 454 SGS method was able to sequence 95% of the reads correctly in homopolymers up to 4 bases, and up to 6 bases could be sequenced with similar success if the results were carefully, visually inspected. The 454 technology was able to detect mixtures or heteroplasmy of approximately 10%. We detected previously unreported heteroplasmy in the GM9947A component of the NIST human mitochondrial DNA SRM-2392 standard reference material.

Complete genome sequence of the Radiation-Resistant bacterium Rubrobacter radiotolerans RSPS-4

Rubrobacter radiotolerans strain RSPS-4 is a slightly thermophilic member of the phylum “Actinobacteria” isolated from a hot spring in São Pedro do Sul, Portugal. This aerobic and halotolerant bacterium is also extremely resistant to gamma and UV radiation, which are the main reasons for the interest in sequencing its genome. Here, we present the complete genome sequence of strain RSPS-4 as well as its assembly and annotation. We also compare the gene sequence of this organism with that of the type strain of the species R. radiotolerans isolated from a hot spring in Japan. The genome of strain RSPS-4 comprises one circular chromosome of 2,875,491 bp with a G+C content of 66.91%, and 3 circular plasmids of 190,889 bp, 149,806 bp and 51,047 bp, harboring 3,214 predicted protein coding genes, 46 tRNA genes and a single rRNA operon.

Peroxidase gene discovery from the horseradish transcriptome

BACKGROUND

Horseradish peroxidases (HRPs) from Armoracia rusticana have long been utilized as reporters in various diagnostic assays and histochemical stainings. Regardless of their increasing importance in the field of life sciences and suggested uses in medical applications, chemical synthesis and other industrial applications, the HRP isoenzymes, their substrate specificities and enzymatic properties are poorly characterized. Due to lacking sequence information of natural isoenzymes and the low levels of HRP expression in heterologous hosts, commercially available HRP is still extracted as a mixture of isoenzymes from the roots of A. rusticana.

RESULTS

In this study, a normalized, size-selected A. rusticana transcriptome library was sequenced using 454 Titanium technology. The resulting reads were assembled into 14871 isotigs with an average length of 1133 bp. Sequence databases, ORF finding and ORF characterization were utilized to identify peroxidase genes from the 14871 isotigs generated by de novo assembly. The sequences were manually reviewed and verified with Sanger sequencing of PCR amplified genomic fragments, resulting in the discovery of 28 secretory peroxidases, 23 of them previously unknown. A total of 22 isoenzymes including allelic variants were successfully expressed in Pichia pastoris and showed peroxidase activity with at least one of the substrates tested, thus enabling their development into commercial pure isoenzymes.

CONCLUSIONS

This study demonstrates that transcriptome sequencing combined with sequence motif search is a powerful concept for the discovery and quick supply of new enzymes and isoenzymes from any plant or other eukaryotic organisms. Identification and manual verification of the sequences of 28 HRP isoenzymes do not only contribute a set of peroxidases for industrial, biological and biomedical applications, but also provide valuable information on the reliability of the approach in identifying and characterizing a large group of isoenzymes.

Development and characterization of 10 microsatellite markers in Sagina nodosa (Caryophyllaceae)

PREMISE OF THE STUDY

We developed 10 novel microsatellite loci for Sagina nodosa, a diploid perennial arctic-alpine herb. To our knowledge, these are the first microsatellite loci for a Sagina species. •

METHODS AND RESULTS

We performed a low-coverage 454 next-generation sequencing of enriched genomic fragments derived from one individual to generate a massive library of contigs containing potential polymorphic microsatellites. We present data for 10 novel polymorphic microsatellite loci containing di-, tri-, tetra-, and hexanucleotide repeats with two to nine alleles per locus assessed in 29 individuals. •

CONCLUSIONS

These polymorphic microsatellite loci in S. nodosa will provide insights on the population structure and life history of S. nodosa in Isle Royale and other North American populations.

Mathematical modeling of primary succession of murine intestinal microbiota

Understanding the nature of interpopulation interactions in host-associated microbial communities is critical to understanding gut colonization, responses to perturbations, and transitions between health and disease. Characterizing these interactions is complicated by the complexity of these communities and the observation that even if populations can be cultured, their in vitro and in vivo phenotypes differ significantly. Dynamic models are the cornerstone of computational systems biology and a key objective of computational systems biologists is the reconstruction of biological networks (i.e., network inference) from high-throughput data. When such computational models reflect biology, they provide an opportunity to generate testable hypotheses as well as to perform experiments that are impractical or not feasible in vivo or in vitro. We modeled time-series data for murine microbial communities using statistical approaches and systems of ordinary differential equations. To obtain the dense time-series data, we sequenced the 16S ribosomal RNA (rRNA) gene from DNA isolated from the fecal material of germfree mice colonized with cecal contents of conventionally raised animals. The modeling results suggested a lack of mutualistic interactions within the community. Among the members of the Bacteroidetes, there was evidence for closely related pairs of populations to exhibit parasitic interactions. Among the Firmicutes, the interactions were all competitive. These results suggest future animal and in silico experiments. Our modeling approach can be applied to other systems to provide a greater understanding of the dynamics of communities associated with health and disease.

Deep sequencing and Circos analyses of antibody libraries reveal antigen-driven selection of Ig VH genes during HIV-1 infection

The vast diversity of antibody repertoires is largely attributed to heavy chain (V(H)) recombination of variable (V), diversity (D) and joining (J) gene segments. We used 454 sequencing information of the variable domains of the antibody heavy chain repertoires from neonates, normal adults and an HIV-1-infected individual, to analyze, with Circos software, the VDJ pairing patterns at birth, adulthood and a time-dependent response to HIV-1 infection. Our comparative analyses of the Ig VDJ repertoires from these libraries indicated that, from birth to adulthood, VDJ recombination patterns remain the same with some slight changes, whereas some V(H) families are selected and preferentially expressed after long-term infection with HIV-1. We also demonstrated that the immune system responds to HIV-1 chronic infection by selectively expanding certain HV families in an attempt to combat infection. Our findings may have implications for understanding immune responses in pathology as well as for development of new therapeutics and vaccines.

Colorado potato beetle (Coleoptera) gut transcriptome analysis: expression of RNA interference-related genes

In the search for new methods of pest control, the potential of RNA interference (RNAi) is being explored. Because the gut is the first barrier for the uptake of double-stranded (ds)RNA, pyrosequencing of the gut transcriptome is a powerful tool for obtaining the necessary sequences for specific dsRNA-mediated pest control. In the present study, a dataset representing the gut transcriptome of the Colorado potato beetle (CPB; Leptinotarsa decemlineata) was generated and analysed for the presence of RNAi-related genes. Almost all selected genes that were implicated in silencing efficiency at different levels in the RNAi pathway (core machinery, associated intracellular factors, dsRNA uptake, antiviral RNAi, nucleases), which uses different types of small RNA (small interfering RNA, microRNA and piwi-RNA), were expressed in the CPB gut. Although the database is of lower quality, the majority of the RNAi genes are also found to be present in the gut transcriptome of the tobacco hornworm [TH; Manduca sexta (19 out of 35 genes analysed)]. The high quality of the CPB transcriptome database will lay the foundation for future gene expression and functional studies regarding the gut and RNAi.

Characterization of microsatellite markers for the vulnerable grassland forb Senecio macrocarpus (Asteraceae)

PREMISE OF THE STUDY

Development of microsatellite markers for the vulnerable forb Senecio macrocarpus was performed to begin an assessment of its population structure and breeding method to aid in the conservation of the species in Victoria, Australia. •

METHODS AND RESULTS

Fifteen microsatellite markers were developed for S. macrocarpus from 454 pyrosequencing. The markers were tested on 104 individuals from four populations. The markers produced between two and seven alleles per locus while the expected heterozygosity ranged from 0.20 to 0.67 and the observed heterozygosity ranged from 0.00 to 1.00. The observed heterozygosity is suggestive that the populations may be apomictic. •

CONCLUSIONS

The microsatellite markers developed for S. macrocarpus are intended to be used on future studies that aim to assess the population genetics and local breeding dynamics of the species with an emphasis on conservation.

Development and characterization of microsatellite primers in Geranium carolinianum (Geraniaceae) with 454 sequencing

PREMISE OF THE STUDY

Microsatellite primers were developed for Geranium carolinanum, a North American winter annual herb, for use in population genetic analyses. •

METHODS AND RESULTS

Genomic DNA enriched for repeat-containing fragments was sequenced on the Roche 454 Titanium platform, resulting in 470 primer pairs developed from 1115 microsatellite-containing sequences. A subset of 37 primer pairs was screened for polymorphism across three native and three invasive populations. We identified four monomorphic and eight polymorphic loci. Polymorphic loci contained between two and seven alleles per locus, and mean within-population expected heterozygosity ranged from 0.100 to 0.290. Within populations, observed heterozygosity for individual loci ranged from zero to 0.857, and expected heterozygosity ranged from 0.046 to 0.559. •

CONCLUSIONS

These microsatellite markers will be useful for future studies of genetic diversity, structure, and mating systems across the geographic range of G. carolinianum, and may be transferable to other closely related species.

High-throughput microsatellite marker development for the distylous herb Primula mistassinica (Primulaceae)

PREMISE OF THE STUDY

Twelve microsatellite markers were developed for Primula mistassinica, a distylous, diploid arctic-alpine plant. The markers will be used to investigate the landscape genetics of a disjunct population on Isle Royale, Michigan, and the phylogeographic patterns of the species. •

METHODS AND RESULTS

We used Roche/454 high-throughput technology to sequence microsatellite-enriched regions in the P. mistassinica genome. We developed 12 polymorphic microsatellite primer sets. These loci contained di-, tri-, and tetranucleotide repeats with two to nine alleles per locus when assessed in 23 individuals. •

CONCLUSIONS

Understanding the historical movements of P. mistassinica will provide insight to the survival prospects of current Arctic plant populations, which face the pressures of global, anthropogenic climate change.

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host-pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1'665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen-mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host-pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context. Using next-generation sequencing, the present manuscript identifies the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of a cyprinid fish: the European minnow (Phoxinus phoxinus). We highlight that the relative importance of neutral versus adaptive processes in shaping immune competence may differ between duplicates as a consequence of alternative selective regimes or different genomic contexts.

Isolation and characterization of microsatellite loci in the invasive herb Solidago altissima (Asteraceae)

PREMISE OF THE STUDY

Microsatellite markers were developed in the perennial herb Solidago altissima from populations within its introduced range in Japan to assess its population structure and to facilitate tracking of invasion expansion. •

METHODS AND RESULTS

Using 454 pyrosequencing, 16 microsatellite primer sets were developed for S. altissima. The primer sets were tested on 70 individuals sampled from three populations in Japan. The primers amplified di- and trinucleotide repeats with five to 25 alleles per locus, and the expected heterozygosity ranged from 0.46 to 0.92. •

CONCLUSIONS

These results indicate the utility of primers in S. altissima for future research on a wide range of applications, including tracking of invasion dynamics and investigating population genetics of the species.

Development and molecular characterization of 55 novel polymorphic cDNA-SSR markers in faba bean (Vicia faba L.) using 454 pyrosequencing

Faba bean (Vicia faba L.) is a major food source and fodder legume, popularly known for its high content of seed-protein. Its role is critical in crop rotation, and for fixing nitrogen effectively. Polymorphic simple sequence repeat markers from transcript sequences (cDNA; simple sequence repeat [SSR]) were developed for faba bean (Vicia faba). We found that 1,729 SSR loci from 81,333 individual sequence reads and 240 primer pairs were designed and synthesized. In total, 55 primer pairs were found to be polymorphic and scorable consistently when screened in 32 accessions. The number of alleles ranged from 2 to 15, frequency of major alleles per locus varied from 0.17 to 0.91, the genotypes number ranged from 2 to 17, observed and expected heterozycosity values ranged from 0.00 to 0.44 and 0.17 to 0.89 and overall PIC values ranged from 0.16 to 0.88 respectively. These markers will be a useful tool for assessing the genetic diversity, understanding the population structure, and breeding patterns of faba bean.

Comparative analysis of genome composition in Triticeae reveals strong variation in transposable element dynamics and nucleotide diversity

A 454 sequencing snapshot was utilised to investigate the genome composition and nucleotide diversity of transposable elements (TEs) for several Triticeae taxa, including Triticum aestivum, Hordeum vulgare, Hordeum spontaneum and Secale cereale together with relatives of the A, B and D genome donors of wheat, Triticum urartu (A), Aegilops speltoides (S) and Aegilops tauschii (D). Additional taxa containing the A genome, Triticum monococcum and its wild relative Triticum boeoticum, were also included. The main focus of the analysis was on the genomic composition of TEs as these make up at least 80% of the overall genome content. Although more than 200 TE families were identified in each species, approximately 50% of the overall genome comprised 12-15 TE families. The BARE1 element was the largest contributor to all genomes, contributing more than 10% to the overall genome. We also found that several TE families differ strongly in their abundance between species, indicating that TE families can thrive extremely successfully in one species while going virtually extinct in another. Additionally, the nucleotide diversity of BARE1 populations within individual genomes was measured. Interestingly, the nucleotide diversity in the domesticated barley H. vulgare cv. Barke was found to be twice as high as in its wild progenitor H. spontaneum, suggesting that the domesticated barley gained nucleotide diversity from the addition of different genotypes during the domestication and breeding process. In the rye/wheat lineage, sequence diversity of BARE1 elements was generally higher, suggesting that factors such as geographical distribution and mating systems might play a role in intragenomic TE diversity.