Reduced Representation Libraries Research Articles

Abstract 7022: Whole genome and reduced representation enzymatic methyl-seq enable cost effective methylomes

Abstract DNA methylation is an epigenetic regulator of gene expression with important functions in development and diseases such as cancer. Traditionally, sodium bisulfite conversion was used to distinguish 5-methylcytosines (5mC) and 5-hydroxymethylcytosines (5hmC) from cytosines. However, this method damages DNA and introduces significant sequencing bias. NEBNext® Enzymatic Methyl-seq (EM-seq™) is an enzymatic conversion approach that minimizes DNA damage therefore enabling longer insert sizes, lower duplication rates and a more accurate quantification of methylation in DNA samples with inputs ranging from 10 ng to 200 ng. An improved enzymatic conversion method, EM-seq v2, has a more streamlined workflow and an expanded DNA input range of 100 pg to 200 ng. Libraries made using NA12878 genomic DNA show improved yields, consistent global methylation and an even GC bias representation. The increased library yields subsequently reduces sequencing duplications and costs. Sequencing costs are an important consideration when performing whole genome methylome analysis as this often needs high sequencing depth to attain meaningful coverage to accurately call methylation. Reduced representation libraries, a commonly used method to investigate DNA methylation, relies on enrichment of CpGs within libraries digested using MspI. These CpGs are commonly found in regulatory regions. Although lacking comprehensive CpG coverage, for many researchers this method provides enough cost effective CpG methylation coverage for their studies. Reduced Representation EM-seq libraries (RREM-seq) and Reduced Representation Bisulfite libraries (RRBS) were made using 200 ng to 100 pg of NA12878 DNA. EM-seq adaptors were ligated to MspI digested DNA followed by conversion using either EM-seq v2 or sodium bisulfite. Libraries were PCR amplified and sequenced on an Illumina NovaSeq platform. RREM-seq resulted in higher yields and reproducible results regardless of input amounts compared to RRBS libraries. RREM-seq libraries also identified a higher number of CpGs with even coverage facilitating accurate methylation calling. RREM-seq libraries have superior sequencing metrics resulting in robust methylation profiling. Focusing on a subset of the genome generates higher coverage DNA methylation data at a lower DNA sequencing cost. Citation Format: Vaishnavi Panchapakesa, Chaithanya Ponnaluri, Daniel Evanich, Ariel Erijman, Bradley Langhorst, Louise Williams. Whole genome and reduced representation enzymatic methyl-seq enable cost effective methylomes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7022.

Genotyping by sequencing for the construction of oil palm (Elaeis guineensis Jacq.) genetic linkage map and mapping of yield related quantitative trait loci.

Oil palm (Elaeis guineensis Jacq.) is one of the major oil-producing crops. Improving the quality and increasing the production yield of oil palm have been the primary focuses of both conventional and modern breeding approaches. However, the conventional breeding approach for oil palm is very challenging due to its longevity, which results in a long breeding cycle. Thus, the establishment of marker assisted selection (MAS) for oil palm breeding programs would speed up the breeding pipeline by generating new oil palm varieties that possess high commercial traits. With the decreasing cost of sequencing, Genotyping-by-sequencing (GBS) is currently feasible to many researchers and it provides a platform to accelerate the discovery of single nucleotide polymorphism (SNP) as well as insertion and deletion (InDel) markers for the construction of a genetic linkage map. A genetic linkage map facilitates the identification of significant DNA regions associated with the trait of interest via quantitative trait loci (QTL) analysis. A mapping population of 112 F1 individuals from a cross of Deli dura and Serdang pisifera was used in this study. GBS libraries were constructed using the double digestion method with HindIII and TaqI enzymes. Reduced representation libraries (RRL) of 112 F1 progeny and their parents were sequenced and the reads were mapped against the E. guineensis reference genome. To construct the oil palm genetic linkage map, informative SNP and InDel markers were used to discover significant DNA regions associated with the traits of interest. The nine traits of interest in this study were fresh fruit bunch (FFB) yield, oil yield (OY), oil to bunch ratio (O/B), oil to dry mesocarp ratio (O/DM) ratio, oil to wet mesocarp ratio (O/WM), mesocarp to fruit ratio (M/F), kernel to fruit ratio (K/F), shell to fruit ratio (S/F), and fruit to bunch ratio (F/B). A total of 2.5 million SNP and 153,547 InDel markers were identified. However, only a subset of 5,278 markers comprising of 4,838 SNPs and 440 InDels were informative for the construction of a genetic linkage map. Sixteen linkage groups were produced, spanning 2,737.6 cM for the maternal map and 4,571.6 cM for the paternal map, with average marker densities of one marker per 2.9 cM and one per 2.0 cM respectively, were produced. A QTL analysis was performed on nine traits; however, only QTL regions linked to M/F, K/F and S/F were declared to be significant. Of those QTLs were detected: two for M/F, four for K/F and one for S/F. These QTLs explained 18.1-25.6% of the phenotypic variance and were located near putative genes, such as casein kinase II and the zinc finger CCCH domain, which are involved in seed germination and growth. The identified QTL regions for M/F, K/F and S/F from this study could be applied in an oil palm breeding program and used to screen palms with desired traits via marker assisted selection (MAS).

Environment predicts the maintenance of reproductive isolation in a mosaic hybrid zone of rubber rabbitbrush.

Widely distributed plants of western North America experience divergent selection across environmental gradients, have complex histories shaped by biogeographic barriers and distributional shifts and often illustrate continuums of reproductive isolation. Rubber rabbitbrush (Ericameria nauseosa) is a foundational shrub species that occurs across diverse environments of western North America. Its remarkable phenotypic diversity is currently ascribed to two subspecies-Ericameria nauseosa nauseosa and Ericameria nauseosa consimilis-and 22 named varieties. To understand how genetic variation is partitioned across subspecies, varieties, and environments, we used high throughput sequencing of reduced representation libraries. We found clear evidence for divergence between the two subspecies, despite largely sympatric distributions. Numerous locations exhibiting admixed ancestry were not geographically localized but were widely distributed across a mosaic hybrid zone. The occurrence of hybrid and subspecific ancestries was strongly predicted by environmental variables as well as the proximity to major ecotones between ecoregions. Although this repeatability illustrates the importance of environmental factors in shaping reproductive isolation, variability in the prevalence of hybridization also indicates these factors likely differ across ecological contexts. There was mixed evidence for the evolutionary cohesiveness of varieties, but several genetically distinct and narrow endemic varieties exhibited admixed subspecific ancestries, hinting at the possibility for transgressive hybridization to contribute to phenotypic novelty and the colonization of new environments in E. nauseosa.

Reduced-representation libraries in insect genetics.

Genotyping-by-sequencing of reduced-representation libraries has ushered in an era where genome-wide data can be gotten for any species. Here, I review research on this topic during the last two years, report meta-analysis of the results, and discuss analysis methods and issues. Scanning the literature from 2021 to 2022 identified 21 papers, the majority of which were on population differences, including local adaptation and migration, but several papers were on genetic maps and their use in assembly scaffolding or analysis of quantitative trait loci, on the origin of incursions of pest insects, or on infection rates of a pathogen in a disease vector. The research reviewed includes 33 species from 25 families and 11 orders. Meta-analysis showed that less than 16%, and most often, less than 1% of the genome was implicated in local adaptation and that the number of adaptive loci correlated with genetic divergence among populations.

Development and application of Single Primer Enrichment Technology (SPET) SNP assay for population genomics analysis and candidate gene discovery in lettuce.

Single primer enrichment technology (SPET) is a novel high-throughput genotyping method based on short-read sequencing of specific genomic regions harboring polymorphisms. SPET provides an efficient and reproducible method for genotyping target loci, overcoming the limits associated with other reduced representation library sequencing methods that are based on a random sampling of genomic loci. The possibility to sequence regions surrounding a target SNP allows the discovery of thousands of closely linked, novel SNPs. In this work, we report the design and application of the first SPET panel in lettuce, consisting of 41,547 probes spanning the whole genome and designed to target both coding (~96%) and intergenic (~4%) regions. A total of 81,531 SNPs were surveyed in 160 lettuce accessions originating from a total of 10 countries in Europe, America, and Asia and representing 10 horticultural types. Model ancestry population structure clearly separated the cultivated accessions (Lactuca sativa) from accessions of its presumed wild progenitor (L. serriola), revealing a total of six genetic subgroups that reflected a differentiation based on cultivar typology. Phylogenetic relationships and principal component analysis revealed a clustering of butterhead types and a general differentiation between germplasm originating from Western and Eastern Europe. To determine the potentiality of SPET for gene discovery, we performed genome-wide association analysis for main agricultural traits in L. sativa using six models (GLM naive, MLM, MLMM, CMLM, FarmCPU, and BLINK) to compare their strength and power for association detection. Robust associations were detected for seed color on chromosome 7 at 50 Mbp. Colocalization of association signals was found for outer leaf color and leaf anthocyanin content on chromosome 9 at 152 Mbp and on chromosome 5 at 86 Mbp. The association for bolting time was detected with the GLM, BLINK, and FarmCPU models on chromosome 7 at 164 Mbp. Associations were detected in chromosomal regions previously reported to harbor candidate genes for these traits, thus confirming the effectiveness of SPET for GWAS. Our findings illustrated the strength of SPET for discovering thousands of variable sites toward the dissection of the genomic diversity of germplasm collections, thus allowing a better characterization of lettuce collections.

High‐throughput genotyping in estimating genetic resources and detecting pathogens in aquaculture

High‐throughput genotyping in estimating genetic resources and detecting pathogens in aquaculture

Next Generation Sequencing Revolutionizes Organismal Biology Research in Bats.

The advent of next generation sequencing technologies (NGS) has greatly accelerated our understanding of critical aspects of organismal biology from non-model organisms. Bats form a particularly interesting group in this regard, as genomic data have helped unearth a vast spectrum of idiosyncrasies in bat genomes associated with bat biology, physiology, and evolution. Bats are important bioindicators and are keystone species to many eco-systems. They often live in proximity to humans and are frequently associated with emerging infectious diseases, including the COVID-19 pandemic. Nearly four dozen bat genomes have been published to date, ranging from drafts to chromosomal level assemblies. Genomic investigations in bats have also become critical towards our understanding of disease biology and host-pathogen coevolution. In addition to whole genome sequencing, low coverage genomic data like reduced representation libraries, resequencing data, etc. have contributed significantly towards our understanding of the evolution of natural populations, and their responses to climatic and anthropogenic perturbations. In this review, we discuss how genomic data have enhanced our understanding of physiological adaptations in bats (particularly related to ageing, immunity, diet, etc.), pathogen discovery, and host pathogen co-evolution. In comparison, the application of NGS towards population genomics, conservation, biodiversity assessment, and functional genomics has been appreciably slower. We reviewed the current areas of focus, identifying emerging topical research directions and providing a roadmap for future genomic studies in bats.

SAG-RAD: A Method for Single-Cell Population Genomics of Unicellular Eukaryotes

Sequencing of reduced representation libraries enables genotyping of many individuals for population genomic studies. However, high amounts of DNA are required, and the method cannot be applied directly on single cells, preventing its use on most microbes. We developed and implemented the analysis of single amplified genomes followed by restriction-site-associated DNA sequencing to bypass labor-intensive culturing and to avoid culturing bias in population genomic studies of unicellular eukaryotes. This method thus opens the way for addressing important questions about the genetic diversity, gene flow, adaptation, dispersal, and biogeography of hitherto unexplored species.

Molecular ecology of the sleeper shark subgenus Somniosus (Somniosus) reveals genetic homogeneity within species and lack of support for S. antarcticus.

Inferences made from molecular data support regional stock assessment goals by providing insights into the genetic population dynamics of enigmatic species. Population genomics metrics, such as genetic diversity and population connectivity, serve as useful proxies for species health and stability. Sleeper sharks (genus Somniosus) are ecologically important deep-sea predators, estimated to reach ages of 250 to 300 yr and taking decades to reach sexual maturity. The subgenus Somniosus (Somniosus) is comprised of 3 species: S. pacificus, S. microcephalus, and S. antarcticus. Given the life history strategy of somniosids, they are vulnerable to overfishing and population declines. Further, data to assess the stocks of these species are limited. To address this deficiency, we used the reduced representation library method Restriction-site Associated DNA sequencing (RADseq) to conduct phylogenomic and population genomics analyses, providing novel information for use in stock assessments. Our results strongly support the species status of S. microcephalus (N = 79), but recover S. antarcticus (N = 2) intermixed within the S. pacificus (N = 170) clade. Population genomics analyses reveal genetic homogeneity within S. pacificus and S. microcephalus, and estimates of effective population size were in the hundreds for both species. Kinship analysis identified 2 first-degree relative pairs within our dataset (1 within each species). Our results contribute new information for stock assessments of these uniquely long-lived species by providing the strongest molecular evidence to date for the synonymization of S. antarcticus and S. pacificus, as well as estimating population genomic metrics for each supported species within the Somniosus (Somniosus) subgenus.

Predictors of sequence capture in a large-scale anchored phylogenomics project

Next-generation sequencing (NGS) technologies have revolutionized phylogenomics by decreasing the cost and time required to generate sequence data from multiple markers or whole genomes. Further, the fragmented DNA of biological specimens collected decades ago can be sequenced with NGS, reducing the need for collecting fresh specimens. Sequence capture, also known as anchored hybrid enrichment, is a method to produce reduced representation libraries for NGS sequencing. The technique uses single-stranded oligonucleotide probes that hybridize with pre-selected regions of the genome that are sequenced via NGS, culminating in a dataset of numerous orthologous loci from multiple taxa. Phylogenetic analyses using these sequences have the potential to resolve deep and shallow phylogenetic relationships. Identifying the factors that affect sequence capture success could save time, money, and valuable specimens that might be destructively sampled despite low likelihood of sequencing success. We investigated the impacts of specimen age, preservation method, and DNA concentration on sequence capture (number of captured sequences and sequence quality) while accounting for taxonomy and extracted tissue type in a large-scale butterfly phylogenomics project. This project used two probe sets to extract 391 loci or a subset of 13 loci from over 6,000 butterfly specimens. We found that sequence capture is a resilient method capable of amplifying loci in samples of varying age (0–111 years), preservation method (alcohol, papered, pinned), and DNA concentration (0.020 ng/μl - 316 ng/ul). Regression analyses demonstrate that sequence capture is positively correlated with DNA concentration. However, sequence capture and DNA concentration are negatively correlated with sample age and preservation method. Our findings suggest that sequence capture projects should prioritize the use of alcohol-preserved samples younger than 20 years old when available. In the absence of such specimens, dried samples of any age can yield sequence data, albeit with returns that diminish with increasing age.

Identification of late blight resistance quantitative trait loci in Solanum pimpinellifolium accession PI 270441.

Late blight (LB), caused by the oomycete Phytophthora infestans, is one of the most destructive diseases of the cultivated tomato (Solanum lycopersicum L.) and potato (Solanum tuberosum L.) worldwide. Genetic changes in the pathogen have resulted in the emergence of new genotypes, overcoming formerly effective fungicides or host resistance genes. We previously reported the identification of a LB-resistant accession (PI 270441) of the wild tomato species S. pimpinellifolium L. and the high heritability of its resistance. In the present study, an F2 population (n = 1,209), derived from a cross between PI 270441 and a LB-susceptible tomato breeding line (Fla. 8059), was screened for response to LB infection. Extreme resistant (n = 44) and susceptible (n = 39) F2 individuals were selected and used in a trait-based marker analysis (TBA; a.k.a selective genotyping) to identify and map quantitative trait loci (QTLs) conferring LB resistance. Reduced representation libraries (RRLs) of Fla. 8059 and PI 270441 were constructed, sequenced, and mapped to the tomato genome. A total of 13,054 single-nucleotide polymorphisms (SNPs) were identified, of which, 200 were used to construct a genetic linkage map and locate QTLs. Four LB resistance QTLs were identified on chromosomes 1, 10, and 11 of PI 270441. The markers associated with these QTLs can be used to transfer LB resistance from PI 270441 into new tomato cultivars and to develop near-isogenic lines for fine mapping of the QTL.

Editorial 2022

Editorial 2022

Genome Wide SNP Identification in <i>Eucalyptus</i> and its Application for Linkage Mapping and Adventitious Rooting Related QTL Analysis

Eucalyptus is a versatile genus with a particular significance for paper pulp, and it is extensively cultivated in India's marginal lands. Eucalyptus camaldulensis and E. tereticornis are highly preferred due to its optimal pulp yield, excellent fibre quality and good abiotic stress tolerance properties. Marker-assisted selections have played a major role in eucalyptus breeding programs. The use of genotyping-by-sequencing (GBS) approach for a largescale single nucleotide polymorphism (SNP) discovery and genotyping of an inter-specific mapping population, Eucalyptus tereticornis × E. camaldulensis are reported in this study. Reduced representation libraries of 81 F1 progeny were sequenced and number of SNPs considered for linkage analysis was 96,377 single nucleotide substitutions, most of which (54.03%) represented transition events. Totally, 1999 and 4206 SNP markers were mapped to female ( Eucalyptus tereticornis ) and male ( Eucalyptus camaldulensis ) map, respectively. Consensus map consisting of 4844 markers were placed on a linkage map, which spanned 1367.6 cM and had an average of one marker every 0.29 cM corresponding to a physical distance of about 0.66 Mb. The map was utilized to identify six quantitative trait loci (QTLs) associated with adventitious rooting traits of stem cuttings accounting for phenotypic variations ranging from 7.18% to 9.17%. Analysis of the genomic sequence corresponding to six QTLs led to the identification of 201 putative candidate genes and 4 key genes were related to adventitious rooting specific expression. The integrated strategy utilizing the identified QTLs along with the available genome could serve as a platform for candidate gene identification for molecular breeding of eucalypts.

How challenging RADseq data turned out to favor coalescent-based species tree inference. A case study in Aichryson (Crassulaceae)

Analysing multiple genomic regions while incorporating detection and qualification of discordance among regions has become standard for understanding phylogenetic relationships. In plants, which usually have comparatively large genomes, this is feasible by the combination of reduced-representation library (RRL) methods and high-throughput sequencing enabling the cost effective acquisition of genomic data for thousands of loci from hundreds of samples. One popular RRL method is RADseq. A major disadvantage of established RADseq approaches is the rather short fragment and sequencing range, leading to loci of little individual phylogenetic information. This issue hampers the application of coalescent-based species tree inference. The modified RADseq protocol presented here targets ca. 5,000 loci of 300-600nt length, sequenced with the latest short-read-sequencing (SRS) technology, has the potential to overcome this drawback. To illustrate the advantages of this approach we use the study group Aichryson Webb & Berthelott (Crassulaceae), a plant genus that diversified on the Canary Islands. The data analysis approach used here aims at a careful quality control of the long loci dataset. It involves an informed selection of thresholds for accurate clustering, a thorough exploration of locus properties, such as locus length, coverage and variability, to identify potential biased data and a comparative phylogenetic inference of filtered datasets, accompanied by an evaluation of resulting BS support, gene and site concordance factor values, to improve overall resolution of the resulting phylogenetic trees. The final dataset contains variable loci with an average length of 373nt and facilitates species tree estimation using a coalescent-based summary approach. Additional improvements brought by the approach are critically discussed.

Hierarchical genetic structure and implications for conservation of the world\u2019s largest salmonid, Hucho taimen

Population genetic analyses can evaluate how evolutionary processes shape diversity and inform conservation and management of imperiled species. Taimen (Hucho taimen), the world’s largest freshwater salmonid, is threatened, endangered, or extirpated across much of its range due to anthropogenic activity including overfishing and habitat degradation. We generated genetic data using high throughput sequencing of reduced representation libraries for taimen from multiple drainages in Mongolia and Russia. Nucleotide diversity estimates were within the range documented in other salmonids, suggesting moderate diversity despite widespread population declines. Similar to other recent studies, our analyses revealed pronounced differentiation among the Arctic (Selenge) and Pacific (Amur and Tugur) drainages, suggesting historical isolation among these systems. However, we found evidence for finer-scale structure within the Pacific drainages, including unexpected differentiation between tributaries and the mainstem of the Tugur River. Differentiation across the Amur and Tugur basins together with coalescent-based demographic modeling suggests the ancestors of Tugur tributary taimen likely diverged in the eastern Amur basin, prior to eventual colonization of the Tugur basin. Our results suggest the potential for differentiation of taimen at different geographic scales, and suggest more thorough geographic and genomic sampling may be needed to inform conservation and management of this iconic salmonid.

Population Genomics of the Commercially Important Gulf of Mexico Pink Shrimp Farfantepenaeus duorarum (Burkenroad, 1939) Support Models of Juvenile Transport Around the Florida Peninsula

The Gulf of Mexico pink shrimp, Farfantepenaeus duorarum, supports large fisheries in the United States and Mexico, with nearly 7,000 tons harvested from the region in 2016. Given the commercial importance of this species, management is critical: in 1997, the southern Gulf of Mexico pink shrimp fishery was declared collapsed and mitigation strategies went into effect, with recovery efforts lasting over a decade. Fisheries management can be informed and improved through a better understanding of how factors associated with early life history impact genetic diversity and population structure in the recruited population. Farfantepenaeus duorarum are short-lived, but highly fecund, and display high variability in recruitment patterns. To date, modeling the impacts of ecological, physical, and behavioral factors on juvenile settlement has focused on recruitment of larval individuals of F. duorarum to nursery grounds in Florida Bay. Here, we articulate testable hypotheses stemming from a recent model of larval transport and evaluate support for each with a population genomics approach, generating reduced representation library sequencing data for F. duorarum collected from seven regions around the Florida Peninsula. Our research represents the first and most molecular data-rich study of population structure in F. duorarum in the Gulf and reveals evidence of a differentiated population in the Dry Tortugas. Our approach largely validates a model of larval transport, allowing us to make management-informative inferences about the impacts of spawning location and recruitment patterns on intraspecific genetic diversity. Such inferences improve our understanding of the roles of non-genetic factors in generating and maintaining genetic diversity in a commercially important penaeid shrimp species.

Contrasting genetic signal of recolonization after rainforest fragmentation in African trees with different dispersal abilities

Although today the forest cover is continuous in Central Africa, this may have not always been the case, as the scarce fossil record in this region suggests that arid conditions might have significantly reduced tree density during the ice ages. Our aim was to investigate whether the dry ice age periods left a genetic signature on tree species that can be used to infer the date of the past fragmentation of the rainforest. We sequenced reduced representation libraries of 182 samples representing five widespread legume trees and seven outgroups. Phylogenetic analyses identified an early divergent lineage for all species in West Africa (Upper Guinea) and two clades in Central Africa: Lower Guinea-North and Lower Guinea-South. As the structure separating the Northern and Southern clades-congruent across species-cannot be explained by geographic barriers, we tested other hypotheses with demographic model testing using δαδι. The best estimates indicate that the two clades split between the Upper Pliocene and the Pleistocene, a date compatible with forest fragmentation driven by ice age climatic oscillations. Furthermore, we found remarkably older split dates for the shade-tolerant tree species with nonassisted seed dispersal than for light-demanding species with long-distance wind-dispersed seeds. Different recolonization abilities after recurrent cycles of forest fragmentation seem to explain why species with long-distance dispersal show more recent genetic admixture between the two clades than species with limited seed dispersal. Despite their old history, our results depict the African rainforests as a dynamic biome where tree species have expanded relatively recently after the last glaciation.

SambaR: An R package for fast, easy and reproducible population-genetic analyses of biallelic SNP data sets.

SNP data sets can be used to infer a wealth of information about natural populations, including information about their structure, genetic diversity, and the presence of loci under selection. However, SNP data analysis can be a time-consuming and challenging process, not in the least because at present many different software packages are needed to execute and depict the wide variety of mainstream population-genetic analyses. Here, we present SambaR, an integrative and user-friendly R package which automates and simplifies quality control and population-genetic analyses of biallelic SNP data sets. SambaR allows users to perform mainstream population-genetic analyses and to generate a wide variety of ready to publish graphs with a minimum number of commands (less than 10). These wrapper commands call functions of existing packages (including adegenet, ape, LEA, poppr, pcadapt and StAMPP) as well as new tools uniquely implemented in SambaR. We tested SambaR on online available SNP data sets and found that SambaR can process data sets of over 100,000 SNPs and hundreds of individuals within hours, given sufficient computing power. Newly developed tools implemented in SambaR facilitate optimization of filter settings, objective interpretation of ordination analyses, enhance comparability of diversity estimates from reduced representation library SNP data sets, and generate reduced SNP panels and structure-like plots with Bayesian population assignment probabilities. SambaR facilitates rapid population genetic analyses on biallelic SNP data sets by removing three major time sinks: file handling, software learning, and data plotting. In addition, SambaR provides a convenient platform for SNP data storage and management, as well as several new utilities, including guidance in setting appropriate data filters. The SambaR source script, manual and example data set are distributed through GitHub: https://github.com/mennodejong1986/SambaR.

Development of a 76k Alpaca (Vicugna pacos) Single Nucleotide Polymorphisms (SNPs) Microarray.

Small farm producers’ sustenance depends on their alpaca herds and the production of fiber. Genetic improvement of fiber characteristics would increase their economic benefits and quality of life. The incorporation of molecular marker technology could overcome current limitations for the implementation of genetic improvement programs. Hence, the aim of this project was the generation of an alpaca single nucleotide polymorphism (SNP) microarray. A sample of 150 Huacaya alpacas from four farms, two each in Puno and Cerro de Pasco were used for SNP discovery by genotyping by sequencing (GBS). Reduced representation libraries, two per animal, were produced after DNA digestion with ApeK1 and double digestion with Pst1-Msp1. Ten alpaca genomes, sequenced at depths between 12× to 30×, and the VicPac3.1 reference genome were used for read alignments. Bioinformatics analysis discovered 76,508 SNPs included in the microarray. Candidate genes SNPs (302) for fiber quality and color are also included. The microarray SNPs cover 90.5% of the genome length with a density of about 39 ± 2.51 SNPs/Mb of DNA at an average interval of 26.45 ± 18.57 kbp. The performance was evaluated by genotyping 30 family trios and comparing them to their pedigrees, as well as comparing microarray to GBS genotypes. Concordance values of 0.93 and 0.94 for ApeK1 and Pst1-Msp1 generated SNPs were observed. Similarly, 290 fiber quality and color candidate gene SNPs were validated. Availability of this microarray will facilitate genome-wide association studies, marker-assisted selection and, in time, genomic selection.

Persistent panmixia despite extreme habitat loss and population decline in the threatened tricolored blackbird (Agelaius tricolor).

Habitat loss and alteration has driven many species into decline, often to the point of requiring protection and intervention to avert extinction. Genomic data provide the opportunity to inform conservation and recovery efforts with details about vital evolutionary processes with a resolution far beyond that of traditional genetic approaches. The tricolored blackbird (Agelaius tricolor) has suffered severe losses during the previous century largely due to anthropogenic impacts on their habitat. Using a dataset composed of a whole genome paired with reduced representation libraries (RAD‐Seq) from samples collected across the species’ range, we find evidence for panmixia using multiple methods, including PCA (no geographic clustering), admixture analyses (ADMIXTURE and TESS conclude K = 1), and comparisons of genetic differentiation (average FST = 0.029). Demographic modeling approaches recovered an ancient decline that had a strong impact on genetic diversity but did not detect any effect from the known recent decline. We also did not detect any evidence for selection, and hence adaptive variation, at any site, either geographic or genomic. These results indicate that species continues to have high vagility across its range despite population decline and habitat loss and should be managed as a single unit.