Genome scans and elusive candidate genes: detecting the variation that matters for speciation

Next‐generation sequencing is providing us with vast amounts of genetic data, yet we are currently struggling in our attempts to make sense of them. In particular, it has proven difficult to link phenotypic divergence and speciation to genome level divergence. In the current issue of Molecular Ecology, Ruegg et al. ( 2014 ) present new empirical results from two closely related bird taxa. They use a promising approach combining genome scan and candidate gene analysis. Their results suggest that we may have been looking in vain for candidate speciation genes by focusing only on genes found within genomic islands of divergence. This is because genes important in divergence and speciation may not be detected by genome scans and because features of the genomic architecture per se may have a large effect on the pattern of genome divergence.     

Fitness Is Strongly Influenced by Rare Mutations of Large Effect in a Microbial Mutation Accumulation Experiment

Our understanding of the evolutionary consequences of mutation relies heavily on estimates of the rate and fitness effect of spontaneous mutations generated by mutation accumulation (MA) experiments. We performed a classic MA experiment in which frequent sampling of MA lines was combined with whole genome resequencing to develop a high-resolution picture of the effect of spontaneous mutations in a hypermutator (ΔmutS) strain of the bacterium Pseudomonas aeruginosa. After ∼644 generations of mutation accumulation, MA lines had accumulated an average of 118 mutations, and we found that average fitness across all lines decayed linearly over time. Detailed analyses of the dynamics of fitness change in individual lines revealed that a large fraction of the total decay in fitness (42.3%) was attributable to the fixation of rare, highly deleterious mutations (comprising only 0.5% of fixed mutations). Furthermore, we found that at least 0.64% of mutations were beneficial and probably fixed due to positive selection. The majority of mutations that fixed (82.4%) were base substitutions and we failed to find any signatures of selection on nonsynonymous or intergenic mutations. Short indels made up a much smaller fraction of the mutations that were fixed (17.4%), but we found evidence of strong selection against indels that caused frameshift mutations in coding regions. These results help to quantify the amount of natural selection present in microbial MA experiments and demonstrate that changes in fitness are strongly influenced by rare mutations of large effect.     

Toward genomic prediction from whole-genome sequence data: impact of sequencing design on genotype imputation and accuracy of predictions

Genomic prediction from whole-genome sequence data is attractive, as the accuracy of genomic prediction is no longer bounded by extent of linkage disequilibrium between DNA markers and causal mutations affecting the trait, given the causal mutations are in the data set. A cost-effective strategy could be to sequence a small proportion of the population, and impute sequence data to the rest of the reference population. Here, we describe strategies for selecting individuals for sequencing, based on either pedigree relationships or haplotype diversity. Performance of these strategies (number of variants detected and accuracy of imputation) were evaluated in sequence data simulated through a real Belgian Blue cattle pedigree. A strategy (AHAP), which selected a subset of individuals for sequencing that maximized the number of unique haplotypes (from single-nucleotide polymorphism panel data) sequenced gave good performance across a range of variant minor allele frequencies. We then investigated the optimum number of individuals to sequence by fold coverage given a maximum total sequencing effort. At 600 total fold coverage (x 600), the optimum strategy was to sequence 75 individuals at eightfold coverage. Finally, we investigated the accuracy of genomic predictions that could be achieved. The advantage of using imputed sequence data compared with dense SNP array genotypes was highly dependent on the allele frequency spectrum of the causative mutations affecting the trait. When this followed a neutral distribution, the advantage of the imputed sequence data was small; however, when the causal mutations all had low minor allele frequencies, using the sequence data improved the accuracy of genomic prediction by up to 30%.     

The considerable genome size variation of Hordeum species (poaceae) is linked to phylogeny, life form, ecology, and speciation rates

Genome size variation in plants is thought to be correlatedwith cytological, physiological, or ecological characters. However,conclusions drawn in several studies were often contradictory.To analyze nuclear genome size evolution in a phylogenetic framework,DNA contents of 134 accessions, representing all but one speciesof the barley genus Hordeum L., were measured by flow cytometry.The 2C DNA contents were in a range from 6.85 to 10.67 pg indiploids (2n = 14) and reached up to 29.85 pg in hexaploid species(2n = 42). The smallest genomes were found in taxa from theNew World, which became secondarily annual, whereas the largestdiploid genomes occur in Eurasian annuals. Genome sizes of polyploidtaxa equaled mostly the added sizes of their proposed progenitorsor were slightly (1% to 5%) smaller. The analysis of ancestralgenome sizes on the base of the phylogeny of the genus revealedlineages with decreasing and with increasing genome sizes. Correlationsof intraspecific genome size variation with the length of vegetationperiod were found in H. marinum populations from Western Europebut were not significant within two species from South America.On a higher taxonomical level (i.e., for species groups or theentire genus), environmental correlations were absent. Thiscould mostly be attributed to the superimposition of life-formchanges and phylogenetic constraints, which conceal ecogeographicalcorrelations.     

Deep genomic analysis of the Chlorella sorokiniana SAG 211-8k chloroplast

The chloroplast genome contains information that is applicable in many scientific fields, such as plant systematics, phylogenetic reconstruction and biotechnology, because its features are highly conserved among species. To date, several complete green algal chloroplast genomes have been sequenced and assembled. In this study, the nucleotide sequence of the chloroplast genome (cpDNA) of Chlorella sorokiniana SAG 211-8k is reported and compared for the first time to the chloroplast genomes of 10 Chlorellaceae. The recently updated Chlorella sorokiniana cpDNA sequence, assembled as a circular map of 109?811 bp, encodes 113 genes. Similar to other Chlorella strains, this chloroplast genome does not show a quadripartite structure and lacks the large rRNA operon-encoding Inverted Repeat (IR). The Chlorella sorokiniana plastid encodes the tRNA(Ile)-lysidine synthetase (tilS), which is responsible for modifying the CAU anticodon of a unique tRNA. Gene ordering and clustering highlight the close relationships among Chlorella clade members and the preservation of crucial gene clusters in photosynthetic strains. The features of Chlorella sorokiniana presented here reinforce the monophyletic character of Chlorellaceae and provide important information that sheds light on chloroplast genome evolution among species of Chlorella.     

Rare symmetric and asymmetric Nicotiana tabacum (+) N. megalosiphon somatic hybrids recovered by selection for nuclear-encoded resistance genes and in the absence of genome inactivation

Following protoplast fusion between Nicotiana tabacum (dhfr) and N. megalosiphon (nptII) somatic hybrids were selected on the basis of dual resistance to kanamycin and methotrexate. Despite strong selection for parental nuclear-encoded resistances, only nine N. tabacum (+) N. megalosiphon somatic hybrids were obtained. A preferential loss of the parental N. tabacum nuclear and organelle genome was apparent in some plants in spite of the lack of genomic inactivation by the irradiation or chemical treatment of the parental protoplasts. Only six of the nine hybrids recovered possessed both parental profiles of nuclear RFLPs and isoenzymes. The remaining three hybrids were highly asymmetric with two being identical to N. megalosiphon except for minor morphological differences and rearranged or recombined mitochondrial DNAs (mtDNA), while the other one was distinguishable only by the presence of a rearranged or recombined mtDNA, and was therefore possibly a cybrid. Overall, eight somatic hybrids possessed rearranged or recombined mtDNAs and chloroplast inheritance was non-random since eight possessed N. megalosiphon-type chloroplasts and only one had N. tabacum chloroplasts. In contrast, using the same selection approach, numerous morphologically similar symmetric somatic hybrids with nuclear RFLPs and isozymes of both the parental species were recovered from control fusions between N. tabacum and the more closely related N. sylvestris. In spite of the low frequency of recovery of symmetric N. tabacum (+) N. megalosiphon hybrids in this study, one of these hybrids displayed a significant degree of self-fertility allowing for back-crosses to transfer N. megalosiphon disease-resistance traits to N. tabacum. Plant Research Centre Contribution No. 1579     

A short note on predicted motifs in the Thaparocleidus wallagonius genome

Thaparocleidus wallagonius is a monogenean parasite and a fish-borne pathogen with a worldwide distribution. The genome for Thaparocleidus wallagonius is known. Therefore, it is of interest to report the DNA motif analysis data in the 18S rDNA of Thaparocleidus wallagonius collected from the fish Wallago attu in India. This data forms a framework for an in-depth analysis of the parasite biology and development, immune evasion strategies, virulence and long-term survival within the definitive host.     

Progress in fold recognition

The prediction experiment reveals that fold recognition has become a powerful tool in structural biology. We applied our fold recognition technique to 13 target sequences. In two cases, replication terminating protein and prosequence of subtilisin, the predicted structures are very similar to the experimentally determined folds. For the first time, in a public blind test, the unknown structures of proteins have been predicted ahead of experiment to an accuracy approaching molecular detail. In two other cases the approximate folds have been predicted correctly. According to the assessors there were 12 recognizable folds among the target proteins. In our postprediction analysis we find that in 7 cases our fold recognition technique is successful. In several of the remaining cases the predicted folds have interesting features in common with the experimental results. We present our procedure, discuss the results, and comment on several fundamental and technical problems encountered in fold recognition. © 1995 Wiley-Liss, Inc.     

Landscape genomics of Colorado potato beetle provides evidence of polygenic adaptation to insecticides

The ability of insect pests to rapidly and repeatedly adapt to insecticides has long challenged entomologists and evolutionary biologists. Since Crow's seminal paper on insecticide resistance in 1957, new data and insights continue to emerge that are relevant to the old questions about how insecticide resistance evolves: such as whether it is predominantly mono‐ or polygenic, and evolving from standing vs. de novo genetic variation. Many studies support the monogenic hypothesis, and current management recommendations assume single‐ or two‐locus models. But inferences could be improved by integrating data from a broader sample of pest populations and genomes. Here, we generate evidence relevant to these questions by applying a landscape genomics framework to the study of insecticide resistance in a major agricultural pest, Colorado potato beetle, Leptinotarsa decemlineata (Say). Genome–environment association tests using genomic variation from 16 populations spanning gradients of landscape variables associated with insecticide exposure over time revealed 42 strong candidate insecticide resistance genes, with potentially overlapping roles in multiple resistance mechanisms. Measurements of resistance to a widely used insecticide, imidacloprid, among 47 L. decemlineata populations revealed heterogeneity at a small (2 km) scale and no spatial signature of origin or spread throughout the landscape. Analysis of nucleotide diversity suggested candidate resistance loci have undergone varying degrees of selective sweeps, often maintaining similar levels of nucleotide diversity to neutral loci. This study suggests that many genes are involved in insecticide resistance in L. decemlineata and that resistance likely evolves from both de novo and standing genetic variation.     

Polyploidy does not control all: Lineage‐specific average chromosome length constrains genome size evolution in ferns

Recent studies investigating the evolution of genome size diversity in ferns have shown that they have a distinctive genome profile compared with other land plants. Ferns are typically characterized by possessing medium‐sized genomes, although a few lineages have evolved very large genomes. Ferns are different from other vascular plant lineages as they are the only group to show evidence for a correlation between genome size and chromosome number. In this study, we aim to explore whether the evolution of fern genome sizes is not only shaped by chromosome number changes arising from polyploidy but also by constraints on the average amount of DNA per chromosome. We selected the genus Asplenium L. as a model genus to study the question because of the unique combination of a highly conserved base chromosome number and a high frequency of polyploidy. New genome size data for Asplenium taxa were combined with existing data and analyzed within a phylogenetic framework. Genome size varied substantially between diploid species, resulting in overlapping genome sizes among diploid and tetraploid spleenworts. The observed additive pattern indicates the absence of genome downsizing following polyploidy. The genome size of diploids varied non‐randomly and we found evidence for clade‐specific trends towards larger or smaller genomes. The 578‐fold range of fern genome sizes have arisen not only from repeated cycles of polyploidy but also through clade‐specific constraints governing accumulation and/or elimination of DNA.