Evaluating the role of reference-genome phylogenetic distance on evolutionary inference

When a high-quality genome assembly of a target species is unavailable, an option to avoid the costly de novo assembly process is a mapping-based assembly. However, mapping shotgun data to a distant relative may lead to biased or erroneous evolutionary inference. Here, we used short-read data from a mammal (beluga whale) and a bird species (rowi kiwi) to evaluate whether reference genome phylogenetic distance can impact downstream demographic (Pairwise Sequentially Markovian Coalescent) and genetic diversity (heterozygosity, runs of homozygosity) analyses. We mapped to assemblies of species of varying phylogenetic distance (from conspecific to genome-wide divergence of >7%), and de novo assemblies created using cross-species scaffolding. We show that while reference genome phylogenetic distance has an impact on demographic analyses, it is not pronounced until using a reference genome with >3% divergence from the target species. When mapping to cross-species scaffolded assemblies, we are unable to replicate the original beluga demographic results, but are able with the rowi kiwi, presumably reflecting the more fragmented nature of the beluga assemblies. We find that increased phylogenetic distance has a pronounced impact on genetic diversity estimates; heterozygosity estimates deviate incrementally with increasing phylogenetic distance. Moreover, runs of homozygosity are largely undetectable when mapping to any nonconspecific assembly. However, these biases can be reduced when mapping to a cross-species scaffolded assembly. Taken together, our results show that caution should be exercised when selecting reference genomes. Cross-species scaffolding may offer a way to avoid a costly, traditional de novo assembly, while still producing robust, evolutionary inference.     

Illuminating the mysteries of wolf history

One of the most enduring surprises about the genetic history of Late Pleistocene populations is that continuity is often disturbed by upheaval. In fact, studies that support population continuity are increasingly rare in humans, a variety of vertebrate taxa, and vascular plants (Hofreiter & Stewart 2009; Burbrink et al. 2016). Perhaps such continuity should not be expected as the Pleistocene is marked by episodes of climate change, glaciation and the invasions of humans into previously isolated areas. Although fossils are one of the primary sources for inferring population continuity, a problem with fossil material is that, even if similar morphological forms might exist in a place over time, they may not be from the same genetic lineage. There are now readily available methods to assess genetic continuity solely from DNA found in fossil material, provided the record is fairly continuous. In a From the Cover article in this issue of Molecular Ecology, Loog et al. (2020) apply some of these readily available methods to analyse mitochondrial genomes and model the demography of wolves over the last 50,000 years.     

The discovery of Antarctic RNA viruses: a new game changer

Antarctic ecosystems are dominated by micro‐organisms, and viruses play particularly important roles in the food webs. Since the first report in 2009 (López‐Bueno et al. 2009 ), ‘omic’‐based studies have greatly enlightened our understanding of Antarctic aquatic microbial diversity and ecosystem function (Wilkins et al. 2013 ; Cavicchioli 2015 ). This has included the discovery of many new eukaryotic viruses (López‐Bueno et al. 2009 ), virophage predators of algal viruses (Yau et al. 2011 ), bacteria with resistance to phage (Lauro et al. 2011 ) and mechanisms of haloarchaeal evasion, defence and adaptation to viruses (Tschitschko et al. 2015 ). In this issue of Molecular Ecology, López‐Bueno et al. ( 2015 ) report the first discovery of RNA viruses from an Antarctic aquatic environment. High sequence coverage enabled genome variation to be assessed for four positive‐sense single‐stranded RNA viruses from the order Picornavirales. By examining the populations present in the water column and in the lake's catchment area, populations of ‘quasispecies’ were able to be linked to local environmental factors. In view of the importance of viruses in Antarctic ecosystems but lack of data describing them, this study represents a significant advance in the field.     

Don't throw the baby out with the bathwater: identifying and mapping paralogs in salmonids

Many eukaryotic genomes contain a large fraction of gene duplicates (or paralogs) as a result of ancient or recent whole‐genome duplications (Ohno 1970 ; Jaillon et al. 2004 ; Kellis et al. 2004 ). Identifying paralogs with NGS data is a pervasive problem in both ancient polyploids and neopolyploids. Likewise, paralogs are often treated as a nuisance that has to be detected and removed (Everett et al. 2012 ). In this issue of Molecular Ecology Resources, Waples et al. ( 2015 ) show that exclusion might not be necessary and how we may miss out on important genomic information in doing so. They present a novel statistical approach to detect paralogs based on the segregation of RAD loci in haploid offspring and test their method by constructing linkage maps with and without these duplicated loci in chum salmon, Oncorhynchus keta (Fig.  1 ). Their linkage map including the resolved paralogs shows that these are mostly located in the distal regions of several linkage groups. Particularly intriguing is their finding that these homoeologous regions appear impoverished in transposable elements (TE). Given the role that TE play in genome remodelling, it is noteworthy that these elements are of low abundance in regions showing residual tetrasomic inheritance. This raises the question whether re‐diploidization is constrained in these regions and whether they might have a role to play in salmonid speciation. This study provides an original approach to identifying duplicated loci in species with a pedigree, as well as providing a dense linkage map for chum salmon, and interesting insights into the retention of gene duplicates in an ancient polyploid.     

New gene-derived simple sequence repeat markers for common bean (Phaseolus vulgaris L.)

Common bean is an important and diverse crop legume with several wild relatives that are all part of the Phaseoleae tribe of tropical crop legumes. Sequence databases have been a good source of sequences to mine for simple sequence repeats (SSRs). The objective of this research was to evaluate 14 sequence collections from common bean for SSRs and to evaluate the diversity of the polymorphic microsatellites derived from these collections. SSRs were found in 10 of the GenBank sequence collections with an average of 11.3% of sequences containing microsatellite motifs. The most common motifs were based on tri- and dinucleotides. In a marker development programme, primers were designed for 125 microsatellites which were tested on a panel of 18 common bean genotypes. The markers were named as part of the bean microsatellite-database (BMd) series, and the average polymorphism information content was 0.404 for polymorphic markers and predicted well the genepool structure of common beans and the status of the wild and cultivated accessions that were included in the study. Therefore, the BMd series of microsatellites is useful for multiple studies of genetic relatedness and as anchor markers in future mapping of wide crosses in the species.     

No signature of Y chromosomal resemblance between possible descendants of the Cimbri in Denmark and Northern Italy

Two European populations are believed to be related to the ancient Germanic tribe Cimbri: one living in Northern Italy, the other living in Jutland, Denmark. The people called Cimbri are documented in the ancient Roman historical record. Arriving from the far north their movements can be tracked from successive battles with the Romans. The Cimbri finally entered Italy from the northeast and were defeated at Vercellae (present day Vercelli) in 101 BC by Gaius Marius and his professional legions. Classical sources from the first centuries AD relate the homeland of the Cimbri to the coasts around the Elb estuary (northern Germany) or specifically towards the north (Himmerland in northern Jutland). In the alpine parts of Veneto, northeast of the historical battlefield, local traditions dating back to late medieval time, identify a local population as Cimbri living in Terra dei Cimbri. They are considered the descendents of the Germanic combatants that fled the battlefield at Vercelli. As the defeated Cimbri that possibly fled to the mountains of Northern Italy most likely would have been male (warriors), the present study investigated the possible Y chromosomal diversity of the two present populations using microsatellite markers and single nucleotide polymorphisms. While Cimbri from Himmerland resembled their geographical neighbors from Denmark for the Y-chromosome markers, Cimbri from Italy were significantly differentiated both from Cimbri from Himmerland and from Danes. Therefore, we were not able to show any biological relationship for uniparentally transmitted markers.     

Critical review of NGS analyses for de novo genotyping multigene families

The genotyping of highly polymorphic multigene families across many individuals used to be a particularly challenging task because of methodological limitations associated with traditional approaches. Next‐generation sequencing (NGS) can overcome most of these limitations, and it is increasingly being applied in population genetic studies of multigene families. Here, we critically review NGS bioinformatic approaches that have been used to genotype the major histocompatibility complex (MHC) immune genes, and we discuss how the significant advances made in this field are applicable to population genetic studies of gene families. Increasingly, approaches are introduced that apply thresholds of sequencing depth and sequence similarity to separate alleles from methodological artefacts. We explain why these approaches are particularly sensitive to methodological biases by violating fundamental genotyping assumptions. An alternative strategy that utilizes ultra‐deep sequencing (hundreds to thousands of sequences per amplicon) to reconstruct genotypes and applies statistical methods on the sequencing depth to separate alleles from artefacts appears to be more robust. Importantly, the ‘degree of change’ (DOC) method avoids using arbitrary cut‐off thresholds by looking for statistical boundaries between the sequencing depth for alleles and artefacts, and hence, it is entirely repeatable across studies. Although the advances made in generating NGS data are still far ahead of our ability to perform reliable processing, analysis and interpretation, the community is developing statistically rigorous protocols that will allow us to address novel questions in evolution, ecology and genetics of multigene families. Future developments in third‐generation single molecule sequencing may potentially help overcome problems that still persist in de novo multigene amplicon genotyping when using current second‐generation sequencing approaches.     

Exploring linkage disequilibrium

Linkage disequilibrium (LD, association of allelic states across loci) is poorly understood by many evolutionary biologists, but as technology for multilocus sampling improves, we ignore LD at our peril. If we sample variation at 10 loci in an organism with 20 chromosomes, we can reasonably treat them as 10 ‘independent witnesses’ of the evolutionary process. If instead, we sample variation at 1000 loci, many are bound to be close together on a chromosome. With only one or two crossovers per meiosis, associations between close neighbours decay so slowly that even LD created far in the past will not have dissipated, so we cannot treat the 1000 loci as independent witnesses (Barton 2011 ). This means that as marker density on genomes increases classic analyses assuming independent loci become mired in the problem of overconfidence: if 1000 independent witnesses are assumed, and that number should be much lower, any conclusion will be overconfident. This is of special concern because our literature suffers from a strong publication bias towards confident answers, even when they turn out to be wrong (Knowles 2008 ). In contrast, analyses that take into account associations across loci both control for overconfidence and can inform us about LD generating events far in the past, for example human/Neanderthal admixture (Fu et al. 2014 ). With increased marker density, biologists must increase their awareness of LD and, in this issue of Molecular Ecology Resources, Kemppainen et al. ( 2015 ) make software available that can only help in this process: LDna allows patterns of LD in a data set to be explored using tools borrowed from network analysis. This has great potential, but realizing that potential requires understanding LD.     

Indo-Pacific population structure and evolutionary history of the coconut crab Birgus latro

Mitochondrial DNA variation was used to examine population structure in a widespread, marine-dispersed species, Birgus latro . Crabs were collected from eight locations throughout the species' Indo-Pacific distribution. Purified mtDNA from 160 individuals was cut with five restriction enzymes, revealing high haplotype diversity (0.96) and moderate nucleotide diversity (0.75%). Island populations from the Indian Ocean (Christmas I.) and Pacific Ocean were significantly different ( G _ST= 0.37) and had distinct mtDNA lineages with a net sequence divergence of 1.4%. Pacific island populations had diverged in a manner consistent with isolation by distance, with only the most peripheral populations being significantly different. The results for mtDNA are largely concordant with those from allozymes, although estimates of gene flow between the Indian and Pacific Oceans were much lower when based on mtDNA. The mtDNA phylogeny also permitted a deeper examination of the evolutionary and demographic history of Birgus latro . Long-term separation of populations is evident in the complete phylogenetic subdivision of mtDNA lineages between the Indian and Pacific Ocean populations sampled. The starlike phylogeny of alleles from the Pacific suggests a rapid population expansion in the Pacific during the Pleistocene. Including information about allele phylogeny, as well as distribution and frequency, obscured contemporary population structure, but provided unique insights into the evolutionary history of the species.     

EST‐SSR markers from five sequenced cDNA libraries of common bean (Phaseolus vulgaris L.) comparing three bioinformatic algorithms

Expressed sequence tags (ESTs) are a rich source of SSR sequences, but the proportion of long Class I microsatellites with many repeats vs. short Class II microsatellites with few repeats is an important factor to consider. Class I microsatellites, with more than 20 bp of repeats, tend to make better markers with higher polymorphism. The goal of this study was to determine the frequency of Class I and Class II microsatellites in a collection of over 21 000 ESTs from a single study of five different tissues of common bean: two types of leaves, nodules, pods and roots. For this objective, we used three different bioinformatics pipelines: Automated Microsatellite Marker Development (AMMD), Batchprimer3 and SSRLocator. In addition, we determined the frequency of single or multiple SSRs in the assembled ESTs, the frequency of perfect and compound repeats and whether Class I microsatellites were mainly di‐nucleotide or tri‐nucleotide motifs with each of the search engines. Primers were designed for a total of 175 microsatellites concentrating on class I microsatellites identified with SSR locator. A few other microsatellites were included from the other search engines, AMMD and Batchprimer3 programs so as to have a representative set of class II markers for comparison sake. The comparison of 95 class I vs. 80 class II markers confirmed that the Class I were more polymorphic and therefore more useful.     

A novel assessment of population structure and gene flow in grey wolf populations of the Northern Rocky Mountains of the United States

The successful re‐introduction of grey wolves to the western United States is an impressive accomplishment for conservation science. However, the degree to which subpopulations are genetically structured and connected, along with the preservation of genetic variation, is an important concern for the continued viability of the metapopulation. We analysed DNA samples from 555 Northern Rocky Mountain wolves from the three recovery areas (Greater Yellowstone Area, Montana, and Idaho), including all 66 re‐introduced founders, for variation in 26 microsatellite loci over the initial 10‐year recovery period (1995–2004). The population maintained high levels of variation (H_O = 0.64–0.72; allelic diversity k = 7.0–10.3) with low levels of inbreeding (F_IS < 0.03) and throughout this period, the population expanded rapidly (n₁₉₉₅ = 101; n₂₀₀₄ = 846). Individual‐based Bayesian analyses revealed significant population genetic structure and identified three subpopulations coinciding with designated recovery areas. Population assignment and migrant detection were difficult because of the presence of related founders among different recovery areas and required a novel approach to determine genetically effective migration and admixture. However, by combining assignment tests, private alleles, sibship reconstruction, and field observations, we detected genetically effective dispersal among the three recovery areas. Successful conservation of Northern Rocky Mountain wolves will rely on management decisions that promote natural dispersal dynamics and minimize anthropogenic factors that reduce genetic connectivity.     

Nonconcordant evolutionary history of maternal and paternal lineages in Adriatic sturgeon

Although analyses of intraspecific variability are an important prerequisite for species identification assays, only a few studies have focused on population genetics and historical biogeography of sturgeon species. Here we present the first study on genetic variability of the last remaining Adriatic sturgeon, Acipenser naccarii, derived from mitochondrial and nuclear DNA. Our mitochondrial DNA analyses arranged individuals into three distinguished mitochondrial DNA haplogroups (Po1, Po2 and Buna). Two haplogroups (Po1 and Buna) were correlated to geographical distribution, whereas the third (Po2) was not. It was, however, very closely related to one lineage of its Ponto-Caspian sister species, A. gueldenstaedtii. The distribution of nuclear markers (microsatellites and amplified fragment length polymorphism) was strongly correlated to geographical distribution. An assignment test based on nuclear data placed no specimen of A. naccarii to A. gueldenstaedtii and vice versa. Therefore, the presence of gueldenstaedtii-like haplotypes within the Po population is either the result of a postglacial introgression or an ancestral polymorphism and does not indicate a hybrid population. The most valuable tool for forensic species identification purposes is one diagnostic deletion separating all A. naccarii from A. gueldenstaedtii. As both A. naccarii populations are genetically differentiated, stocking of sturgeon from the Po River in Italy into waters of the Buna River would jeopardize the genetic differences between both populations and should thus be avoided.     

Distribution of Y chromosomes among native North Americans: a study of Athapaskan population history

In this study, 231 Y chromosomes from 12 populations were typed for four diagnostic single nucleotide polymorphisms (SNPs) to determine haplogroup membership and 43 Y chromosomes from three of these populations were typed for eight short tandem repeats (STRs) to determine haplotypes. These data were combined with previously published data, amounting to 724 Y chromosomes from 26 populations in North America, and analyzed to investigate the geographic distribution of Y chromosomes among native North Americans and to test the Southern Athapaskan migration hypothesis. The results suggest that European admixture has significantly altered the distribution of Y chromosomes in North America and because of this caution should be taken when inferring prehistoric population events in North America using Y chromosome data alone. However, consistent with studies of other genetic systems, we are still able to identify close relationships among Y chromosomes in Athapaskans from the Subarctic and the Southwest, suggesting that a small number of proto-Apachean migrants from the Subarctic founded the Southwest Athapaskan populations.     

An improved algorithm for the detection of genomic variation using short oligonucleotide expression microarrays

High‐throughput microarray experiments often generate far more biological information than is required to test the experimental hypotheses. Many microarray analyses are considered finished after differential expression and additional analyses are typically not performed, leaving untapped biological information left undiscovered. This is especially true if the microarray experiment is from an ecological study of multiple populations. Comparisons across populations may also contain important genomic polymorphisms, and a subset of these polymorphisms may be identified with microarrays using techniques for the detection of single feature polymorphisms (SFP). SFPs are differences in microarray probe level intensities caused by genetic polymorphisms such as single‐nucleotide polymorphisms and small insertions/deletions and not expression differences. In this study, we provide a new algorithm for the detection of SFPs, evaluate the algorithm using existing data from two publicly available Affymetrix Barley (Hordeum vulgare) microarray data sets and compare them to two previously published SFP detection algorithms. Results show that our algorithm provides more consistent and sensitive calling of SFPs with a lower false discovery rate. Simultaneous analysis of SFPs and differential expression is a low‐cost method for the enhanced analysis of microarray data, enabling additional biological inferences to be made.     

Loss of chromosome Y in primary tumors

1. Download : Download high-res image (162KB)
2. Download : Download full-size image

     

Mitochondrial DNA phylogeography and population history of the grey wolf canis lupus

The grey wolf (Canis lupus) and coyote (C. latrans) are highly mobile carnivores that disperse over great distances in search of territories and mates. Previous genetic studies have shown little geographical structure in either species. However, population genetic structure is also influenced by past isolation events and population fluctuations during glacial periods. In this study, control region sequence data from a worldwide sample of grey wolves and a more limited sample of coyotes were analysed. The results suggest that fluctuating population sizes during the late Pleistocene have left a genetic signature on levels of variation in both species. Genealogical measures of nucleotide diversity suggest that historical population sizes were much larger in both species and grey wolves were more numerous than coyotes. Currently, about 300 000 wolves and 7 million coyotes exist. In grey wolves, genetic diversity is greater than that predicted from census population size, reflecting recent historical population declines. By contrast, nucleotide diversity in coyotes is smaller than that predicted by census population size, reflecting a recent population expansion following the extirpation of wolves from much of North America. Both species show little partitioning of haplotypes on continental or regional scales. However, a statistical parsimony analysis indicates local genetic structure that suggests recent restricted gene flow.