Search

Abstract Background Local adaptation is a key evolutionary process that enhances the growth of plants in their native habitat compared to non-native habitats, resulting in patterns of adaptive genetic variation across the entire geographic range of the species. The study of population adaptation to local environments and predicting their response to future climate change is important because of climate change. Results Here, we explored the genetic diversity of candidate genes associated with bud burst in pedunculate oak individuals sampled from 6 populations in Poland. Single nucleotide polymorphism (SNP) diversity was assessed in 720 candidate genes using the sequence capture technique, yielding 18,799 SNPs. Using landscape genomic approaches, we identified 8 F ST outliers and 781 unique SNPs in 389 genes associated with geography, climate, and phenotypic variables (individual/family spring and autumn phenology, family diameter at breast height (DBH), height, and survival) that are potentially involved in local adaptation. Then, using a nonlinear multivariate model, Gradient Forests, we identified vulnerable areas of the pedunculate oak distribution in Poland that are at risk from climate change. Conclusions The model revealed that pedunculate oak populations in the eastern part of the analyzed geographical region are the most sensitive to climate change. Our results might offer an initial evaluation of a potential management strategy for preserving the genetic diversity of pedunculate oak.

Germline DNA testing using the next-generation sequencing (NGS) technology has become the analytical standard for the diagnostics of hereditary diseases, including cancer. Its increasing use places high demands on correct sample identification, independent confirmation of prioritized variants, and their functional and clinical interpretation. To streamline these processes, we introduced parallel DNA and RNA capture-based NGS using identical capture panel CZECANCA, which is routinely used for DNA analysis of hereditary cancer predisposition. Here, we present the analytical workflow for RNA sample processing and its analytical and diagnostic performance. Parallel DNA/RNA analysis allowed credible sample identification by calculating the kinship coefficient. The RNA capture-based approach enriched transcriptional targets for the majority of clinically relevant cancer predisposition genes to a degree that allowed analysis of the effect of identified DNA variants on mRNA processing. By comparing the panel and whole-exome RNA enrichment, we demonstrated that the tissue-specific gene expression pattern is independent of the capture panel. Moreover, technical replicates confirmed high reproducibility of the tested RNA analysis. We concluded that parallel DNA/RNA NGS using the identical gene panel is a robust and cost-effective diagnostic strategy. In our setting, it allows routine analysis of 48 DNA/RNA pairs using NextSeq 500/550 Mid Output Kit v2.5 (150 cycles) in a single run with sufficient coverage to analyse 226 cancer predisposition and candidate genes. This approach can replace laborious Sanger confirmatory sequencing, increase testing turnaround, reduce analysis costs, and improve interpretation of the impact of variants by analysing their effect on mRNA processing. [ABSTRACT FROM AUTHOR]

Despite the economic, ecological, and scientific importance of the genera Salix L. (willows) and Populus L. (poplars, cottonwoods, and aspens) Salicaceae, we know little about the sources of differences in species diversity between the genera and of the phylogenetic conflict that often confounds estimating phylogenetic trees. Salix subgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a "spectacular failure" due to a speculated radiation of the subgenera Vetrix and Chamaetia. Here, we use targeted sequence capture to understand the evolutionary history of this portion of the Salicaceae plant family. Our phylogenetic hypothesis was based on 787 gene regions and identified extensive phylogenetic conflict among genes. Our analysis supported some previously described subgeneric relationships and confirmed the polyphyly of others. Using an f branch analysis, we identified several cases of hybridization in deep branches of the phylogeny, which likely contributed to discordance among gene trees. In addition, we identified a rapid increase in diversification rate near the origination of the Vetrix–Chamaetia clade in Salix. This region of the tree coincided with several nodes that lacked strong statistical support, indicating a possible increase in incomplete lineage sorting due to rapid diversification. The extraordinary level of both recent and ancient hybridization in both Salix and Populus have played important roles in the diversification and diversity in these two genera. [ABSTRACT FROM AUTHOR]

The alpine butterfly genus Parnassius is a popular model group for studying biogeography, evolution, conservation biology, and ecology. Despite its scientific importance, a comprehensive and robust phylogeny of this group is still lacking. In this study, we used an amplicon capture strategy to sequence 144 nuclear protein‐coding genes and complete mitochondrial genomes for 60 Parnassius specimens covering 42 species and all eight subgenera of Parnassius. Our results strongly support the monophyly of the genus and eight subgenera. The relationships among subgenera are robustly resolved as (Sachaia, (Kreizbergia, (Driopa, (Parnassius, (Tadumia, Lingamius), (Kailasius, Koramius))))), which is different from all previous results. Biogeographic and divergence time analyses indicate that the ancestor of Parnassius originated in an area including the Himalayas and Tibetan Plateau (HTP) and Mongolian steppes in the middle Miocene approximately 13.19 Mya. The middle Miocene global cooling event (starting from ~13.9 Mya) probably provided climatic opportunities for the diversification of cold‐adapted Parnassius. The ancestral state reconstruction analyses suggest that the ancestor of Parnassius butterflies most likely lived in a medium elevational area (2000–4000 m) and fed on Papaveraceae plants. The host shift from Papaveraceae to Crassulaceae in the subgenus Parnassius increases the species diversity of this subgenus, concurring with the "escape and radiate" hypothesis. Overall, our work provides valuable nuclear gene and mitochondrial genome data and a robust phylogenetic framework of Parnassius for future studies of the taxonomy, evolution, and ecology of this group. [ABSTRACT FROM AUTHOR]

Described in 2009, the Patch-nosed Salamander (Urspelerpes brucei) is a miniature species of lungless salamander with a geographic range of only ~ 45 km2. This species is endemic to the foothills of the Appalachian Mountains in extreme northeastern Georgia and northwestern South Carolina. The Tugaloo River—a waterway of some 50 m in width that forms the political boundary between the two states—bisects the tiny range of U. brucei and likely acts as a barrier to gene flow. Using RADcap data and a suite of complementary population genomic analyses, we evaluated the role that this river and its tributaries may play in enabling and/or interrupting gene flow among populations of U. brucei, and we investigated patterns of within-population and between-population genetic variation. Our results revealed a general pattern of isolation-by-stream distance and indicated that a population separated by the Tugaloo River is moderately more differentiated than what is explainable by stream distance alone. Unique in both its physiography and geologic history, this region in which U. brucei lives also harbors more than a dozen other species of lungless salamanders. Therefore, the genetic patterns that we have elucidated may have larger implications for differentiation among populations of other species with similar dispersal abilities. [ABSTRACT FROM AUTHOR]

Abstract Effective management of protected species requires information on appropriate evolutionary and geographic population boundaries and knowledge of how the physical environment and life‐history traits combine to shape the population structure and connectivity. Saltwater crocodiles (Crocodylus porosus) are the largest and most widely distributed of living crocodilians, extending from Sri Lanka to Southeast Asia and down to northern Australia. Given the long‐distance movement capabilities reported for C. porosus, management units are hypothesised to be highly connected by migration. However, the magnitude, scale, and consistency of connection across managed populations are not fully understood. Here we used an efficient genotyping method that combines DArTseq and sequence capture to survey ≈3000 high‐quality genome‐wide single nucleotide polymorphisms from 1176 C. porosus sampled across nearly the entire range of the species in Queensland, Australia. We investigated historical and present‐day connectivity patterns using fixation and diversity indices coupled with clustering methods and the spatial distribution of kin pairs. We inferred kinship using forward simulation coupled with a kinship estimation method that is robust to unspecified population structure. The results demonstrated that the C. porosus population has substantial genetic structure with six broad populations correlated with geographical location. The rate of gene flow was highly correlated with spatial distance, with greater differentiation along the east coast compared to the west. Kinship analyses revealed evidence of reproductive philopatry and limited dispersal, with approximately 90% of reported first and second‐degree relatives showing a pairwise distance of

As whole‐genome sequencing has become pervasive, some have suggested that reduced genomic representation approaches, for example, sequence capture, are becoming obsolete. In the present study, we argue that these techniques still provide excellent tools in terms of price and quality of data as well as in their ability to provide markers with specific features, as required, for example, in phylogenomics. A potential drawback of the wide‐scale application of reduced representation approaches could be their drop in efficiency with increasing phylogenetic distance from the reference species. While some studies have focused on the degree and performance of reduced representation techniques in such situations, to our knowledge, none of them evaluated their applicability to inter‐specific hybrids and polyploids. This highlights a significant gap in current knowledge since there is increasing evidence for the frequent occurrence of natural hybrids and polyploids, as well as for the major importance of both phenomena in evolution. The main aim of the present study was to carry out a thorough validation of SEQcap applicability to (1) a set of non‐model taxa with a wide range of phylogenetic relatedness and (2) inter‐specific hybrids of various ploidies and genomic compositions. Considering the latter point, we especially focused on mechanisms causing allelic bias and consequent allelic dropout, as these could have confounding effects with respect to the evolutionary genomic dynamics of hybrids, especially in asexuals, which virtually reproduce as a frozen F1 generation. [ABSTRACT FROM AUTHOR]

Effective management of protected species requires information on appropriate evolutionary and geographic population boundaries and knowledge of how the physical environment and life‐history traits combine to shape the population structure and connectivity. Saltwater crocodiles (Crocodylus porosus) are the largest and most widely distributed of living crocodilians, extending from Sri Lanka to Southeast Asia and down to northern Australia. Given the long‐distance movement capabilities reported for C. porosus, management units are hypothesised to be highly connected by migration. However, the magnitude, scale, and consistency of connection across managed populations are not fully understood. Here we used an efficient genotyping method that combines DArTseq and sequence capture to survey ≈3000 high‐quality genome‐wide single nucleotide polymorphisms from 1176 C. porosus sampled across nearly the entire range of the species in Queensland, Australia. We investigated historical and present‐day connectivity patterns using fixation and diversity indices coupled with clustering methods and the spatial distribution of kin pairs. We inferred kinship using forward simulation coupled with a kinship estimation method that is robust to unspecified population structure. The results demonstrated that the C. porosus population has substantial genetic structure with six broad populations correlated with geographical location. The rate of gene flow was highly correlated with spatial distance, with greater differentiation along the east coast compared to the west. Kinship analyses revealed evidence of reproductive philopatry and limited dispersal, with approximately 90% of reported first and second‐degree relatives showing a pairwise distance of <50 km between sampling locations. Given the limited dispersal, lack of suitable habitat, low densities of crocodiles and the high proportion of immature animals in the population, future management and conservation interventions should be considered at regional and state‐wide scales. [ABSTRACT FROM AUTHOR]

BackgroundPopulus trichocarpa is an important forest tree species for the generation of lignocellulosic ethanol. Understanding the genomic basis of biomass production and chemical composition of wood is fundamental in supporting genetic improvement programs. Considerable variation has been observed in this species for complex traits related to growth, phenology, ecophysiology and wood chemistry. Those traits are influenced by both polygenic control and environmental effects, and their genome architecture and regulation are only partially understood. Genome wide association studies (GWAS) represent an approach to advance that aim using thousands of single nucleotide polymorphisms (SNPs). Genotyping using exome capture methodologies represent an efficient approach to identify specific functional regions of genomes underlying phenotypic variation.ResultsWe identified 813 K SNPs, which were utilized for genotyping 461 P. trichocarpa clones, representing 101 provenances collected from Oregon and Washington, and established in California. A GWAS performed on 20 traits, considering single SNP-marker tests identified a variable number of significant SNPs (p-value

Intra- and interspecific communication is crucial to fitness via its role in facilitating mating, territoriality and defence. Yet, the evolution of animal communication systems is puzzling—how do they originate and change over time? Studying stridulatory morphology provides a tractable opportunity to deduce the origin and diversification of a communication mechanism. Stridulation occurs when two sclerotized structures rub together to produce vibratory and acoustic (vibroacoustic) signals, such as a cricket ‘chirp’. We investigated the evolution of stridulatory mechanisms in the superfamily Coreoidea (Hemiptera: Heteroptera), a group of insects known for elaborate male fighting behaviours and enlarged hindlegs. We surveyed a large sampling of taxa and used a phylogenomic dataset to investigate the evolution of stridulatory mechanisms. We identified four mechanisms, with at least five evolutionary gains. One mechanism, occurring only in male Harmostini (Rhopalidae), is described for the first time. Some stridulatory mechanisms appear to be non-homoplastic apomorphies within Rhopalidae, while others are homoplastic or potentially homoplastic within Coreidae and Alydidae, respectively. We detected no losses of these mechanisms once evolved, suggesting they are adaptive. Our work sets the stage for further behavioural, evolutionary and ecological studies to better understand the context in which these traits evolve and change.

Sapindales is an angiosperm order of high economic and ecological value comprising nine families, c. 479 genera, and c. 6570 species. However, family and subfamily relationships in Sapindales remain unclear, making reconstruction of the order's spatio-temporal and morphological evolution difficult. In this study, we used Angiosperms353 target capture data to generate the most densely sampled phylogenetic trees of Sapindales to date, with 448 samples and c. 85% of genera represented. The percentage of paralogous loci and allele divergence was characterized across the phylogeny, which was time-calibrated using 29 rigorously assessed fossil calibrations. All families were supported as monophyletic. Two core family clades subdivide the order, the first comprising Kirkiaceae, Burseraceae, and Anacardiaceae, the second comprising Simaroubaceae, Meliaceae, and Rutaceae. Kirkiaceae is sister to Burseraceae and Anacardiaceae, and, contrary to current understanding, Simaroubaceae is sister to Meliaceae and Rutaceae. Sapindaceae is placed with Nitrariaceae and Biebersteiniaceae as sister to the core Sapindales families, but the relationships between these families remain unclear, likely due to their rapid and ancient diversification. Sapindales families emerged in rapid succession, coincident with the climatic change of the Mid-Cretaceous Hothouse event. Subfamily and tribal relationships within the major families need revision, particularly in Sapindaceae, Rutaceae and Meliaceae. Much of the difficulty in reconstructing relationships at this level may be caused by the prevalence of paralogous loci, particularly in Meliaceae and Rutaceae, that are likely indicative of ancient gene duplication events such as hybridization and polyploidization playing a role in the evolutionary history of these families. This study provides key insights into factors that may affect phylogenetic reconstructions in Sapindales across multiple scales, and provides a state-of-the-art phylogenetic framework for further research. [ABSTRACT FROM AUTHOR]

Sapindales is an angiosperm order of high economic and ecological value comprising nine families, c. 479 genera, and c. 6570 species. However, family and subfamily relationships in Sapindales remain unclear, making reconstruction of the order’s spatio-temporal and morphological evolution difficult. In this study, we used Angiosperms353 target capture data to generate the most densely sampled phylogenetic trees of Sapindales to date, with 448 samples and c. 85% of genera represented. The percentage of paralogous loci and allele divergence was characterized across the phylogeny, which was time-calibrated using 29 rigorously assessed fossil calibrations. All families were supported as monophyletic. Two core family clades subdivide the order, the first comprising Kirkiaceae, Burseraceae, and Anacardiaceae, the second comprising Simaroubaceae, Meliaceae, and Rutaceae. Kirkiaceae is sister to Burseraceae and Anacardiaceae, and, contrary to current understanding, Simaroubaceae is sister to Meliaceae and Rutaceae. Sapindaceae is placed with Nitrariaceae and Biebersteiniaceae as sister to the core Sapindales families, but the relationships between these families remain unclear, likely due to their rapid and ancient diversification. Sapindales families emerged in rapid succession, coincident with the climatic change of the Mid-Cretaceous Hothouse event. Subfamily and tribal relationships within the major families need revision, particularly in Sapindaceae, Rutaceae and Meliaceae. Much of the difficulty in reconstructing relationships at this level may be caused by the prevalence of paralogous loci, particularly in Meliaceae and Rutaceae, that are likely indicative of ancient gene duplication events such as hybridization and polyploidization playing a role in the evolutionary history of these families. This study provides key insights into factors that may affect phylogenetic reconstructions in Sapindales across multiple scales, and provides a state-of-the-art phylogenetic framework for further research.

The plant family Bignoniaceae is a conspicuous and charismatic element of the tropical flora. The family has a complex taxonomic history, with substantial changes in the classification of the group during the past two centuries. Recent re-classifications at the tribal and generic levels have been largely possible by the availability of molecular phylogenies reconstructed using Sanger sequencing data. However, our complete understanding of the systematics, evolution, and biogeography of the family remains incomplete, especially due to the low resolution and support of different portions of the Bignoniaceae phylogeny. To overcome these limitations and increase the amount of molecular data available for phylogeny reconstruction within this plant family, we developed a bait kit targeting 762 nuclear genes, including 329 genes selected specifically for the Bignoniaceae; 348 genes obtained from the Angiosperms353 with baits designed specifically for the family; and, 85 low-copy genes of known function. On average, 77.4% of the reads mapped to the targets, and 755 genes were obtained per species. After removing genes with putative paralogs, 677 loci were used for phylogenetic analyses. On-target genes were compared and combined in the Exon-Only dataset, and on-target + off-target regions were combined in the Supercontig dataset. We tested the performance of the bait kit at different taxonomic levels, from family to species-level, using 38 specimens of 36 different species of Bignoniaceae, representing: 1) six (out of eight) tribal level-clades (e.g., Bignonieae, Oroxyleae, Tabebuia Alliance, Paleotropical Clade, Tecomeae, and Jacarandeae), only Tourrettieae and Catalpeae were not sampled; 2) all 20 genera of Bignonieae; 3) seven (out of nine) species of Dolichandra (e.g., D. chodatii, D. cynanchoides, D. dentata, D. hispida, D. quadrivalvis, D. uncata, and D. uniguis-cati), only D. steyermarkii and D. unguiculata were not sampled; and 4) three individuals of Dolichandra unguis-cati. Our data reconstructed a well-supported phylogeny of the Bignoniaceae at different taxonomic scales, opening new perspectives for a comprehensive phylogenetic framework for the family as a whole.

Resolving the order of speciation events that occurred in rapid succession is inherently hard and typically requires a phylogenomic approach. A case in point concerns the previously unresolved phylogeny of the three species of banded newt (genus Ommatotriton). We obtain c. 7k nuclear DNA markers using target enrichment by sequence capture and analyze the dataset using maximum likelihood inference of concatenated data with RAxML, summary multi-species coalescent analysis with ASTRAL and Bayesian species tree inference using a diffusion model with SNAPPER, and use TreeMix and PhyloNet to test for interspecific gene flow. All analyses recover three distinct species with no evidence of interspecific gene flow. All analyses retrieved the topology (O. nesterovi, (O. ophryticus, O. vittatus)), with high support. SNAPPER did show the tendency to get stuck in a local optimum, resulting in a different but still highly supported topology. Furthermore, we notice that fewer SNAPPER runs get stuck in a local optimum when we include an outgroup. Therefore, we recommend the exploration of multiple independent runs and the use of an outgroup with this approach. The banded newt radiation illustrates the use of genome-wide data to tackle formerly unresolved phylogenies., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Pelagic seabirds are highly mobile, reducing opportunities for population isolation that might promote differentiation and speciation. At the same time, many seabirds are philopatric, and their tendency to return to their natal islands to breed might reduce gene flow sufficiently to permit local adaptation and differentiation. To test the net impact of these competing processes, estimates of differentiation and gene flow based on comprehensive geographical sampling are required. We leveraged diverse source material to achieve comprehensive geographical sampling in a widespread seabird, the Wedge-tailed Shearwater (Ardenna pacifica). Using data from sequence capture and high-throughput sequencing of 2402 loci containing 20 780 single nucleotide polymorphisms, we tested for population differentiation and gene flow among breeding areas. We found little evidence of deep divergences within A. pacifica but were able to resolve fine-scale differentiation across island groups. This differentiation was sufficient to assign individuals sampled away from breeding areas to their likely source populations. Estimated effective migration surfaces revealed reduced migration between the Indian Ocean and Pacific Ocean, presumably owing to land barriers, and across the equatorial Pacific Ocean, perhaps associated with differences in breeding schedule. Our results reveal that, despite their mobility, pelagic seabirds can exhibit fine-scale population differentiation and reduced gene flow among ocean basins. [ABSTRACT FROM AUTHOR]

Recent phylogenomic analyses within the insect superfamily Coreoidea (Hemiptera: Heteroptera) have begun to challenge previous phylogenetic hypotheses of the Coreidae and Alydidae based on more traditional cladistic and non-cladistic studies. Phylogenomic studies have found the coreid subfamilies Hydarinae and Pseudophloeinae to be more closely related to a potentially paraphyletic Alydidae (an "AHP" clade) in contrast with traditional cladistic studies. However, taxon sampling within these higher-level groups has remained sparse in current phylogenetic analyses, and the taxonomic positions and monophyly of some of these taxa continue to be unclear. Here, we expand upon previous phylogenomic studies using ultraconserved element loci by increasing taxon sampling within the AHP clade. Using concatenation and summary coalescent approaches, we specifically tested previous support for an AHP clade, the paraphyly of Alydidae, the phylogenetic position of Hydarinae, and the monophyly of the two tribes of Pseudophloeinae. Our results robustly support an AHP clade and resolved the position of Hydarinae as the sister group to a clade consisting of a paraphyletic Alydidae and Pseudophloeinae, regardless of analytical method and locus/gene tree filtering strategies we employed. We also found support for the monophyly of the pseudophloeine tribes Clavigrallini and Pseudophloeini, but generic relationships within each of these tribes varied across analyses. We discuss past non-cladistic morphological studies that have suggested the potential for an AHP clade in light of our results, and we highlight further systematic work needed to discern the AHP clade as a morphologically diagnosable group for future re-classification of the Alydidae and Coreidae. [ABSTRACT FROM AUTHOR]

Variation in breeding systems between species of the same taxonomic group complicates the consistent application of species concepts, and perhaps even the logically consistent circumscription of species. Several genera of arid-zone ephemerals in the Angianthus clade (Asteraceae: Gnaphalieae: Gnaphaliinae) contain both outcrossing and non-outcrossing species. The latter are recognised by producing an order of magnitude fewer pollen grains per anther and an often reduced number of corolla lobes, and they are frequently more widespread than are the former. In its current taxonomy, the genus Pogonolepis comprises an otherwise morphologically indistinguishable pair of one outcrossing and one non-outcrossing species. I generated sequence capture data to test the genetic segregation of P. stricta and P. muelleriana and the utility of sequence capture data for species circumscription and diagnostics. Phylogenetic analysis showed the two species to form two specimen clades, supporting the current taxonomy. Contrary to expectations, non-outcrossing P. muelleriana exhibited lower gene concordance, in line with values expected from recombination, as well as higher heterozygosity than its outcrossing sister species. More research on the breeding system and population structure of the two species may be required to explain these results. The genus Pogonolepis comprises an otherwise morphologically indistinguishable pair of one outcrossing and one non-outcrossing species. I generated sequence capture data to test the genetic segregation of P. stricta and P. muelleriana and the utility of sequence capture data for species circumscription and diagnostics. Phylogenetic analysis showed the two species to form two specimen clades, supporting the current taxonomy, but contrary to expectations, non-outcrossing P. muelleriana exhibited lower gene concordance than its outcrossing sister species. [ABSTRACT FROM AUTHOR]

Grasslands are widespread and economically relevant ecosystems at the basis of sustainable roughage production. Plant genetic diversity (PGD; i.e., within‐species diversity) is related to many beneficial effects on the ecosystem functioning of grasslands. The monitoring of PGD in temperate grasslands is complicated by the multiplicity of species present and by a shortage of methods for large‐scale assessments. However, the continuous advancement of high‐throughput DNA sequencing approaches has improved the prospects of broad, multispecies PGD monitoring. Among them, amplicon sequencing stands out as a robust and cost‐effective method. Here, we report a set of 12 multispecies primer pairs that can be used for high‐throughput PGD assessments in multiple grassland plant species. The target loci were selected and tested in two phases: a "discovery phase" based on a sequence capture assay (611 nuclear loci assessed in 16 grassland plant species), which resulted in the selection of 11 loci; and a "validation phase", in which the selected loci were targeted and sequenced using multispecies primers in test populations of Dactylis glomerata L., Lolium perenne L., Festuca pratensis Huds., Trifolium pratense L. and T. repens L. The multispecies amplicons had nucleotide diversities per species from 5.19 × 10−3 to 1.29 × 10−2, which is in the range of flowering‐related genes but slightly lower than pathogen resistance genes. We conclude that the methodology, the DNA sequence resources, and the primer pairs reported in this study provide the basis for large‐scale, multispecies PGD monitoring in grassland plants. [ABSTRACT FROM AUTHOR]

Schistosomatidae Stiles and Hassall 1898 is a medically significant family of digenetic trematodes (Trematoda: Digenea), members of which infect mammals or birds as definitive hosts and aquatic or amphibious gastropods as intermediate hosts. Currently, there are 17 named genera, for many of which evolutionary interrelationships remain unresolved. The lack of a resolved phylogeny has encumbered our understanding of schistosomatid evolution, specifically patterns of host-use and the role of host-switching in diversification. Here, we used targeted sequence capture of ultra-conserved elements (UCEs) from representatives of 13 of the 17 named genera and 11 undescribed lineages that are presumed to represent either novel genera or species to generate a phylogenomic dataset for the estimation of schistosomatid interrelationships. This study represents the largest phylogenetic effort within the Schistosomatidae in both the number of loci and breadth of taxon sampling. We present a near-comprehensive family-level phylogeny providing resolution to several clades of long-standing uncertainty within Schistosomatidae, including resolution for the placement of the North American mammalian schistosomes, implying a second separate capture of mammalian hosts. Additionally, we present evidence for the placement of Macrobilharzia at the base of the Schistosoma + Bivitellobilharzia radiation. Patterns of definitive and intermediate host use and a strong role for intermediate host-switching are discussed relative to schistosomatid diversification. [ABSTRACT FROM AUTHOR]

Abstract Background Diversity among phenology-related genes is predicted to be a contributing factor in local adaptations seen in widely distributed plant species that grow in climatically variable geographic areas, such as forest trees. European beech (Fagus sylvatica L.) is widespread, and is one of the most important broadleaved tree species in Europe; however, its potential for adaptation to climate change is a matter of uncertainty, and little is known about the molecular basis of climate change-relevant traits like bud burst. Results We explored single nucleotide polymorphisms (SNP) at candidate genes related to bud burst in beech individuals sampled across 47 populations from Europe. SNP diversity was monitored for 380 candidate genes using a sequence capture approach, providing 2909 unlinked SNP loci. We used two complementary analytical methods to find loci significantly associated with geographic variables, climatic variables (expressed as principal components), or phenotypic variables (spring and autumn phenology, height, survival). Redundancy analysis (RDA) was used to detect candidate markers across two spatial scales (entire study area and within subregions). We revealed 201 candidate SNPs at the broadest scale, 53.2% of which were associated with phenotypic variables. Additive polygenic scores, which provide a measure of the cumulative signal across significant candidate SNPs, were correlated with a climate variable (first principal component, PC1) related to temperature and precipitation availability, and spring phenology. However, different genotype-environment associations were identified within Southeastern Europe as compared to the entire geographic range of European beech. Conclusions Environmental conditions play important roles as drivers of genetic diversity of phenology-related genes that could influence local adaptation in European beech. Selection in beech favors genotypes with earlier bud burst under warmer and wetter habitats within its range; however, selection pressures may differ across spatial scales.

Resolving intractable phylogenetic relationships often requires simultaneously analyzing a large number of coding and non-coding orthologous loci. To gather both coding and non-coding data, traditional sequence capture methods require custom-designed commercial probes. Here, we present a cost-effective sequence capture method based on homemade probes, to capture thousands of coding and non-coding orthologous loci simultaneously, suitable for all organisms. This approach, called “FLc-Capture,” synthesizes biotinylated full-length cDNAs from mRNA as capture probes, eliminates the need for costly commercial probe design and synthesis. To demonstrate the utility of FLc-Capture, we prepared full-length cDNA probes from mRNA extracted from a common colubrid snake. We performed capture experiments with these homemade cDNA probes and successfully obtained thousands of coding and non-coding genomic loci from 24 Colubridae species and 12 distantly related snake species of other families. The average capture specificity of FLc-Capture across all tested snake species is 35%, similar to the previously published EecSeq method. We constructed two phylogenomic data sets, one including 1,075 coding loci (∼817,000 bp) and the other including 1,948 non-coding loci (∼1,114,000 bp), to study the phylogeny of Colubridae. Both data sets yielded highly similar and well-resolved trees, with 85% of nodes having >95% bootstrap support. Our experimental tests show that FLc-Capture is a flexible, fast, and cost-effective sequence capture approach for simultaneously gathering coding and non-coding phylogenomic data sets to study intractable phylogenetic questions. We hope that this method will serve as a new data collection tool for evolutionary biologists working in the era of phylogenomics.

Fungal dye-decolorizing peroxidases (DyPs) have found applications in the treatment of dye-contaminated industrial wastes or to improve biomass digestibility. Their roles in fungal biology are uncertain, although it has been repeatedly suggested that they could participate in lignin degradation and/or modification. Using a comprehensive set of 162 fully sequenced fungal species, we defined seven distinct fungal DyP clades on basis of a sequence similarity network. Sequences from one of these clades clearly diverged from all others, having on average the lower isoelectric points and hydropathy indices, the highest number of N-glycosylation sites, and N-terminal sequence peptides for secretion. Putative proteins from this clade are absent from brown-rot and ectomycorrhizal species that have lost the capability of degrading lignin enzymatically. They are almost exclusively present in white-rot and other saprotrophic Basidiomycota that digest lignin enzymatically, thus lending support for a specific role of DyPs from this clade in biochemical lignin modification. Additional nearly full-length fungal DyP genes were isolated from the environment by sequence capture by hybridization; they all belonged to the clade of the presumably secreted DyPs and to another related clade. We suggest focusing our attention on the presumably intracellular DyPs from the other clades, which have not been characterized thus far and could represent enzyme proteins with novel catalytic properties. Key points: • A fungal DyP phylogeny delineates seven main sequence clades. • Putative extracellular DyPs form a single clade of Basidiomycota sequences. • Extracellular DyPs are associated to white-rot fungi. [ABSTRACT FROM AUTHOR]

Aim: Specialized life histories can restrict the way that organisms move through the landscape. Some fishes, such as mudskippers in the family Gobiidae, are obligately amphibious and require a terrestrial phase in which adults build burrows and lay eggs in mudflats. To evaluate whether this restricted, out‐of‐water phase might influence population structure in one such mudskipper, we tested whether populations separated by water are more differentiated than those connected along a coastline. Location: South and East China Seas. Taxon: Periophthalmus modestus (family: Gobiidae, subfamily: Oxudercinae) (Cantor 1842). Methods: We collected 236 individuals from 10 locations around the East and South China Seas. We used a targeted capture sequencing method, RADcap, to call SNPs from 365 loci. We used multiple spatial and population genetics analyses and coalescent models to address the contemporary and historic biogeography of the species. Results: The observed levels of genetic variation and population structure show populations of P. modestus are connected across large distances of continuous coastline and small spans of water. Restricted long‐distance dispersal across deep water has resulted in fragmentation of populations. We recovered no signal of isolation by distance, but instead, found genetic differentiation associated with regional groupings separated by large expanses of water. Main conclusions: A limited dispersal phase significantly affects population structure in this amphibious fish. The restrictive amphibious life history may be a key factor in the present‐day distribution and population structure of P. modestus. Interactions between organismal life history and geography are critical for explaining the population genetic structure of a species. Significance Statement: Life histories that restrict individuals to particular habitats during specific life stages can have a drastic effect on population structure. Some fishes, such as mudskippers, spend a large amount of time on terrestrial habitats. We propose that this behaviour might be important in shaping the population dynamics of the species, which is in line with our findings that populations are genetically connected along coastlines, but differentiated across open expanses of water. [ABSTRACT FROM AUTHOR]

The anomaly zone, defined by the presence of gene tree topologies that are more probable than the true species tree, presents a major challenge to the accurate resolution of many parts of the Tree of Life. This discrepancy can result from consecutive rapid speciation events in the species tree. Similar to the problem of long-branch attraction, including more data via loci concatenation will only reinforce the support for the incorrect species tree. Empirical phylogenetic studies often employ coalescent-based species tree methods to avoid the anomaly zone, but to this point these studies have not had a method for providing any direct evidence that the species tree is actually in the anomaly zone. In this study, we use 16 species of lizards in the family Scincidae to investigate whether nodes that are difficult to resolve place the species tree within the anomaly zone. We analyze new phylogenomic data (429 loci), using both concatenation and coalescent-based species tree estimation, to locate conflicting topological signal. We then use the unifying principle of the anomaly zone, together with estimates of ancestral population sizes and species persistence times, to determine whether the observed phylogenetic conflict is a result of the anomaly zone. We identify at least three regions of the Scincidae phylogeny that provide demographic signatures consistent with the anomaly zone, and this new information helps reconcile the phylogenetic conflict in previously published studies on these lizards. The anomaly zone presents a real problem in phylogenetics, and our new framework for identifying anomalous relationships will help empiricists leverage their resources appropriately for investigating and overcoming this challenge.

Massively parallel sequencing has revolutionized many areas of biology, but sequencing large amounts of DNA in many individuals is cost-prohibitive and unnecessary for many studies. Genomic complexity reduction techniques such as sequence capture and restriction enzyme-based methods enable the analysis of many more individuals per unit cost. Despite their utility, current complexity reduction methods have limitations, especially when large numbers of individuals are analyzed. Here we develop a much improved restriction site-associated DNA (RAD) sequencing protocol and a new method called Rapture ( R: AD c APTURE: ). The new RAD protocol improves versatility by separating RAD tag isolation and sequencing library preparation into two distinct steps. This protocol also recovers more unique (nonclonal) RAD fragments, which improves both standard RAD and Rapture analysis. Rapture then uses an in-solution capture of chosen RAD tags to target sequencing reads to desired loci. Rapture combines the benefits of both RAD and sequence capture, i.e., very inexpensive and rapid library preparation for many individuals as well as high specificity in the number and location of genomic loci analyzed. Our results demonstrate that Rapture is a rapid and flexible technology capable of analyzing a very large number of individuals with minimal sequencing and library preparation cost. The methods presented here should improve the efficiency of genetic analysis for many aspects of agricultural, environmental, and biomedical science.

The diversity of avian visual phenotypes provides a framework for studying mechanisms of trait diversification generally, and the evolution of vertebrate vision, specifically. Previous research has focused on opsins, but to fully understand visual adaptation, we must study the complete phototransduction cascade (PTC). Here, we developed a probe set that captures exonic regions of 46 genes representing the PTC and other light responses. For a subset of species, we directly compared gene capture between our probe set and low‐coverage whole genome sequencing (WGS), and we discuss considerations for choosing between these methods. Finally, we developed a unique strategy to avoid chimeric assembly by using "decoy" reference sequences. We successfully captured an average of 64% of our targeted exome in 46 species across 14 orders using the probe set and had similar recovery using the WGS data. Compared to WGS or transcriptomes, our probe set: (1) reduces sequencing requirements by efficiently capturing vision genes, (2) employs a simpler bioinformatic pipeline by limiting required assembly and negating annotation, and (3) eliminates the need for fresh tissues, enabling researchers to leverage existing museum collections. We then utilized our vision exome data to identify positively selected genes in two evolutionary scenarios—evolution of night vision in nocturnal birds and evolution of high‐speed vision specific to manakins (Pipridae). We found parallel positive selection of SLC24A1 in both scenarios, implicating the alteration of rod response kinetics, which could improve color discrimination in dim light conditions and/or facilitate higher temporal resolution. [ABSTRACT FROM AUTHOR]

Background: Madin–Darby canine kidney (MDCK) cells are widely used for vaccine production, however, the safety of MDCK cells needs to be considered seriously because of high tumorigenicity. Micro RNAs (miRNAs) that are involved in the tumorigenicity of MDCK cells have been never been reported. Objective: To reveal the role of miRNA in the tumorigenic phenotype of MDCK cell line. Methods: The miRNA expression profiles of two monoclonal MDCK cells (M09CL and M35CL) with low tumorigenicity and one MDCK cell line (M73P) with high tumorigenicity were characterized and investigated by using small RNA-seq technology. Results: A total of 5 known miRNAs and 5 novel miRNAs were highly expressed in M73P. In addition, 4 known miRNAs and 4 novel miRNAs were highly expressed in M09CL and M35CL. The target genes of the differentially expressed miRNAs were significantly enriched in several biological processes, and the majority of these genes were involved in pathways in cancer and the MAPK signaling pathway. Through interaction analysis, 4 up-regulated miRNAs (cfa-miR-452, cfa-miR-8826, cfa-miR-224, and cfa-miR-2387) and their crucial target genes related to the tumor regulation network were identified. Results indicated these 4 miRNAs might play crucial roles in the tumorigenesis of MDCK cells. Conclusion: Our findings, which were based on the functional prediction of miRNAs and target genes, suggested that miRNAs might influence the tumorigenicity of MDCK cells by regulating target genes. Moreover, the results provided important data for understanding the miRNA-mediated regulatory networks that control the tumorigenicities of MDCK cells. [ABSTRACT FROM AUTHOR]

Sequence capture studies result in rich data sets comprising hundreds to thousands of targeted genomic regions that are superseding Sanger‐based data sets comprised of a few well‐known loci with historical uses in phylogenetics ('legacy loci'). However, integrating sequence capture and Sanger‐based data sets is of interest as legacy loci can include different types of loci (e.g. mitochondrial and nuclear) across a potentially larger sample of species from past studies. Sequence capture data sets include nontargeted sequences, and there has been recent interest in extracting legacy loci from invertebrate data sets. Here, we use published legacy data from leaf‐footed bugs (Hemiptera: Coreoidea) to recover 15 mitochondrial and seven nuclear legacy loci from off‐target sequences in a sequence capture data set, explore approaches to improve legacy locus recovery, and combine these loci with sequence capture data for phylogenetic analysis. Two nuclear loci were determined to already be targeted by sequence capture baits. Most of the remaining loci were successfully recovered from off‐target sequences, but this recovery varied greatly. Additionally, complementing complete mitogenomes with additional reference mitochondrial sequences from a genetic depository did not offer improvement for most of our taxa; however, supplementing these reference sequences with extracted legacy loci offered ≥6% improvement across taxa for a given mitochondrial locus (negligible improvement for nuclear loci). Phylogenetic analysis of legacy and sequence capture data produced a topology generally congruent with recent studies, but support was lower. Thus, future studies may employ the approaches used in this study to integrate legacy data with newly generated sequence capture data sets without added expenses. [ABSTRACT FROM AUTHOR]

Abstract For population genetic studies in nonmodel organisms, it is important to use every single source of genomic information. This paper presents EXFI, a Python pipeline that predicts the splice graph and exon sequences using an assembled transcriptome and raw whole‐genome sequencing reads. The main algorithm uses Bloom filters to remove reads that are not part of the transcriptome, to predict the intron–exon boundaries, to then proceed to call exons from the assembly, and to generate the underlying splice graph. The results are returned in GFA1 format, which encodes both the predicted exon sequences and how they are connected to form transcripts. EXFI is written in Python, tested on Linux platforms, and the source code is available under the MIT License at https://github.com/jlanga/exfi.

An overlooked aspect of current microbiome studies is the role of viruses in human health. Compared to bacterial studies, laboratory and analytical methods to study the entirety of viral communities in clinical samples are rudimentary and need further refinement. In order to address this need, we developed Virobiome-Seq, a sequence capture method and an accompanying bioinformatics analysis pipeline, that identifies viral reads in human samples. Virobiome-Seq is able to enrich for and detect multiple types of viruses in human samples, including novel subtypes that diverge at the sequence level. In addition, Virobiome-Seq is able to detect RNA transcripts from DNA viruses and may provide a sensitive method for detecting viral activity in vivo. Since Virobiome-Seq also yields the viral sequence, it makes it possible to investigate associations between viral genotype and psychiatric illness. In this proof of concept study, we detected HIV1, Torque Teno, Pegi, Herpes and Papilloma virus sequences in Peripheral Blood Mononuclear Cells, plasma and stool samples collected from individuals with psychiatric disorders. We also detected the presence of numerous novel circular RNA viruses but were unable to determine whether these viruses originate from the sample or represent contaminants. Despite this challenge, we demonstrate that our knowledge of viral diversity is incomplete and opportunities for novel virus discovery exist. Virobiome-Seq will enable a more sophisticated analysis of the virome and has the potential of uncovering complex interactions between viral activity and psychiatric disease. [ABSTRACT FROM AUTHOR]

Background: Diversity among phenology-related genes is predicted to be a contributing factor in local adaptations seen in widely distributed plant species that grow in climatically variable geographic areas, such as forest trees. European beech (Fagus sylvatica L.) is widespread, and is one of the most important broadleaved tree species in Europe; however, its potential for adaptation to climate change is a matter of uncertainty, and little is known about the molecular basis of climate change-relevant traits like bud burst. Results: We explored single nucleotide polymorphisms (SNP) at candidate genes related to bud burst in beech individuals sampled across 47 populations from Europe. SNP diversity was monitored for 380 candidate genes using a sequence capture approach, providing 2909 unlinked SNP loci. We used two complementary analytical methods to find loci significantly associated with geographic variables, climatic variables (expressed as principal components), or phenotypic variables (spring and autumn phenology, height, survival). Redundancy analysis (RDA) was used to detect candidate markers across two spatial scales (entire study area and within subregions). We revealed 201 candidate SNPs at the broadest scale, 53.2% of which were associated with phenotypic variables. Additive polygenic scores, which provide a measure of the cumulative signal across significant candidate SNPs, were correlated with a climate variable (first principal component, PC1) related to temperature and precipitation availability, and spring phenology. However, different genotype-environment associations were identified within Southeastern Europe as compared to the entire geographic range of European beech. Conclusions: Environmental conditions play important roles as drivers of genetic diversity of phenology-related genes that could influence local adaptation in European beech. Selection in beech favors genotypes with earlier bud burst under warmer and wetter habitats within its range; however, selection pressures may differ across spatial scales. [ABSTRACT FROM AUTHOR]

[Display omitted] • The >250 species of Delias live throughout Asia, Australia, and Melanesia. • Frequent inter-island dispersal during the Neogene triggered diversification. • Lineages on New Guinea began diversifying ca. 5 Mya as the island's central highlands were forming. • We use our large taxon sample to categorize all described species into species groups. The world's largest butterfly genus Delias , commonly known as Jezebels, comprises ca. 251 species found throughout Asia, Australia, and Melanesia. Most species are endemic to islands in the Indo-Australian Archipelago or to New Guinea and nearby islands in Melanesia, and many species are restricted to montane habitats over 1200 m. We inferred an extensively sampled and well-supported molecular phylogeny of the group to better understand the spatial and temporal dimensions of its diversification. The remarkable diversity of Delias evolved in just ca. 15–16 Myr (crown age). The most recent common ancestor of a clade with most of the species dispersed out of New Guinea ca. 14 Mya, but at least six subsequently diverging lineages dispersed back to the island. Diversification was associated with frequent dispersal of lineages among the islands of the Indo-Australian Archipelago, and the divergence of sister taxa on a single landmass was rare and occurred only on the largest islands, most notably on New Guinea. We conclude that frequent inter-island dispersal during the Neogene—likely facilitated by frequent sea level change—sparked much diversification during that period. Many extant New Guinea lineages started diversifying 5 Mya, suggesting that orogeny facilitated their diversification. Our results largely agree with the most recently proposed species group classification system, and we use our large taxon sample to extend this system to all described species. Finally, we summarize recent insights to speculate how wing pattern evolution, mimicry, and sexual selection might also contribute to these butterflies' rapid speciation and diversification. [ABSTRACT FROM AUTHOR]

Abstract Background Populus trichocarpa is an important forest tree species for the generation of lignocellulosic ethanol. Understanding the genomic basis of biomass production and chemical composition of wood is fundamental in supporting genetic improvement programs. Considerable variation has been observed in this species for complex traits related to growth, phenology, ecophysiology and wood chemistry. Those traits are influenced by both polygenic control and environmental effects, and their genome architecture and regulation are only partially understood. Genome wide association studies (GWAS) represent an approach to advance that aim using thousands of single nucleotide polymorphisms (SNPs). Genotyping using exome capture methodologies represent an efficient approach to identify specific functional regions of genomes underlying phenotypic variation. Results We identified 813 K SNPs, which were utilized for genotyping 461 P. trichocarpa clones, representing 101 provenances collected from Oregon and Washington, and established in California. A GWAS performed on 20 traits, considering single SNP-marker tests identified a variable number of significant SNPs (p-value

Abstract Objective Hybrid zones are geographic regions where genetically distinct taxa interbreed, resulting in offspring of mixed ancestry. California Scrub-Jays (Aphelocoma californica) and Woodhouse’s Scrub-Jays (A. woodhouseii) come into secondary contact and hybridize in western Nevada. Although previous work investigated divergence and gene flow between these species using a handful of microsatellite markers, the hybrid zone has not been studied using genome-scale markers, providing an opportunity to assess genome-wide introgression, test for a genetic basis for ecomorphological traits, and compare these estimates to those derived from microsatellites. Results Using variant sites flanking ultraconserved elements (UCEs), we performed population assignment and quantified hybrid ancestry for 16 individuals across the zone of secondary contact. Our study included 2468 SNPs distributed throughout the genome, allowing discrimination of genetic affinities of hybrid individuals that were similar to estimates from microsatellites. We show a relationship between bill and wing length and the genetic composition of individuals that was not found in prior work using microsatellites, suggesting a genetic basis for these traits. Our analyses demonstrate the utility of UCEs for the analysis of hybrid zones and provide a basis for future studies to identify the genomic architecture of speciation and phenotypic differences between these incipient species.

Abstract Target sequence capture is an efficient technique to enrich specific genomic regions for high‐throughput sequencing in ecological and evolutionary studies. In recent years, many sequence capture approaches have been proposed, but most of them rely on commercial synthetic baits which make the experiment expensive. Here, we present a novel sequence capture approach called AFLP‐based genome sequence capture (AFLP Capture). This method uses the AFLP (amplified fragment length polymorphism) technique to generate homemade capture baits without the need for prior genome information, thus is applicable to any organisms. In this approach, biotinylated AFLP fragments representing a random fraction of the genome are used as baits to capture the homologous fragments from genomic shotgun sequencing libraries. In a trial study, by using AFLP Capture, we successfully obtained 511 orthologous loci (>700,000 bp in total length) from 11 Odorrana species and more than 100,000 single nucleotide polymorphisms (SNPs) in four analyzed individuals of an Odorrana species. This result shows that our method can be used to address questions of various evolutionary depths (from interspecies level to intraspecies level). We also discuss the flexibility in bait preparation and how the sequencing data are analyzed. In summary, AFLP Capture is a rapid and flexible tool and can significantly reduce the experimental cost for phylogenetic studies that require analyzing genome‐scale data (hundreds or thousands of loci).

Norway spruce (Picea abies (L.) Karst.) is a conifer species of substanital economic and ecological importance. In common with most conifers, the P. abies genome is very large (∼20 Gbp) and contains a high fraction of repetitive DNA. The current P. abies genome assembly (v1.0) covers approximately 60% of the total genome size but is highly fragmented, consisting of >10 million scaffolds. The genome annotation contains 66,632 gene models that are at least partially validated (www.congenie.org), however, the fragmented nature of the assembly means that there is currently little information available on how these genes are physically distributed over the 12 P. abies chromosomes. By creating an ultra-dense genetic linkage map, we anchored and ordered scaffolds into linkage groups, which complements the fine-scale information available in assembly contigs. Our ultra-dense haploid consensus genetic map consists of 21,056 markers derived from 14,336 scaffolds that contain 17,079 gene models (25.6% of the validated gene models) that we have anchored to the 12 linkage groups. We used data from three independent component maps, as well as comparisons with previously published Picea maps to evaluate the accuracy and marker ordering of the linkage groups. We demonstrate that approximately 3.8% of the anchored scaffolds and 1.6% of the gene models covered by the consensus map have likely assembly errors as they contain genetic markers that map to different regions within or between linkage groups. We further evaluate the utility of the genetic map for the conifer research community by using an independent data set of unrelated individuals to assess genome-wide variation in genetic diversity using the genomic regions anchored to linkage groups. The results show that our map is sufficiently dense to enable detailed evolutionary analyses across the P. abies genome.

Schistosomatidae Stiles and Hassall 1898 is a medically significant family of digenetic trematodes (Trematoda: Digenea), members of which infect mammals or birds as definitive hosts and aquatic or amphibious gastropods as intermediate hosts. Currently, there are 17 named genera, for many of which evolutionary interrelationships remain unresolved. The lack of a resolved phylogeny has encumbered our understanding of schistosomatid evolution, specifically patterns of host-use and the role of host-switching in diversification. Here, we used targeted sequence capture of ultra-conserved elements (UCEs) from representatives of 13 of the 17 named genera and 11 undescribed lineages that are presumed to represent either novel genera or species to generate a phylogenomic dataset for the estimation of schistosomatid interrelationships. This study represents the largest phylogenetic effort within the Schistosomatidae in both the number of loci and breadth of taxon sampling. We present a near-comprehensive family-level phylogeny providing resolution to several clades of long-standing uncertainty within Schistosomatidae, including resolution for the placement of the North American mammalian schistosomes, implying a second separate capture of mammalian hosts. Additionally, we present evidence for the placement of Macrobilharzia at the base of the Schistosoma + Bivitellobilharzia radiation. Patterns of definitive and intermediate host use and a strong role for intermediate host-switching are discussed relative to schistosomatid diversification.

Genomic resources are under-developed for rodents, including well-studied species such as prairie dogs (Cynomys spp.). We conducted whole-genome resequencing on 10 Gunnison's prairie dogs (C. gunnisoni, GUPD) and identified 12,842,055 high-quality SNPs, from which four sets of bait sequences were created. We designed two sets each (containing either 20k or 60k baits) for projects using contemporary (120 bp baits) or historical (100 bp baits) DNA. These bait sets can be used to study a variety of ecological and evolutionary questions in prairie dogs and other ground squirrel taxa. [ABSTRACT FROM AUTHOR]

Cis-regulatory sequences play a crucial role in regulating gene expression and are evolutionary hot spots that drive phenotypic divergence among organisms. Sequencing some cis-regulatory regions of interest in many different species is common in comparative genetic studies. For nonmodel organisms lacking genomic data, genome walking is often the preferred method for this type of application. However, applying genome walking will be laborious and time-consuming when the number of cis-regulatory regions and species to be analyzed is large. In this study, we propose a novel method called amplicon genome fishing (AGF), which can isolate and sequence cis-regulatory regions of interest for any organism. The main idea of the AGF method is to use fragments amplified from the target cis-regulatory regions as enrichment baits to capture and sequence the whole target cis-regulatory regions from genomic library pools. Unlike genome walking, the AGF method is based on hybridization capture and high-throughput sequencing, which makes this method rapid and efficient for projects where some cis-regulatory regions have to be sequenced for many species. We used human amplicons as capture baits and successfully sequenced five target enhancer regions of Homo sapiens, Mus musculus, Gallus gallus, and Xenopus tropicalis, proving the feasibility and repeatability of AGF. To show the utility of the AGF method in real studies, we used it to sequence the ZRS enhancer, a cis-regulatory region associated with the limb loss of snakes, for twenty-three vertebrate species (includes many limbless species never sequenced before). The newly obtained ZRS sequences provide new perspectives into the relationship between the ZRS enhancer's evolution and limb loss in major tetrapod lineages. [ABSTRACT FROM AUTHOR]

Range fragmentation poses challenges for species persistence over time and can be caused by both historical and contemporary processes. We combined genomic data, phylogeographical model testing and palaeoclimatic niche modelling to infer the evolutionary history of the Pine Barrens tree frog (Hyla andersonii), a seepage bog specialist, in eastern North America to gain a better understanding of the historical context of its fragmented distribution. We sampled H. andersonii populations across the three disjunct regions of the species range: Alabama/Florida (AF), the Carolinas (CL) and New Jersey (NJ). Phylogenetic relationships within H. andersonii were consistent between the nuclear species tree and mitochondrial analyses, indicating divergence between AF and CL/NJ (Atlantic clade) ~0.9 Mya and divergence of the NJ clade ~0.15 Mya. Several predictions of north-eastern expansion along the Atlantic coast were supported by phylogeographical analyses. Model testing using genome-wide single nucleotide polymorphism data and species distribution models both provided evidence for multiple disjunct refugia. This comprehensive phylogeographical study of H. andersonii demonstrates a long history of range fragmentation within an endemic coastal plain species and highlights the influence of historical climate change on the current distribution of species and their genetic diversity. [ABSTRACT FROM AUTHOR]

PREMISE: During plant domestication, traits can be subject to a variety of types of selection, ranging from strong directional selection for traits such as seed or fruit size to diversifying selection for traits like color or flavor. These types of selection interact with other evolutionary processes including genetic bottlenecks and interspecific gene flow to generate different levels of genetic diversity across the genome and at target genes in domesticated lineages, but little is known about the impacts of these processes in perennial fruit crops. METHODS: We used sequence capture by hybridization to examine patterns of diversity at a suite of candidate domestication and anonymous background genes in domesticated apple (Malus ×domestica) in comparison to its wild relatives Malus sieversii and Malus orientalis. RESULTS: We found no change in average diversity at these candidate domestication genes across the three species. However, a subset of the genes did exhibit patterns of very high or very low diversity in M. ×domestica compared to its progenitor, M. sieversii. Of the genes with characterized function, the low‐diversity genes mainly contributed to fruit quality traits like color and flavor, predicted to be under conscious, directional selection relatively late in the domestication process, while the high‐diversity genes included a variety of functions. CONCLUSIONS: Overall, these results are consistent with predictions based on the likely timing and nature of selection during domestication and open new avenues for understanding genes with high diversity in a perennial crop compared to its wild relatives. [ABSTRACT FROM AUTHOR]

Many wheat varieties have the potential to develop unacceptably high levels of α-amylase in the grains if exposed to a cool temperature shock or simply cool temperature during the early to middle stages of grain filling. This phenomenon is referred to as late maturity α-amylase (LMA). The enzyme persists in the grain until harvest and may result in wheat with a low Falling Number that does not meet receival and export specifications. Resistance to LMA is therefore a valuable target for wheat breeders and wheat industries in general. Genetic evidence implicating a locus on the long arm of chromosome 7B in variation in LMA phenotype was confirmed in this investigation. Through intensive fine-mapping an ent-copalyl diphosphate synthase (CPS), hitherto named LMA-1, was identified as the likely candidate gene associated with variation in LMA phenotype. Single Nucleotide Polymorphisms (SNPs) within the LMA-1 coding sequence of Chinese Spring, Maringa and Halberd result in either prematurely terminated or functionally altered proteins that are associated with useful levels of resistance to LMA. LMA-1 transcripts detected in de-embryonated grain tissue from around 15 days after anthesis, several days before the synthesis of α-amylase, were low in the resistant varieties Chinese Spring and Maringa compared with LMA susceptible genotype Spica. This was associated with a dramatic reduction in the concentrations of intermediates in the gibberellin biosynthesis pathway such as GA19, evidence that LMA-1 was functioning as CPS in the gibberellin biosynthesis pathway. A survey of a large collection of Australian and international wheat varieties distinguished 9 major haplotypes at the LMA-1 locus. Generally, within classes, there was notable variation for LMA phenotype and evidence for genotypes whose resistance is presumed to be due to genetic loci located elsewhere on the wheat genome. Further investigation is required to characterize the sequence of steps between LMA-1 and α-amylase synthesis as well as to gain a better understanding of the role and potential impact of other genetic loci. Diagnostic markers for sources of resistance and SNP variation reported in this study should assist breeders to deploy resistance associated with LMA-1 variants in breeding programs.

Many wheat varieties have the potential to develop unacceptably high levels of α-amylase in the grains if exposed to a cool temperature shock or simply cool temperature during the early to middle stages of grain filling. This phenomenon is referred to as late maturity α-amylase (LMA). The enzyme persists in the grain until harvest and may result in wheat with a low Falling Number that does not meet receival and export specifications. Resistance to LMA is therefore a valuable target for wheat breeders and wheat industries in general. Genetic evidence implicating a locus on the long arm of chromosome 7B in variation in LMA phenotype was confirmed in this investigation. Through intensive fine-mapping an ent -copalyl diphosphate synthase (CPS), hitherto named LMA-1 , was identified as the likely candidate gene associated with variation in LMA phenotype. Single Nucleotide Polymorphisms (SNPs) within the LMA-1 coding sequence of Chinese Spring, Maringa and Halberd result in either prematurely terminated or functionally altered proteins that are associated with useful levels of resistance to LMA. LMA-1 transcripts detected in de-embryonated grain tissue from around 15 days after anthesis, several days before the synthesis of α-amylase, were low in the resistant varieties Chinese Spring and Maringa compared with LMA susceptible genotype Spica. This was associated with a dramatic reduction in the concentrations of intermediates in the gibberellin biosynthesis pathway such as GA19, evidence that LMA-1 was functioning as CPS in the gibberellin biosynthesis pathway. A survey of a large collection of Australian and international wheat varieties distinguished 9 major haplotypes at the LMA-1 locus. Generally, within classes, there was notable variation for LMA phenotype and evidence for genotypes whose resistance is presumed to be due to genetic loci located elsewhere on the wheat genome. Further investigation is required to characterize the sequence of steps between LMA-1 and α-amylase synthesis as well as to gain a better understanding of the role and potential impact of other genetic loci. Diagnostic markers for sources of resistance and SNP variation reported in this study should assist breeders to deploy resistance associated with LMA-1 variants in breeding programs. [ABSTRACT FROM AUTHOR]

With its large proportion of endemic taxa, complex geological past, and location at the confluence of the highly diverse Malesian and Australian floristic regions, Papuasia – the floristic region comprising the Bismarck Archipelago, New Guinea, and the Solomon Islands – represents an ideal natural experiment in plant biogeography. However, scattered knowledge of its flora and limited representation in herbaria have hindered our understanding of the drivers of its diversity. Focusing on the woody angiosperm genus Schefflera (Araliaceae), we ask whether its morphologically defined infrageneric groupings are monophyletic, when these lineages diverged, and where (within Papuasia or elsewhere) they diversified. To address these questions, we use a high-throughput sequencing approach (Hyb-Seq) which combines target capture (with an angiosperm-wide bait kit targeting 353 single-copy nuclear loci) and genome shotgun sequencing (which allows retrieval of regions in high-copy number, e.g., organellar DNA) of historical herbarium collections. To reconstruct the evolutionary history of the genus we reconstruct molecular phylogenies with Bayesian inference, maximum likelihood, and pseudo-coalescent approaches, and co-estimate divergence times and ancestral areas in a Bayesian framework. We find strong support for most infrageneric morphological groupings, as currently circumscribed, and we show the efficacy of the Angiosperms-353 probe kit in resolving both deep and shallow phylogenetic relationships. We infer a sequence of colonization to explain the present-day distribution of Schefflera in Papuasia: from the Sunda Shelf, Schefflera arrived to the Woodlark plate (present-day eastern New Guinea) in the late Oligocene (when most of New Guinea was submerged) and, subsequently (throughout the Miocene), it migrated westwards (to the Maoke and Bird’s Head Plates and thereon) and further diversified, in agreement with previous reconstructions.

Abstract Background Bread wheat has a large complex genome that makes whole genome resequencing costly. Therefore, genome complexity reduction techniques such as sequence capture make re-sequencing cost effective. With a high-quality draft wheat genome now available it is possible to design capture probe sets and to use them to accurately genotype and anchor SNPs to the genome. Furthermore, in addition to genetic variation, epigenetic variation provides a source of natural variation contributing to changes in gene expression and phenotype that can be profiled at the base pair level using sequence capture coupled with bisulphite treatment. Here, we present a new 12 Mbp wheat capture probe set, that allows both the profiling of genotype and methylation from the same DNA sample. Furthermore, we present a method, based on Agilent SureSelect Methyl-Seq, that will use a single capture assay as a starting point to allow both DNA sequencing and methyl-seq. Results Our method uses a single capture assay that is sequentially split and used for both DNA sequencing and methyl-seq. The resultant genotype and epi-type data is highly comparable in terms of coverage and SNP/methylation site identification to that generated from separate captures for DNA sequencing and methyl-seq. Furthermore, by defining SNP frequencies in a diverse landrace from the Watkins collection we highlight the importance of having genotype data to prevent false positive methylation calls. Finally, we present the design of a new 12 Mbp wheat capture and demonstrate its successful application to re-sequence wheat. Conclusions We present a cost-effective method for performing both DNA sequencing and methyl-seq from a single capture reaction thus reducing reagent costs, sample preparation time and DNA requirements for these complementary analyses.

Abstract Pathogens are recognized as major drivers of local adaptation in wildlife systems. By determining which gene variants are favored in local interactions among populations with and without disease, spatially explicit adaptive responses to pathogens can be elucidated. Much of our current understanding of host responses to disease comes from a small number of genes associated with an immune response. High‐throughput sequencing (HTS) technologies, such as genotype‐by‐sequencing (GBS), facilitate expanded explorations of genomic variation among populations. Hybridization‐based GBS techniques can be leveraged in systems not well characterized for specific variants associated with disease outcome to “capture” specific genes and regulatory regions known to influence expression and disease outcome. We developed a multiplexed, sequence capture assay for red foxes to simultaneously assess ~300‐kbp of genomic sequence from 116 adaptive, intrinsic, and innate immunity genes of predicted adaptive significance and their putative upstream regulatory regions along with 23 neutral microsatellite regions to control for demographic effects. The assay was applied to 45 fox DNA samples from Alaska, where three arctic rabies strains are geographically restricted and endemic to coastal tundra regions, yet absent from the boreal interior. The assay provided 61.5% on‐target enrichment with relatively even sequence coverage across all targeted loci and samples (mean = 50×), which allowed us to elucidate genetic variation across introns, exons, and potential regulatory regions (4,819 SNPs). Challenges remained in accurately describing microsatellite variation using this technique; however, longer‐read HTS technologies should overcome these issues. We used these data to conduct preliminary analyses and detected genetic structure in a subset of red fox immune‐related genes between regions with and without endemic arctic rabies. This assay provides a template to assess immunogenetic variation in wildlife disease systems.