Genomic Differentiation Research Articles

Genome size of grass Festuca mountain species from the southwestern European Pyrenees: variation, evolution, and new assessments

Festuca is one of the most ecologically and economically important genera of temperate grasses. Species of its type subgenus, Festuca, are common components of cool-seasonal pastures and are especially frequent in mountainous areas, where there are often several morphologically similar species that grow in the same or comparable habitats and sometimes live in sympatry. Nuclear DNA assessments by flow cytometry can be used to identify taxa and uncover new ploidy levels in species complexes for which new chromosome data are provided or previous chromosome counts and genome sizes are known. Holoploid (2C) values of newly studied Pyrenean Festuca subgen. Festuca sects. Eskia, Festuca and Aulaxyper species fall within the expected ranges for these taxonomic groups and include 2x, 4x, 6x and 8x ploidy levels. Monoploid (1Cx) genome sizes of diploids and polyploids are larger in the species of the more ancestral F. sect. Eskia group showing a decreasing trend in the species of the more recently evolved F. sects. Festuca and Aulaxyper lineages. 1Cx values of high polyploid Aulaxyper taxa are among the smallest of the three Festuca sections, corroborating previous findings. Our analysis provides new genome size values and inferred ploidy levels for hexaploid F.heteromalla and octoploid F.trichophylla and highlights the genomic and ecological differentiation of tetraploid F.gautieri susbsp. gautieri from diploid F.gautieri subsp. scoparia.

Generation of a human haploid neural stem cell line for genome-wide genetic screening.

Haploid embryonic stem cells (haESCs) have been established in many species. Differentiated haploid cell line types in mammals are lacking due to spontaneous diploidization during differentiation that compromises lineage-specific screens. To derive human haploid neural stem cells (haNSCs) to carry out lineage-specific screens. Human haNSCs were differentiated from human extended haESCs with the help of Y27632 (ROCK signaling pathway inhibitor) and a series of cytokines to reduce diploidization. Neuronal differentiation of haNSCs was performed to examine their neural differentiation potency. Global gene expression analysis was con-ducted to compare haNSCs with diploid NSCs and haESCs. Fluorescence activated cell sorting was performed to assess the diploidization rate of extended haESCs and haNSCs. Genetic manipulation and screening were utilized to evaluate the significance of human haNSCs as genetic screening tools. Human haESCs in extended pluripotent culture medium showed more compact and smaller colonies, a higher efficiency in neural differentiation, a higher cell survival ratio and higher stability in haploidy maintenance. These characteristics effectively facilitated the derivation of human haNSCs. These human haNSCs can be generated by differentiation and maintain haploidy and multipotency to neurons and glia in the long term in vitro. After PiggyBac transfection, there were multiple insertion sites in the human haNSCs' genome, and the insertion sites were evenly spread across all chromosomes. In addition, after the cells were treated with manganese, we were able to generate a list of manganese-induced toxicity genes, demonstrating their utility as genetic screening tools. This is the first report of a generated human haploid somatic cell line with a complete genome, proliferative ability and neural differentiation potential that provides cell resources for recessive inheritance and drug targeted screening.

Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss).

New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae. We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally. Using sliding-window population genetic analyses we find that differentiation is concentrated within 'islands' of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap. We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth's climate rapidly changes.

Whole-genome sequencing reveals fine-scale environment-associated divergence near the range limits of a temperate reef fish.

Environmental variation is increasingly recognized as an important driver of diversity in marine species despite the lack of physical barriers to dispersal and the presence of pelagic stages in many taxa. A robust understanding of the genomic and ecological processes involved in structuring populations is lacking for most marine species, often hindering management and conservation action. Cunner (Tautogolabrus adspersus) is a temperate reef fish with both pelagic early life-history stages and strong site-associated homing as adults; the species is also of interest for use as a cleaner fish in salmonid aquaculture in Atlantic Canada. We aimed to characterize genomic and geographic differentiation of cunner in the Northwest Atlantic. To achieve this, a chromosome-level genome assembly for cunner was produced and used to characterize spatial population structure throughout Atlantic Canada using whole-genome sequencing. The genome assembly spanned 0.72 Gbp and 24 chromosomes; whole-genome sequencing of 803 individuals from 20 locations from Newfoundland to New Jersey identified approximately 11 million genetic variants. Principal component analysis revealed four regional Atlantic Canadian groups. Pairwise FST and selection scans revealed signals of differentiation and selection at discrete genomic regions, including adjacent peaks on chromosome 10 across multiple pairwise comparisons (i.e. FST 0.5-0.75). Redundancy analysis suggested association of environmental variables related to benthic temperature and oxygen range with genomic structure. Results suggest regional scale diversity in this temperate reef fish and can directly inform the collection and translocation of cunner for aquaculture applications and the conservation of wild populations throughout the Northwest Atlantic.

The Metasequoia genome and evolutionary relationships among redwoods

Redwood trees (Sequoioideae), including Metasequoia glyptostroboides (dawn redwood), Sequoiadendron giganteum (giant sequoia), and Sequoia sempervirens (coast redwood), are threatened and widely recognized iconic tree species. Genomic resources for redwood trees could provide clues to their evolutionary relationships. Here, we report the 8-Gb reference genome of M. glyptostroboides and a comparative analysis with two related species. More than 62% of the M. glyptostroboides genome is composed of repetitive sequences. Clade-specific bursts of long terminal repeat retrotransposons may have contributed to genomic differentiation in the three species. The chromosomal synteny between M. glyptostroboides and S. giganteum is extremely high, whereas there has been significant chromosome reorganization in S. sempervirens. Phylogenetic analysis of marker genes indicates that S. sempervirens is an autopolyploid, and more than 48% of the gene trees are incongruent with the species tree. Results of multiple analyses suggest that incomplete lineage sorting (ILS) rather than hybridization explains the inconsistent phylogeny, indicating that genetic variation among redwoods may be due to random retention of polymorphisms in ancestral populations. Functional analysis of ortholog groups indicates that gene families of ion channels, tannin biosynthesis enzymes, and transcription factors for meristem maintenance have expanded in S. giganteum and S. sempervirens, which is consistent with their extreme height. As a wetland-tolerant species, M. glyptostroboides shows a transcriptional response to flooding stress that is conserved with that of analyzed angiosperm species. Our study offers insights into redwood evolution and adaptation and provides genomic resources to aid in their conservation and management.

O-107 Sperm chromosome chaos induced by MEIKIN gene mutations leads to blastocyst formation failure

Abstract Study question What are the effect and corresponding mechanisms of chromosome chaos of sperm caused by mutation of the MEIKIN gene on preimplantation embryo development? Summary answer Sperm chromosome chaos induced by MEIKIN gene mutations directly causes blastocyst formation failure by affecting the activation of the embryonic genome and lineage differentiation. What is known already In mice, Meikin is a meiosis-specific kinetochore factor that plays a crucial role in mono-orientation and protection of centromeric cohesion in meiosis I and the regulation of chromosome alignment during meiosis II. Both male and female Meikin gene knockout mice are completely infertile. Infertility in male Meikin-/- mice originates from defects in meiotic chromosome segregation and subsequent spermatogenesis. To our knowledge, no mutation of the MEIKIN gene has been reported in humans. Study design, size, duration In this study, semen from three patients with MEIKIN mutations and arrested embryos from one of them were analyzed to determine the effect of MEIKIN mutations on sperm and early embryonic development. Participants/materials, setting, methods Three male patients with MEIKIN mutations with a specific clinical phenotype were recruited. We assessed mutant sperm morphology using Papanicolaou staining. Sperm FISH and single sperm chromosome aneuploidy analysis（CNV-seq） were applied to identify mutant sperm Chromosomal abnormality. We also investigated the chromosome constitution and gene expression of the arrested embryos from one patient by parallel single-cell genome and transcriptome sequencing. Main results and the role of chance Novel mutations in the MEIKIN gene were identified in three infertility males with a specific clinical phenotype（They had good-quality embryos on day 3, but these embryos repeatedly failed to implant or develop to the blastocyst stage）. Routine semen assessments were conducted based on the WHO’s guidelines, and all these patients were diagnosed with oligoasthenoteratozoospermia. Sperm FISH showed that the frequencies of aneuploidies in spermatozoa were significantly higher than the normal control（89.05%,70.9%,93.24% vs. 5.03%）, and the CNV-seq analysis of single spermatozoa confirmed mutant sperm had severe chromosome chaos. Each cell of the arrested embryos from one patient had complicated chromosomal abnormalities, and the abnormal chromosomes were mainly derived from sperm. At the gene expression level, these embryos were arrested at the 8-cell to morula stage transition, and defects activation of a large ZGA and lineage specification genes led to blastocyst formation failure. Limitations, reasons for caution In this study, we discovered that mutations in MEIKIN caused male infertility manifesting as sperm chromosome chaos and blastocyst formation failure. However, the molecular mechanism of the sperm chromosome chaos involved in blastocyst formation needs to be further illuminated by using Meikin knockout mouse models. Wider implications of the findings We first reported that male with MEIKIN mutations is related to severe sperm chromosome abnormalities and blastocyst formation failure. Therefore, In the process of ICSI, if there are good-quality cleavage embryos that repeatedly fail to implant or develop to the blastocyst stage, it is recommended to examine male factors concurrently Trial registration number not applicable

Genetic differentiation following recent domestication events: A study of farmed Nile tilapia (Oreochromis niloticus) populations.

Nile tilapia (Oreochromis niloticus) is among the most farmed finfish worldwide, distributed across different environmental conditions. Its wide distribution has mainly been facilitated by several breeding programs and widespread dissemination of genetically improved strains. In the first Nile tilapia study exploiting a whole-genome pooled sequencing (Poolseq) approach, we identified the genetic structure and signatures of selection in diverse, farmed Nile tilapia populations, with a particular focus on the GIFT strain, developed in the 1980s, and currently managed by WorldFish (GIFTw). We also investigated important farmed strains from The Philippines and Africa. Using both SNP array data and Poolseq SNPs, we characterized the population structure of these samples. We observed the greatest separation between the Asian and African populations and greater admixture in the Asian populations than in the African ones. We also established that the SNP array data were able to successfully resolve relationships between these diverse Nile tilapia populations. The Poolseq data identified genomic regions with high levels of differentiation (F ST) between GIFTw and the other populations. Gene ontology terms associated with mesoderm development were significantly enriched in the genes located in these regions. A region on chromosome Oni06 was genetically differentiated in pairwise comparisons between GIFTw and all other populations. This region contains genes associated with muscle-related traits and overlaps with a previously published QTL for fillet yield, suggesting that these traits may have been direct targets for selection on GIFT. A nearby region was also identified using XP-EHH to detect genomic differentiation using the SNP array data. Genomic regions with high or extended homozygosity within each population were also identified. This study provides putative genomic landmarks associated with the recent domestication process in several Nile tilapia populations, which could help to inform their genetic management and improvement.

Comparative genomics of two Empidonax flycatchers reveal candidate genes for bird song production.

Whole-genome level comparisons of sister taxa that vary in phenotype against a background of high genomic similarity can be used to identify the genomic regions that might underlie their phenotypic differences. In wild birds, this exploratory approach has detected markers associated with plumage coloration, beak and wing morphology, and complex behavioral traits like migration. Here, we use genomic comparisons of two closely related suboscine flycatchers (Empidonax difficilis and E. occidentalis) and their hybrids to search for candidate genes underlying their variation in innate vocal signals. We sequenced the genomes of 20 flycatchers that sang one of two species-specific parental song types and 14 putative hybrid individuals with intermediate song types. In the resulting genomic comparisons, we found six areas of high differentiation that may be associated with variation in non-learned songs. These narrow regions of genomic differentiation contain a total of 67 described genes, three of which have been previously associated with forms of language impairment and dyslexia in humans, and 18 of which are known to be differentially expressed in the song nuclei regions of the avian brain compared to adjacent parts of the avian brain. This 'natural experiment' therefore may help identify loci associated with song differences that merit further study across bird lineages with both learned and innate vocalizations.

Semipermeable species boundaries create opportunities for gene flow and adaptive potential.

Hybridisation and gene flow can have both deleterious and adaptive consequences for natural populations and species. To better understand the extent of hybridisation in nature and the balance between its beneficial and deleterious outcomes in a changing environment, information on naturally hybridising nonmodel organisms is needed. This requires the characterisation of the structure and extent of natural hybrid zones. Here, we study natural populations of five keystone mound-building wood ant species in the Formica rufa group across Finland. No genomic studies across the species group exist, and the extent of hybridisation and genomic differentiation in sympatry is unknown. Combining genome-wide and morphological data, we demonstrate more extensive hybridisation than was previously detected between all five species in Finland. Specifically, we reveal a mosaic hybrid zone between Formica aquilonia, F. rufa and F. polyctena, comprising further generation hybrid populations. Despite this, we find that F. rufa, F. aquilonia, F. lugubris and F. pratensis form distinct gene pools in Finland. We also find that hybrids occupy warmer microhabitats than the nonadmixed populations of cold-adapted F. aquilonia, and suggest that warm winters and springs, in particular, may benefit hybrids over F. aquilonia, the most abundant F. rufa group species in Finland. In summary, our results indicate that extensive hybridisation may create adaptive potential that could promote wood ant persistence in a changing climate. Additionally, they highlight the potentially significant ecological and evolutionary consequences of extensive mosaic hybrid zones, within which independent hybrid populations face an array of ecological and intrinsic selection pressures.

Incomplete reproductive barriers and genomic differentiation impact the spread of resistance mutations between green- and red-colour morphs of a cosmopolitan mite pest.

Pesticide resistance represents a clear and trackable case of adaptive evolution with a strong societal impact. Understanding the factors associated with the evolution and spread of resistance is imperative to develop sustainable crop management strategies. The two-spotted spider mite Tetranychus urticae, a major crop pest with worldwide distribution and a polyphagous lifestyle, has evolved resistance to most classes of pesticides. Tetranychus urticae exists as either a green- or a red-coloured morph. However, the extent of genetic divergence and reproductive compatibility vary across populations of these colour morphs, complicating their taxonomic resolution at the species level. Here, we studied patterns of genetic differentiation and barriers to gene flow within and between morphs of T. urticae in order to understand the factors that influence the spread of resistance mutations across its populations. We derived multiple iso-female lines from Tetranychus populations collected from agricultural crops. We generated genomic and morphological data, characterized their bacterial communities and performed controlled crosses. Despite morphological similarities, we found large genomic differentiation between the morphs. This pattern was reflected in the incomplete, but strong postzygotic incompatibility in crosses between colour morphs, while crosses within morphs from different geographical locations were largely compatible. In addition, our results suggest recent/on-going gene flow between green-coloured T. urticae and T. turkestani. By screening the sequences of 10 resistance genes, we found evidence for multiple independent origins and for single evolutionary origins of target-site resistance mutations. Our results indicate that target-site mutations mostly evolve independently in populations on different geographical locations, and that these mutations can spread due to incomplete barriers to gene flow within and between populations.

Genome diversity and highland-adaptative variation in Tibet barley landrace population of China.

Barley landraces accumulated variation in adapting to extreme highland environments during long-term domestication in Tibet, but little is known about their population structure and genomic selection traces. In this study, tGBS (tunable genotyping by sequencing) sequencing, molecular marker and phenotypic analyses were conducted on 1,308 highland and 58 inland barley landraces in China. The accessions were divided into six sub-populations and clearly distinguished most six-rowed, naked barley accessions (Qingke in Tibet) from inland barley. Genome-wide differentiation was observed in all five sub-populations of Qingke and inland barley accessions. High genetic differentiation in the pericentric regions of chromosomes 2H and 3H contributed to formation of five types of Qingke. Ten haplotypes of the pericentric regions of 2H, 3H, 6H and 7H were further identified as associated with ecological diversification of these sub-populations. There was genetic exchange between eastern and western Qingke but they shared the same progenitor. The identification of 20 inland barley types indicated multiple origins of Qingke in Tibet. The distribution of the five types of Qingke corresponded to specific environments. Two predominant highland-adaptative variations were identified for low temperature tolerance and grain color. Our results provide new insights into the origin, genome differentiation, population structure and highland adaptation in highland barley which will benefit both germplasm enhancement and breeding of naked barley.

Giant Sex Chromosomes in Omophoita Species (Coleoptera, Chrysomelidae): Structural and Evolutionary Relationships Revealed by Zoo-FISH and Comparative Genomic Hybridization (CGH)

The beetles of the subtribe Oedionychina (Chrysomelidae, Alticinae) are the only ones that have the atypical giant and achiasmatic sex chromosomes, which are substantially larger than the autosomes. Previous cytogenetic analyses suggest a large accumulation of repetitive DNA in the sex chromosomes. In this study, we examined the similarity of X and Y chromosomes in four Omophoita species and compared genomic differentiation to better understand the evolutionary process and the giant sex chromosomes origin. Intraspecific genomic comparation using male and female genomes of O. octoguttata and interespecific analyses using genomic DNA of O. octoguttata, O. sexnotata, O. magniguttis, and O. personata were performed. In addition, whole chromosome painting (WCP) experiments were performed with X and Y chromosome probes of O. octogutatta. CGH analysis revealed great genomic similarity between the sexes and a sex-specific region on the Y chromosome, and interspecific analysis revealed a genomic divergence between species. In contrast, WCP results revealed that the sex chromosomes of O. octoguttata have high intra- and interspecific similarity with the studied species. Our data support a common origin under the canonical evolution of the sex chromosomes in this group, as they have high genomic similarity between them.

Highly differentiated genomic properties underpin the different cell walls of Poaceae and eudicots.

Plant cell walls of Poaceae and eudicots differ substantially, both in the content and composition of their components. However, the genomic and genetic basis underlying these differences is not fully resolved. In this research, we analyzed multiple genomic properties of 150 cell wall gene families across 169 angiosperm genomes. The properties analyzed include gene presence/absence, copy number, synteny, occurrence of tandem gene clusters, and phylogenetic gene diversity. Results revealed a profound genomic differentiation of cell wall genes between Poaceae and eudicots, often associated with the cell wall diversity between these plant groups. For example, overall patterns of gene copy number variation and synteny were clearly divergent between Poaceae and eudicot species. Moreover, differential Poaceae-eudicot copy number and genomic contexts were observed for all the genes within the BEL1-like HOMEODOMAIN 6 regulatory pathway, which respectively induces and represses secondary cell wall synthesis in Poaceae and eudicots. Similarly, divergent synteny, copy number, and phylogenetic gene diversification were observed for the major biosynthetic genes of xyloglucans, mannans, and xylans, potentially contributing to the differences in content and types of hemicellulosic polysaccharides differences in Poaceae and eudicot cell walls. Additionally, the Poaceae-specific tandem clusters and/or higher copy number of PHENYLALANINE AMMONIA-LYASE, CAFFEIC ACID O-METHYLTRANSFERASE, or PEROXIDASE genes may underly the higher content and larger variety of phenylpropanoid compounds observed in Poaceae cell walls. All these patterns are discussed in detail in this study, along with their evolutionary and biological relevance for cell wall (genomic) diversification between Poaceae and eudicots.

Biogeographic provinces and genomically delineated stocks are congruent in snapper (Chrysophrys auratus) from southeastern Australia

AbstractMarine species often exhibit genetic discontinuities concordant with biogeographic boundaries, frequently occurring due to changes in ocean circulation, bathymetry, coastline topography, and temperature. Here, we used 10 916 single nucleotide polymorphisms (SNPs) to assess the concordance between population genomic differentiation and coastal biogeography in the fishery-important snapper (Chrysophrys auratus) across southeastern Australia. Additionally, we investigated whether spatial scales of assessment and management of snapper align with evidence from population genomics. We detected genomic structure concordant with the region's three biogeographic provinces across snapper from 11 localities (n = 488) between the west coast of South Australia and the south coast of New South Wales. We also detected fine-scale genetic structuring relating to spatial variation in spawning and recruitment dynamics, as well as temporal stability in the genomic signal associated with two important spawning areas. The current management boundaries in the region coincided with either the genetic breaks at bioregional boundaries or with local-scale variation. Our study highlights the value of population genomic surveys in species with high dispersal potential for uncovering stock boundaries and demographic variation related to spawning and recruitment. It also illustrates the importance of marine biogeography in shaping population structure in commercial species with high dispersal potential.

Whole Genome Resequencing Reveals Genetic Diversity and Selection Signatures of Ethiopian Indigenous Cattle Adapted to Local Environments

Cattle are among the most important domesticated bovid species in the world, of which Ethiopia possesses large populations adapted to different agro-ecologies and production systems. Though several molecular population genetic studies have been done on Ethiopian indigenous cattle, genomic diversity and selection signatures pertinent to adaptation to the different local environments have yet to be comprehensively characterized. Here, the whole genome sequences of 151 samples from 14 Ethiopian indigenous cattle populations were analyzed to assess genomic diversity and differentiation as well as signatures of positive selection (using Hp, iHS, FST, and XP-CLR) in comparison to Sudanese zebu, Asian zebu, Ankole, and African and European taurine cattle. High genomic differentiation was observed between Ethiopian and non-Ethiopian cattle populations, while low genomic differentiation and inbreeding were present between and within Ethiopian cattle populations. Sixteen genome regions overlapping with 40 candidate genes were commonly identified by at least three genome scan methods. High frequencies of missense SNPs in ITPR2, CHADL, GNAS, STING1, and KIT genes with high haplotype differentiations were observed in Ethiopian cattle compared to non-Ethiopian cattle. The candidate genes were significantly associated with several biological functions and molecular pathways responsible for nutrient metabolism, skeletal development, immune response, reproduction, water balance, coat color pigmentation, and circulatory homeostasis. Our results provide new insights into the adaptation of the Ethiopian indigenous cattle to the country’s diverse environments.

Isolation by resistance explains genetic diversity in the Arremon brushfinches of northern Mesoamerica.

Genetic differentiation between and within natural populations is the result of the joint effects of neutral and adaptative processes. In addition, the spatial arrangement of the landscape promotes connectivity or creates barriers to gene flow, directly affecting speciation processes. In this study, we carried out a landscape genomics analysis using NextRAD data from a montane forest specialist bird complex, the Mesoamerican Chestnut-capped/Green-striped Brushfinch of the genus Arremon. Specifically, we examined population genomic structure using different assignment methods and genomic differentiation and diversity, and we tested alternative genetic isolation hypotheses at the individual level (e.g., isolation by barrier, IBB; isolation by environment, IBE; isolation by resistance, IBR). We found well-delimited genomic structuring (K = 5) across Mesoamerican montane forests in the studied group. Individual-level genetic distances among major montane ranges were mainly explained by IBR hypotheses in this sedentary Neotropical taxon. Our results uncover genetic distances/differentiation and patterns of gene flow in allopatric species that support the role of tropical mountains as spatial landscape drivers of biodiversity. IBR clearly supports a pattern of conserved niche-tracking of suitable habitat conditions and topographic complexity throughout glacial-interglacial dynamics.

Recurrent selection and reduction in recombination shape the genomic landscape of divergence across multiple population pairs of Green-backed Tit.

Speciation is fundamental for building and maintaining biodiversity. The formation of the highly differentiated genomic regions between diverging taxa has been interpreted as a result of divergence with gene flow, linked selection, and reduction in recombination. It is challenging to unravel these nonexclusive processes in shaping genomic divergence. Here, we investigate the relative roles of these processes in shaping genomic differentiation in a montane bird, the Green-backed Tit (Parus monticolus). Our genetic structure and demographic analyses identify that four genetic lineages diverge between 838 and 113 thousand years ago and there is evidence of secondary gene flow. The highly divergent genomic regions do not increase with the divergence time, as we found that the old lineages show relatively fewer numbers and smaller sizes of highly differentiated regions than the young divergent lineages (numbers, 118-138 vs. 156-289; sizes, 5.9-6.9 vs. 7.8-14.5 megabase). Across the genome, the outlier windows show a reduction in nucleotide diversity, absolute genetic divergence, and recombination rate, suggesting recurrent selection in regions with low recombination being the major driver of genomic divergence. Finally, we show that secondary gene flow tends to affect the highly differentiated genomic regions if these regions are less likely to be the targets of selection. Altogether, our study shows how common ancestry, recurrent selection, low recombination rate, and gene flow have contributed to the emergence of genomic islands at different stages of speciation.

Geographic-genomic and geographic-phenotypic differentiation of the Aquilegia viridiflora complex.

How species diverge into different lineages is a central issue in evolutionary biology. Despite the increasing evidence indicating that such divergences do not need geographic isolation, the correlation between lineage divergence and the adaptive ecological divergence of phenotype corresponding to distribution is still unknown. In addition, gene flow has been widely detected during and through such diverging processes. We used one widely distributed Aquilegia viridiflora complex as a model system to examine genomic differentiation and corresponding phenotypic variations along geographic gradients. Our phenotypic analyses of 20 populations from northwest to northeast China identified two phenotypic groups along the geographic cline. All examined traits are distinct from each other, although a few intermediate individuals occur in their contacting regions. We further sequenced the genomes of representative individuals of each population. However, four distinct genetic lineages were detected based on nuclear genomes. In particular, we recovered numerous genetic hybrids in the contact regions of four lineages. Gene flow is widespread and continuous between four lineages but much higher between contacting lineages than geographically isolated lineages. Gene flow and natural selection might result in inconsistency between heredity and phenotype. Moreover, many genes with fast lineage-specific mutations were identified to be involved in local adaptation. Our results suggest that both geographic isolation and local selection exerted by the environment and pollinators may together create geographic distributions of phenotypic variations as well as the underlying genomic divergences in numerous lineages.

Where to draw the line? Expanding the delineation of conservation units to highly mobile taxa.

Conservation units (CUs) are an essential tool for maximizing evolutionary potential and prioritizing areas across a species' range for protection when implementing conservation and management measures. However, current workflows for identifying CUs on the basis of neutral and adaptive genomic variation largely ignore information contained in patterns of isolation by distance (IBD), frequently the primary signal of population structure in highly mobile taxa, such as birds, bats, and marine organisms with pelagic larval stages. While individuals located on either end of a species' distribution may exhibit clear genetic, phenotypic, and ecological differences, IBD produces subtle changes in allele frequencies across space, making it difficult to draw clear boundaries for conservation purposes in the absence of discrete population structure. Here, we highlight potential pitfalls that arise when applying common methods for delineating CUs to continuously distributed organisms and review existing methods for detecting subtle breakpoints in patterns of IBD that can indicate barriers to gene flow in highly mobile taxa. In addition, we propose a new framework for identifying CUs in all organisms, including those characterized by continuous genomic differentiation, and suggest several possible ways to harness the information contained in patterns of IBD to guide conservation and management decisions.

Putative Signals of Generalist Plant Species Adaptation to Local Pollinator Communities and Abiotic Factors.

The reproductive success of flowering plants with generalized pollination systems is influenced by interactions with a diverse pollinator community and abiotic factors. However, knowledge about the adaptative potential of plants to complex ecological networks and the underlying genetic mechanisms is still limited. Based on a pool-sequencing approach of 21 natural populations of Brassica incana in Southern Italy, we combined a genome-environmental association analysis with a genome scan for signals of population genomic differentiation to discover genetic variants associated with the ecological variation. We identified genomic regions putatively involved in the adaptation of B. incana to the identity of local pollinator functional categories and pollinator community composition. Interestingly, we observed several shared candidate genes associated with long-tongue bees, soil texture, and temperature variation. We established a genomic map of potential generalist flowering plant local adaptation to complex biotic interactions, and the importance of considering multiple environmental factors to describe the adaptive landscape of plant populations.