Search

The genetic diversity found in natural populations is the result of the evolutionary forces in response to historical and contemporary factors. The environmental characteristics and geological history of Mexico promoted the evolution and diversification of plant species, including wild relatives of crops such as the wild pumpkins (Cucurbita). Wild pumpkin species are found in a variety of habitats, evidencing their capability to adapt to different environments. Despite the potential value of wild Cucurbita as a genetic reservoir for crops, there is a lack of studies on their genetic diversity. Cucurbita radicans is an endangered species threatened by habitat destruction leading to low densities in small and isolated populations. Here, we analyze Genotype by Sequencing genomic data of the wild pumpkin C. radicans to evaluate the influence of factors like isolation, demographic history, and the environment shaping the amount and distribution of its genetic variation. We analyzed 91 individuals from 14 localities along its reported distribution. We obtained 5,107 SNPs and found medium-high levels of genetic diversity and genetic structure distributed in four main geographic areas with different environmental conditions. Moreover, we found signals of demographic growth related to historical climatic shifts. Outlier loci analysis showed significant association with the environment, principally with precipitation variables. Also, the outlier loci displayed differential changes in their frequencies in response to future global climate change scenarios. Using the results of genetic structure, outlier loci and multivariate analyses of the environmental conditions, we propose priority localities for conservation that encompass most of the genetic diversity of C. radicans., (© 2024. The Author(s).)

Past sea level changes and geological instability along watershed boundaries have largely influenced fish distribution across coastal basins, either by dispersal via palaeodrainages now submerged or by headwater captures, respectively. Accordingly, the South American Atlantic coast encompasses several small and isolated drainages that share a similar species composition, representing a suitable model to infer historical processes. Leporinus bahiensis is a freshwater fish species widespread along adjacent coastal basins over narrow continental shelf with no evidence of palaeodrainage connections at low sea level periods. Therefore, this study aimed to reconstruct its evolutionary history to infer the role of headwater captures in the dispersal process. To accomplish this, we employed molecular-level phylogenetic and population structure analyses based on Sanger sequences (5 genes) and genome-wide SNP data. Phylogenetic trees based on Sanger data were inconclusive, but SNPs data did support the monophyletic status of L. bahiensis. Both COI and SNP data revealed structured populations according to each hydrographic basin. Species delimitation analyses revealed from 3 (COI) to 5 (multilocus approach) MOTUs, corresponding to the sampled basins. An intricate biogeographic scenario was inferred and supported by Approximate Bayesian Computation (ABC) analysis. Specifically, a staggered pattern was revealed and characterized by sequential headwater captures from basins adjacent to upland drainages into small coastal basins at different periods. These headwater captures resulted in dispersal throughout contiguous coastal basins, followed by deep genetic divergence among lineages. To decipher such recent divergences, as herein represented by L. bahiensis populations, we used genome-wide SNPs data. Indeed, the combined use of genome-wide SNPs data and ABC method allowed us to reconstruct the evolutionary history and speciation of L. bahiensis. This framework might be useful in disentangling the diversification process in other neotropical fishes subject to a reticulate geological history., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Society of Systematic Biologists.)

Rice is an important crop that is consumed by approximately half of the world's population on a regular basis. Plant height is an important characteristic with shorter rice often having higher lodging resistance and better soil nutrient utilization allowing for lower fertilizer use. We used a Chromosome Segment Substitution Line (CSSL) population generated by introgressing segments of CT9993 and IR62266 into KDML 105. We identified height QTLs on chromosomes 1 and 4. We performed whole genome sequencing of the parental lines and found that IR62266 has the deletion in Gibberellin 20-oxidase 2 corresponding to the semi-dwarf 1 locus. However, short height on chromosome 1 came from CT9993 with no mutation in Gibberellin 20-oxidase 2, or any known height genes. The height QTL on chromosome 4 contains mutations in Peroxisome biogenesis protein 6, which has been linked to a reduced growth phenotype in A. thaliana, making this a good candidate height gene., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Use of genomic prediction (GP) in tetraploid is becoming more common. Therefore, we think it is the right time for a comparison of GP models for tetraploid potato. GP models were compared that contrasted shrinkage with variable selection, parametric vs. non-parametric models and different ways of accounting for non-additive genetic effects. As a complement to GP, association studies were carried out in an attempt to understand the differences in prediction accuracy. We compared our GP models on a data set consisting of 147 cultivars, representing worldwide diversity, with over 39 k GBS markers and measurements on four tuber traits collected in six trials at three locations during 2 years. GP accuracies ranged from 0.32 for tuber count to 0.77 for dry matter content. For all traits, differences between GP models that utilised shrinkage penalties and those that performed variable selection were negligible. This was surprising for dry matter, as only a few additive markers explained over 50% of phenotypic variation. Accuracy for tuber count increased from 0.35 to 0.41, when dominance was included in the model. This result is supported by Genome Wide Association Study (GWAS) that found additive and dominance effects accounted for 37% of phenotypic variation, while significant additive effects alone accounted for 14%. For tuber weight, the Reproducing Kernel Hilbert Space (RKHS) model gave a larger improvement in prediction accuracy than explicitly modelling epistatic effects. This is an indication that capturing the between locus epistatic effects of tuber weight can be done more effectively using the semi-parametric RKHS model. Our results show good opportunities for GP in 4x potato., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Wilson, Zheng, Maliepaard, Mulder, Visser, van der Burgt and van Eeuwijk.)

Premise: Long-distance dispersal has been important in explaining the present distributions of many plant species. Despite being infrequent, such dispersal events have considerable evolutionary consequences, because bottlenecks during colonization can result in reduced genetic diversity. We examined the phylogeographic history of Lycium carolinianum, a widespread taxon that ranges from southeastern North America to several Pacific islands, with intraspecific diversity in sexual and mating systems., Methods: We used Bayesian, likelihood, and coalescent approaches with nuclear and plastid sequence data and genome-wide single nucleotide polymorphisms to reconstruct the dispersal history of this species. We also compared patterns of genetic variation in mainland and island populations using single nucleotide polymorphisms and allelic diversity at the S-RNase mating system gene., Results: Lycium carolinianum is monophyletic and dispersed once from the North American mainland, colonizing the Pacific islands ca. 40,100 years ago. This dispersal was accompanied by a loss of genetic diversity in SNPs and the S-RNase locus due to a colonization bottleneck and the loss of self-incompatibility. Additionally, we documented at least two independent transitions to gynodioecy: once following the colonization of the Hawaiian Islands and loss of self-incompatibility, and a second time associated with polyploidy in the Yucatán region of Mexico., Conclusions: Long-distance dispersal via fleshy, bird dispersed fruits best explains the unusually widespread distribution of L. carolinianum. The collapse of diversity at the S-RNase locus in island populations suggests that self-fertilization may have facilitated the subsequent colonization of Pacific islands following a single dispersal from mainland North America., (© 2021 Botanical Society of America.)

Uncovering the genetic basis of local adaptation is a major goal of evolutionary biology and conservation science alike. In an era of climate change, an understanding of how environmental factors shape adaptive diversity is crucial to predicting species response and directing management. Here, we investigate patterns of genomic variation in giant sequoia, an iconic and ecologically important tree species, using 1,364 bi-allelic single nucleotide polymorphisms (SNPs). We use an F ST outlier test and two genotype-environment association methods, latent factor mixed models (LFMMs) and redundancy analysis (RDA), to detect complex signatures of local adaptation. Results indicate 79 genomic regions of potential adaptive importance, with limited overlap between the detection methods. Of the 58 loci detected by LFMM, 51 showed strong correlations to a precipitation-driven composite variable and seven to a temperature-related variable. RDA revealed 24 outlier loci with association to climate variables, all of which showed strongest relationship to summer precipitation. Nine candidate loci were indicated by two methods. After correcting for geographic distance, RDA models using climate predictors accounted for 49% of the explained variance and showed significant correlations between SNPs and climatic factors. Here, we present evidence of local adaptation in giant sequoia along gradients of precipitation and provide a first step toward identifying genomic regions of adaptive significance. The results of this study will provide information to guide management strategies that seek to maximize adaptive potential in the face of climate change., Competing Interests: The authors declare no conflicts of interest., (© 2020 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Background: The development of genome-wide genotyping resources has provided terrestrial livestock and crop industries with the unique ability to accurately assess genomic relationships between individuals, uncover the genetic architecture of commercial traits, as well as identify superior individuals for selection based on their specific genetic profile. Utilising recent advancements in de-novo genome-wide genotyping technologies, it is now possible to provide aquaculture industries with these same important genotyping resources, even in the absence of existing genome assemblies. Here, we present the development of a genome-wide SNP assay for the Black Tiger shrimp (Penaeus monodon) through utilisation of a reduced-representation whole-genome genotyping approach (DArTseq)., Results: Based on a single reduced-representation library, 31,262 polymorphic SNPs were identified across 650 individuals obtained from Australian wild stocks and commercial aquaculture populations. After filtering to remove SNPs with low read depth, low MAF, low call rate, deviation from HWE, and non-Mendelian inheritance, 7542 high-quality SNPs were retained. From these, 4236 high-quality genome-wide loci were selected for baits-probe development and 4194 SNPs were included within a finalized target-capture genotype-by-sequence assay (DArTcap). This assay was designed for routine and cost effective commercial application in large scale breeding programs, and demonstrates higher confidence in genotype calls through increased call rate (from 80.2 ± 14.7 to 93.0% ± 3.5%), increased read depth (from 20.4 ± 15.6 to 80.0 ± 88.7), as well as a 3-fold reduction in cost over traditional genotype-by-sequencing approaches., Conclusion: Importantly, this assay equips the P. monodon industry with the ability to simultaneously assign parentage of communally reared animals, undertake genomic relationship analysis, manage mate pairings between cryptic family lines, as well as undertake advance studies of genome and trait architecture. Critically this assay can be cost effectively applied as P. monodon breeding programs transition to undertaking genomic selection.

• Herein we propose a framework for assembling and analyzing Genotype by Sequencing (GBS) data to better understand evolutionary relationships within a group of closely related species using the mastiff bats ( Molossus ) as our model system. Many species within this genus have low-levels of genetic variation within and between morphologically distinct species, and the relationships among them remain unresolved using traditional Sanger sequencing methods. Given that both de novo and reference genome pipelines can be used to assemble next generation sequences, and that several tree inference methodologies have been proposed for single nucleotide polymorphism (SNP) data, we test whether different alignments and phylogenetic approaches produce similar results. We also examined how the process of SNP identification and mapping can affect the consistency of the analyses. Different alignments and phylogenetic inferences produced consistent results, supporting the GBS approach for answering evolutionary questions on a macroevolutionary scale when the genetic distance among phenotypically identifiable clades is low. We highlight the importance of exploring the relationships among groups using different assembly assumptions and also distinct phylogenetic inference methods, particularly when addressing phylogenetic questions in genetic and morphologically conservative taxa. • The method uses the comparison of several filter settings, alignments, and tree inference approaches on Genotype by Sequencing data. • Consistent results were found among several approaches. • The methodology successfully recovered well supported species boundaries and phylogenetic relationships among species of mastiff bats not hypothesized by previous methods., (© 2020 The Author(s). Published by Elsevier B.V.)

The broad variability of Cucumis melo (melon, Cucurbitaceae) presents a challenge to conventional classification and organization within the species. To shed further light on the infraspecific relationships within C. melo , we compared genotypic and metabolomic similarities among 44 accessions representative of most of the cultivar-groups. Genotyping-by-sequencing (GBS) provided over 20,000 single-nucleotide polymorphisms (SNPs). Metabolomics data of the mature fruit flesh and rind provided over 80,000 metabolomic and elemental features via an orchestra of six complementary metabolomic platforms. These technologies probed polar, semi-polar, and non-polar metabolite fractions as well as a set of mineral elements and included both flavor- and taste-relevant volatile and non-volatile metabolites. Together these results enabled an estimate of "metabolomic/elemental distance" and its correlation with the genetic GBS distance of melon accessions. This study indicates that extensive and non-targeted metabolomics/elemental characterization produced classifications that strongly, but not completely, reflect the current and extensive genetic classification. Certain melon Groups, such as Inodorous, clustered in parallel with the genetic classifications while other genome to metabolome/element associations proved less clear. We suggest that the combined genomic, metabolic, and element data reflect the extensive sexual compatibility among melon accessions and the breeding history that has, for example, targeted metabolic quality traits, such as taste and flavor.

Mammals are one of the better known groups of animals, and in the Neotropics bats typically comprise about half of the mammalian species diversity. But, well resolved species-level phylogenies are still lacking for most taxa of bats. One broadly distributed genus is the mastiff bats, Molossus. Species within this genus are morphologically very similar, which results in a confusing and unstable taxonomy. In addition, low levels of genetic divergence among some clades make resolution of phylogenetic relationships difficult. Most authors recognize Molossus as being monophyletic, however, phylogenetic relationships within the genus remain poorly understood based on traditional Sanger sequencing of individual genes. We propose a more comprehensive framework based on large-scale genomic data derived from Next Generation Sequencing techniques to better understand evolutionary relationships within a group of closely related species with a rich taxonomic history. In this study, we utilized the NGS method of Genotype by Sequencing (GBS) to test the monophyly of the genus, understand evolutionary relationships within Molossus and investigate the genetic integrity of currently recognized species. Given that both de novo and reference genome pipelines are often used in the assembly of Single Nucleotide Polymorphism data from GBS, and that several tree inference methodologies have been proposed for SNP data, we test whether different alignments and phylogenetic approaches produce similar results. We also examined how the process of SNP identification and mapping can affect the consistency of the analyses. Our data provide the first high resolution phylogeny for the genus Molossus, bringing new insights into recognition of species boundaries and relationships among taxa. This study clarifies the taxonomy of Molossus and elevates the number of species in the genus from 11 to 14. We suggest the revalidation of the names M. nigricans, and M. fluminensis, which were synonymized under the name M. rufus; and M. bondae, previously synonymized under the name M. currentium. Different alignments and phylogenetic inferences produce consistent results, supporting use of SNP approach in addressing evolutionary questions on a macroevolutionary scale where the genetic distance among clades is low., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The complete or partial loss of shattering ability occurred independently during the domestication of several crops. Therefore, the study of this trait can provide an understanding of the link between phenotypic and molecular convergent evolution. The genetic dissection of 'pod shattering' in Phaseolus vulgaris is achieved here using a population of introgression lines and next-generation sequencing techniques. The 'occurrence' of the indehiscent phenotype (indehiscent versus dehiscent) depends on a major locus on chromosome 5. Furthermore, at least two additional genes are associated with the 'level' of shattering (number of shattering pods per plant: low versus high) and the 'mode' of shattering (non-twisting versus twisting pods), with all of these loci contributing to the phenotype by epistatic interactions. Comparative mapping indicates that the major gene identified on common bean chromosome 5 corresponds to one of the four quantitative trait loci for pod shattering in Vigna unguiculata. None of the loci identified comprised genes that are homologs of the known shattering genes in Glycine max. Therefore, although convergent domestication can be determined by mutations at orthologous loci, this was only partially true for P. vulgaris and V. unguiculata, which are two phylogenetically closely related crop species, and this was not the case for the more distant P. vulgaris and G. max. Conversely, comparative mapping suggests that the convergent evolution of the indehiscent phenotype arose through mutations in different genes from the same underlying gene networks that are involved in secondary cell-wall biosynthesis and lignin deposition patterning at the pod level., (© 2018 The Authors The Plant Journal © 2018 John Wiley & Sons Ltd.)

Interactions between environmental factors and complex life-history characteristics of marine organisms produce the genetic diversity and structure observed within species. Our main goal was to test for genetic differentiation among eastern oyster populations from the coastal region of Canadian Maritimes against expected genetic homogeneity caused by historical events, taking into account spatial and environmental (temperature, salinity, turbidity) variation. This was achieved by genotyping 486 individuals originating from 13 locations using RADSeq. A total of 11,321 filtered SNPs were used in a combination of population genomics and environmental association analyses. We revealed significant neutral genetic differentiation (mean F ST  = 0.009) between sampling locations, and the occurrence of six major genetic clusters within the studied system. Redundancy analyses (RDAs) revealed that spatial and environmental variables explained 3.1% and 4.9% of the neutral genetic variation and 38.6% and 12.2% of the putatively adaptive genetic variation, respectively. These results indicate that these environmental factors play a role in the distribution of both neutral and putatively adaptive genetic diversity in the system. Moreover, polygenic selection was suggested by genotype-environment association analysis and significant correlations between additive polygenic scores and temperature and salinity. We discuss our results in the context of their conservation and management implications for the eastern oyster.

Fiber quality and yield improvement are crucial for cotton domestication and breeding. With the transformation in spinning techniques and multiplicity needs, the development of cotton fiber quality and yield is of great importance. A genetic map of 5178 Single Nucleotide Polymorphism (SNP) markers were generated using 277 F 2:3 population, from an intra-specific cross between two upland cotton accessions, CCRI35 a high fiber quality as female and Nan Dan Ba Di Da Hua (NH), with good yield properties as male parent. The map spanned 4768.098 cM with an average distance of 0.92 cM. A total of 110 Quantitative Traits Loci (QTLs) were identified for 11 traits, but only 30 QTLs were consistent in at least two environments. The highest percentage of phenotypic variance explained by a single QTL was 15.45%. Two major cluster regions were found, cluster 1 (chromosome17-D03) and cluster 2 (chromosome26-D12). Five candidate genes were identified in the two QTL cluster regions. Based on GO functional annotation, all the genes were highly correlated with fiber development, with functions such as protein kinase and phosphorylation. The five genes were associated with various fiber traits as follows: Gh_D03G0889 linked to qFM-D03_cb, Gh_D12G0093 , Gh_D12G0410 , Gh_D12G0435 associated with qFS-D12_cb and Gh_D12G0969 linked to qFY-D12_cb. Further structural annotation and fine mapping is needed to determine the specific role played by the five identified genes in fiber quality and yield related pathway., Competing Interests: The authors declare no conflict of interest.

The necrotrophic fungal pathogen Pyrenophora teres f. teres causes the foliar disease net form net blotch (NFNB) on barley. To investigate the genetics of virulence in the barley- P. teres f. teres pathosystem, we evaluated 118 progeny derived from a cross between the California isolates 15A and 6A on the barley lines Rika and Kombar, chosen based on their differential reactions to isolates 15A and 6A for NFNB disease. Genetic maps generated with SNP, SSR, and AFLP markers were scanned for quantitative trait loci (QTL) associated with virulence in P. teres f. teres. Loci underlying two major QTL, VR1 and VR2, were associated with virulence on Rika barley, accounting for 35% and 20% of the disease reaction type variation, respectively. Two different loci, VK1 and VK2, were shown to underlie two major QTL associated with virulence on Kombar barley accounting for 26% and 19% of the disease reaction type variation, respectively. Progeny isolates harboring VK1, VK2, or VR2 alone were inoculated onto a Rika×Kombar recombinant inbred line mapping population and the susceptibility induced by each pathogen genotype corresponded to the same region on barley chromosome 6H as that identified for the parental isolates 15A and 6A. The data presented here indicate that the P. teres f. teres - barley interaction can at least partially be explained by pathogen-produced necrotrophic effectors (NEs) that interact with dominant barley susceptibility genes resulting in NE triggered susceptibility (NETS)., (Published by Elsevier Inc.)