Restriction Site-associated DNA Research Articles

A genome-wide association study of coat color in Chinese Rex rabbits.

Coat color is an important phenotypic characteristic of the domestic rabbit (Oryctolagus cuniculus) and has specific economic importance in the Rex rabbit industry. Coat color varies considerably among different populations of rabbits, and several causal genes for this variation have been thoroughly studied. Nevertheless, the candidate genes affecting coat color variation in Chinese Rex rabbits remained to be investigated. In this study, we collected blood samples from 250 Chinese Rex rabbits with six different coat colors. We performed genome sequencing using a restriction site-associated DNA sequencing approach. A total of 91,546 single nucleotide polymorphisms (SNPs), evenly distributed among 21 autosomes, were identified. Genome-wide association studies (GWAS) were performed using a mixed linear model, in which the individual polygenic effect was fitted as a random effect. We detected a total of 24 significant SNPs that were located within a genomic region on chromosome 4 (OCU4). After re-fitting the most significant SNP (OCU4:13,434,448, p = 1.31e-12) as a covariate, another near-significant SNP (OCU4:11,344,946, p = 7.03e-07) was still present. Hence, we conclude that the 2.1-Mb genomic region located between these two significant SNPs is significantly associated with coat color in Chinese Rex rabbits. The well-studied coat-color-associated agouti signaling protein (ASIP) gene is located within this region. Furthermore, low genetic differentiation was also observed among the six coat color varieties. In conclusion, our results confirmed that ASIP is a putative causal gene affecting coat color variation in Chinese Rex rabbits.

Extensive gene flow suggests lack of reproductive barriers between the two subspecies of the green whip snake, Hierophis viridiflavus (Squamata: Colubridae)

Abstract Speciation is the fundamental process that underlies the diversification of species. Different steps of speciation, such as reproductive isolation or ecological niche separation, can evolve in different orders and over different time scales. The two subspecies of the green whip snake, Hierophis viridiflavus viridiflavus and Hierophis viridiflavus carbonarius, sometimes considered as separate species, are ideal to study speciation because they diverged in isolated glacial refugia. So far, little is known about whether these two subspecies are reproductively isolated from each other. Furthermore, it is unknown whether the common melanism of the subspecies H. v. carbonarius might be genetically encoded. We quantified gene flow along two transects across the contact zone of the two subspecies using double-digest restriction site-associated DNA. Our analyses implicated the existence of two well-defined clades that experience substantial gene flow along a contact zone that is ~300 km wide. We detected a clear association between the genetic assignment and the quantitative extent of melanism within the contact zone, which is likely to be genetically polygenic. Given the level of gene flow, it seems that strong postzygotic reproductive barriers have not evolved. We conclude that the two subspecies represent evolutionarily significant units for conservation and suggest treating them taxonomically as subspecies.

Warming appears as the main risk of non-adaptedness for western Mediterranean relict fir forests under expected climate change scenarios.

Circum-Mediterranean firs are considered among the most drought-sensitive species to climate change. Understanding the genetic basis of trees' adaptive capacity and intra-specific variability to drought avoidance is mandatory to define conservation measures, thus potentially preventing their extinction. We focus here on Abies pinsapo and Abies marocana, both relict tree species, endemic from south Spain and north Morocco, respectively. A total of 607 samples were collected from eight nuclei: six from Spanish fir and two from Moroccan fir. A genotyping by sequencing technique called double digestion restriction site-associated DNA sequencing (ddRAD-seq) was performed to obtain a genetic matrix based on single-nucleotide polymorphisms (SNPs). This matrix was utilized to study the genetic structure of A. pinsapo populations and to carry out selection signature studies. In order to understand how Spanish fir and Moroccan fir cope with climate change, genotype-environment associations (GEAs) were identified. Further, the vulnerability of these species to climate variations was estimated by the risk of non-adaptedness (RONA). The filtering of the de novo assembly of A. pinsapo provided 3,982 SNPs from 504 out of 509 trees sequenced. Principal component analysis (PCA) genetically separated Grazalema from the rest of the Spanish populations. However, FST values showed significant differences among the sampling points. We found 51 loci potentially under selection. Homolog sequences were found for some proteins related to abiotic stress response, such as dehydration-responsive element binding transcription factor, regulation of abscisic acid signaling, and methylation pathway. A total of 15 associations with 11 different loci were observed in the GEA studies, with the maximum temperature of the warmest month being the variable with the highest number of associated loci. This temperature sensitivity was also supported by the risk of non-adaptedness, which yielded a higher risk for both A. pinsapo and A. marocana under the high emission scenario (Representative Concentration Pathway (RCP) 8.5). This study sheds light on the response to climate change of these two endemic species.

Ice holes microrefugia harbor genetically and functionally distinct populations of Vaccinium vitis-idaea (Ericaceae)

In the mountain terrain, ice holes are little depressions between rock boulders that are characterized by the exit of cold air able to cool down the rock surface even in summer. This cold air creates cold microrefugia in warmer surroundings that preserve plant species probably over thousands of years under extra-zonal climatic conditions. We hypothesized that ice hole populations of the model species Vaccinium vitis-idaea (Ericaceae) show genetic differentiation from nearby zonal subalpine populations, and high functional trait distinctiveness, in agreement with genetic patterns. We genotyped almost 30,000 single nucleotide polymorphisms using restriction site-associated DNA sequencing and measured eight functional traits indicative of individual performance and ecological strategies. Genetic results showed high differentiation among the six populations suggesting isolation. On siliceous bedrock, ice hole individuals exhibited higher levels of admixture than those from subalpine populations which could have experienced more bottlenecks during demographic fluctuations related to glacial cycles. Ice hole and subalpine calcareous populations clearly separated from siliceous populations, indicating a possible effect of bedrock in shaping genetic patterns. Trait analysis reflected the bedrock effect on populations’ differentiation. The significant correlation between trait and genetic distances suggests the genetic contribution in shaping intraspecific functional differentiation. In conclusion, extra-zonal populations reveal a prominent genetic and phenotypic differentiation determined by history and ecological contingency. Therefore, microrefugia populations can contribute to the overall variability of the species and lead to intraspecific-driven responses to upcoming environmental changes.

Mapping of a major QTL for increased robustness and detection of genome assembly errors in Asian seabass (Lates calcarifer)

BackgroundFor Asian seabass (Lates calcarifer, Bloch 1790) cultured at sea cages various aquatic pathogens, complex environmental and stress factors are considered as leading causes of disease, causing tens of millions of dollars of annual economic losses. Over the years, we conducted farm-based challenges by exposing Asian seabass juveniles to complex natural environmental conditions. In one of these challenges, we collected a total of 1,250 fish classified as either ‘sensitive’ or ‘robust’ individuals during the 28-day observation period.ResultsWe constructed a high-resolution linkage map with 3,089 SNPs for Asian seabass using the double digest Restriction-site Associated DNA (ddRAD) technology and a performed a search for Quantitative Trait Loci (QTL) associated with robustness. The search detected a major genome-wide significant QTL for increased robustness in pathogen-infected marine environment on linkage group 11 (ASB_LG11; 88.9 cM to 93.6 cM) with phenotypic variation explained of 81.0%. The QTL was positioned within a > 800 kb genomic region located at the tip of chromosome ASB_LG11 with two Single Nucleotide Polymorphism markers, R1-38468 and R1-61252, located near to the two ends of the QTL. When the R1-61252 marker was validated experimentally in a different mass cross population, it showed a statistically significant association with increased robustness. The majority of thirty-six potential candidate genes located within the QTL have known functions related to innate immunity, stress response or disease. By utilizing this ddRAD-based map, we detected five mis-assemblies corresponding to four chromosomes, namely ASB_LG8, ASB_LG9, ASB_LG15 and ASB_LG20, in the current Asian seabass reference genome assembly.ConclusionAccording to our knowledge, the QTL associated with increased robustness is the first such finding from a tropical fish species. Depending on further validation in other stocks and populations, it might be potentially useful for selecting robust Asian seabass lines in selection programs.

Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study.

Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.

Source tracing of Anthonomus grandis captured in areas of the USA where the species had previously been eradicated

Anthonomus grandis Boheman (Coleoptera: Curculionidae) is comprised of two subspecies: the boll weevil (A. g. grandis) and the thurberia weevil (A. g. thurberiae). The thurberia weevil is generally regarded as a non-pest variant of the species, whereas the boll weevil is an important economic pest of commercial cotton. In the USA, a nationwide eradication program has been largely successful against the boll weevil, but reintroductions of eradicated zones remain a threat. In 2015, A. grandis was trapped in the functionally eradicated Winter Garden growing region in Texas, where it had not been found since 2012. Weevils established in the area and triggered a multi-year, multi-million-dollar mitigation effort. In 2017, the species was trapped in Hidalgo County, New Mexico, another area where it had been previously eradicated. In this study, we obtained specimens from both trapping events and used double digest restriction site-associated DNA sequencing to obtain genetic data in homology with previously published A. grandis population genetic data. We conducted source tracing to match the captured weevils to previously described populations to determine a likely source population. We also compared the accuracy of the approach to other diagnostic methods currently used by the United States Department of Agriculture. The results suggested that a small founder population of weevils dispersed from the Lower Rio Grande Valley population in southern Texas and northern Tamaulipas and established in the Winter Garden. In Hidalgo Co., NM, the weevils were found to be genetically indistinct from non-pest populations of thurberia weevil associated with Arizona wild cotton (Gossypium thurberi). Thurberia weevils in Arizona remain unlikely to threaten commercial growers. Improved rapid diagnostic tools may be desirable to enable better assessment of trap captures and allow for a reduction in the application of mitigation efforts when the non-pest variant is identified.

Genome-Wide RADseq Reveals Genetic Differentiation of Wild and Cultured Populations of Large Yellow Croaker.

Larimichthys crocea (also known as the large yellow croaker) is one of the most economically important marine fishes in China, and research on the ecology and genetics of this species is of immense significance. In this study, we performed restriction site-associated DNA sequencing (RAD-seq) of 54 individuals collected from four sites in China to analyze the genetic structure and diversity of large yellow croaker at the genome level. It revealed that the large yellow croaker populations in the Ningde and Zhoushan coastal waters can be clearly distinguished. Different genetic diversity indices were used to analyze the genetic diversity of the large yellow croaker, which showed that there was a differentiation trend between the wild and farmed populations in Ningde. Moreover, we identified genetically differentiated genomic regions between the populations. GO gene enrichment analysis identified genes that are related to fatty acid metabolism and growth. These findings enhance our understanding of genetic differentiation and adaptation to different living environments, providing a theoretical basis for the preservation and restoration of the genetic resources of the large yellow croaker.

Parsing variance by marker type: Testing biogeographic hypotheses and differential contribution of historical processes to population structure in a desert lizard.

A fundamental goal of population genetic studies is to identify historical biogeographic patterns and understand the processes that generate them. However, localized demographic events can skew population genetic inference. Assessing populations with multiple types of genetic markers, each with unique mutation rates and responses to changes in population size, can help to identify potentially confounding population-specific demographic processes. Here, we compared population structure and connectivity inferred from microsatellites and restriction site-associated DNA loci among 17 populations of an arid-specialist lizard, the desert night lizard, Xantusia vigilis, in central California to test among historical processes structuring population genetic diversity. We found that both marker types yielded generally concordant insights into population genetic structure including a major phylogenetic break maintained between two populations separated by less than 10 km, suggesting that either marker type could be used to understand generalized demographic patterns across the region for management purposes. However, we also found that the effects of demography on marker discordance could be used to elucidate population histories and distinguish among competing biogeographic hypotheses. Our results suggest that comparisons of within-population diversity across marker types provide powerful opportunities for leveraging marker discordance, particularly for understanding the creation and maintenance of contact zones among clades.

Chromosome-level genome assembly and population genomic analysis reveal evolution and local adaptation in common hairfin anchovy (Setipinna tenuifilis).

Understanding the genetic structure and the factors associated with adaptive diversity has significant implications for the effective management of wild populations under threat from overfishing and climate change. The common hairfin anchovy (Setipinna tenuifilis) is an economically and ecologically important pelagic fish species, spanning a broad latitudinal gradient along marginal seas of the Northwest Pacific. In this study, we constructed the first reference genome of S. tenuifilis using PacBio long reads and high-resolution chromosome conformation capture (Hi-C) technology. The assembled genome was 798.38 Mb with a contig N50 of 1.43 Mb and a scaffold N50 of 32.42 Mb, which were anchored onto 24 pseudochromosomes. A total of 22,019 genes were functionally annotated, which accounted for 95.27% of the predicted protein-coding genes. Chromosomal collinearity analysis revealed chromosome fusion or fission events in Clupeiformes species. Three genetic groups of S. tenuifilis were revealed along the Chinese coast using restriction site-associated DNA sequencing (RADseq). We investigated the influence of four bioclimatic variables as potential drivers of adaptive divergence in S. tenuifilis, suggesting that these environmental variables, especially sea surface temperature, may play important roles as drivers of spatially varying selection for S. tenuifilis. We also identified candidate functional genes underlying adaptive mechanisms and ecological tradeoffs using redundancy analysis (RDA) and BayeScan analysis. In summary, this study sheds light on the evolution and spatial patterns of genetic variation of S. tenuifilis, providing a valuable genomic resource for further biological and genetic studies on this species and other closely related Clupeiformes.

Clonality, local population structure and gametophyte sex ratios in cryptic species of the Sphagnum magellanicum complex.

Sphagnum (peatmoss) comprises a moss (Bryophyta) clade with ~300-500 species. The genus has unparalleled ecological importance because Sphagnum-dominated peatlands store almost a third of the terrestrial carbon pool and peatmosses engineer the formation and microtopography of peatlands. Genomic resources for Sphagnum are being actively expanded, but many aspects of their biology are still poorly known. Among these are the degree to which Sphagnum species reproduce asexually, and the relative frequencies of male and female gametophytes in these haploid-dominant plants. We assess clonality and gametophyte sex ratios and test hypotheses about the local-scale distribution of clones and sexes in four North American species of the S. magellanicum complex. These four species are difficult to distinguish morphologically and are very closely related. We also assess microbial communities associated with Sphagnum host plant clones and sexes at two sites. Four hundred and five samples of the four species, representing 57 populations, were subjected to restriction site-associated DNA sequencing (RADseq). Analyses of population structure and clonality based on the molecular data utilized both phylogenetic and phenetic approaches. Multi-locus genotypes (genets) were identified using the RADseq data. Sexes of sampled ramets were determined using a molecular approach that utilized coverage of loci on the sex chromosomes after the method was validated using a sample of plants that expressed sex phenotypically. Sex ratios were estimated for each species, and populations within species. Difference in fitness between genets was estimated as the numbers of ramets each genet comprised. Degrees of clonality [numbers of genets/numbers of ramets (samples)] within species, among sites, and between gametophyte sexes were estimated. Sex ratios were estimated for each species, and populations within species. Sphagnum-associated microbial communities were assessed at two sites in relation to Sphagnum clonality and sex. All four species appear to engage in a mixture of sexual and asexual (clonal) reproduction. A single ramet represents most genets but two to eight ramets were dsumbers ansd text etected for some genets. Only one genet is represented by ramets in multiple populations; all other genets are restricted to a single population. Within populations ramets of individual genets are spatially clustered, suggesting limited dispersal even within peatlands. Sex ratios are male-biased in S. diabolicum but female-biased in the other three species, although significantly so only in S. divinum. Neither species nor males/females differ in levels of clonal propagation. At St Regis Lake (NY) and Franklin Bog (VT), microbial community composition is strongly differentiated between the sites, but differences between species, genets and sexes were not detected. Within S. divinum, however, female gametophytes harboured two to three times the number of microbial taxa as males. These four Sphagnum species all exhibit similar reproductive patterns that result from a mixture of sexual and asexual reproduction. The spatial patterns of clonally replicated ramets of genets suggest that these species fall between the so-called phalanx patterns, where genets abut one another but do not extensively mix because of limited ramet fragmentation, and the guerrilla patterns, where extensive genet fragmentation and dispersal result in greater mixing of different genets. Although sex ratios in bryophytes are most often female-biased, both male and female biases occur in this complex of closely related species. The association of far greater microbial diversity for female gametophytes in S. divinum, which has a female-biased sex ratio, suggests additional research to determine if levels of microbial diversity are consistently correlated with differing patterns of sex ratio biases.

Adaptive divergence and genetic vulnerability of relict species under climate change: a case study of Pterocarya macroptera.

Understanding adaptive genetic variation and whether it can keep pace with predicted future climate change is critical in assessing the genetic vulnerability of species and developing conservation management strategies. The lack of information on adaptive genetic variation in relict species carrying abundant genetic resources hinders the assessment of genetic vulnerability. Using a landscape genomics approach, this study aimed to determine how adaptive genetic variation shapes population divergence and to predict the adaptive potential of Pterocarya macroptera (a vulnerable relict species in China) under future climate scenarios. We applied restriction site-associated DNA sequencing (RAD-seq) to obtain 8244 single-nucleotide polymorphisms (SNPs) from 160 individuals across 28 populations. We examined the pattern of genetic diversity and divergence, and then identified outliers by genetic differentiation (FST) and genotype-environment association (GEA) methods. We further dissected the effect of geographical/environmental gradients on genetic variation. Finally, we predicted genetic vulnerability and adaptive risk under future climate scenarios. We identified three genetic lineages within P. macroptera: the Qinling-Daba-Tianmu Mountains (QDT), Western Sichuan (WS) and Northwest Yunnan (NWY) lineages, which showed significant signals of isolation by distance (IBD) and isolation by environment (IBE). IBD and IBE explained 3.7-5.7 and 8.6-12.8 % of the genetic structure, respectively. The identified GEA SNP-related genes were involved in chemical defence and gene regulation and may exhibit higher genetic variation to adapt to the environment. Gradient forest analysis revealed that the genetic variation was mainly shaped by temperature-related variables, indicating its adaptation to local thermal environments. A limited adaptive potential was suggested by the high levels of genetic vulnerability in marginal populations. Environmental gradient mainly shaped the population differentiation of P. macroptera. Marginal populations may be at high risk of extinction, and thus proactive management measures, such as assisted gene flow, are required to ensure the survival of these populations.

Diversification processes in Gerp's mouse lemur demonstrate the importance of rivers and altitude as biogeographic barriers in Madagascar's humid rainforests.

Madagascar exhibits exceptionally high levels of biodiversity and endemism. Models to explain the diversification and distribution of species in Madagascar stress the importance of historical variability in climate conditions which may have led to the formation of geographic barriers by changing water and habitat availability. The relative importance of these models for the diversification of the various forest-adapted taxa of Madagascar has yet to be understood. Here, we reconstructed the phylogeographic history of Gerp's mouse lemur (Microcebus gerpi) to identify relevant mechanisms and drivers of diversification in Madagascar's humid rainforests. We used restriction site associated DNA (RAD) markers and applied population genomic and coalescent-based techniques to estimate genetic diversity, population structure, gene flow and divergence times among M. gerpi populations and its two sister species M. jollyae and M. marohita. Genomic results were complemented with ecological niche models to better understand the relative barrier function of rivers and altitude. We show that M. gerpi diversified during the late Pleistocene. The inferred ecological niche, patterns of gene flow and genetic differentiation in M. gerpi suggest that the potential for rivers to act as biogeographic barriers depended on both size and elevation of headwaters. Populations on opposite sides of the largest river in the area with headwaters that extend far into the highlands show particularly high genetic differentiation, whereas rivers with lower elevation headwaters have weaker barrier functions, indicated by higher migration rates and admixture. We conclude that M. gerpi likely diversified through repeated cycles of dispersal punctuated by isolation to refugia as a result of paleoclimatic fluctuations during the Pleistocene. We argue that this diversification scenario serves as a model of diversification for other rainforest taxa that are similarly limited by geographic factors. In addition, we highlight conservation implications for this critically endangered species, which faces extreme habitat loss and fragmentation.

Retained chromosomal integrity following CRISPR-Cas9-based mutational correction in human embryos

Human germline gene correction by targeted nucleases holds great promise for reducing mutation transmission. However, recent studies have reported concerning observations in CRISPR-Cas9-targeted human embryos, including mosaicism and loss of heterozygosity (LOH). The latter has been associated with either gene conversion or (partial) chromosome loss events. In this study, we aimed to correct a heterozygous basepair substitution in PLCZ1, related to infertility. In 36% of the targeted embryos that originated from mutant sperm, only wild-type alleles were observed. By performing genome-wide double-digest restriction site-associated DNA sequencing, integrity of the targeted chromosome (i.e., no deletions larger than 3 Mb or chromosome loss) was confirmed in all seven targetedGENType-analyzed embryos (mutant editing and absence of mutation), while short-range LOH events (shorter than 10 Mb) were clearly observed by single-nucleotide polymorphism assessment in two of these embryos. These results fuel the currently ongoing discussion on double-strand break repair in early human embryos, making a case for the occurrence of gene conversion events or partial template-based homology-directed repair.

Microsatellite markers for Monitipora digitata designed using restriction-site associated DNA sequencing

Montipora digitata is a species belonging to the Acroporidae. In the Indo-Pacific region, M. digitata is widely distributed and is the dominant species of scleractinian coral in the South China Sea, however, there are currently no molecular markers suitable for assessing the species genetic diversity. Here, restriction-site associated DNA sequencing (RAD-seq) was used to isolate and characterize polymorphic microsatellite loci. A total of 317,361 RAD-tags were obtained using RAD-seq, including 6,778 microsatellite loci. Primer pairs for 106 loci were ordered and twenty-one polymorphic loci, that amplified reliably were identified. The number of alleles per locus were 2-7, observed heterozygosity was 0.111-0.556 with an average value of 0.285, and expected heterozygosity was 0.105- 0.802, with an average value of 0.536. Before Bonferroni correction 13 loci deviated significantly from the expectations of Hardy-Weinberg equilibrium (P < 0.05), after correction, two microsatellite loci deviated significantly (P < 0.0002). The polymorphic information content (PIC) ranged from 0.100-0.778, with 12 loci highly polymorphic (PIC > 0.5), six moderately polymorphic (0.25 < PIC < 0.5), and three loci with low polymorphism (PIC < 0.25). The microsatellite loci developed here will be effective tools for conservation genetic research on M. digitata.

Combination of GWAS and FST-based approaches identified loci associated with economic traits in sugarcane.

Sugarcane is a globally important plant for both sugar and biofuel production. Although conventional breeding has played an important role in increasing the productivity of sugarcane, it takes a long time to achieve breeding goals such as high yield and resistant to diseases. Molecular breeding, including marker-assisted breeding and genomic selection, can accelerate genetic improvement by selecting elites at the seedling stage with DNA markers. However, only a few DNA markers associated with important traits were identified in sugarcane. The purpose of this study was to identify DNA markers associated with sugar content, stalk diameter, and sugarcane top borer resistance. The sugarcane samples with trait records were genotyped using the restriction site-associated DNA sequencing (RADseq) technology. Using FST analysis and genome-wide association study (GWAS), a total of 9, 23 and 9 DNA variants (single nucleotide polymorphisms (SNPs)/insertions and deletions (indels)) were associated with sugar content, stalk diameter, and sugarcane top borer resistance, respectively. The identified genetic variants were on different chromosomes, suggesting that these traits are complex and determined by multiple genetic factors. These DNA markers identified by both approaches have the potential to be used in selecting elite clones at the seeding stage in our sugarcane breeding program to accelerate genetic improvement. Certainly, it is essential to verify the reliability of the identified DNA markers associated with traits before they are used in molecular breeding in other populations.

A pilot-scale comparison between single and double-digest RAD markers generated using GBS strategy in sesame (Sesamum indicum L.).

To reduce the genome sequence representation, restriction site-associated DNA sequencing (RAD-seq) protocols is being widely used either with single-digest or double-digest methods. In this study, we genotyped the sesame population (48 sample size) in a pilot scale to compare single and double-digest RAD-seq (sd and ddRAD-seq) methods. We analysed the resulting short-read data generated from both protocols and assessed their performance impacting the downstream analysis using various parameters. The distinct k-mer count and gene presence absence variation (PAV) showed a significant difference between the sesame samples studied. Additionally, the variant calling from both datasets (sdRAD-seq and ddRAD-seq) exhibits a significant difference between them. The combined variants from both datasets helped in identifying the most diverse samples and possible sub-groups in the sesame population. The most diverse samples identified from each analysis (k-mer, gene PAV, SNP count, Heterozygosity, NJ and PCA) can possibly be representative samples holding major diversity of the small sesame population used in this study. The best possible strategies with suggested inputs for modifications to utilize the RAD-seq strategy efficiently on a large dataset containing thousands of samples to be subjected to molecular analysis like diversity, population structure and core development studies were discussed.

Betula mcallisteri sp. nov. (sect. Acuminatae, Betulaceae), a new diploid species overlooked in the wild and in cultivation, and its relation to the widespread B. luminifera.

Taxa are traditionally identified using morphological proxies for groups of evolutionarily isolated populations. These proxies are common characters deemed by taxonomists as significant. However, there is no general rule on which character or sets of characters are appropriate to circumscribe taxa, leading to discussions and uncertainty. Birch species are notoriously hard to identify due to strong morphological variability and factors such as hybridization and the existence of several ploidy levels. Here, we present evidence for an evolutionarily isolated line of birches from China that are not distinguishable by traditionally assumed taxon recognition proxies, such as fruit or leaf characters. We have discovered that some wild material in China and some cultivated in the Royal Botanic Gardens Edinburgh, formerly recognized as Betula luminifera, differ from other individuals by having a peeling bark and a lack of cambial fragrance. We use restriction site-associated DNA sequencing and flow cytometry to study the evolutionary status of the unidentified Betula samples to assess the extent of hybridization between the unidentified Betula samples and typical B. luminifera in natural populations. Molecular analyses show the unidentified Betula samples as a distinct lineage and reveal very little genetic admixture between the unidentified samples and B. luminifera. This may also be facilitated by the finding that B. luminifera is tetraploid, while the unidentified samples turned out to be diploid. We therefore conclude that the samples represent a yet unrecognized species, which is here described as Betula mcallisteri.

Restriction site-associated DNA sequencing technologies as an alternative to low-density SNP chips for genomic selection: a simulation study in layer chickens

BackgroundTo reduce the cost of genomic selection, a low-density (LD) single nucleotide polymorphism (SNP) chip can be used in combination with imputation for genotyping selection candidates instead of using a high-density (HD) SNP chip. Next-generation sequencing (NGS) techniques have been increasingly used in livestock species but remain expensive for routine use for genomic selection. An alternative and cost-efficient solution is to use restriction site-associated DNA sequencing (RADseq) techniques to sequence only a fraction of the genome using restriction enzymes. From this perspective, use of RADseq techniques followed by an imputation step on HD chip as alternatives to LD chips for genomic selection was studied in a pure layer line.ResultsGenome reduction and sequencing fragments were identified on reference genome using four restriction enzymes (EcoRI, TaqI, AvaII and PstI) and a double-digest RADseq (ddRADseq) method (TaqI-PstI). The SNPs contained in these fragments were detected from the 20X sequence data of the individuals in our population. Imputation accuracy on HD chip with these genotypes was assessed as the mean correlation between true and imputed genotypes. Several production traits were evaluated using single-step GBLUP methodology. The impact of imputation errors on the ranking of the selection candidates was assessed by comparing a genomic evaluation based on ancestry using true HD or imputed HD genotyping. The relative accuracy of genomic estimated breeding values (GEBVs) was investigated by considering the GEBVs estimated on offspring as a reference. With AvaII or PstI and ddRADseq with TaqI and PstI, more than 10 K SNPs were detected in common with the HD SNP chip, resulting in an imputation accuracy greater than 0.97. The impact of imputation errors on genomic evaluation of the breeders was reduced, with a Spearman correlation greater than 0.99. Finally, the relative accuracy of GEBVs was equivalent.ConclusionsRADseq approaches can be interesting alternatives to low-density SNP chips for genomic selection. With more than 10 K SNPs in common with the SNPs of the HD SNP chip, good imputation and genomic evaluation results can be obtained. However, with real data, heterogeneity between individuals with missing data must be considered.

Genomic insights into genetic diversity and local adaptation of a dominant desert steppe feather grass, Stipa breviflora Griseb.

Investigating the genetic mechanisms of local adaptation is critical to understanding how species adapt to heterogeneous environments. In the present study, we analyzed restriction site-associated DNA sequencing (RADseq) data in order to explore genetic diversity, genetic structure, genetic differentiation, and local adaptation of Stipa breviflora. In total, 135 individual plants were sequenced and 25,786 polymorphic loci were obtained. We found low genetic diversity (He=0.1284) within populations of S. breviflora. Four genetic clusters were identified along its distribution range. The Mantel test, partial Mantel test, and multiple matrix regression with randomization (MMRR) indicate that population differentiation was caused by both geographic distance and environmental factors. Through the FST outlier test and environmental association analysis (EAA), 113 candidate loci were identified as putatively adaptive loci. RPK2 and CPRF1, which are associated with meristem maintenance and light responsiveness, respectively, were annotated. To explore the effects of climatic factors on genetic differentiation and local adaptation of S. breviflora, gradient forest (GF) analysis was applied to 25,786 single nucleotide polymorphisms (SNPs) and 113 candidate loci, respectively. The results showed that both temperature and precipitation affected the genetic differentiation of S. breviflora, and precipitation was strongly related to local adaptation. Our study provides a theoretical basis for understanding the local adaptation of S. breviflora.