Forces driving transposable element load variation during Arabidopsis range expansion.

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

Arabidopsis EXECUTER1 interacts with WRKY transcription factors to mediate plastid-to-nucleus singlet oxygen signaling.

Chloroplasts produce singlet oxygen (1O2), which causes changes in nuclear gene expression through plastid-to-nucleus retrograde signaling to increase plant fitness. However, the identity of this 1O2-triggered pathway remains unclear. Here, we identify mutations in GENOMES UNCOUPLED4 (GUN4) and GUN5 as suppressors of phytochrome-interacting factor1 (pif1) pif3 in regulating the photo-oxidative response in Arabidopsis thaliana. GUN4 and GUN5 specifically interact with EXECUTER1 (EX1) and EX2 in plastids, and this interaction is alleviated by treatment with Rose Bengal (RB) or white light. Impaired expression of GUN4, GUN5, EX1, or EX2 leads to insensitivity to excess light and overexpression of EX1 triggers photo-oxidative responses. Strikingly, upon light irradiation or RB treatment, EX1 transiently accumulates in the nucleus and the nuclear fraction of EX1 shows a similar molecular weight as the plastid-located protein. Point mutagenesis analysis indicated that nuclear localization of EX1 is required for its function. EX1 acts as a transcriptional co-activator and interacts with the transcription factors WRKY18 and WRKY40 to promote the expression of 1O2-responsive genes. This study suggests that EX1 may act in plastid-to-nucleus signaling and establishes a 1O2-triggered retrograde signaling pathway that allows plants adapt to changing light environments during chloroplast development.

Signatures of Adaptation and Purifying Selection in Highland Populations of Dasiphora fruticosa.

High mountains harbor a considerable proportion of biodiversity, but we know little about how diverse plants adapt to the harsh environment. Here we finished a high-quality genome assembly for Dasiphora fruticosa, an ecologically important plant distributed in the Qinghai-Tibetan Plateau and lowland of the Northern Hemisphere, and resequenced 592 natural individuals to address how this horticulture plant adapts to highland. Demographic analysis revealed D. fruticosa underwent a bottleneck after Naynayxungla Glaciation. Selective sweep analysis of two pairs of lowland and highland populations identified 63 shared genes related to cell wall organization or biogenesis, cellular component organization, and dwarfism, suggesting parallel adaptation to highland habitats. Most importantly, we found that stronger purging of estimated genetic load due to inbreeding in highland populations apparently contributed to their adaptation to the highest mountain. Our results revealed how plants could tolerate the extreme plateau, which could provide potential insights for species conservation and crop breeding.

Polygenic adaptation of rosette growth in Arabidopsis thaliana.

The rate at which plants grow is a major functional trait in plant ecology. However, little is known about its evolution in natural populations. Here, we investigate evolutionary and environmental factors shaping variation in the growth rate of Arabidopsis thaliana. We used plant diameter as a proxy to monitor plant growth over time in environments that mimicked latitudinal differences in the intensity of natural light radiation, across a set of 278 genotypes sampled within four broad regions, including an outgroup set of genotypes from China. A field experiment conducted under natural conditions confirmed the ecological relevance of the observed variation. All genotypes markedly expanded their rosette diameter when the light supply was decreased, demonstrating that environmental plasticity is a predominant source of variation to adapt plant size to prevailing light conditions. Yet, we detected significant levels of genetic variation both in growth rate and growth plasticity. Genome-wide association studies revealed that only 2 single nucleotide polymorphisms associate with genetic variation for growth above Bonferroni confidence levels. However, marginally associated variants were significantly enriched among genes with an annotated role in growth and stress reactions. Polygenic scores computed from marginally associated variants confirmed the polygenic basis of growth variation. For both light regimes, phenotypic divergence between the most distantly related population (China) and the various regions in Europe is smaller than the variation observed within Europe, indicating that the evolution of growth rate is likely to be constrained by stabilizing selection. We observed that Spanish genotypes, however, reach a significantly larger size than Northern European genotypes. Tests of adaptive divergence and analysis of the individual burden of deleterious mutations reveal that adaptive processes have played a more important role in shaping regional differences in rosette growth than maladaptive evolution.

Transcriptional repression specifies the central cell for double fertilization.

Double fertilization is a key innovation for the evolutionary success of angiosperms by which the two fertilized female gametes, the egg cell and central cell, generate the embryo and endosperm, respectively. The female gametophyte (embryo sac) enclosed in the sporophyte is derived from a one-celled haploid cell lineage. It undergoes successive events of mitotic divisions, cellularization, and cell specification to give rise to the mature embryo sac, which contains the two female gametes accompanied by two types of accessory cells, namely synergids and antipodals. How the cell fate of the central cell is specified has long been equivocal and is further complicated by the structural diversity of female gametophyte across plant taxa. Here, MADS-box protein AGL80 was verified as a transcriptional repressor that directly suppresses the expression of accessory cell-specific genes to specify the central cell. Further genetic rescue and phylogenetic assay of the AGL80 orthologs revealed a possible conserved mechanism in the Brassicaceae family. Results from this study provide insight into the molecular determination of the second female gamete cell in Brassicaceae.

Scaling-up and proteomic analysis reveals photosynthetic and metabolic insights toward prolonged H2 photoproduction in Chlamydomonas hpm91 mutant lacking proton gradient regulation 5 (PGR5).

Clean and sustainable H2 production is crucial to a carbon-neutral world. H2 generation by Chlamydomonas reinhardtii is an attractive approach for solar-H2 from H2O. However, it is currently not large-scalable because of lacking desirable strains with both optimal H2 productivity and sufficient knowledge of underlying molecular mechanism. We hereby carried out extensive and in-depth investigations of H2 photoproduction of hpm91 mutant lacking PGR5 (Proton Gradient Regulation 5) toward its up-scaling and fundamental mechanism issues. We show that hpm91 is at least 100-fold scalable (up to 10 L) with continuous H2 collection of 7287 ml H2/10L-HPBR in averagely 26 days under sulfur deprivation. Also, we show that hpm91 is robust and active during sustained H2 photoproduction, most likely due to decreased intracellular ROS relative to wild type. Moreover, we obtained quantitative proteomic profiles of wild type and hpm91 at four representing time points of H2 evolution, leading to 2229 and 1350 differentially expressed proteins, respectively. Compared to wild type, major proteome alterations of hpm91 include not only core subunits of photosystems and those related to anti-oxidative responses but also essential proteins in photosynthetic antenna, C/N metabolic balance, and sulfur assimilation toward both cysteine biosynthesis and sulfation of metabolites during sulfur-deprived H2 production. These results reveal not only new insights of cellular and molecular basis of enhanced H2 production in hpm91 but also provide additional candidate gene targets and modules for further genetic modifications and/or in artificial photosynthesis mimics toward basic and applied research aiming at advancing solar-H2 technology.

Adaptation and Phenotypic Diversification in Arabidopsis through Loss-of-Function Mutations in Protein-Coding Genes.

According to the less-is-more hypothesis, gene loss is an engine for evolutionary change. Loss-of-function (LoF) mutations resulting in the natural knockout of protein-coding genes not only provide information about gene function but also play important roles in adaptation and phenotypic diversification. Although the less-is-more hypothesis was proposed two decades ago, it remains to be explored on a large scale. In this study, we identified 60,819 LoF variants in 1071 Arabidopsis (Arabidopsis thaliana) genomes and found that 34% of Arabidopsis protein-coding genes annotated in the Columbia-0 genome do not have any LoF variants. We found that nucleotide diversity, transposable element density, and gene family size are strongly correlated with the presence of LoF variants. Intriguingly, 0.9% of LoF variants with minor allele frequency larger than 0.5% are associated with climate change. In addition, in the Yangtze River basin population, 1% of genes with LoF mutations were under positive selection, providing important insights into the contribution of LoF mutations to adaptation. In particular, our results demonstrate that LoF mutations shape diverse phenotypic traits. Overall, our results highlight the importance of the LoF variants for the adaptation and phenotypic diversification of plants.

Transposable elements drive rapid phenotypic variation in Capsella rubella.

Rapid phenotypic changes in traits of adaptive significance are crucial for organisms to thrive in changing environments. How such phenotypic variation is achieved rapidly, despite limited genetic variation in species that experience a genetic bottleneck is unknown. Capsella rubella, an annual and inbreeding forb (Brassicaceae), is a great system for studying this basic question. Its distribution is wider than those of its congeneric species, despite an extreme genetic bottleneck event that severely diminished its genetic variation. Here, we demonstrate that transposable elements (TEs) are an important source of genetic variation that could account for its high phenotypic diversity. TEs are (i) highly enriched in C. rubella compared with its outcrossing sister species Capsella grandiflora, and (ii) 4.2% of polymorphic TEs in C. rubella are associated with variation in the expression levels of their adjacent genes. Furthermore, we show that frequent TE insertions at FLOWERING LOCUS C (FLC) in natural populations of C. rubella could explain 12.5% of the natural variation in flowering time, a key life history trait correlated with fitness and adaptation. In particular, we show that a recent TE insertion at the 3' UTR of FLC affects mRNA stability, which results in reducing its steady-state expression levels, to promote the onset of flowering. Our results highlight that TE insertions can drive rapid phenotypic variation, which could potentially help with adaptation to changing environments in a species with limited standing genetic variation.

Natural variation in the HAN1 gene confers chilling tolerance in rice and allowed adaptation to a temperate climate.

Rice (Oryza sativa L.) is a chilling-sensitive staple crop that originated in subtropical regions of Asia. Introduction of the chilling tolerance trait enables the expansion of rice cultivation to temperate regions. Here we report the cloning and characterization of HAN1, a quantitative trait locus (QTL) that confers chilling tolerance on temperate japonica rice. HAN1 encodes an oxidase that catalyzes the conversion of biologically active jasmonoyl-L-isoleucine (JA-Ile) to the inactive form 12-hydroxy-JA-Ile (12OH-JA-Ile) and fine-tunes the JA-mediated chilling response. Natural variants in HAN1 diverged between indica and japonica rice during domestication. A specific allele from temperate japonica rice, which gained a putative MYB cis-element in the promoter of HAN1 during the divergence of the two japonica ecotypes, enhances the chilling tolerance of temperate japonica rice and allows it to adapt to a temperate climate. The results of this study extend our understanding of the northward expansion of rice cultivation and provide a target gene for the improvement of chilling tolerance in rice.

Identification of long noncoding natural antisense transcripts (lncNATs) correlated with drought stress response in wild rice (Oryza nivara).

BACKGROUND: Wild rice, including Oryza nivara and Oryza rufipogon, which are considered as the ancestors of Asian cultivated rice (Oryza sativa), possess high genetic diversity and serve as a crucial resource for breeding novel cultivars of cultivated rice. Although rice domestication related traits, such as seed shattering and plant architecture, have been intensively studied at the phenotypic and genomic levels, further investigation is needed to understand the molecular basis of phenotypic differences between cultivated and wild rice. Drought stress is one of the most severe abiotic stresses affecting rice growth and production. Adaptation to drought stress involves a cascade of genes and regulatory factors that form complex networks. O. nivara inhabits swampy areas with a seasonally dry climate, which is an ideal material to discover drought tolerance alleles. Long noncoding natural antisense transcripts (lncNATs), a class of long noncoding RNAs (lncRNAs), regulate the corresponding sense transcripts and play an important role in plant growth and development. However, the contribution of lncNATs to drought stress response in wild rice remains largely unknown. RESULTS: Here, we conducted strand-specific RNA sequencing (ssRNA-seq) analysis of Nipponbare (O. sativa) and two O. nivara accessions (BJ89 and BJ278) to determine the role of lncNATs in drought stress response in wild rice. A total of 1246 lncRNAs were identified, including 1091 coding-noncoding NAT pairs, of which 50 were expressed only in Nipponbare, and 77 were expressed only in BJ89 and/or BJ278. Of the 1091 coding-noncoding NAT pairs, 240 were differentially expressed between control and drought stress conditions. Among these 240 NAT pairs, 12 were detected only in Nipponbare, and 187 were detected uniquely in O. nivara. Furthermore, 10 of the 240 coding-noncoding NAT pairs were correlated with genes enriched in stress responsive GO terms; among these, nine pairs were uniquely found in O. nivara, and one pair was shared between O. nivara and Nipponbare. CONCLUSION: We identified lncNATs associated with drought stress response in cultivated rice and O. nivara. These results will improve our understanding of the function of lncNATs in drought tolerance and accelerate rice breeding.

Parallel Evolution of Common Allelic Variants Confers Flowering Diversity in Capsella rubella.

Flowering time is an adaptive life history trait. Capsella rubella, a close relative of Arabidopsis thaliana and a young species, displays extensive variation for flowering time but low standing genetic variation due to an extreme bottleneck event, providing an excellent opportunity to understand how phenotypic diversity can occur with a limited initial gene pool. Here, we demonstrate that common allelic variation and parallel evolution at the FLC locus confer variation in flowering time in C. rubella. We show that two overlapping deletions in the 5' untranslated region (UTR) of C. rubella FLC, which are associated with local changes in chromatin conformation and histone modifications, reduce its expression levels and promote flowering. We further show that these two pervasive variants originated independently in natural C. rubella populations after speciation and spread to an intermediate frequency, suggesting a role of this parallel cis-regulatory change in adaptive evolution. Our results provide an example of how parallel mutations in the same 5' UTR region can shape phenotypic evolution in plants.

Parallel Speciation of Wild Rice Associated with Habitat Shifts.

The occurrence of parallel speciation strongly implies the action of natural selection. However, it is unclear how general a phenomena parallel speciation is since it was only shown in a small number of animal species. In particular, the adaptive process and mechanisms underlying the process of parallel speciation remain elusive. Here, we used an integrative approach incorporating population genomics, common garden, and crossing experiments to investigate parallel speciation of the wild rice species Oryza nivara from O. rufipogon. We demonstrated that O. nivara originated multiple times from different O. rufipogon populations and revealed that different O. nivara populations have evolved similar phenotypes under divergent selection, a reflection of recurrent local adaptation of ancient O. rufipogon populations to dry habitats. Almost completed premating isolation was detected between O. nivara and O. rufipogon in the absence of any postmating barriers between and within these species. These results suggest that flowering time is a "magic" trait that contributes to both local adaptation and reproductive isolation in the origin of wild rice species. Our study thus demonstrates a convincing case of parallel ecological speciation as a consequence of adaptation to new environments.

Rice FLUORESCENT1 Is Involved in the Regulation of Chlorophyll.

Chlorophyll biosynthesis plays essential roles in photosynthesis and plant growth in response to environmental conditions. The accumulation of excess chlorophyll biosynthesis intermediates under light results in the production of reactive oxygen species and oxidative stress. In this study, we identified a rice (Oryza sativa) mutant, oxidation under photoperiod (oxp), that displayed photobleached lesions on its leaves, reduced growth and decreased chlorophyll content during light/dark cycles or following a dark-to-light transition. The oxp mutant accumulated more chlorophyll precursors (5-aminolevulinic acid and protochlorophyllide) than the wild type in the dark, and more singlet oxygen following light exposure. Several singlet-oxygen-responsive genes were greatly upregulated in oxp, whereas the expression patterns of OsPORA and OsPORB, two genes encoding the chlorophyll biosynthesis enzyme NADPH:protochlorop hyllide oxidoreductase, were altered in de-etiolated oxp seedlings. Molecular and complementation studies revealed that oxp is a loss-of-function mutant in LOC_Os01g32730, a homolog of FLUORESCENT (FLU) in Arabidopsis thaliana. Rice PHYTOCHROME-INTERACTING FACTOR-LIKE14 (OsPIL14) transcription factor directly bound to the OsFLU1 promoter and activated its expression. Dark-grown transgenic rice seedlings overexpressing OsPIL14 accumulated more chlorophyll and turned green faster than the wild type upon light illumination. Thus, OsFLU1 is an important regulator of chlorophyll biosynthesis in rice.

Gene family evolution in green plants with emphasis on the origination and evolution of Arabidopsis thaliana genes.

Gene family size variation is an important mechanism that shapes the natural variation for adaptation in various species. Despite its importance, the pattern of gene family size variation in green plants is still not well understood. In particular, the evolutionary pattern of genes and gene families remains unknown in the model plant Arabidopsis thaliana in the context of green plants. In this study, eight representative genomes of green plants are sampled to study gene family evolution and characterize the origination of A. thaliana genes, respectively. Four important insights gained are that: (i) the rate of gene gains and losses is about 0.001359 per gene every million years, similar to the rate in yeast, Drosophila, and mammals; (ii) some gene families evolved rapidly with extreme expansions or contractions, and 2745 gene families present in all the eight species represent the 'core' proteome of green plants; (iii) 70% of A. thaliana genes could be traced back to 450 million years ago; and (iv) intriguingly, A. thaliana genes with early origination are under stronger purifying selection and more conserved. In summary, the present study provides genome-wide insights into evolutionary history and mechanisms of genes and gene families in green plants and especially in A. thaliana.

Transposable elements and small RNAs contribute to gene expression divergence between Arabidopsis thaliana and Arabidopsis lyrata.

Transposable elements (TEs) are often the primary determinant of genome size differences among eukaryotes. In plants, the proliferation of TEs is countered through epigenetic silencing mechanisms that prevent mobility. Recent studies using the model plant Arabidopsis thaliana have revealed that methylated TE insertions are often associated with reduced expression of nearby genes, and these insertions may be subject to purifying selection due to this effect. Less is known about the genome-wide patterns of epigenetic silencing of TEs in other plant species. Here, we compare the 24-nt siRNA complement from A. thaliana and a closely related congener with a two- to threefold higher TE copy number, Arabidopsis lyrata. We show that TEs--particularly siRNA-targeted TEs--are associated with reduced gene expression within both species and also with gene expression differences between orthologs. In addition, A. lyrata TEs are targeted by a lower fraction of uniquely matching siRNAs, which are associated with more effective silencing of TE expression. Our results suggest that the efficacy of RNA-directed DNA methylation silencing is lower in A. lyrata, a finding that may shed light on the causes of differential TE proliferation among species.

[Improved unilateral puncture PVP based on 3D printing technology for the treatment of osteoporotic vertebral compression fracture].

OBJECTIVE: To investigate the clinical effect of unilateral percutaneous vertebroplasty (PVP) combined with 3D printing technology for the treatment of thoracolumbar osteoporotic compression fracture. METHODS: A total of 77 patients with thoracolumbar osteoporotic compression fractures from October 2020 to April 2022 were included in the study, all of which were vertebral body compression fractures caused by trauma. According to different treatment methods, they were divided into experimental group and control group. Thirty-two patients used 3D printing technology to improve unilateral transpedicle puncture vertebroplasty in the experimental group, there were 5 males and 27 females, aged from 63 to 91 years old with an average of (77.59±8.75) years old. Forty-five patients were treated with traditional bilateral pedicle puncture vertebroplasty, including 7 males and 38 females, aged from 60 to 88 years old with an average of(74.89±7.37) years old. Operation time, intraoperative C-arm X-ray times, anesthetic dosage, bone cement injection amount, bone cement diffusion good and good rate, complications, vertebral height, kyphotic angle (Cobb angle), visual analogue scale(VAS), Oswestry disability index (ODI) and other indicators were recorded before and after surgery, and statistically analyzed. RESULTS: All patients were followed up for 6 to 23 months, with preoperative imaging studies, confirmed for thoracolumbar osteoporosis compression fractures, two groups of patients with postoperative complications, no special two groups of patients' age, gender, body mass index (BMI), time were injured, the injured vertebral distribution had no statistical difference(P>0.05), comparable data. Two groups of patients with bone cement injection, bone cement dispersion rate, preoperative and postoperative vertebral body height, protruding after spine angle(Cobb angle), VAS, ODI had no statistical difference(P>0.05). The operative time, intraoperative fluoroscopy times and anesthetic dosage were statistically different between the two groups(P<0.05). Compared with the traditional bilateral puncture group, the modified unilateral puncture group combined with 3D printing technology had shorter operation time, fewer intraoperative fluoroscopy times and less anesthetic dosage. The height of anterior vertebral edge, kyphosis angle (Cobb angle), VAS score and ODI of the affected vertebrae were statistically different between two groups at each time point after surgery(P<0.05). CONCLUSION: In the treatment of thoracolumbar osteoporotic compression fractures, 3D printing technology is used to improve unilateral puncture PVP, which is convenient and simple, less trauma, short operation time, fewer fluoroscopy times, satisfactory distribution of bone cement, vertebral height recovery and kyphotic Angle correction, and good functional improvement.

The plant immune receptor SNC1 monitors helper NLRs targeted by a bacterial effector.

Plants deploy intracellular receptors to counteract pathogen effectors that suppress cell-surface-receptor-mediated immunity. To what extent pathogens manipulate intracellular receptor-mediated immunity, and how plants tackle such manipulation, remains unknown. Arabidopsis thaliana encodes three similar ADR1 class helper nucleotide-binding domain leucine-rich repeat receptors (ADR1, ADR1-L1, and ADR1-L2), which are crucial in plant immunity initiated by intracellular receptors. Here, we report that Pseudomonas syringae effector AvrPtoB suppresses ADR1-L1- and ADR1-L2-mediated cell death. ADR1, however, evades such suppression by diversifying into two ubiquitination sites targeted by AvrPtoB. The intracellular sensor SNC1 interacts with and guards the CCR domains of ADR1-L1/L2. Removal of ADR1-L1/L2 or delivery of AvrPtoB activates SNC1, which then signals through ADR1 to trigger immunity. Our work elucidates the long-sought-after function of SNC1 in defense, and also how plants can use dual strategies, sequence diversification, and a multi-layered guard-guardee system, to counteract pathogen's attack on core immunity functions.

Evolution of the S-locus region in Arabidopsis relatives.

The S locus, a single polymorphic locus, is responsible for self-incompatibility (SI) in the Brassicaceae family and many related plant families. Despite its importance, our knowledge of S-locus evolution is largely restricted to the causal genes encoding the S-locus receptor kinase (SRK) receptor and S-locus cysteine-rich protein (SCR) ligand of the SI system. Here, we present high-quality sequences of the genomic region of six S-locus haplotypes: Arabidopsis (Arabidopsis thaliana; one haplotype), Arabidopsis lyrata (four haplotypes), and Capsella rubella (one haplotype). We compared these with reference S-locus haplotypes of the self-compatible Arabidopsis and its SI congener A. lyrata. We subsequently reconstructed the likely genomic organization of the S locus in the most recent common ancestor of Arabidopsis and Capsella. As previously reported, the two SI-determining genes, SCR and SRK, showed a pattern of coevolution. In addition, consistent with previous studies, we found that duplication, gene conversion, and positive selection have been important factors in the evolution of these two genes and appear to contribute to the generation of new recognition specificities. Intriguingly, the inactive pseudo-S-locus haplotype in the self-compatible species C. rubella is likely to be an old S-locus haplotype that only very recently became fixed when C. rubella split off from its SI ancestor, Capsella grandiflora.

Genome-wide comparison of nucleotide-binding site-leucine-rich repeat-encoding genes in Arabidopsis.

Plants, like animals, use several lines of defense against pathogen attack. Prominent among genes that confer disease resistance are those encoding nucleotide-binding site-leucine-rich repeat (NB-LRR) proteins. Likely due to selection pressures caused by pathogens, NB-LRR genes are the most variable gene family in plants, but there appear to be species-specific limits to the number of NB-LRR genes in a genome. Allelic diversity within an individual is also increased by obligatory outcrossing, which leads to genome-wide heterozygosity. In this study, we compared the NB-LRR gene complement of the selfer Arabidopsis thaliana and its outcrossing close relative Arabidopsis lyrata. We then complemented and contrasted the interspecific patterns with studies of NB-LRR diversity within A. thaliana. Three important insights are as follows: (1) that both species have similar numbers of NB-LRR genes; (2) that loci with single NB-LRR genes are less variable than tandem arrays; and (3) that presence-absence polymorphisms within A. thaliana are not strongly correlated with the presence or absence of orthologs in A. lyrata. Although A. thaliana individuals are mostly homozygous and thus potentially less likely to suffer from aberrant interaction of NB-LRR proteins with newly introduced alleles, the number of NB-LRR genes is similar to that in A. lyrata. In intraspecific and interspecific comparisons, NB-LRR genes are also more variable than receptor-like protein genes. Finally, in contrast to Drosophila, there is a clearly positive relationship between interspecific divergence and intraspecific polymorphisms.

Recent speciation of Capsella rubella from Capsella grandiflora, associated with loss of self-incompatibility and an extreme bottleneck.

Flowering plants often prevent selfing through mechanisms of self-incompatibility (S.I.). The loss of S.I. has occurred many times independently, because it provides short-term advantages in situations where pollinators or mates are rare. The genus Capsella, which is closely related to Arabidopsis, contains a pair of closely related diploid species, the self-incompatible Capsella grandiflora and the self-compatible Capsella rubella. To elucidate the transition to selfing and its relationship to speciation of C. rubella, we have made use of comparative sequence information. Our analyses indicate that C. rubella separated from C. grandiflora recently ( approximately 30,000-50,000 years ago) and that breakdown of S.I. occurred at approximately the same time. Contrasting the nucleotide diversity patterns of the 2 species, we found that C. rubella has only 1 or 2 alleles at most loci, suggesting that it originated through an extreme population bottleneck. Our data are consistent with diploid speciation by a single, selfing individual, most likely living in Greece. The new species subsequently colonized the Mediterranean by Northern and Southern routes, at a time that also saw the spread of agriculture. The presence of phenotypic diversity within modern C. rubella suggests that this species will be an interesting model to understand divergence and adaptation, starting from very limited standing genetic variation.