Nitrate transporter protein NPF5.12 and major latex-like protein MLP6 are important defense factors against Verticillium longisporum.

Plant defense responses to the soil-borne fungus Verticillium longisporum causing stem stripe disease on oilseed rape (Brassica napus) are poorly understood. In this study, a population of recombinant inbred lines (RILs) using the Arabidopsis thaliana accessions Sei-0 and Can-0 was established. Composite interval mapping, transcriptome data and T-DNA mutant screening identified the NITRATE/PEPTIDE TRANSPORTER FAMILY 5.12 gene (AtNPF5.12) as being associated with disease susceptibility in Can-0. Coimmunoprecipitation revealed interaction between AtNPF5.12 and the MAJOR LATEX PROTEIN family member AtMLP6, and fluorescence microscopy confirmed interaction in the plasma membrane and endoplasmic reticulum. CRISPR/Cas9 technology was applied to mutate the NPF5.12 and MLP6 genes in B. napus. Elevated fungal growth in the npf5.12 mlp6 double mutant of both oilseed rape and Arabidopsis demonstrated their importance in defense against V. longisporum. Colonization of this fungus depends also on available nitrates in the host root. Accordingly, the negative effect of nitrate depletion on fungal growth was less pronounced in Atnpf5.12 plants with impaired nitrate transport. In addition, suberin staining revealed involvement of the NPF5.12 and MLP6 genes in suberin barrier formation. Together, these results demonstrate a dependency of multiple plant factors that lead to successful V. longisporum root infection.

Major latex protein-like encoding genes contribute to Rhizoctonia solani defense responses in sugar beet.

Sugar beets are attacked by several pathogens that cause root damages. Rhizoctonia (Greek for "root killer") is one of them. Rhizoctonia root rot has become an increasing problem for sugar beet production and to decrease yield losses agronomical measures are adopted. Here, two partially resistant and two susceptible sugar beet genotypes were used for transcriptome analysis to discover new defense genes to this fungal disease, information to be implemented in molecular resistance breeding. Among 217 transcripts with increased expression at 2 days post-infection (dpi), three resistance-like genes were found. These genes were not significantly elevated at 5 dpi, a time point when increased expression of three Bet v I/Major latex protein (MLP) homologous genes BvMLP1, BvMLP2 and BvML3 was observed in the partially resistant genotypes. Quantitative RT-PCR analysis on diseased sugar beet seedlings validated the activity of BvMLP1 and BvMLP3 observed in the transcriptome during challenge by R. solani. The three BvMLP genes were cloned and overexpressed in Arabidopsis thaliana to further dissect their individual contribution. Transgenic plants were also compared to T-DNA mutants of orthologous MLP genes. Plants overexpressing BvMLP1 and BvMLP3 showed significantly less infection whereas additive effects were seen on Atmlp1/Atmlp3 double mutants. The data suggest that BvMLP1 and BvMLP3 may contribute to the reduction of the Rhizoctonia root rot disease in sugar beet. Impact on the defense reaction from other differential expressed genes observed in the study is discussed.

Analysis of the hybrid genomes of two field isolates of the soil-borne fungal species Verticillium longisporum.

BACKGROUND: Brassica plant species are attacked by a number of pathogens; among them, the ones with a soil-borne lifestyle have become increasingly important. Verticillium stem stripe caused by Verticillium longisporum is one example. This fungal species is thought to be of a hybrid origin, having a genome composed of combinations of lineages denominated A and D. In this study we report the draft genomes of 2 V. longisporum field isolates sequenced using the Illumina technology. Genomic characterization and lineage composition, followed by selected gene analysis to facilitate the comprehension of its genomic features and potential effector categories were performed. RESULTS: The draft genomes of 2 Verticillium longisporum single spore isolates (VL1 and VL2) have an estimated ungapped size of about 70 Mb. The total number of protein encoding genes identified in VL1 was 20,793, whereas 21,072 gene models were predicted in VL2. The predicted genome size, gene contents, including the gene families coding for carbohydrate active enzymes were almost double the numbers found in V. dahliae and V. albo-atrum. Single nucleotide polymorphisms (SNPs) were frequently distributed in the two genomes but the distribution of heterozygosity and depth was not independent. Further analysis of potential parental lineages suggests that the V. longisporum genome is composed of two parts, A1 and D1, where A1 is more ancient than the parental lineage genome D1, the latter being more closer related to V. dahliae. Presence of the mating-type genes MAT1-1-1 and MAT1-2-1 in the V. longisporum genomes were confirmed. However, the MAT genes in V. dahliae, V. albo-atrum and V. longisporum have experienced extensive nucleotide changes at least partly explaining the present asexual nature of these fungal species. CONCLUSIONS: The established draft genome of V. longisporum is comparatively large compared to other studied ascomycete fungi. Consequently, high numbers of genes were predicted in the two V. longisporum genomes, among them many secreted proteins and carbohydrate active enzyme (CAZy) encoding genes. The genome is composed of two parts, where one lineage is more ancient than the part being more closely related to V. dahliae. Dissimilar mating-type sequences were identified indicating possible ancient hybridization events.

Evolutionary Origins of Rhizarian Parasites.

The SAR group (Stramenopila, Alveolata, Rhizaria) is one of the largest clades in the tree of eukaryotes and includes a great number of parasitic lineages. Rhizarian parasites are obligate and have devastating effects on commercially important plants and animals but despite this fact, our knowledge of their biology and evolution is limited. Here, we present rhizarian transcriptomes from all major parasitic lineages in order to elucidate their evolutionary relationships using a phylogenomic approach. Our results suggest that Ascetosporea, parasites of marine invertebrates, are sister to the novel clade Apofilosa. The phytomyxean plant parasites branch sister to the vampyrellid algal ectoparasites in the novel clade Phytorhiza. They also show that Ascetosporea + Apofilosa + Retaria + Filosa + Phytorhiza form a monophyletic clade, although the branching pattern within this clade is difficult to resolve and appears to be model-dependent. Our study does not support the monophyly of the rhizarian parasitic lineages (Endomyxa), suggesting independent origins for rhizarian animal and plant parasites.

Detection of Verticillium species in Swedish soils using real-time PCR.

Verticillium species are soilborne plant pathogens, responsible for big yield losses worldwide. Here, we report improved procedures to generate DNA from Verticillium species imbedded in farm soils. Using new genomic sequence information, primers for V. dahliae, V. albo-atrum, V. tricorpus, and V. longisporum were designed. In a survey of 429 samples from intensively farmed soil of two Swedish regions, only V. dahliae and V. longisporum were identified. A bias towards V. longisporum (40%) was seen in the south, whereas V. dahliae was more frequent in the western region (19%). Analyses of soil and leaf samples from 20 sugar beet fields, where foliar wilting had been observed, revealed V. dahliae DNA in all leaf and soil samples and V. longisporum in 18 soil samples, illustrating host choice and longevity of the V. longisporum microsclerotia. This study demonstrates the applicability of new molecular diagnostic tools that are important for growers of variable crops.

Genome analysis of the sugar beet pathogen Rhizoctonia solani AG2-2IIIB revealed high numbers in secreted proteins and cell wall degrading enzymes.

BACKGROUND: Sugar beet (Beta vulgaris) is a crop cultivated for its high content in sugar, but it is vulnerable to many soil-borne pathogens. One of them is the basidiomycete Rhizoctonia solani. This fungal species has a compatibility system regulating hyphal fusions (anastomosis). Consequently, R. solani species are categorized in anastomosis groups (AGs). AG2-2IIIB isolates are most aggressive on sugar beet. In the present study, we report on the draft genome of R. solani AG2-2IIIB using the Illumina technology. Genome analysis, interpretation and comparative genomics of five sequenced R. solani isolates were carried out. RESULTS: The draft genome of R. solani AG2-2IIIB has an estimated size of 56.02 Mb. In addition, two normalized EST libraries were sequenced. In total 20,790 of 21,980 AG2-2IIIB isotigs (transcript isoforms) were mapped on the genome with more than 95 % sequence identity. The genome of R. solani AG2-2IIIB was predicted to harbor 11,897 genes and 4908 were found to be isolate-specific. R. solani AG2-2IIIB was predicted to contain 1142 putatively secreted proteins and 473 of them were found to be unique for this isolate. The R. solani AG2-2IIIB genome encodes a high number of carbohydrate active enzymes. The highest numbers were observed for the polysaccharide lyases family 1 (PL-1), glycoside hydrolase family 43 (GH-43) and carbohydrate estarase family 12 (CE-12). Transcription analysis of selected genes representing different enzyme clades revealed a mixed pattern of up- and down-regulation six days after infection on sugar beets featuring variable levels of resistance compared to mycelia of the fungus grown in vitro. CONCLUSIONS: The established R. solani AG2-2IIIB genome and EST sequences provide important information on the gene content, gene structure and transcriptional activity for this sugar beet pathogen. The enriched genomic platform provides an important platform to enhance our understanding of R. solani biology.

Phytophthora infestans Argonaute 1 binds microRNA and small RNAs from effector genes and transposable elements.

Phytophthora spp. encode large sets of effector proteins and distinct populations of small RNAs (sRNAs). Recent evidence has suggested that pathogen-derived sRNAs can modulate the expression of plant defense genes. Here, we studied the sRNA classes and functions associated with Phytophthora infestans Argonaute (Ago) proteins. sRNAs were co-immunoprecipitated with three PiAgo proteins and deep sequenced. Twenty- to twenty-two-nucleotide (nt) sRNAs were identified as the main interaction partners of PiAgo1 and high enrichment of 24-26-nt sRNAs was seen in the PiAgo4-bound sample. The frequencies and sizes of transposable element (TE)-derived sRNAs in the different PiAgo libraries suggested diversified roles of the PiAgo proteins in the control of different TE classes. We further provide evidence for the involvement of PiAgo1 in the P. infestans microRNA (miRNA) pathway. Protein-coding genes are probably regulated by the shared action of PiAgo1 and PiAgo5, as demonstrated by analysis of differential expression. An abundance of sRNAs from genes encoding host cell death-inducing Crinkler (CRN) effectors was bound to PiAgo1, implicating this protein in the regulation of the expanded CRN gene family. The data suggest that PiAgo1 plays an essential role in gene regulation and that at least two RNA silencing pathways regulate TEs in the plant-pathogenic oomycete P. infestans.

Genetic and morphological evidence for introgression between three species of willows.

BACKGROUND: Hybridization and introgression are said to occur relatively frequently in plants, and in particular among different species of willows. However, data on the actual frequency of natural hybridization and introgression is rare. Here, we report the first fine-scale genetic analysis of a contact zone shared between the three basket willow species, Salix dasyclados, S. schwerinii and S. viminalis in the vicinity of the Lake Baikal in Southern Siberia. Individuals were sampled in fourteen populations and classified as pure species or hybrids based on a set of morphological characters. They were then genotyped at 384 nuclear SNP and four chloroplast SSR loci. The STRUCTURE and NewHybrids softwares were used to estimate the frequency and direction of hybridization using genotypic data at the nuclear SNP loci. RESULTS: As many as 19 % of the genotyped individuals were classified as introgressed individuals and these were mainly encountered in the centre of the contact zone. All introgressed individuals were backcrosses to S. viminalis or S. schwerinii and no F1 or F2 hybrids were found. The rest of the genotyped individuals were classified as pure species and formed two clusters, one with S. schwerinii individuals and the other with S. viminalis and S. dasyclados individuals. The two clusters were significantly genetically differentiated, with F ST = 0.333 (0.282-0.382, p < 0.001). In contrast, for the chloroplast haplotypes, no genetic differentiation was observed as they were completely shared between the species. Based on morphological classification only 5 % of the individuals were classified as introgressed individuals, which was much less than what was detected using genotypic data. CONCLUSIONS: We have discovered a new willow hybrid zone with relatively high frequency of introgressed individuals. The low frequency of F1 hybrids indicates that ongoing hybridization is limited, which could be because of the presence of reproductive barriers or simply because the conditions are not favorable for hybridization. We further conclude that in order to get a complete picture of the species composition of a hybrid zone it is necessary to use a combination of morphological characters and genetic data from both nuclear and chloroplast markers.

Plant-mediated gene silencing restricts growth of the potato late blight pathogen Phytophthora infestans.

Phytophthora infestans is an oomycete that causes severe damage to potato, and is well known for its ability to evolve rapidly in order to overcome resistant potato varieties. An RNA silencing strategy was evaluated here to clarify if small interfering RNA homologous to selected genes in P. infestans could be targeted from the plant host to reduce the magnitude of the infection. As a proof-of-concept, a hairpin RNA (hp-RNA) construct using the GFP marker gene was designed and introduced in potato. At 72 hpi, a 55-fold reduction of the signal intensity of a corresponding GFP expressing P. infestans strain on leaf samples of transgenic plants, compared with wild-type potato, was detected. This suggests that an RNA interference construct in the potato host could be processed and target a transcript of the pathogen. Three genes important in the infection process of P. infestans, PiGPB1, PiCESA2, and PiPEC, together with PiGAPDH taking part in basic cell maintenance were subsequently tested using an analogous transgenic strategy. Out of these gene candidates, the hp-PiGPB1 targeting the G protein ß-subunit (PiGPB1) important for pathogenicity resulted in most restricted disease progress. Further, Illumina sequencing of inoculated transgenic potato leaves revealed sRNAs of 24/25 nt size homologous to the PiGPB1 gene in the transgenic plants indicating post-transcriptional silencing of the target gene. The work demonstrates that a host-induced gene-silencing approach is functional against P. infestans but is highly dependent on target gene for a successful outcome. This finding broadens the arsenal of control strategies to this important plant disease.

Loss and retention of resistance genes in five species of the Brassicaceae family.

BACKGROUND: Plants have evolved disease resistance (R) genes encoding for nucleotide-binding site (NB) and leucine-rich repeat (LRR) proteins with N-terminals represented by either Toll/Interleukin-1 receptor (TIR) or coiled-coil (CC) domains. Here, a genome-wide study of presence and diversification of CC-NB-LRR and TIR-NB-LRR encoding genes, and shorter domain combinations in 19 Arabidopsis thaliana accessions and Arabidopsis lyrata, Capsella rubella, Brassica rapa and Eutrema salsugineum are presented. RESULTS: Out of 528 R genes analyzed, 12 CC-NB-LRR and 17 TIR-NB-LRR genes were conserved among the 19 A. thaliana genotypes, while only two CC-NB-LRRs, including ZAR1, and three TIR-NB-LRRs were conserved when comparing the five species. The RESISTANCE TO LEPTOSPHAERIA MACULANS 1 (RLM1) locus confers resistance to the Brassica pathogen L. maculans the causal agent of blackleg disease and has undergone conservation and diversification events particularly in B. rapa. On the contrary, the RLM3 locus important in the immune response towards Botrytis cinerea and Alternaria spp. has recently evolved in the Arabidopsis genus. CONCLUSION: Our genome-wide analysis of the R gene repertoire revealed a large sequence variation in the 23 cruciferous genomes. The data provides further insights into evolutionary processes impacting this important gene family.

Fragmentation of tRNA in Phytophthora infestans asexual life cycle stages and during host plant infection.

BACKGROUND: The oomycete Phytophthora infestans possesses active RNA silencing pathways, which presumably enable this plant pathogen to control the large numbers of transposable elements present in its 240 Mb genome. Small RNAs (sRNAs), central molecules in RNA silencing, are known to also play key roles in this organism, notably in regulation of critical effector genes needed for infection of its potato host. RESULTS: To identify additional classes of sRNAs in oomycetes, we mapped deep sequencing reads to transfer RNAs (tRNAs) thereby revealing the presence of 19-40 nt tRNA-derived RNA fragments (tRFs). Northern blot analysis identified abundant tRFs corresponding to half tRNA molecules. Some tRFs accumulated differentially during infection, as seen by examining sRNAs sequenced from P. infestans-potato interaction libraries. The putative connection between tRF biogenesis and the canonical RNA silencing pathways was investigated by employing hairpin RNA-mediated RNAi to silence the genes encoding P. infestans Argonaute (PiAgo) and Dicer (PiDcl) endoribonucleases. By sRNA sequencing we show that tRF accumulation is PiDcl1-independent, while Northern hybridizations detected reduced levels of specific tRNA-derived species in the PiAgo1 knockdown line. CONCLUSIONS: Our findings extend the sRNA diversity in oomycetes to include fragments derived from non-protein-coding RNA transcripts and identify tRFs with elevated levels during infection of potato by P. infestans.

Dominance of Mating Type A1 and Indication of Epigenetic Effects During Early Stages of Mating in Phytophthora infestans.

The potato late blight pathogen Phytophthora infestans has both an asexual and a sexual mode of reproduction. In Scandinavia, the pathogen is reproducing sexually on a regular basis, whereas clonal lineages dominate in other geographical regions. This study aimed at elucidating events or key genes underlying this difference in sexual behavior. First, the transcriptomes of eight strains, known as either clonal or sexual, were compared during early stages of mating. Principal component analysis (PCA) divided the samples in two clusters A and B and a clear grouping of the mating samples together with the A1 mating type parents was observed. Induction of genes encoding DNA adenine N6-methylation (6mA) methyl-transferases clearly showed a bias toward the cluster A. In contrast, the Avrblb2 effector gene family was highly induced in most of the mating samples and was associated with cluster B in the PCA, similarly to genes coding for acetyl-transferases, which play an important role in RXLR modification prior to secretion. Avrblb2 knock-down strains displayed a reduction in virulence and oospore formation, suggesting a role during the mating process. In conclusion, a number of gene candidates important for the reproductive processes were revealed. The results suggest a possible epigenetic influence and involvement of specific RXLR effectors in mating-related processes.

The architecture of the Plasmodiophora brassicae nuclear and mitochondrial genomes.

Plasmodiophora brassicae is a soil-borne pathogen that attacks roots of cruciferous plants causing clubroot disease. The pathogen belongs to the Plasmodiophorida order in Phytomyxea. Here we used long-read SMRT technology to clarify the P. brassicae e3 genomic constituents along with comparative and phylogenetic analyses. Twenty contigs representing the nuclear genome and one mitochondrial (mt) contig were generated, together comprising 25.1 Mbp. Thirteen of the 20 nuclear contigs represented chromosomes from telomere to telomere characterized by [TTTTAGGG] sequences. Seven active gene candidates encoding synaptonemal complex-associated and meiotic-related protein homologs were identified, a finding that argues for possible genetic recombination events. The circular mt genome is large (114,663 bp), gene dense and intron rich. It shares high synteny with the mt genome of Spongospora subterranea, except in a unique 12 kb region delimited by shifts in GC content and containing tandem minisatellite- and microsatellite repeats with partially palindromic sequences. De novo annotation identified 32 protein-coding genes, 28 structural RNA genes and 19 ORFs. ORFs predicted in the repeat-rich region showed similarities to diverse organisms suggesting possible evolutionary connections. The data generated here form a refined platform for the next step involving functional analysis, all to clarify the complex biology of P. brassicae.

Association mapping in Salix viminalis L. (Salicaceae) - identification of candidate genes associated with growth and phenology.

Willow species (Salix) are important as short-rotation biomass crops for bioenergy, which creates a demand for faster genetic improvement and breeding through deployment of molecular marker-assisted selection (MAS). To find markers associated with important adaptive traits, such as growth and phenology, for use in MAS, we genetically dissected the trait variation of a Salix viminalis (L.) population of 323 accessions. The accessions were sampled throughout northern Europe and were established at two field sites in Pustnäs, Sweden, and at Woburn, UK, offering the opportunity to assess the impact of genotype-by-environment interactions (G × E) on trait-marker associations. Field measurements were recorded for growth and phenology traits. The accessions were genotyped using 1536 SNP markers developed from phenology candidate genes and from genes previously observed to be differentially expressed in contrasting environments. Association mapping between 1233 of these SNPs and the measured traits was performed taking into account population structure and threshold selection bias. At a false discovery rate (FDR) of 0.2, 29 SNPs were associated with bud burst, leaf senescence, number of shoots or shoot diameter. The percentage of accession variation (Radj2) explained by these associations ranged from 0.3% to 4.4%, suggesting that the studied traits are controlled by many loci of limited individual impact. Despite this, a SNP in the EARLY FLOWERING 3 gene was repeatedly associated (FDR < 0.2) with bud burst. The rare homozygous genotype exhibited 0.4-1.0 lower bud burst scores than the other genotype classes on a five-grade scale. Consequently, this marker could be promising for use in MAS and the gene deserves further study. Otherwise, associations were less consistent across sites, likely due to their small Radj2 estimates and to considerable G × E interactions indicated by multivariate association analyses and modest trait accession correlations across sites (0.32-0.61).

The Plasmodiophora brassicae genome reveals insights in its life cycle and ancestry of chitin synthases.

Plasmodiophora brassicae causes clubroot, a major disease of Brassica oil and vegetable crops worldwide. P. brassicae is a Plasmodiophorid, obligate biotrophic protist in the eukaryotic kingdom of Rhizaria. Here we present the 25.5 Mb genome draft of P. brassicae, developmental stage-specific transcriptomes and a transcriptome of Spongospora subterranea, the Plasmodiophorid causing powdery scab on potato. Like other biotrophic pathogens both Plasmodiophorids are reduced in metabolic pathways. Phytohormones contribute to the gall phenotypes of infected roots. We report a protein (PbGH3) that can modify auxin and jasmonic acid. Plasmodiophorids contain chitin in cell walls of the resilient resting spores. If recognized, chitin can trigger defense responses in plants. Interestingly, chitin-related enzymes of Plasmodiophorids built specific families and the carbohydrate/chitin binding (CBM18) domain is enriched in the Plasmodiophorid secretome. Plasmodiophorids chitin synthases belong to two families, which were present before the split of the eukaryotic Stramenopiles/Alveolates/Rhizaria/Plantae and Metazoa/Fungi/Amoebozoa megagroups, suggesting chitin synthesis to be an ancient feature of eukaryotes. This exemplifies the importance of genomic data from unexplored eukaryotic groups, such as the Plasmodiophorids, to decipher evolutionary relationships and gene diversification of early eukaryotes.

High rates of gene flow by pollen and seed in oak populations across Europe.

Gene flow is a key factor in the evolution of species, influencing effective population size, hybridisation and local adaptation. We analysed local gene flow in eight stands of white oak (mostly Quercus petraea and Q. robur, but also Q. pubescens and Q. faginea) distributed across Europe. Adult trees within a given area in each stand were exhaustively sampled (range [239, 754], mean 423), mapped, and acorns were collected ([17,147], 51) from several mother trees ([3], [47], 23). Seedlings ([65,387], 178) were harvested and geo-referenced in six of the eight stands. Genetic information was obtained from screening distinct molecular markers spread across the genome, genotyping each tree, acorn or seedling. All samples were thus genotyped at 5-8 nuclear microsatellite loci. Fathers/parents were assigned to acorns and seedlings using likelihood methods. Mating success of male and female parents, pollen and seed dispersal curves, and also hybridisation rates were estimated in each stand and compared on a continental scale. On average, the percentage of the wind-borne pollen from outside the stand was 60%, with large variation among stands (21-88%). Mean seed immigration into the stand was 40%, a high value for oaks that are generally considered to have limited seed dispersal. However, this estimate varied greatly among stands (20-66%). Gene flow was mostly intraspecific, with large variation, as some trees and stands showed particularly high rates of hybridisation. Our results show that mating success was unevenly distributed among trees. The high levels of gene flow suggest that geographically remote oak stands are unlikely to be genetically isolated, questioning the static definition of gene reserves and seed stands.

Evidence for small RNAs homologous to effector-encoding genes and transposable elements in the oomycete Phytophthora infestans.

Phytophthora infestans is the oomycete pathogen responsible for the devastating late blight disease on potato and tomato. There is presently an intense research focus on the role(s) of effectors in promoting late blight disease development. However, little is known about how they are regulated, or how diversity in their expression may be generated among different isolates. Here we present data from investigation of RNA silencing processes, characterized by non-coding small RNA molecules (sRNA) of 19-40 nt. From deep sequencing of sRNAs we have identified sRNAs matching numerous RxLR and Crinkler (CRN) effector protein genes in two isolates differing in pathogenicity. Effector gene-derived sRNAs were present in both isolates, but exhibited marked differences in abundance, especially for CRN effectors. Small RNAs in P. infestans grouped into three clear size classes of 21, 25/26 and 32 nt. Small RNAs from all size classes mapped to RxLR effector genes, but notably 21 nt sRNAs were the predominant size class mapping to CRN effector genes. Some effector genes, such as PiAvr3a, to which sRNAs were found, also exhibited differences in transcript accumulation between the two isolates. The P. infestans genome is rich in transposable elements, and the majority of sRNAs of all size classes mapped to these sequences, predominantly to long terminal repeat (LTR) retrotransposons. RNA silencing of Dicer and Argonaute genes provided evidence that generation of 21 nt sRNAs is Dicer-dependent, while accumulation of longer sRNAs was impacted by silencing of Argonaute genes. Additionally, we identified six microRNA (miRNA) candidates from our sequencing data, their precursor sequences from the genome sequence, and target mRNAs. These miRNA candidates have features characteristic of both plant and metazoan miRNAs.

Polymorphism and Divergence in Two Willow Species, Salix viminalis L. and Salix schwerinii E. Wolf.

We investigated species divergence, present and past gene flow, levels of nucleotide polymorphism, and linkage disequilibrium in two willows from the plant genus Salix. Salix belongs together with Populus to the Salicaceae family; however, most population genetic studies of Salicaceae have been performed in Populus, the model genus in forest biology. Here we present a study on two closely related willow species Salix viminalis and S. schwerinii, in which we have resequenced 33 and 32 nuclear gene segments representing parts of 18 nuclear loci in 24 individuals for each species. We used coalescent simulations and estimated the split time to around 600,000 years ago and found that there is currently limited gene flow between the species. Mean intronic nucleotide diversity across gene segments was slightly higher in S. schwerinii (π(i) = 0.00849) than in S. viminalis (π(i) = 0.00655). Compared with other angiosperm trees, the two willows harbor intermediate levels of silent polymorphisms. The decay of linkage disequilibrium was slower in S. viminalis compared with S. schwerinii, and we speculate that this is due to different demographic histories as S. viminalis has been partly domesticated in Europe.

Cryptic population genetic structure: the number of inferred clusters depends on sample size.

Clustering methods have been used extensively to unravel cryptic population genetic structure. We investigated the effect of the number of individuals sampled in each location on the resulting number of clusters. Our study was motivated by recent results in Arabidopsis thaliana: studies in which more than one individual was sampled per location apparently have led to a much higher number of clusters than studies where only one individual was sampled in each location, as is generally done in this species. We show, using computer simulations and microsatellite data in A. thaliana, that the number of sampled individuals indeed has a strong impact on the number of resulting clusters. This effect is smaller if the sampled populations have a hierarchical structure. In most cases, sampling 5-10 individuals per population should be enough. The results argue for abandoning the concept of 'accessions' in partially selfing organisms.