Single nucleotide polymorphism SNP19140160 A > C is a potential breeding locus for fast-growth largemouth bass (Micropterus salmoides).

BACKGROUND: Largemouth bass (Micropterus salmoides) has significant economic value as a high-yielding fish species in China's freshwater aquaculture industry. Determining the major genes related to growth traits and identifying molecular markers associated with these traits serve as the foundation for breeding strategies involving gene pyramiding. In this study, we screened restriction-site associated DNA sequencing (RAD-seq) data to identify single nucleotide polymorphism (SNP) loci potentially associated with extreme growth differences between fast-growth and slow-growth groups in the F1 generation of a largemouth bass population. RESULTS: We subsequently identified associations between these loci and specific candidate genes related to four key growth traits (body weight, body length, body height, and body thickness) based on SNP genotyping. In total, 4,196,486 high-quality SNPs were distributed across 23 chromosomes. Using a population-specific genotype frequency threshold of 0.7, we identified 30 potential SNPs associated with growth traits. Among the 30 SNPs, SNP19140160, SNP9639603, SNP9639605, and SNP23355498 showed significant associations; three of them (SNP9639603, SNP9639605, and SNP23355498) were significantly associated with one trait, body length, in the F1 generation, and one (SNP19140160) was significantly linked with four traits (body weight, height, length, and thickness) in the F1 generation. The markers SNP19140160 and SNP23355498 were located near two growth candidate genes, fam174b and ppip5k1b, respectively, and these candidate genes were closely linked with growth, development, and feeding. The average body weight of the group with four dominant genotypes at these SNP loci in the F1 generation population (703.86 g) was 19.63% higher than that of the group without dominant genotypes at these loci (588.36 g). CONCLUSIONS: Thus, these four markers could be used to construct a population with dominant genotypes at loci related to fast growth. These findings demonstrate how markers can be used to identify genes related to fast growth, and will be useful for molecular marker-assisted selection in the breeding of high-quality largemouth bass.

Which contributes more to the relict flora distribution pattern in East Asia, geographical processes or climate change? New evidence from the phylogeography of Rehderodendron kwangtungense.

BACKGROUND: Relict species are important for enhancing the understanding of modern biogeographic distribution patterns. Although both geological and climatic changes since the Cenozoic have affected the relict flora in East Asia, the contributions of geographical processes remain unclear. In this study, we employed restriction-site associated DNA sequencing (RAD-seq) and shallow genome sequencing data, in conjunction with ecological niche modeling (ENM), to investigate the spatial genetic patterns and population differentiation history of the relict species Rehderodendron kwangtungense Chun. RESULTS: A total of 138 individuals from 16 populations were collected, largely covering the natural distribution of R. kwangtungense. The genetic diversity within the R. kwangtungense populations was extremely low (HO = 0.048 ± 0.019; HE = 0.033 ± 0.011). Mantel tests revealed isolation-by-distance pattern (R2 = 0.38, P < 0.001), and AMOVA analysis showed that the genetic variation of R. kwangtungense occurs mainly between populations (86.88%, K = 7). Between 23 and 21 Ma, R. kwangtungense underwent a period of rapid differentiation that coincided with the rise of the Himalayas and the establishment of the East Asian monsoon. According to ENM and population demographic history, the suitable area and effective population size of R. kwangtungense decreased sharply during the glacial period and expanded after the last glacial maximum (LGM). CONCLUSION: Our study shows that the distribution pattern of southern China mountain relict flora may have developed during the panplain stage between the middle Oligocene and the early Miocene. Then, the flora later fragmented under the force of orogenesis, including intermittent uplift during the Cenozoic Himalayan orogeny and the formation of abundant rainfall associated with the East Asian monsoon. The findings emphasized the predominant role of geographical processes in shaping relict plant distribution patterns.

Phylogenetic insights into the Salicaceae: The evolution of willows and beyond.

The Salicaceae includes approximately 54 genera and over 1,400 species with a cosmopolitan distribution. Members of the family are well-known for their diverse secondary plant metabolites, and they play crucial roles in tropical and temperate forest ecosystems. Phylogenetic reconstruction of the Salicaceae has been historically challenging due to the limitations of molecular markers and the extensive history of hybridization and polyploidy within the family. Our study employs whole-genome sequencing of 74 species to generate an extensive phylogeny of the Salicaceae. We generated two RAD-Seq enriched whole-genome sequence datasets and extracted two additional gene sets corresponding to the universal Angiosperms353 and Salicaceae-specific targeted-capture arrays. We reconstructed maximum likelihood-based molecular phylogenies using supermatrix and coalescent-based supertree approaches. Our fossil-calibrated phylogeny estimates that the Salicaceae originated around 128 million years ago and unravels the complex taxonomic relationships within the family. Our findings confirm the non-monophyly of the subgenus Salix s.l. and further support the merging of subgenera Chamaetia and Vetrix, both of which exhibit intricate patterns within and among different sections. Overall, our study not only enhances our understanding of the evolution of the Salicaceae, but also provides valuable insights into the complex relationships within the family.

The genomic diversity of the Eliurus genus in northern Madagascar with a putative new species.

Madagascar exhibits extraordinarily high level of species richness and endemism, while being severely threatened by habitat loss and fragmentation (HL&F). In front of these threats to biodiversity, conservation effort can be directed, for instance, in the documentation of species that are still unknown to science, or in investigating how species respond to HL&F. The tufted-tail rats genus (Eliurus spp.) is the most speciose genus of endemic rodents in Madagascar, with 13 described species, which occupy two major habitat types: dry or humid forests. The large species diversity and association to specific habitat types make the Eliurus genus a suitable model for investigating species adaptation to new environments, as well as response to HL&F (dry vs humid). In the present study, we investigated Eliurus spp. genomic diversity across northern Madagascar, a region covered by both dry and humid fragmented forests. From the mitochondrial DNA (mtDNA) and nuclear genomic (RAD-seq) data of 124 Eliurus individuals sampled in poorly studied forests of northern Madagascar, we identified an undescribed Eliurus taxon (Eliurus sp. nova). We tested the hypothesis of a new Eliurus species using several approaches: i) DNA barcoding; ii) phylogenetic inferences; iii) species delimitation tests based on the Multi-Species Coalescent (MSC) model, iv) genealogical divergence index (gdi); v) an ad-hoc test of isolation-by-distance within versus between sister-taxa, vi) comparisons of %GC content patterns and vii) morphological analyses. All analyses support the recognition of the undescribed lineage as a putative distinct species. In addition, we show that Eliurus myoxinus, a species known from the dry forests of western Madagascar, is, surprisingly, found mostly in humid forests in northern Madagascar. In conclusion, we discuss the implications of such findings in the context of Eliurus species evolution and diversification, and use the distribution of northern Eliurus species as a proxy for reconstructing past changes in forest cover and vegetation type in northern Madagascar.

Parallel patterns of genetic diversity and structure in circumboreal species of the Sphagnum capillifolium complex.

PREMISE: Shared geographical patterns of population genetic variation among related species is a powerful means to identify the historical events that drive diversification. The Sphagnum capillifolium complex is a group of closely related peat mosses within the Sphagnum subgenus Acutifolia and contains several circumboreal species whose ranges encompass both glaciated and unglaciated regions across the northern hemisphere. In this paper, we (1) inferred the phylogeny of subg. Acutifolia and (2) investigated patterns of population structure and genetic diversity among five circumboreal species within the S. capillifolium complex. METHODS: We generated RAD sequencing data from most species of the subg. Acutifolia and samples from across the distribution ranges of circumboreal species within the S. capillifolium complex. RESULTS: We resolved at least 14 phylogenetic clusters within the S. capillifolium complex. Five circumboreal species show some common patterns: One population system comprises plants in eastern North America and Europe, and another comprises plants in the Pacific Northwest or around the Beringian and Arctic regions. Alaska appears to be a hotspot for genetic admixture, genetic diversity, and sometimes endemic subclades. CONCLUSIONS: Our results support the hypothesis that populations of five circumboreal species within the S. capillifolium complex survived in multiple refugia during the last glacial maximum. Long-distance dispersal out of refugia, population bottlenecks, and possible adaptations to conditions unique to each refugium could have contributed to current geographic patterns. These results indicate the important role of historical events in shaping the complex population structure of plants with broad distribution ranges.

LsMybW-encoding R2R3-MYB transcription factor is responsible for a shift from black to white in lettuce seed.

KEY MESSAGE: We identified LsMybW as the allele responsible for the shift in color from black to white seeds in wild ancestors of lettuce to modern cultivars. Successfully selected white seeds are a key agronomic trait for lettuce cultivation and breeding; however, the mechanism underlying the shift from black-in its wild ancestor-to white seeds remains uncertain. We aimed to identify the gene/s responsible for white seed trait in lettuce. White seeds accumulated less proanthocyanidins than black seeds, similar to the phenotype observed in Arabidopsis TT2 mutants. Genetic mapping of a candidate gene was performed with double-digest RAD sequencing using an F2 population derived from a cross between "ShinanoPower" (white) and "Escort" (black). The white seed trait was controlled by a single recessive locus (48.055-50.197 Mbp) in linkage group 7. Using five PCR-based markers and numerous cultivars, eight candidate genes were mapped in the locus. Only the LG7_v8_49.251Mbp_HinfI marker, employing a single-nucleotide mutation in the stop codon of Lsat_1_v5_gn_7_35020.1, was completely linked to seed color phenotype. In addition, the coding region sequences for other candidate genes were identical in the resequence analysis of "ShinanoPower" and "Escort." Therefore, we proposed Lsat_1_v5_gn_7_35020.1 as the candidate gene and designated it as LsMybW (Lactuca sativa Myb White seeds), an ortholog encoding the R2R3-MYB transcription factor in Arabidopsis. When we validated the role of LsMybW through genome editing, LsMybW knockout mutants harboring an early termination codon showed a change in seed color from black to white. Therefore, LsMybW was the allele responsible for the shift in seed color. The development of a robust marker for marker-assisted selection and identification of the gene responsible for white seeds have implications for future breeding technology and physiological analysis.

Genetic analysis of pungency deficiency in Japanese chili pepper 'Shishito' (Capsicum annuum) revealed its unique heredity and brought the discovery of two genetic loci involved with the reduction of pungency.

The sensation of pungency generated by capsaicinoids is a characteristic trait of chili peppers (Capsicum spp.), and the presence or absence of pungency is central in determining its usage as a spice or a vegetable. In the present study, we aimed to clarify the heredity and genetic factors involved in the deficiency of pungency (quite low pungency) that is uniquely observed in the Japanese chili pepper 'Shishito' (Capsicum annuum). First, the F2 population ('Shishito' × pungent variety 'Takanotsume') was used for segregation analysis, and pungency level was investigated using capsaicinoid quantification with high-performance liquid chromatography. Also, restriction site associated DNA sequencing of the F2 population was performed, and genetic map construction and quantitative trait locus (QTL) mapping were implemented. The results indicated that the F2 population showed varying capsaicinoid content and two major QTLs were detected, Shql3 and Shql7, which explained 39.8 and 19.7% of the genetic variance, respectively. According to these results, the quite low pungency of 'Shishito' was a quantitative trait that involved at least the two loci. Further, this trait was completely separate from general non-pungent traits controlled by individual recessive genes, as described in previous studies. The present study is the first report to investigate the genetic mechanism of pungency deficiency in Japanese chili peppers, and our results provide new insights into the genetic regulation of pungency in chili pepper.

Genomic hotspots of chromosome rearrangements explain conserved synteny despite high rates of chromosome evolution in a holocentric lineage.

Holocentric organisms, unlike typical monocentric organisms, have kinetochore activity distributed along almost the whole length of the chromosome. Because of this, chromosome rearrangements through fission and fusion are more likely to become fixed in holocentric species, which may account for the extraordinary rates of chromosome evolution that many holocentric lineages exhibit. Long blocks of genome synteny have been reported in animals with holocentric chromosomes despite high rates of chromosome rearrangements. Nothing is known from plants, however, despite the fact that holocentricity appears to have played a key role in the diversification of one of the largest angiosperm genera, Carex (Cyperaceae). In the current study, we compared genomes of Carex species and a distantly related Cyperaceae species to characterize conserved and rearranged genome regions. Our analyses span divergence times ranging between 2 and 50 million years. We also compared a C. scoparia chromosome-level genome assembly with a linkage map of the same species to study rearrangements at a population level and suppression of recombination patterns. We found longer genome synteny blocks than expected under a null model of random rearrangement breakpoints, even between very distantly related species. We also found repetitive DNA to be non-randomly associated with holocentromeres and rearranged regions of the genome. The evidence of conserved synteny in sedges despite high rates of chromosome fission and fusion suggests that conserved genomic hotspots of chromosome evolution related to repetitive DNA shape the evolution of recombination, gene order and crossability in sedges. This finding may help explain why sedges are able to maintain species cohesion even in the face of high interspecific chromosome rearrangements.

Transitioning from microsatellites to SNP-based microhaplotypes in genetic monitoring programmes: Lessons from paired data spanning 20 years.

Many long-term genetic monitoring programmes began before next-generation sequencing became widely available. Older programmes can now transition to new marker systems usually consisting of 1000s of SNP loci, but there are still important questions about comparability, precision, and accuracy of key metrics estimated using SNPs. Ideally, transitioned programmes should capitalize on new information without sacrificing continuity of inference across the time series. We combined existing microsatellite-based genetic monitoring information with SNP-based microhaplotypes obtained from archived samples of Rio Grande silvery minnow (Hybognathus amarus) across a 20-year time series to evaluate point estimates and trajectories of key genetic metrics. Demographic and genetic monitoring bracketed multiple collapses of the wild population and included cases where captive-born repatriates comprised the majority of spawners in the wild. Even with smaller sample sizes, microhaplotypes yielded comparable and in some cases more precise estimates of variance genetic effective population size, multilocus heterozygosity and inbreeding compared to microsatellites because many more microhaplotype loci were available. Microhaplotypes also recorded shifts in allele frequencies associated with population bottlenecks. Trends in microhaplotype-based inbreeding metrics were associated with the fraction of hatchery-reared repatriates to the wild and should be incorporated into future genomic monitoring. Although differences in accuracy and precision of some metrics were observed between marker types, biological inferences and management recommendations were consistent.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes).

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide a comprehensive phylogeographic framework for the A. obstetricans complex, as well as a fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

Genome-wide data reveal cryptic diversity and hybridization in a group of tree ferns.

Discovery of cryptic diversity is essential to understanding both the process of speciation and the conservation of species. Determining species boundaries in fern lineages represents a major challenge due to lack of morphologically diagnostic characters and frequent hybridization. Genomic data has substantially enhanced our understanding of the speciation process, increased the resolution of species delimitation studies, and led to the discovery of cryptic diversity. Here, we employed restriction-site-associated DNA sequencing (RAD-seq) and integrated phylogenomic and population genomic analyses to investigate phylogenetic relationships and evolutionary history of 16 tree ferns with marginate scales (Cyatheaceae) from China and Vietnam. We conducted multiple species delimitation analyses using the multispecies coalescent (MSC) model and novel approaches based on genealogical divergence index (gdi) and isolation by distance (IBD). In addition, we inferred species trees using concatenation and several coalescent-based methods, and assessed hybridization patterns and rate of gene flow across the phylogeny. We obtained highly supported and generally congruent phylogenies inferred from concatenated and summary-coalescent methods, and the monophyly of all currently recognized species were strongly supported. Our results revealed substantial evidence of cryptic diversity in three widely distributed Gymnosphaera species, each of which was composite of two highly structure lineages that may correspond to cryptic species. We found that hybridization was fairly common between not only closely related species, but also distantly related species. Collectively, it appears that scaly tree ferns may contain cryptic diversity and hybridization has played an important role throughout the evolutionary history of this group.

Interspecific gene flow obscures phylogenetic relationships in an important insect pest species complex.

As genomic data proliferates, the prevalence of post-speciation gene flow is making species boundaries and relationships increasingly ambiguous. Although current approaches inferring fully bifurcating phylogenies based on concatenated datasets provide simple and robust answers to many species relationships, they may be inaccurate because the models ignore inter-specific gene flow and incomplete lineage sorting. To examine the potential error resulting from ignoring gene flow, we generated both a RAD-seq and a 500 protein-coding loci highly multiplexed amplicon (HiMAP) dataset for a monophyletic group of 12 species defined as the Bactrocera dorsalis sensu lato clade. With some of the world's worst agricultural pests, the taxonomy of the B. dorsalis s.l. clade is important for trade and quarantines. However, taxonomic confusion confounds resolution due to intra- and interspecific phenotypic variation and convergence, mitochondrial introgression across half of the species, and viable hybrids. We compared the topological convergence of our datasets using concatenated phylogenetic and various multispecies coalescent approaches, some of which account for gene flow. All analyses agreed on species delimitation, but there was incongruence between species relationships. Under concatenation, both datasets suggest identical species relationships with mostly high statistical support. However, multispecies coalescent and multispecies network approaches suggest markedly different hypotheses and detected significant gene flow. We suggest that the network approaches are likely more accurate because gene flow violates the assumptions of the concatenated phylogenetic analyses, but the data-reductive requirements of network approaches resulted in reduced statistical support and could not unambiguously resolve gene flow directions. Our study highlights the importance of testing for gene flow, particularly with phylogenomic datasets, even when concatenated approaches receive high statistical support.

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

BACKGROUND AND AIMS: 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. METHODS: 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. KEY RESULTS: 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. CONCLUSIONS: 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.

Phylogenomic and morphological data reveal hidden patterns of diversity in the national tree of Brazil, Paubrasilia echinata.

PREMISE: Paubrasilia echinata (common names, pau brasil, brazilwood) is the national tree of Brazil and an endangered species endemic to the Brazilian Atlantic Forest. Over its wide distribution of 2000 km, its leaflets morphology exhibits extensive plasticity. Three morphotypes are commonly identified based on leaf size, but it is unclear if they represent distinct taxa or a single polymorphic species. This study aims to clarify the taxonomic position of the three morphotypes to inform conservation decisions. METHODS: A morphometric study of leaf characters of herbarium specimens was coupled with genetic analyses using genotype-by-sequencing data. We used maximum-likelihood and coalescent methods to evaluate the phylogenetic and population structure of the species. We compared these with a morphological dendrogram built from hierarchical clustering. RESULTS: Two of the three morphotypes formed separately evolving lineages, the third morphotype formed two geographically separate lineages, and northern trees with intermediate leaf morphology formed a separate fifth lineage. Leaflet size varied by over 35-fold, and although morphological clustering generally matched the genetic patterns, there were some overlaps, highlighting the cryptic diversity within this group. CONCLUSIONS: Our genetic and morphological results provide some evidence that cultivated trees from different states in Brazil seem to have a limited genetic origin and do not reflect the broader genetic and geographical diversity of the species. As a result, more care is likely needed to preserve the overall genomic diversity of this endangered and iconic species.

An updated phylogeny, biogeography, and PhyloCode-based classification of Cornaceae based on three sets of genomic data.

PREMISE: A major goal of systematic biology is to uncover the evolutionary history of organisms and translate that knowledge into stable classification systems. Here, we integrate three sets of genome-wide data to resolve phylogenetic relationships in Cornaceae (containing only Cornus s.l.), reconstruct the biogeographic history of the clade, and provide a revised classification using the PhyloCode to stabilize names for this taxonomically controversial group. METHODS: We conducted phylogenetic analyses using 312 single-copy nuclear genes and 70 plastid genes from Angiosperms353 Hyb-Seq, plus numerous loci from RAD-Seq. We integrated fossils using morphological data and produced a dated phylogeny for biogeographical analysis. RESULTS: A well-resolved, strongly supported, comprehensive phylogeny was obtained. Biogeographic analyses support an origin and rapid diversification of Cornus into four morphologically distinct major clades in the Northern Hemisphere (with an eastern Asian ancestor) during the late Cretaceous. Dispersal into Africa from eastern Asia likely occurred along the Tethys Seaway during the Paleogene, whereas dispersal into South America likely occurred during the Neogene. Diversification within the northern hemisphere likely involved repeated independent colonization of new areas during the Paleogene and Neogene along the Bering Land Bridge, the North Atlantic Land Bridge, and the Tethys Seaway. Thirteen strongly supported clades were named following rules of the PhyloCode. CONCLUSIONS: Our study provides an example of integrating genomic and morphological data to produce a robust, explicit species phylogeny that includes fossil taxa, which we translate into an updated classification scheme using the PhyloCode to stabilize names.

Genome-wide analysis facilitates estimation of the amount of male contribution in meiotic gynogenetic three-spined stickleback (Gasterosteus aculeatus).

Gynogenetic embryos - those inheriting only maternal DNA - can be experimentally created by fertilizing eggs with radiation-treated sperm containing inactivated paternal chromosomes. Diploidy in the zygotes can be maintained through prevention of the second meiosis or restored by preventing the first mitosis after the maternal chromosome complement has been replicated. These gynogenetic organisms are useful in many fields including aquaculture, evolutionary biology and genomics. Although gynogenetic organisms have been created in numerous species, the completeness of uni-parental inheritance has often been assumed rather than thoroughly quantified across the genome. Instead, when tests of uni-parental inheritance occur, they typically rely on well-studied genetically determined phenotypes that represent a very small sub-set of the genome. Only assessing small genomic regions for paternal inheritance leaves the question of whether some paternal contributions to offspring might still have occurred. In this study, the authors quantify the efficacy of creating gynogenetic diploid three-spined stickleback fish (Gasterosteus aculeatus). To this end, the authors mirrored previous assessments of paternal contribution using well-studied genetically determined phenotypes including sex and genetically dominant morphological traits but expanded on previous studies using dense restriction site-associated DNA sequencing (RAD-seq) markers in parents and offspring to assess paternal inheritance genome-wide. In the gynogenetic diploids, the authors found no male genotypes underlying their phenotypes of interest - sex and dominant phenotypic traits. Using genome-wide assessments of paternal contribution, nevertheless, the authors found evidence of a small, yet potentially important, amount of paternally "leaked" genetic material. The application of this genome-wide approach identifies the need for more widespread assessment of paternal contributions to gynogenetic animals and promises benefits for many aspects of aquaculture, evolutionary biology and genomics.

Migration through a Major Andean Ecogeographic Disruption as a Driver of Genetic and Phenotypic Diversity in a Wild Tomato Species.

Evolutionary dynamics at the population level play a central role in creating the diversity of life on our planet. In this study, we sought to understand the origins of such population-level variation in mating systems and defensive acylsugar chemistry in Solanum habrochaites-a wild tomato species found in diverse Andean habitats in Ecuador and Peru. Using Restriction-site-Associated-DNA-Sequencing (RAD-seq) of 50 S. habrochaites accessions, we identified eight population clusters generated via isolation and hybridization dynamics of 4-6 ancestral populations. Detailed characterization of mating systems of these clusters revealed emergence of multiple self-compatible (SC) groups from progenitor self-incompatible populations in the northern part of the species range. Emergence of these SC groups was also associated with fixation of deleterious alleles inactivating acylsugar acetylation. The Amotape-Huancabamba Zone-a geographical landmark in the Andes with high endemism and isolated microhabitats-was identified as a major driver of differentiation in the northern species range, whereas large geographical distances contributed to population structure and evolution of a novel SC group in the central and southern parts of the range, where the species was also inferred to have originated. Findings presented here highlight the role of the diverse ecogeography of Peru and Ecuador in generating population differentiation, and enhance our understanding of the microevolutionary processes that create biological diversity.

RAD-seq as an effective strategy for heterogenous variety identification in plants-a case study in Italian Ryegrass (Lolium multiflorum).

The primary approach for variety distinction in Italian ryegrass is currently the DUS (distinctness, uniformity and stability) test based on phenotypic traits. Considering the diverse genetic background within the population and the complexity of the environment, however, it is challenging to accurately distinguish varieties based on DUS criteria alone. In this study, we proposed the application of high-throughput RAD-seq to distinguish 11 Italian ryegrass varieties with three bulks of 50 individuals per variety. Our findings revealed significant differences among the 11 tested varieties. The PCA, DAPC and STRUCTURE analysis indicated a heterogeneous genetic background for all of them, and the AMOVA analysis also showed large genetic variance among these varieties (ΦST = 0.373), which were clearly distinguished based on phylogenetic analysis. Further nucleotide diversity (Pi) analysis showed that the variety 'Changjiang No.2' had the best intra-variety consistency among 11 tested varieties. Our findings suggest that the RAD-seq could be an effectively alternative method for the variety distinction of Italian ryegrass, as well as a potential tool for open-pollinated varieties (OPVs) of other allogamous species.

Untying Gordian knots: unraveling reticulate polyploid plant evolution by genomic data using the large Ranunculus auricomus species complex.

Speciation via hybridization and polyploidization is a major evolutionary force in plant evolution but is still poorly understood for neopolyploid groups. Challenges are attributed to high heterozygosity, low genetic divergence, and missing information on progenitors, ploidy, and reproduction. We study the large Eurasian Ranunculus auricomus species complex and use a comprehensive workflow integrating reduced-representation sequencing (RRS) genomic data to unravel reticulate evolution, genome diversity and composition of polyploids. We rely on 97 312 restriction site-associated DNA sequencing (RAD-Seq) loci, 576 targeted nuclear genes (48 phased), and 71 plastid regions derived from 78 polyploid apomictic taxa and four diploid and one tetraploid putative sexual progenitor species. We applied (phylo)genomic structure, network, and single nucleotide polymorphism (SNP)-origin analyses. Results consistently showed only 3-5 supported and geographically structured polyploid genetic groups, each containing extant sexual and one unknown progenitor species. Combined analyses demonstrated predominantly allopolyploid origins, each involving 2-3 different diploid sexual progenitor species. Young allotetraploids were characterized by subgenome dominance and nonhybrid SNPs, suggesting substantial post-origin but little lineage-specific evolution. The biodiversity of neopolyploid complexes can result from multiple hybrid origins involving different progenitors and substantial post-origin evolution (e.g. homoeologous exchanges, hybrid segregation, gene flow). Reduced-representation sequencing genomic data including multi-approach information is efficient to delimit shallow reticulate relationships.

Ecological adaptation shaped the genetic structure of homoploid ferns against strong dispersal capacity.

The formation of spatial genetic structure with the presence of extensive gene flow, an evolutionary force which is generally expected to eliminate population-specific variation and maintain genetic homogeneity, remains poorly understood. Homosporous ferns, which spread by spores through wind and possess long distance dispersal capacity, provide an ideal system to investigate such a process. Here, using a homoploid fern lineage, the Athyrium sinense complex, we used reduced-representation genomic data to examine spatial genetic structure and explored potential driving forces including geographical distance, environment, climatic history and external dispersal constraints. Our findings showed a clear north-south divergence at the genetic, morphological and ecological levels between both sides of 35°N in East Asia. Fluctuant and heterogeneous climatic condition was demonstrated to play a crucial role during the formation of the divergence. Our results suggested that this lineage was able to migrate southward and colonize new habitat as a result of the Quaternary climatic fluctuation. Furthermore, the present genetic structure is attributed to adaptation to heterogeneous environments, especially temperature difference. In addition to ecological adaptation, we found clues showing that canopy density, wind direction as well as habitat continuity were all likely to constrain the effect of gene flow. These results demonstrated a diversification process without ploidy changes in ferns providing new insights for our present knowledge on ferns' spatio-temporal evolutionary pattern. In particular, our study highlights the influence of environmental heterogeneity in driving genetic divergence against strong dispersal capacity.