Different filtering strategies of genotyping‐by‐sequencing data provide complementary resolutions of species boundaries and relationships in a clade of sexually deceptive orchids

Ongoing hybridization and retained ancestral polymorphism in rapidly radiating lineages could mask recent cladogenetic events. This presents a challenge for the application of molecular phylogenetic methods to resolve differences between closely related taxa. We reanalyzed published genotyping‐by‐sequencing (GBS) data to infer the phylogeny of four species within the Ophrys sphegodes complex, a recently radiated clade of orchids. We used different data filtering approaches to detect different signals contained in the dataset generated by GBS and estimated their effects on maximum likelihood trees, global F_ST and bootstrap support values. We obtained a maximum likelihood tree with high bootstrap support, separating the species by using a large dataset based on loci shared by at least 30% of accessions. Bootstrap and F_ST values progressively decreased when filtering for loci shared by a higher number of accessions. However, when filtering more stringently to retain homozygous and organellar loci, we identified two main clades. These clades group individuals independently from their a priori species assignment, but were associated with two organellar haplotype clusters. We infer that a less stringent filtering preferentially selects for rapidly evolving lineage‐specific loci, which might better delimit lineages. In contrast, when using homozygous/organellar DNA loci the signature of a putative hybridization event in the lineage prevails over the most recent phylogenetic signal. These results show that using differing filtering strategies on GBS data could dissect the organellar and nuclear DNA phylogenetic signal and yield novel insights into relationships between closely related species.     

Estimating and accounting for genotyping errors in RAD‐seq experiments

In non‐model organisms, evolutionary questions are frequently addressed using reduced representation sequencing techniques due to their low cost, ease of use, and because they do not require genomic resources such as a reference genome. However, evidence is accumulating that such techniques may be affected by specific biases, questioning the accuracy of obtained genotypes, and as a consequence, their usefulness in evolutionary studies. Here, we introduce three strategies to estimate genotyping error rates from such data: through the comparison to high quality genotypes obtained with a different technique, from individual replicates, or from a population sample when assuming Hardy‐Weinberg equilibrium. Applying these strategies to data obtained with Restriction site Associated DNA sequencing (RAD‐seq), arguably the most popular reduced representation sequencing technique, revealed per‐allele genotyping error rates that were much higher than sequencing error rates, particularly at heterozygous sites that were wrongly inferred as homozygous. As we exemplify through the inference of genome‐wide and local ancestry of well characterized hybrids of two Eurasian poplar (Populus) species, such high error rates may lead to wrong biological conclusions. By properly accounting for these error rates in downstream analyses, either by incorporating genotyping errors directly or by recalibrating genotype likelihoods, we were nevertheless able to use the RAD‐seq data to support biologically meaningful and robust inferences of ancestry among Populus hybrids. Based on these findings, we strongly recommend carefully assessing genotyping error rates in reduced representation sequencing experiments, and to properly account for these in downstream analyses, for instance using the tools presented here.     

Genotyping‐by‐sequencing approaches to characterize crop genomes: choosing the right tool for the right application

In the last decade, the revolution in sequencing technologies has deeply impacted crop genotyping practice. New methods allowing rapid, high‐throughput genotyping of entire crop populations have proliferated and opened the door to wider use of molecular tools in plant breeding. These new genotyping‐by‐sequencing (GBS) methods include over a dozen reduced‐representation sequencing (RRS) approaches and at least four whole‐genome resequencing (WGR) approaches. The diversity of methods available, each often producing different types of data at different cost, can make selection of the best‐suited method seem a daunting task. We review the most common genotyping methods used today and compare their suitability for linkage mapping, genomewide association studies (GWAS), marker‐assisted and genomic selection and genome assembly and improvement in crops with various genome sizes and complexity. Furthermore, we give an outline of bioinformatics tools for analysis of genotyping data. WGR is well suited to genotyping biparental cross populations with complex, small‐ to moderate‐sized genomes and provides the lowest cost per marker data point. RRS approaches differ in their suitability for various tasks, but demonstrate similar costs per marker data point. These approaches are generally better suited for de novo applications and more cost‐effective when genotyping populations with large genomes or high heterozygosity. We expect that although RRS approaches will remain the most cost‐effective for some time, WGR will become more widespread for crop genotyping as sequencing costs continue to decrease.     

Estimating admixture pedigrees of recent hybrids without a contiguous reference genome

The genome of recently admixed individuals or hybrids has characteristic genetic patterns that can be used to learn about their recent admixture history. One of these are patterns of interancestry heterozygosity, which can be inferred from SNP data from either called genotypes or genotype likelihoods, without the need for information on genomic location. This makes them applicable to a wide range of data that are often used in evolutionary and conservation genomic studies, such as low-depth sequencing mapped to scaffolds and reduced representation sequencing. Here we implement maximum likelihood estimation of interancestry heterozygosity patterns using two complementary models. We furthermore develop apoh (Admixture Pedigrees of Hybrids), a software that uses estimates of paired ancestry proportions to detect recently admixed individuals or hybrids, and to suggest possible admixture pedigrees. It furthermore calculates several hybrid indices that make it easier to identify and rank possible admixture pedigrees that could give rise to the estimated patterns. We implemented apoh both as a command line tool and as a Graphical User Interface that allows the user to automatically and interactively explore, rank and visualize compatible recent admixture pedigrees, and calculate the different summary indices. We validate the performance of the method using admixed family trios from the 1000 Genomes Project. In addition, we show its applicability on identifying recent hybrids from RAD-seq data of Grant's gazelle (Nanger granti and Nanger petersii) and whole genome low-depth data of waterbuck (Kobus ellipsiprymnus) which shows complex admixture of up to four populations.     

Evidence for natural hybridization between Primula beesiana and P. bulleyana,two heterostylous primroses in NW Yunnan,China

Natural hybridization was assumed to play an essential role for the diversification of Primula; however, only one study of hybridization in the region of the Himalayas has been undertaken. In the present study, we examined another natural hybrid zone where morphologically putative hybrids as well as P. beesiana Forrest, P. bulleyana Forrest, and P. poissonii Franch. co-occurred. We used molecular data to confirm the parental species of putative hybrids and the unidirectional hybridization pattern between P. beesiana and P. bulleyana. Moreover, with reference to synthetic F₁s, most hybrids examined are possibly advanced generations, although the possibility of F₁ hybrids currently examined could not be completely excluded. In addition, pollinator observations on experimental arrays of transplanted parental species showed interspecific pollen flows during visitations of shared pollinators, indicating an incomplete pre-zygotic barrier between P. beesiana and P. bulleyana. Seed productions from both flower morphs of putative hybrids were significantly lower than parental species, suggesting lower reproductive success in these hybrids. Combined with the evidence of recent habitat disturbance in the study area, we might witness the early process of hybridization between P. beesianaand P. bulleyana.     

Filtering for Compound Heterozygous Sequence Variants in Non-Consanguineous Pedigrees

The identification of disease-causing mutations in next-generation sequencing (NGS) data requires efficient filtering techniques. In patients with rare recessive diseases, compound heterozygosity of pathogenic mutations is the most likely inheritance model if the parents are non-consanguineous. We developed a web-based compound heterozygous filter that is suited for data from NGS projects and that is easy to use for non-bioinformaticians. We analyzed the power of compound heterozygous mutation filtering by deriving background distributions for healthy individuals from different ethnicities and studied the effectiveness in trios as well as more complex pedigree structures. While usually more then 30 genes harbor potential compound heterozygotes in single exomes, this number can be markedly reduced with every additional member of the pedigree that is included in the analysis. In a real data set with exomes of four family members, two sisters affected by Mabry syndrome and their healthy parents, the disease-causing gene PIGO, which harbors the pathogenic compound heterozygous variants, could be readily identified. Compound heterozygous filtering is an efficient means to reduce the number of candidate mutations in studies aiming at identifying recessive disease genes in non-consanguineous families. A web-server is provided to make this filtering strategy available at www.gene-talk.de.     

Reduced meiotic fitness in hybrids with heterozygosity for heterochromatin in the speciating<Emphasis Type="Italic">Mus terricolor</Emphasis> complex

Mus terricolor I, II and III are the three chromosomal species which differ in stable autosomal short-arm heterochromatin variations established in homozygous condition. Analysis of meiosis in the laboratory-generated F₁ male hybrids from crosses (both ways) betweenM. terricolor I and II and betweenM. terricolor I and III shows high frequencies of pairing abnormalities at pachytene. The backcross (N3 generation) male hybrids betweenM. terricolor I and II have meiotic abnormalities as in the F₁male hybrids, though to a lesser extent. They show difference in pairing abnormalities in the different karyotypic forms; the backcross hybrids heterozygous for the heterochromatic short arms have more anomalies compared to the homokaryotypic hybrids. This suggests a negative influence of the heterochromatin heterozygosity in meiotic pairing. The results indicate a role for heterochromatin variations in the development of a reproductive barrier in the speciatingM. terricolor complex.     

A Novel Approach to Estimating Heterozygosity from Low-Coverage Genome Sequence

High-throughput shotgun sequence data make it possible in principle to accurately estimate population genetic parameters without confounding by SNP ascertainment bias. One such statistic of interest is the proportion of heterozygous sites within an individual’s genome, which is informative about inbreeding and effective population size. However, in many cases, the available sequence data of an individual are limited to low coverage, preventing the confident calling of genotypes necessary to directly count the proportion of heterozygous sites. Here, we present a method for estimating an individual’s genome-wide rate of heterozygosity from low-coverage sequence data, without an intermediate step that calls genotypes. Our method jointly learns the shared allele distribution between the individual and a panel of other individuals, together with the sequencing error distributions and the reference bias. We show our method works well, first, by its performance on simulated sequence data and, second, on real sequence data where we obtain estimates using low-coverage data consistent with those from higher coverage. We apply our method to obtain estimates of the rate of heterozygosity for 11 humans from diverse worldwide populations and through this analysis reveal the complex dependency of local sequencing coverage on the true underlying heterozygosity, which complicates the estimation of heterozygosity from sequence data. We show how we can use filters to correct for the confounding arising from sequencing depth. We find in practice that ratios of heterozygosity are more interpretable than absolute estimates and show that we obtain excellent conformity of ratios of heterozygosity with previous estimates from higher-coverage data.     

Combining geometric morphometrics with molecular genetics to investigate a putative hybrid complex: a case study with barbels Barbus spp. (Teleostei: Cyprinidae)

This integrative study examined the morphological and genetic affinities of three endemic barbel species from Italy (brook barbel Barbus caninus, Italian barbel Barbus plebejus and horse barbel Barbus tyberinus) and of putative hybrid specimens to their species of origin. Two of the species frequently occur together with the non‐native barbel Barbus barbus. DNA barcoding indicates that mitochondrial (mt) haplotypes often do not match the species expected from morphology. Linear distance measurements and meristics are not informative for discrimination of the species and putative hybrids, but a discriminant analysis of principal components (DAPC) of geometric landmark data produces reassignments largely in congruence with mt and nuclear genetic data. Cyto‐nuclear conflicts confirm the presence of hybridization in B. plebejus and B. tyberinus and identify additional introgressed specimens. A comparison between mixed genotypes and their morphology‐based assignment reveals no predictable pattern. The finding that most individuals of the morphologically similar B. plebejus and B. tyberinus have very high assignment probabilities to their respective species suggests that the presented approach may serve as a valuable tool to distinguish morphologically very similar taxa.     

Management units of the endangered herb Primula sieboldii based on microsatellite variation among and within populations throughout Japan

To promote programs for the conservation and restoration of the endangered species Primula sieboldii, we examined genetic variation at eight microsatellite loci among and within 32 remnant wild populations throughout Japan. Total allelic diversity within a population was higher in larger populations, but not so after rarefaction adjustment. The positive relationship between population size and the inbreeding coefficient may suggest that more heterozygous genets tend to survive the habitat contraction possibly because of the higher fitness associated with heterozygosity. By principal coordinate analysis and Bayesian analysis, we detected four genetic groups (Hokkaido, northern Honshu, central Honshu, and western Japan), which could be recognized as management units of P. sieboldii. If supplementation with plants from other populations were planned, it should be conducted among populations which belong to the same management unit and which are likely to represent the same adaptive variation.     

Complex patterns of differentiation and gene flow underly the divergence of aposematic phenotypes in Oophaga poison frogs

Hybridization and introgression can have complex consequences for both species evolution and conservation. Here, we investigated the origin and characteristics of a putative hybrid zone between two South American poison dart frog species, Oophaga anchicayensis and the critically endangered Oophaga lehmanni, which are heavily sought after on the illegal pet market. Using a combination of phenotypic (49 traits) and genomic (ddRADseq) data, we found that the putative hybrids are morphologically distinct from their parental species and confirmed genomic signatures of admixture in these populations. Several lines of evidence (hybrid indices, interspecific hybrid heterozygosity, genomic clines, comparisons with simulated hybrids and demographic modelling) support the conclusion that these populations are not comprised of early‐generation hybrids and thus, they probably did not arise as a result of illegal translocations associated with wildlife trafficking. Instead, they probably represent an independent lineage which has persisted through isolation and has only relatively recently re‐established gene flow with both parental species. Furthermore, we detected signals of differential introgression from parental species into these hybrid populations which suggest relaxed stabilizing selection on these aposematic colour morphs, potentially via context‐dependent female choice. These populations thus provide a fascinating window into the role of hybridization, isolation and female choice in the diversification of South American poison dart frogs. In addition, our results underline the importance of landscape conservation measures to protect, not only known localities of nominal species, but also the phenotypic and genomic variation harbored by admixed lineages which represent crucial repositories for the impressive diversity in this system.     

崖壁植物太行菊与长裂太行菊全基因组大小及特征分析北大核心CSCD

太行菊(Opisthopappus taihangensis)、长裂太行菊(O.longilobus),为太行山特有多年生崖壁草本植物,菊科(Compositae)重要野生资源,具有较高的经济与生态价值。为确定适合两物种的全基因组测序策略,该研究利用流式细胞法和高通量测序技术,分析两物种基因组大小、杂合率、重复序列及GC含量等信息。结果表明:(1)流式细胞法估算太行菊基因组大小约为2.1 Gb,长裂太行菊基因组大小约为2.4 Gb。(2)高通量测序修正后太行菊基因组大小为3.13 Gb,重复序列比例为84.35%,杂合度为0.99%,GC含量为36.56%;长裂太行菊基因组为3.18 Gb,重复序列比例为83.83%,杂合度为1.17%,GC含量为36.62%。(3)初步组装后GC含量分布及平均深度存在异常,出现分层现象,可能是两物种基因组杂合率较高所致。综上结果表明,太行菊、长裂太行菊均属于高重复、高杂合、大基因组的复杂基因组,建议使用Illumina+PacBio测序组装策略,进行全基因组测序分析。     

An efficient method for screening faecal DNA genotypes and detecting new individuals and hybrids in the red wolf (<Emphasis Type="Italic">Canis rufus</Emphasis>) experimental population area

Previously, sequencing of mitochondrial DNA (mtDNA) from non-invasively collected faecal material (scat) has been used to help manage hybridization in the wild red wolf (Canis rufus) population. This method is limited by the maternal inheritance of mtDNA and the inability to obtain individual identification. Here, we optimize the use of nuclear DNA microsatellite markers on red wolf scat DNA to distinguish between individuals and detect hybrids. We develop a data filtering method in which scat genotypes are compared to known blood genotypes to reduce the number of PCR amplifications needed. We apply our data filtering method and the more conservative maximum likelihood ratio method (MLR) of Miller et al. (2002 Genetics 160:357–366) to a scat dataset previously screened for hybrids by sequencing of mtDNA. Using seven microsatellite loci, we obtained genotypes for 105 scats, which were matched to 17 individuals. The PCR amplification success rate was 50% and genotyping error rates ranged from 6.6% to 52.1% per locus. Our data filtering method produced comparable results to the MLR method, and decreased the time and cost of analysis by 25%. Analysis of this dataset using our data filtering method verified that no hybrid individuals were present in the Alligator River National Wildlife Refuge, North Carolina in 2000. Our results demonstrate that nuclear DNA microsatellite analysis of red wolf scats provides an efficient and accurate approach to screen for new individuals and hybrids.     

Differences in ploidy levels of inter-specific hybrids obtained by reciprocal crosses between Primula sieboldii and P. kisoana

Differences in ploidy level were found in inter-specific hybrids obtained by reciprocal crosses between Primula sieboldii and P. kisoana. When P. sieboldii was used as the maternal parent, the inter-specific hybrids were triploids; when P. kisoana was the maternal parent, the inter-specific hybrids were diploids. The possibility of diploid female gamete formation in P. sieboldii is discussed as a causal factor in the production of triploids occasionally found in crosses between diploids of this species. Received: 28 December 1999 / Accepted: 10 January 2000     

DYNAMICS OF A NOVEL CHROMOSOMAL POLYMORPHISM WITHIN A HYBRID ZONE BETWEEN TWO CHROMOSOME RACES OF THE SCELOPORUS GRAMMICUS COMPLEX (SAURIA,PHRYNOSOMATIDAE)

Several chromosome races of the mesquite lizard, Sceloporus grammicus complex, hybridize at localities in central Mexico. In most cases, the hybridizing populations are delineated by centric fissions at one or more of the macrochromosomes. One notable exception is the Tulancingo hybrid zone between the F5 and FM2 cytotypes. In addition to fission and/or inversion differences at chromosomes 1, 3, 4, and 6, these races differ by a complex rearrangement of chromosome 2, which carries the nucleolus-organizer region in this species. The meiotic consequences of heterozygosity at this chromosome were examined in males to assess the potential for this chromosome to contribute to the dynamics of the hybrid zone. Chromosomal analysis revealed several putative F₁ hybrids and confirmed the production of nonparental chromosomal morphologies through recombination. Pachytene analysis revealed meiotic pairing difficulties associated with chromosome 2 in males heterozygous for the parental chromosomal morphologies. Significant aneuploidy is expected because of random disjunction of the chromosome-2 elements. As a result, these males likely suffer reduced fertiliity and fitness. In contrast, males heterozygous for recombinant chromosomal morphologies displayed low levels of meiotic irregularities and presumably exhibit higher fertility than individuals heterozygous for parental morphologies. It is hypothesized that the recombinant phenotypes facilitate gene flow between the F5 and FM2 cytotypes.     

Hybridization and asymmetric introgression between Cypripedium tibeticum and C. yunnanense in Shangrila County,Yunnan Province,China

Hybridization and introgression are thought to be important for speciation and adaptation in many plants. However, little is known about the hybridization and introgression among Cypripedium species. To investigate the evidence for hybridization and the pattern of introgression between Cypripedium yunnanense and C. tibeticum in Shangrila County, Yunnan Province, China, morphological characters and amplified fragment length polymorphism (AFLP) data for both the species and their putative hybrids were studied. Hand pollination was also performed to verify the crossability of the putative parents. Principal coordinate analysis based on morphological characters and the AFLP data suggested that the putative hybrids were true hybrids of these two Cypripedium species. Analysis with the NewHybrids software indicated that the putative hybrids were F1 generation individuals and backcrosses to C. yunnanese, but no F2 generation was found. Analysis with the Structure software demonstrated asymmetric introgression from C. tibeticum to C. yunnanense. We conclude that natural hybridization and introgression can occur between these two species and that in situ conservation of the parental species is required before fully assessing the evolutionary potential of hybrids.     

Paralogs are revealed by proportion of heterozygotes and deviations in read ratios in genotyping‐by‐sequencing data from natural populations

Whole‐genome duplications have occurred in the recent ancestors of many plants, fish, and amphibians, resulting in a pervasiveness of paralogous loci and the potential for both disomic and tetrasomic inheritance in the same genome. Paralogs can be difficult to reliably genotype and are often excluded from genotyping‐by‐sequencing (GBS) analyses; however, removal requires paralogs to be identified which is difficult without a reference genome. We present a method for identifying paralogs in natural populations by combining two properties of duplicated loci: (i) the expected frequency of heterozygotes exceeds that for singleton loci, and (ii) within heterozygotes, observed read ratios for each allele in GBS data will deviate from the 1:1 expected for singleton (diploid) loci. These deviations are often not apparent within individuals, particularly when sequence coverage is low; but, we postulated that summing allele reads for each locus over all heterozygous individuals in a population would provide sufficient power to detect deviations at those loci. We identified paralogous loci in three species: Chinook salmon (Oncorhynchus tshawytscha) which retains regions with ongoing residual tetrasomy on eight chromosome arms following a recent whole‐genome duplication, mountain barberry (Berberis alpina) which has a large proportion of paralogs that arose through an unknown mechanism, and dusky parrotfish (Scarus niger) which has largely rediploidized following an ancient whole‐genome duplication. Importantly, this approach only requires the genotype and allele‐specific read counts for each individual, information which is readily obtained from most GBS analysis pipelines.     

Polymorphic microsatellite loci in Lithophragma maximum (Saxifragaceae), an endemic plant of San Clemente Island

In this paper, we report the isolation of 11 polymorphic microsatellite loci from Lithophragma maximum (Saxifragaceae), a rare species restricted to San Clemente Island, California. We found moderate levels of allelic variation (mean N_A= 5.5) and high levels of expected heterozygosity (mean H_E= 0.531) across sampled individuals. However, eight loci showed a significant absence of heterozygous individuals, which we attribute to small population size and demographic history of the species.     

Good from far,but far from good: The impact of a reference genome on evolutionary inference

Genomic diversity and past population histories are key considerations in the fields of conservation and evolutionary biology. In this issue of Molecular Ecology Resources, Prasad et al. (Mol. Ecol. Resour., 2021) examine how the quality and phylogenetic divergence of reference genomes influences the outcomes of downstream analyses such as diversity and demographic history inference. Using the beluga whale and rowi kiwi as examples (Figure 1), they systematically estimate heterozygosity, runs of homozygosity (ROH), and demographic history (PSMC) using reference genomes of varying quality and phylogenetic divergence from the target species. They show that demographic history analyses are impacted by phylogenetic distance, although this is not pronounced until divergence exceeds 3% from the target species. Similarly, their results imply that heterozygosity estimates are dependent on phylogenetic distance and the method used to perform the estimates, and ROHs are potentially undetectable when a nonconspecific reference is used. This investigation into the role of divergence and quality of reference genomes highlights the impact and potential biases generated by genome selection on downstream analyses, and provides a possible alternative in cross-species scaffolding in instances where a conspecific reference genome is not available.     

Hybridization between delta smelt and two other species within the family Osmeridae in the San Francisco Bay-Delta

Hybridization among closely related species may pose a threat to species persistence, especially between native and introduced species. We analyzed nine microsatellite loci, mitochondrial sequences and 16 species-specific single nucleotide polymorphisms (SNPs) in two native species (delta smelt and longfin smelt) and one introduced species (wakasagi smelt) in the family Osmeridae to describe the extent of hybridization among these species in the San Francisco Bay-Delta, CA, USA. We identified 29 putative hybrids with a microsatellite-based Bayesian assignment method, and we further screened these putative hybrids with the SNP loci and mitochondrial DNA (mtDNA) sequencing. From the Yolo Bypass, 11 % of morphologically ambiguous individuals were F₁ hybrids and 0.1 % of positively identified delta smelt from throughout the San Francisco Bay-Delta were F₁ hybrids according to their SNP genotypes. mtDNA sequencing revealed wakasagi smelt as the maternal parent for all five delta smelt × wakasagi smelt hybrids and longfin smelt as the maternal parent for the single longfin smelt × delta smelt hybrid. Hybridization among these three species appears to occur at relatively low frequencies and may not be an immediate threat to the persistence of the imperiled native species; however, the presence of hybrid individuals warrants continued monitoring.