The influence of population structure on gene expression and flowering time variation in the ubiquitous weedCapsella bursa-pastoris(Brassicaceae)

Origin and diversity of Capsella bursa-pastoris from the genomic point of view

BACKGROUND

Capsella bursa-pastoris, a cosmopolitan weed of hybrid origin, is an emerging model object for the study of early consequences of polyploidy, being a fast growing annual and a close relative of Arabidopsis thaliana. The development of this model is hampered by the absence of a reference genome sequence.

RESULTS

We present here a subgenome-resolved chromosome-scale assembly and a genetic map of the genome of Capsella bursa-pastoris. It shows that the subgenomes are mostly colinear, with no massive deletions, insertions, or rearrangements in any of them. A subgenome-aware annotation reveals the lack of genome dominance-both subgenomes carry similar number of genes. While most chromosomes can be unambiguously recognized as derived from either paternal or maternal parent, we also found homeologous exchange between two chromosomes. It led to an emergence of two hybrid chromosomes; this event is shared between distant populations of C. bursa-pastoris. The whole-genome analysis of 119 samples belonging to C. bursa-pastoris and its parental species C. grandiflora/rubella and C. orientalis reveals introgression from C. orientalis but not from C. grandiflora/rubella.

CONCLUSIONS

C. bursa-pastoris does not show genome dominance. In the earliest stages of evolution of this species, a homeologous exchange occurred; its presence in all present-day populations of C. bursa-pastoris indicates on a single origin of this species. The evidence coming from whole-genome analysis challenges the current view that C. grandiflora/rubella was a direct progenitor of C. bursa-pastoris; we hypothesize that it was an extinct (or undiscovered) species sister to C. grandiflora/rubella.

Separating phases of allopolyploid evolution with resynthesized and natural Capsella bursa-pastoris

Allopolyploidization is a frequent evolutionary transition in plants that combines whole-genome duplication (WGD) and interspecific hybridization. The genome of an allopolyploid species results from initial interactions between parental genomes and long-term evolution. Distinguishing the contributions of these two phases is essential to understanding the evolutionary trajectory of allopolyploid species. Here, we compared phenotypic and transcriptomic changes in natural and resynthesized Capsella allotetraploids with their diploid parental species. We focused on phenotypic traits associated with the selfing syndrome and on transcription-level phenomena such as expression-level dominance (ELD), transgressive expression (TRE), and homoeolog expression bias (HEB). We found that selfing syndrome, high pollen, and seed quality in natural allotetraploids likely resulted from long-term evolution. Similarly, TRE and most down-regulated ELD were only found in natural allopolyploids. Natural allotetraploids also had more ELD toward the self-fertilizing parental species than resynthesized allotetraploids, mirroring the establishment of the selfing syndrome. However, short-term changes mattered, and 40% of the cases of ELD in natural allotetraploids were already observed in resynthesized allotetraploids. Resynthesized allotetraploids showed striking variation of HEB among chromosomes and individuals. Homoeologous synapsis was its primary source and may still be a source of genetic variation in natural allotetraploids. In conclusion, both short- and long-term mechanisms contributed to transcriptomic and phenotypic changes in natural allotetraploids. However, the initial gene expression changes were largely reshaped during long-term evolution leading to further morphological changes.

Direct inference and control of genetic population structure from RNA sequencing data

RNAseq data can be used to infer genetic variants, yet its use for estimating genetic population structure remains underexplored. Here, we construct a freely available computational tool (RGStraP) to estimate RNAseq-based genetic principal components (RG-PCs) and assess whether RG-PCs can be used to control for population structure in gene expression analyses. Using whole blood samples from understudied Nepalese populations and the Geuvadis study, we show that RG-PCs had comparable results to paired array-based genotypes, with high genotype concordance and high correlations of genetic principal components, capturing subpopulations within the dataset. In differential gene expression analysis, we found that inclusion of RG-PCs as covariates reduced test statistic inflation. Our paper demonstrates that genetic population structure can be directly inferred and controlled for using RNAseq data, thus facilitating improved retrospective and future analyses of transcriptomic data.

Ancestral self-compatibility facilitates the establishment of allopolyploids in Brassicaceae

Self-incompatibility systems based on self-recognition evolved in hermaphroditic plants to maintain genetic variation of offspring and mitigate inbreeding depression. Despite these benefits in diploid plants, for polyploids who often face a scarcity of mating partners, self-incompatibility can thwart reproduction. In contrast, self-compatibility provides an immediate advantage: a route to reproductive viability. Thus, diploid selfing lineages may facilitate the formation of new allopolyploid species. Here, we describe the mechanism of establishment of at least four allopolyploid species in Brassicaceae (Arabidopsis suecica, Arabidopsis kamchatica, Capsella bursa-pastoris, and Brassica napus), in a manner dependent on the prior loss of the self-incompatibility mechanism in one of the ancestors. In each case, the degraded S-locus from one parental lineage was dominant over the functional S-locus of the outcrossing parental lineage. Such dominant loss-of-function mutations promote an immediate transition to selfing in allopolyploids and may facilitate their establishment.

Evolution of weedy giant ragweed (Ambrosia trifida): Multiple origins and gene expression variability facilitates weediness

Agricultural weeds may originate from wild populations, but the origination patterns and genetics underlying this transition remain largely unknown. Analysis of weedy-wild paired populations from independent locations may provide evidence to identify key genetic variation contributing to this adaptive shift. We performed genetic variation and expression analyses on transcriptome data from 67 giant ragweed samples collected from different locations in Ohio, Iowa, and Minnesota and found geographically separated weedy populations likely originated independently from their adjacent wild populations, but subsequent spreading of weedy populations also occurred locally. By using eight closely related weedy-wild paired populations, we identified thousands of unique transcripts in weedy populations that reflect shared or specific functions corresponding, respectively, to both convergently evolved and population-specific weediness processes. In addition, differential expression of specific groups of genes was detected between weedy and wild giant ragweed populations using gene expression diversity and gene co-expression network analyses. Our study suggests an integrated route of weedy giant ragweed origination, consisting of independent origination combined with the subsequent spreading of certain weedy populations, and provides several lines of evidence to support the hypothesis that gene expression variability plays a key role in the evolution of weedy species.

House mouse subspecies do differ in their social structure

It is widely acknowledged that population structure can have a substantial impact on evolutionary trajectories. In social animals, this structure is strongly influenced by relationships among the population members, so studies of differences in social structure between diverging populations or nascent species are of prime interest. Ideal models for such a study are two house mouse subspecies, Mus musculus musculus and M. m. domesticus, meeting in Europe along a secondary contact zone. Though the latter subspecies has usually been supposed to form tighter and more isolated social units than the former, the evidence is still inconclusive. Here, we carried out a series of radiofrequency identification experiments in semi-natural enclosures to gather large longitudinal data sets on individual mouse movements. The data were summarized in the form of uni- and multi-layer social networks. Within them, we could delimit and describe the social units ("modules"). While the number of estimated units was similar in both subspecies, domesticus revealed a more "modular" structure. This subspecies also showed more intramodular social interactions, higher spatial module separation, higher intramodular persistence of parent-offspring contacts, and lower multiple paternity, suggesting more effective control of dominant males over reproduction. We also demonstrate that long-lasting modules can be identified with basic reproductive units or demes. We thus provide the first robust evidence that the two subspecies differ in their social structure and dynamics of the structure formation.

The relative role of plasticity and demographic history in Capsella bursa-pastoris: a common garden experiment in Asia and Europe

The colonization success of a species depends on the interplay between its phenotypic plasticity, adaptive potential and demographic history. Assessing their relative contributions during the different phases of a species range expansion is challenging, and requires large-scale experiments. Here, we investigated the relative contributions of plasticity, performance and demographic history to the worldwide expansion of the shepherd's purse, Capsella bursa-pastoris. We installed two large common gardens of the shepherd's purse, a young, self-fertilizing, allopolyploid weed with a worldwide distribution. One common garden was located in Europe, the other in Asia. We used accessions from three distinct genetic clusters (Middle East, Europe and Asia) that reflect the demographic history of the species. Several life-history traits were measured. To explain the phenotypic variation between and within genetic clusters, we analysed the effects of (i) the genetic clusters, (ii) the phenotypic plasticity and its association to fitness and (iii) the distance in terms of bioclimatic variables between the sampling site of an accession and the common garden, i.e. the environmental distance. Our experiment showed that (i) the performance of C. bursa-pastoris is closely related to its high phenotypic plasticity; (ii) within a common garden, genetic cluster was a main determinant of phenotypic differences; and (iii) at the scale of the experiment, the effect of environmental distance to the common garden could not be distinguished from that of genetic clusters. Phenotypic plasticity and demographic history both play important role at different stages of range expansion. The success of the worldwide expansion of C. bursa-pastoris was undoubtedly influenced by its strong phenotypic plasticity.

Competitive ability depends on mating system and ploidy level across Capsella species

BACKGROUND AND AIMS

Self-fertilization is often associated with ecological traits corresponding to the ruderal strategy, and selfers are expected to be less competitive than outcrossers, either because of a colonization/competition trade-off or because of the deleterious genetic effects of selfing. Range expansion could reduce further competitive ability while polyploidy could mitigate the effects of selfing. If pollinators are not limited, individual fitness is thus expected to be higher in outcrossers than in selfers and, within selfers, in polyploids than in diploids. Although often proposed in the botanical literature and also suggested by meta-analyses, these predictions have not been directly tested yet.

METHODS

In order to compare fitness and the competitive ability of four Capsella species with a different mating system and ploidy level, we combined two complementary experiments. First, we carried out an experiment outdoors in north-west Greece, i.e. within the range of the obligate outcrossing species, C. grandiflora, where several life history traits were measured under two different disturbance treatments, weeded plots vs. unweeded plots. To better control competition and to remove potential effects of local adaptation of the outcrosser, we also performed a similar competition experiment but under growth chamber conditions.

KEY RESULTS

In the outdoor experiment, disturbance of the environment did not affect the phenotype in any of the four species. For most traits, the obligate outcrossing species performed better than all selfing species. In contrast, polyploids did not survive or reproduce better than diploids. Under controlled conditions, as in the field experiment, the outcrosser had a higher fitness than selfing species and was less affected by competition. Finally, contrary to the outdoor experiment where the two behaved identically, polyploid selfers were less affected by competition than diploid selfes.

CONCLUSIONS

In the Capsella genus, selfing induces lower fitness than outcrossing and can also reduce competitive ability. The effect of polyploidy is, however, unclear. These results highlight the possible roles of ecological context in the evolution of selfing species.

Optimizing bat bioacoustic surveys in human-modified Neotropical landscapes

During the last decades, the use of bioacoustics as a non-invasive and cost-effective sampling method has greatly increased worldwide. For bats, acoustic surveys have long been known to complement traditional mist-netting, however, appropriate protocol guidelines are still lacking for tropical regions. Establishing the minimum sampling effort needed to detect ecological changes in bat assemblages (e.g., activity, composition, and richness) is crucial in view of workload and project cost constraints, and because detecting such changes must be reliable enough to support effective conservation management. Using one of the most comprehensive tropical bat acoustic data sets, collected in the Amazon, we assessed the minimum survey effort required to accurately assess the completeness of assemblage inventories and habitat selection in fragmented forest landscapes for aerial insectivorous bats. We evaluated a combination of 20 different temporal sampling schemes, which differed regarding number of hours per night, number of nights per site, and sampling only during the wet or dry season, or both. This was assessed under two different landscape scenarios: in primary forest fragments embedded in a matrix of secondary forest and in the same forest fragments, but after they had been re-isolated through clearing of the secondary forest. We found that the sampling effort required to achieve 90% inventory completeness varied considerably depending on the research aim and the landscape scenario evaluated, averaging ~80 and 10 nights before and after fragment re-isolation, respectively. Recording for more than 4 h per night did not result in a substantial reduction in the required number of sampling nights. Regarding the effects of habitat selection, except for assemblage composition, bat responses in terms of richness, diversity, and activity were similar across all sampling schemes after fragment re-isolation. However, before re-isolation, a minimum of four to six sampling hours per night after dusk and three to five nights of sampling per site were needed to detect significant effects that could otherwise go unnoticed. Based on our results, we propose guidelines that will aid to optimize sampling protocols for bat acoustic surveys in the Neotropics.

Genomic Variation Landscape of the Model Salt Cress Eutrema salsugineum

Eutrema salsugineum has long been used as the model for examining salt and other abiotic stress in plants. In addition to the forward genetics approaches widely used in the lab, natural variations undoubtedly will provide a rich genetic resource for studying molecular mechanisms underlying the stress tolerance and local adaptation of this species. We used 90 resequencing whole genomes of natural populations of this species across its Asian and North American distributions to detect the selection signals for genes involved in salt and other stresses at the species-range level and local distribution. We detected selection signals for genes involved in salt and other abiotic tolerance at the species level. In addition, several cold-induced and defense genes showed selection signals due to local adaptation in North America-NE Russia or northern China, respectively. These variations and findings provide valuable resources for further deciphering genetic mechanisms underlying the stress tolerance and local adaptations of this model species.

Population genomics of selectively neutral genetic structure and herbicide resistance in UK populations of Alopecurus myosuroides

BACKGROUND

Alopecurus myosuroides (blackgrass) is a major weed in Europe with known resistance to multiple herbicide modes of action. In the UK, there is evidence that blackgrass has undergone a range expansion. In this paper, genotyping-by-sequencing and population-level herbicide resistance phenotypes are used to explore spatial patterns of selectively neutral genetic variation and resistance. We also perform a preliminary genome-wide association study (GWAS) and genomic prediction analysis to evaluate the potential of these approaches for investigating nontarget site herbicide resistance.

RESULTS

Blackgrass was collected from 47 fields across the British Isles and up to eight plants per field population (n = 369) were genotyped by Restriction site-associated DNA (RAD)-sequencing. A total of 20 426 polymorphic loci were identified and used for population genetic analyses. Phenotypic assays revealed significant variation in herbicide resistance between populations. Population structure was weak (F_ST = 0.024-0.048), but spatial patterns were consistent with an ongoing westward and northward range expansion. We detected strong and consistent Wahlund effects (F_IS = 0.30). There were no spatial patterns of herbicide resistance or evidence for confounding with population structure. Using a combination of population-level GWAS and genomic prediction we found that the top 20, 200, and 2000 GWAS loci had higher predictive abilities for fenoxaprop resistance compared to all markers.

CONCLUSION

There is likely extensive human-mediated gene flow between field populations of the weed blackgrass at a national scale. The lack of confounding of adaptive and neutral genetic variation can enable future, more extensive GWAS analyses to identify the genetic architecture of evolved herbicide resistance. © 2020 Society of Chemical Industry.

Transcriptional differentiation of UV-B protectant genes in maize landraces spanning an elevational gradient in Chiapas, Mexico

Globally, farmers cultivate and maintain crop landraces (i.e., traditional varieties). Landraces contain unique diversity shaped in part by natural and human-mediated selection and are an indispensable resource for farmers. Since environmental conditions change with elevation, crop landraces grown along elevational gradients have provided ideal locations to explore patterns of local adaptation. To further probe traits underlying this differentiation, transcriptome signatures can help provide a foundation for understanding the ways in which functional genetic diversity may be shaped by environment. In this study, we returned to an elevational gradient in Chiapas, Mexico, to assess transcriptional differentiation of genes underlying UV-B protection in locally adapted maize landraces from multiple elevations. We collected and planted landraces from three elevational zones (lowland, approximately 600 m; midland, approximately 1,550 m; highland approximately 2,100 m) in a common garden at 1,531 m. Using RNA-seq data derived from leaf tissue, we performed differential expression analysis between maize from these distinct elevations. Highland and lowland landraces displayed differential expression in phenylpropanoid and flavonoid biosynthesis genes involved in the production of UV-B protectants and did so at a rate greater than expected based on observed background transcriptional differentiation across the genome. These findings provide evidence for the differentiation of suites of genes involved in complex ecologically relevant pathways. Thus, while neutral evolutionary processes may have played a role in the observed patterns of differentiation, UV-B may have also acted as a selective pressure to differentiate maize landraces in the region. Studies of the distribution of functional crop genetic diversity across variable landscapes can aid us in understanding the response of diversity to abiotic/biotic change and, ultimately, may facilitate its conservation and utilization.

Shift in ecological strategy helps marginal populations of shepherd's purse (Capsella bursa-pastoris) to overcome a high genetic load

The outcome of species range expansion depends on the interplay of demographic, environmental and genetic factors. Self-fertilizing species usually show a higher invasive ability than outcrossers but selfing and bottlenecks during colonization also lead to an increased genetic load. The relationship between genomic and phenotypic characteristics of expanding populations has, hitherto, rarely been tested experimentally. We analysed how accessions of the shepherd's purse, Capsella bursa-pastoris, from the colonization front or from the core of the natural range performed under increasing density of competitors. First, accessions from the front showed a lower fitness than those from the core. Second, for all accessions, competitor density impacted negatively both vegetative growth and fruit production. However, despite their higher genetic load and lower absolute performances, accessions from the front were less affected by competition than accessions from the core. This seems to be due to phenotypic trade-offs and a shift in phenology that allow accessions from the front to avoid competition.

In the presence of population structure: From genomics to candidate genes underlying local adaptation

Understanding the genomic signatures, genes, and traits underlying local adaptation of organisms to heterogeneous environments is of central importance to the field evolutionary biology. To identify loci underlying local adaptation, models that combine allelic and environmental variation while controlling for the effects of population structure have emerged as the method of choice. Despite being evaluated in simulation studies, there has not been a thorough investigation of empirical evidence supporting local adaptation across these alleles. To evaluate these methods, we use 875 Arabidopsis thaliana Eurasian accessions and two mixed models (GEMMA and LFMM) to identify candidate SNPs underlying local adaptation to climate. Subsequently, to assess evidence of local adaptation and function among significant SNPs, we examine allele frequency differentiation and recent selection across Eurasian populations, in addition to their distribution along quantitative trait loci (QTL) explaining fitness variation between Italy and Sweden populations and cis-regulatory/nonsynonymous sites showing significant selective constraint. Our results indicate that significant LFMM/GEMMA SNPs show low allele frequency differentiation and linkage disequilibrium across locally adapted Italy and Sweden populations, in addition to a poor association with fitness QTL peaks (highest logarithm of odds score). Furthermore, when examining derived allele frequencies across the Eurasian range, we find that these SNPs are enriched in low-frequency variants that show very large climatic differentiation but low levels of linkage disequilibrium. These results suggest that their enrichment along putative functional sites most likely represents deleterious variation that is independent of local adaptation. Among all the genomic signatures examined, only SNPs showing high absolute allele frequency differentiation (AFD) and linkage disequilibrium (LD) between Italy and Sweden populations showed a strong association with fitness QTL peaks and were enriched along selectively constrained cis-regulatory/nonsynonymous sites. Using these SNPs, we find strong evidence linking flowering time, freezing tolerance, and the abscisic-acid pathway to local adaptation.

Towards the new normal: Transcriptomic convergence and genomic legacy of the two subgenomes of an allopolyploid weed (Capsella bursa-pastoris)

Allopolyploidy has played a major role in plant evolution but its impact on genome diversity and expression patterns remains to be understood. Some studies found important genomic and transcriptomic changes in allopolyploids, whereas others detected a strong parental legacy and more subtle changes. The allotetraploid C. bursa-pastoris originated around 100,000 years ago and one could expect the genetic polymorphism of the two subgenomes to follow similar trajectories and their transcriptomes to start functioning together. To test this hypothesis, we sequenced the genomes and the transcriptomes (three tissues) of allotetraploid C. bursa-pastoris and its parental species, the outcrossing C. grandiflora and the self-fertilizing C. orientalis. Comparison of the divergence in expression between subgenomes, on the one hand, and divergence in expression between the parental species, on the other hand, indicated a strong parental legacy with a majority of genes exhibiting a conserved pattern and cis-regulation. However, a large proportion of the genes that were differentially expressed between the two subgenomes, were also under trans-regulation reflecting the establishment of a new regulatory pattern. Parental dominance varied among tissues: expression in flowers was closer to that of C. orientalis and expression in root and leaf to that of C. grandiflora. Since deleterious mutations accumulated preferentially on the C. orientalis subgenome, the bias in expression towards C. orientalis observed in flowers indicates that expression changes could be adaptive and related to the selfing syndrome, while biases in the roots and leaves towards the C. grandiflora subgenome may be reflective of the differential genetic load.

Most plant species have a polyploid at some stage of their ancestry. Polyploidy, genome doubling through either multiple copies of a single species or through genomes of different species coming into the same nucleus, is therefore a crucial step in plant evolution. Understanding its impact on basic biological functions is thus a matter of interest. Shepherd’s purse (Capsella bursa-pastoris) is a major weed that appeared about 100,000 years ago through hybridization of two diploid species of the same genus. In the present project, we measured genetic diversity and analyzed gene expression patterns in flowers, roots, and leaves of C. bursa-pastoris individuals as well as in its two parental species, the outcrossing C. grandiflora and the self-fertilizing C. orientalis. Our data shows that, after 100,000 generations of evolution, the origin of the two subgenomes can still be seen: the genome inherited from C. grandiflora still differs from the one inherited from self-fertilizing C. orientalis. However, there are also signs that the two genomes have started to work together and are jointly regulated, and the way expression pattern varied across the three tissues indicates that the evolution of gene expression was adaptive.

Parental legacy, demography, and admixture influenced the evolution of the two subgenomes of the tetraploid Capsella bursa-pastoris (Brassicaceae)

Allopolyploidy is generally perceived as a major source of evolutionary novelties and as an instantaneous way to create isolation barriers. However, we do not have a clear understanding of how two subgenomes evolve and interact once they have fused in an allopolyploid species nor how isolated they are from their relatives. Here, we address these questions by analyzing genomic and transcriptomic data of allotetraploid Capsella bursa-pastoris in three differentiated populations, Asia, Europe, and the Middle East. We phased the two subgenomes, one descended from the outcrossing and highly diverse Capsella grandiflora (CbpCg) and the other one from the selfing and genetically depauperate Capsella orientalis (CbpCo). For each subgenome, we assessed its relationship with the diploid relatives, temporal changes of effective population size (Ne), signatures of positive and negative selection, and gene expression patterns. In all three regions, Ne of the two subgenomes decreased gradually over time and the CbpCo subgenome accumulated more deleterious changes than CbpCg. There were signs of widespread admixture between C. bursa-pastoris and its diploid relatives. The two subgenomes were impacted differentially depending on geographic region suggesting either strong interploidy gene flow or multiple origins of C. bursa-pastoris. Selective sweeps were more common on the CbpCg subgenome in Europe and the Middle East, and on the CbpCo subgenome in Asia. In contrast, differences in expression were limited with the CbpCg subgenome slightly more expressed than CbpCo in Europe and the Middle-East. In summary, after more than 100,000 generations of co-existence, the two subgenomes of C. bursa-pastoris still retained a strong signature of parental legacy but their evolutionary trajectory strongly varied across geographic regions.

Integrating Networks, Phylogenomics, and Population Genomics for the Study of Polyploidy

Duplication events are regarded as sources of evolutionary novelty, but our understanding of general trends for the long-term trajectory of additional genomic material is still lacking. Organisms with a history of whole genome duplication (WGD) offer a unique opportunity to study potential trends in the context of gene retention and/or loss, gene and network dosage, and changes in gene expression. In this review, we discuss the prevalence of polyploidy across the tree of life, followed by an overview of studies investigating genome evolution and gene expression. We then provide an overview of methods in network biology, phylogenomics, and population genomics that are critical for advancing our understanding of evolution post-WGD, highlighting the need for models that can accommodate polyploids. Finally, we close with a brief note on the importance of random processes in the evolution of polyploids with respect to neutral versus selective forces, ancestral polymorphisms, and the formation of autopolyploids versus allopolyploids.

Homeologue-specific expression divergence in the recently formed tetraploid Capsella bursa-pastoris (Brassicaceae)

Following allopolyploid formation, extensive genome evolution occurs, with the eventual loss of many homeologous gene copies. Although this process of diploidization has occurred many times independently, the evolutionary forces determining the probability and rate of gene loss remain poorly understood. Here, we conduct genome and transcriptome sequencing in a broad sample of Chinese accessions of Capsella bursa-pastoris, a recently formed allotetraploid. Our whole genome data reveal three groups of these accessions: an Eastern group from low-altitude regions, a Western group from high-altitude regions, and a much more differentiated Northwestern group. Population differentiation in total expression was limited among closely related populations; by contrast, the relative expression of the two homeologous copies closely mirrors the genome-wide SNP divergence. Consistent with this, we observe a negative correlation between expression changes in the two homeologues. However, genes showing population genomic evidence for adaptive evolution do not show an enrichment for expression divergence between homeologues, providing no clear evidence for adaptive shifts in relative gene expression. Overall, these patterns suggest that neutral drift may contribute to the population differentiation in the expression of the homeologues, and drive eventual gene loss over longer periods of time.

The role of ecotypic variation in driving worldwide colonization by a cosmopolitan plant

For almost 100 years now, ecotypic differentiation of plant species has been a major topic of research. In changing environments, the question needs to be answered as to how long it takes to adapt, and which parameters are subject to this fast adaptation. Short-living colonizing plant species are excellent examples, especially when they are selfing. Shepherd's Purse Capsella bursa-pastoris (Brassicaceae) is one of the most wide-spread flowering species on earth and avoids only the hot and humid tropics. Many studies demonstrated the ecotypic differentiation of C. bursa-pastoris in various regions of the world but ecotypic differentiation regarding adaptability of anatomy and physiology of rosette leaves so far remained less recognized. However, the leaves are relevant for subsequent seed set; in particular, winter-annual accessions require a robust rosette to survive adverse conditions. Leaf-related traits such as the thickness of the mesophyll and epidermis, stomatal density, photosynthetic capacity and the ability to withstand and even use high light conditions were therefore analysed in provenances from various climatic zones. Photosynthetic capacity depends on leaf anatomy and cellular physiological parameters. In particular, the ability to dynamically adjust the photosynthetic capacity to changing environmental conditions results in higher fitness. Here, we attempt to relate these results to the four Mendelian leaf types according to Shull.

Development and Evaluation of a Novel Set of EST-SSR Markers Based on Transcriptome Sequences of Black Locust (Robinia pseudoacacia L.)

Black locust (Robinia pseudoacacia L. of the family Fabaceae) is an ecologically and economically important deciduous tree. However, few genomic resources are available for this forest species, and few effective expressed sequence tag-derived simple sequence repeat (EST-SSR) markers have been developed to date. In this study, paired-end sequencing was used to sequence transcriptomes of R. pseudoacacia by the Illumina HiSeq TM2000 platform, and EST-SSR loci were identified by de novo assembly. Furthermore, a total of 1697 primer pairs were successfully designed, from which 286 primers met the selection screening criteria; 94 pairs were randomly selected and tested for validation using polymerase chain reaction amplification. Forty-five primers were verified as polymorphic, with clear bands. The polymorphism information content values were 0.033–0.765, the number of alleles per locus ranged from 2 to 10, and the observed and expected heterozygosities were 0.000–0.931 and 0.035–0.810, respectively, indicating a high level of informativeness. Subsequently, 45 polymorphic EST-SSR loci were tested for amplification efficiency, using the verified primers, in an additional nine species of Leguminosae, 23 loci were amplified in more than three species, of which two loci were amplified successfully in all species. These EST-SSR markers provide a valuable tool for investigating the genetic diversity and population structure of R. pseudoacacia, constructing a DNA fingerprint database, performing quantitative trait locus mapping, and preserving genetic information.

High-quality genome assembly of Capsella bursa-pastoris reveals asymmetry of regulatory elements at early stages of polyploid genome evolution

Polyploidization and subsequent sub- and neofunctionalization of duplicated genes represent a major mechanism of plant genome evolution. Capsella bursa-pastoris, a widespread ruderal plant, is a recent allotetraploid and, thus, is an ideal model organism for studying early changes following polyploidization. We constructed a high-quality assembly of C. bursa-pastoris genome and a transcriptome atlas covering a broad sample of organs and developmental stages (available online at http://travadb.org/browse/Species=Cbp). We demonstrate that expression of homeologs is mostly symmetric between subgenomes, and identify a set of homeolog pairs with discordant expression. Comparison of promoters within such pairs revealed emerging asymmetry of regulatory elements. Among them there are multiple binding sites for transcription factors controlling the regulation of photosynthesis and plant development by light (PIF3, HY5) and cold stress response (CBF). These results suggest that polyploidization in C. bursa-pastoris enhanced its plasticity of response to light and temperature, and allowed substantial expansion of its distribution range.

Genomic signature of successful colonization of Eurasia by the allopolyploid shepherd's purse (Capsella bursa-pastoris)

Polyploidization is a dominant feature of flowering plant evolution. However, detailed genomic analyses of the interpopulation diversification of polyploids following genome duplication are still in their infancy, mainly because of methodological limits, both in terms of sequencing and computational analyses. The shepherd's purse (Capsella bursa-pastoris) is one of the most common weed species in the world. It is highly self-fertilizing, and recent genomic data indicate that it is an allopolyploid, resulting from hybridization between the ancestors of the diploid species Capsella grandiflora and Capsella orientalis. Here, we investigated the genomic diversity of C. bursa-pastoris, its population structure and demographic history, following allopolyploidization in Eurasia. To that end, we genotyped 261 C. bursa-pastoris accessions spread across Europe, the Middle East and Asia, using genotyping-by-sequencing, leading to a total of 4274 SNPs after quality control. Bayesian clustering analyses revealed three distinct genetic clusters in Eurasia: one cluster grouping samples from Western Europe and Southeastern Siberia, the second one centred on Eastern Asia and the third one in the Middle East. Approximate Bayesian computation (ABC) supported the hypothesis that C. bursa-pastoris underwent a typical colonization history involving low gene flow among colonizing populations, likely starting from the Middle East towards Europe and followed by successive human-mediated expansions into Eastern Asia. Altogether, these findings bring new insights into the recent multistage colonization history of the allotetraploid C. bursa-pastoris and highlight ABC and genotyping-by-sequencing data as promising but still challenging tools to infer demographic histories of selfing allopolyploids.