Search

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Background: Whole genome duplication (WGD) events are common in the evolutionary history of many living organisms. For decades, researchers have been trying to understand the genetic and epigenetic impact of WGD and its underlying molecular mechanisms. Particular attention was given to allopolyploid study systems, species resulting from an hybridization event accompanied by WGD. Investigating the mechanisms behind the survival of a newly formed allopolyploid highlighted the key role of DNA methylation. With the improvement of high-throughput methods, such as whole genome bisulfite sequencing (WGBS), an opportunity opened to further understand the role of DNA methylation at a larger scale and higher resolution. However, only a few studies have applied WGBS to allopolyploids, which might be due to lack of genomic resources combined with a burdensome data analysis process. To overcome these problems, we developed the Automated Reproducible Polyploid EpiGenetic GuIdance workflOw (ARPEGGIO): the first workflow for the analysis of epigenetic data in polyploids. This workflow analyzes WGBS data from allopolyploid species via the genome assemblies of the allopolyploid’s parent species. ARPEGGIO utilizes an updated read classification algorithm (EAGLE-RC), to tackle the challenge of sequence similarity amongst parental genomes. ARPEGGIO offers automation, but more importantly, a complete set of analyses including spot checks starting from raw WGBS data: quality checks, trimming, alignment, methylation extraction, statistical analyses and downstream analyses. A full run of ARPEGGIO outputs a list of genes showing differential methylation. ARPEGGIO was made simple to set up, run and interpret, and its implementation ensures reproducibility by including both package management and containerization. Results: We evaluated ARPEGGIO in two ways. First, we tested EAGLE-RC’s performance with publicly available datasets given a ground truth, and we show that EAGLE-RC decreases the error rate by 3 to 4 times compared to standard approaches. Second, using the same initial dataset, we show agreement between ARPEGGIO’s output and published results. Compared to other similar workflows, ARPEGGIO is the only one supporting polyploid data. Conclusions: The goal of ARPEGGIO is to promote, support and improve polyploid research with a reproducible and automated set of analyses in a convenient implementation. ARPEGGIO is available at https://github.com/supermaxiste/ARPEGGIO., MC genomics, 22 (1), ISSN:1471-2164

Complex cell shapes are generated first by breaking symmetry, and subsequent polar growth. Localized bending of anticlinal walls initiates lobe formation in the epidermal pavement cells of cotyledons and leaves, but how the microtubule cytoskeleton mediates local cell growth, and how plant pavement cells benefit from adopting jigsaw puzzle-like shapes, are poorly understood. In Arabidopsis (Arabidopsis thaliana), the basic Pro-rich protein (BPP) microtubule-associated protein family comprises seven members. We analyzed lobe morphogenesis in cotyledon pavement cells of a BPP1;BPP2;BPP5 triple knockout mutant. New image analysis methods (MtCurv and BQuant) showed that anticlinal microtubule bundles were significantly reduced and cortical microtubules that fan out radially across the periclinal wall did not enrich at the convex side of developing lobes. Despite these microtubule defects, new lobes were initiated at the same frequency as in wild-type cells, but they did not expand into well-defined protrusions. Eventually, mutant cells formed nearly polygonal shapes and adopted concentric microtubule patterns. The mutant periclinal cell wall bulged outward. The radius of the calculated inscribed circle of the pavement cells, a proposed proxy for maximal stress in the cell wall, was consistently larger in the mutant cells during cotyledon development, and correlated with an increase in cell height. These bpp mutant phenotypes provide genetic and cell biological evidence that initiation and growth of lobes are distinct morphogenetic processes, and that interdigitated cell geometry effectively suppresses large outward bulging of pavement cells.

Contemporary speciation provides a unique opportunity to directly observe the traits and environmental responses of a new species. Cardamine insueta is an allotriploid species that appeared within the past 150 years in a Swiss village, Urnerboden. In contrast to its two progenitor species, C. amara and C. rivularis that live in wet and open habitats, respectively, C. insueta is found in-between their habitats with temporal water level fluctuation. This triploid species propagates clonally and serves as a triploid bridge to form higher ploidy species. Although niche separation is observed in field studies, the mechanisms underlying the environmental robustness of C. insueta are not clear. To characterize responses to a fluctuating environment, we performed a time-course analysis of homeolog gene expression in C. insueta in response to submergence treatment. For this purpose, the two parental (C. amara and C. rivularis) genome sequences were assembled with a reference-guided approach, and homeolog-specific gene expression was quantified by using HomeoRoq software. We found that C. insueta and C. rivularis initiated vegetative propagation by forming ectopic meristems on leaves, while C. amara did not. We examined homeolog-specific gene expression of three species at nine time points during the treatment. The genome-wide expression ratio of homeolog pairs was 2:1 over the time-course, consistent with the ploidy number. By searching the genes with high coefficient of variation of expression over time-course transcriptome data, we found many known key transcriptional factors related to meristem development and formation upregulated in both C. rivularis and rivularis-homeolog of C. insueta, but not in C. amara. Moreover, some amara-homeologs of these genes were also upregulated in the triploid, suggesting trans-regulation. In turn, Gene Ontology analysis suggested that the expression pattern of submergence tolerant genes in the triploid was inherited from C. amara. These results suggest that the triploid C. insueta combined advantageous patterns of parental transcriptomes to contribute to its establishment in a new niche along a water-usage gradient.

Population genetics studies of allopolyploid species lag behind those of diploid species because of practical difficulties in analysis of homeologs-duplicated gene copies originating from hybridized parental species. Pool-Seq, i.e. massive parallel sequencing of pooled individuals, has high potential for detecting nucleotide polymorphisms within and among multiple populations; however, its use has been limited to diploid species. We applied Pool-Seq to an allopolyploid species by developing a bioinformatic pipeline that assigns reads to each homeolog as well as to each polymorphic allele within each homeolog. We simultaneously sequenced eight genes from twenty individuals from each of 24 populations, and found over 100 polymorphic sites in each homeolog. For two sites, we estimated allele frequencies using the number of reads and then validated these estimations by making individual-based estimations. Pool-Seq using our bioinformatic pipeline allows efficient evaluation of nucleotide polymorphisms in a large number of individuals, even in allopolyploid species.

Mouse APOBEC3 (mA3) inhibits murine leukemia virus (MuLV) replication by a deamination-independent mechanism in which the reverse transcription is considered the main target process. However, other steps in virus replication that can be targeted by mA3 have not been examined. We have investigated the possible effect of mA3 on MuLV protease-mediated processes and found that mA3 binds both mature viral protease and Pr180gag-pol precursor polyprotein. Using replication-competent MuLVs, we also show that mA3 inhibits the processing of Pr65 Gag precursor. Furthermore, we demonstrate that the autoprocessing of Pr180gag-pol is impeded by mA3, resulting in reduced production of mature viral protease. This reduction appears to link with the above inefficient Pr65gag processing in the presence of mA3. Two major isoforms of mA3, exon 5-containing and -lacking ones, equally exhibit this antiviral activity. Importantly, physiologically expressed levels of mA3 impedes both Pr180gag-pol autocatalysis and Pr65gag processing. This blockade is independent of the deaminase activity and requires the C-terminal region of mA3. These results suggest that the above impairment of Pr180gag-pol autoprocessing may significantly contribute to the deaminase-independent antiretroviral activity exerted by mA3., Author summary Soon after the identification of the polynucleotide cytidine deaminase APOBEC3 as a host restriction factor against vif-deficient HIV, it was noticed that deamination-independent mechanisms are involved in the inhibition of viral replication in addition to the deaminase-dependent mechanism. We previously showed that mouse APOBEC3 (mA3) physiologically restricted mouse retrovirus replication in their natural hosts without causing significant G-to-A hypermutations. Inhibition of reverse transcription is reported to be the most plausible mechanism for the deamination-independent antiretroviral function. However, it remains unknown whether the inhibition of reverse transcription is the only way to explain the whole picture of deamination-independent antiviral activity exerted by APOBEC3. Here we show that mA3 targets the autoprocessing of Pr180gag-pol polyprotein. This activity does not require the deaminase catalytic center and mainly exerted by the C-terminal half of mA3. mA3 physically interacts with murine retroviral protease and its precursor Pr180gag-pol. mA3-induced disruption of the autocatalytic Pr180gag-pol cleavage leads to a significant reduction of mature viral protease, resulting in the inhibition of Pr65gag processing to mature Gag proteins. As the Pr180gag-pol autoprocessing is necessary for the maturation of other viral enzymes including the reverse transcriptase, its inhibition by host APOBEC3 may precede the previously described impairment of reverse transcription. Our discovery may lead to the development of novel antiretroviral drugs through the future identification of detailed molecular interfaces between retroviral Gag-Pol polyprotein and APOBEC3.

The self-incompatible species Arabidopsis halleri is a close relative of the self-compatible model plant Arabidopsis thaliana. The broad European and Asian distribution and heavy metal hyperaccumulation ability make A. halleri a useful model for ecological genomics studies. We used long-insert mate-pair libraries to improve the genome assembly of the A. halleri ssp. gemmifera Tada mine genotype (W302) collected from a site with high contamination by heavy metals in Japan. After five rounds of forced selfing, heterozygosity was reduced to 0.04%, which facilitated subsequent genome assembly. Our assembly now covers 196 Mb or 78% of the estimated genome size and achieved scaffold N50 length of 712 kb. To validate assembly and annotation, we used synteny of A. halleri Tada mine with a previously published high-quality reference assembly of a closely related species, Arabidopsis lyrata. Further validation of the assembly quality comes from synteny and phylogenetic analysis of the HEAVY METAL ATPASE4 (HMA4) and METAL TOLERANCE PROTEIN1 (MTP1) regions using published sequences from European A. halleri for comparison. Three tandemly duplicated copies of HMA4, key gene involved in cadmium and zinc hyperaccumulation, were assembled on a single scaffold. The assembly will enhance the genomewide studies of A. halleri as well as the allopolyploid Arabidopsis kamchatica derived from A. lyrata and A. halleri.

The habitats of polyploid species are generally distinct from their parental species. Stebbins described polyploids as 'general purpose genotypes', which can tolerate a wide range of environmental conditions. However, little is known about its molecular basis because of the complexity of polyploid genomes. We hypothesized that allopolyploid species might utilize the expression patterns of both parents depending on environments (polyploid plasticity hypothesis). We focused on hydrological niche segregation along fine-scale soil moisture and waterlogging gradients. Two diploid species, Cardamine amara and Cardamine hirsuta, grew best in submerged and unsubmerged conditions, respectively, consistent with their natural habitats. Interestingly, the allotetraploid Cardamine flexuosa derived from them grew similarly in fluctuating as well as submerged and unsubmerged conditions, consistent with its wide environmental tolerance. A similar pattern was found in another species trio: allotetraploid Cardamine scutata and its parents. Using the close relatedness of Cardamine and Arabidopsis, we quantified genomewide expression patterns following dry and wet treatments using an Arabidopsis microarray. Hierarchical clustering analysis revealed that the expression pattern of C. flexuosa clustered with C. hirsuta in the dry condition and with C. amara in the wet condition, supporting our hypothesis. Furthermore, the induction levels of most genes in the allopolyploid were lower than in a specialist diploid species. This reflects a disadvantage of being allopolyploid arising from fixed heterozygosity. We propose that recurrent allopolyploid speciation along soil moisture and waterlogging gradients confers niche differentiation and reproductive isolation simultaneously and serves as a model for studying the molecular basis of ecological speciation and adaptive radiation.

Finger millet (Eleusine coracana (L.) Gaertn) is an important crop for food security because of its tolerance to drought, which is expected to be exacerbated by global climate changes. Nevertheless, it is often classified as an orphan/underutilized crop because of the paucity of scientific attention. Among several small millets, finger millet is considered as an excellent source of essential nutrient elements, such as iron and zinc; hence, it has potential as an alternate coarse cereal. However, high-quality genome sequence data of finger millet are currently not available. One of the major problems encountered in the genome assembly of this species was its polyploidy, which hampers genome assembly compared with a diploid genome. To overcome this problem, we sequenced its genome using diverse technologies with sufficient coverage and assembled it via a novel multiple hybrid assembly workflow that combines next-generation with single-molecule sequencing, followed by whole-genome optical mapping using the Bionano Irys® system. The total number of scaffolds was 1,897 with an N50 length >2.6 Mb and detection of 96% of the universal single-copy orthologs. The majority of the homeologs were assembled separately. This indicates that the proposed workflow is applicable to the assembly of other allotetraploid genomes., DNA Research, 25 (1), ISSN:1340-2838, ISSN:1756-1663

Background and Aims The coexistence of hermaphrodites and female-sterile individuals, or androdioecy, has been documented in only a handful of plants and animals. This study reports its existence in the plant species Cardamine amara (Brassicaceae), in which female-sterile individuals have shorter pistils than seed-producing hermaphrodites. Methods Morphological analysis, in situ manual pollination, microsatellite genotyping and differential gene expression analysis using Arabidopsis microarrays were used to delimit variation between female-sterile individuals and hermaphrodites. Key Results Female sterility in C. amara appears to be caused by disrupted ovule development. It was associated with a 2.4- to 2.9-fold increase in clonal propagation. This made the pollen number of female-sterile genets more than double that of hermaphrodite genets, which fulfils a condition of co-existence predicted by simple androdioecy theories. When female-sterile individuals were observed in wild androdioecious populations, their ramet frequencies ranged from 5 to 54 %; however, their genet frequencies ranged from 11 to 29 %, which is consistent with the theoretically predicted upper limit of 50 %. Conclusions The results suggest that a combination of sexual reproduction and increased asexual proliferation by female-sterile individuals probably explains the invasion and maintenance of female sterility in otherwise hermaphroditic populations. To our knowledge, this is the first report of the coexistence of female sterility and hermaphrodites in the Brassicaceae.

This article describes the use of cytogenomic and molecular approaches to explore the origin and evolution of Cardamine schulzii, a textbook example of a recent allopolyploid, in its ;110-year history of human-induced hybridization and allopolyploidy in the Swiss Alps. Triploids are typically viewed as bridges between diploids and tetraploids but rarely as parental genomes of high-level hybrids and polyploids. The genome of the triploid semifertile hybrid Cardamine 3 insueta (2n = 24, RRA) was shown to combine the parental genomes of two diploid (2n =2 x = 16) species, Cardamine amara (AA) and Cardamine rivularis (RR). These parental genomes have remained structurally stable within the triploid genome over the >100 years since its origin. Furthermore, we provide compelling evidence that the alleged recent polyploid C. schulzii is not an autohexaploid derivative of C. 3 insueta. Instead, at least two hybridization events involving C. 3 insueta and the hypotetraploid Cardamine pratensis (PPPP, 2n =4 x22 = 30) have resulted in the origin of the trigenomic hypopentaploid (2n =5 x22 = 38, PPRRA) and hypohexaploid (2n =6 x22 = 46, PPPPRA). These data show that the semifertile triploid hybrid can promote a merger of three different genomes and demonstrate how important it is to reexamine the routinely repeated textbook examples using modern techniques.

Allopolyploidization combines parental genomes and often confers broader species distribution. However, little is known about parentally transmitted gene expression underlying quantitative traits following allopolyploidization because of the complexity of polyploid genomes. The allopolyploid species Arabidopsis kamchatica is a natural hybrid of the zinc hyperaccumulator Arabidopsis halleri and of the nonaccumulator Arabidopsis lyrata. We found that A. kamchatica retained the ability to hyperaccumulate zinc from A. halleri and grows in soils with both low and high metal content. Hyperaccumulation of zinc by A. kamchatica was reduced to about half of A. halleri, but is 10-fold greater than A. lyrata. Homeologs derived from A. halleri had significantly higher levels of expression of genes such as HEAVY METAL ATPASE4 (HMA4), METAL TRANSPORTER PROTEIN1 and other metal ion transporters than those derived from A. lyrata, which suggests cis-regulatory differences. A. kamchatica has on average about half the expression of these genes compared with A. halleri due to fixed heterozygosity inherent in allopolyploids. Zinc treatment significantly changed the ratios of expression of 1% of homeologous pairs, including genes putatively involved in metal homeostasis. Resequencing data showed a significant reduction in genetic diversity over a large genomic region (290 kb) surrounding the HMA4 locus derived from the A. halleri parent compared with the syntenic A. lyrata-derived region, which suggests different evolutionary histories. We also estimated that three A. halleri-derived HMA4 copies are present in A. kamchatica. Our findings support a transcriptomic model in which environment-related transcriptional patterns of both parents are conserved but attenuated in the allopolyploids., Molecular Biology and Evolution, 33 (11), ISSN:0737-4038, ISSN:1537-1719

Polyploidization, or genome duplication, has played a critical role in the diversification of animals, fungi and plants. Little is known about the population structure and multiple origins of polyploid species because of the difficulty in identifying multiple homeologous nuclear genes. The allotetraploid species Arabidopsis kamchatica is closely related to the model species Arabidopsis thaliana and is distributed in a broader climatic niche than its parental species. Here, we performed direct sequencing of homeologous pairs of the low-copy nuclear genes WER and CHS by designing homeolog-specific primers, and obtained also chloroplast and ribosomal internal transcribed spacer sequences. Phylogenetic analysis showed that 50 individuals covering the distribution range including North America are allopolyploids derived from Arabidopsis lyrata and Arabidopsis halleri. Three major clusters within A. kamchatica were detected using Bayesian clustering. One cluster has widespread distribution. The other two are restricted to the southern part of the distribution range including Japan, where the parent A. lyrata is not currently distributed. This suggests that the mountains in Central Honshu and surrounding areas in Japan served as refugia during glacial-interglacial cycles and retained this diversity. We also found that multiple haplotypes of nuclear and chloroplast sequences of A. kamchatica are identical to those of their parental species. This indicates that multiple diploid individuals contributed to the origin of A. kamchatica. The haplotypes of low-copy nuclear genes in Japan suggest independent polyploidization events rather than introgression. Our findings suggest that self-compatibility and gene silencing occurred independently in different origins.

The evolution from outcrossing based on self-incompatibility (SI) to a selfing system is one of the most prevalent transitions in flowering plants. It has been suggested that the loss of SI in Arabidopsis thaliana is associated with pseudogene formation at the SCR male component of the S locus. Recent work, however, suggests that alternative alleles with large deletions at the S locus are also present and may be responsible for the evolution of self-compatibility in this species. We demonstrate that most of these deletion alleles are evolutionarily derived from an S haplotype (haplogroups A) that already possessed the SCR pseudogene. This haplotype and its deletion variants are nearly fixed in Europe. Together with previous transgenic data, these results suggest that the pseudogenization of PsiSCR1 gene changed the SI phenotype in the majority of A. thaliana accessions, and was a critical step in the evolution of selfing in this species. Two other haplogroups (B and C) were also identified, the former of which contains a novel and possibly functional SCR allele. In contrast to haplogroups A, these two haplogroups are found primarily in Africa and Asia. These results suggest that self-compatibility, which appears to be fixed in this species, arose multiple times with different genetic bases, and indicates that a species-specific trait is associated with parallel evolution at the molecular level.

Genomewide markers enable us to study genetic differentiation within a species and the factors underlying it at a much higher resolution than before, which advances our understanding of adaptation in organisms. We investigated genomic divergence in Metrosideros polymorpha, a woody species that occupies a wide range of ecological habitats across the Hawaiian Islands and shows remarkable phenotypic variation. Using 1659 single nucleotide polymorphism (SNP) markers annotated with the genome assembly, we examined the population genetic structure and demographic history of nine populations across five elevations and two ages of substrates on Mauna Loa, the island of Hawaii. The nine populations were differentiated into two genetic clusters distributed on the lower and higher elevations and were largely admixed on the middle elevation. Demographic modelling revealed that the two genetic clusters have been maintained in the face of gene flow, and the effective population size of the high-altitude cluster was much smaller. A FST -based outlier search among the 1659 SNPs revealed that 34 SNPs (2.05%) were likely to be under divergent selection and the allele frequencies of 21 of them were associated with environmental changes along elevations, such as temperature and precipitation. This study shows a genomic mosaic of M. polymorpha, in which contrasting divergence patterns were found. While most genomic polymorphisms were shared among populations, a small fraction of the genome was significantly differentiated between populations in diverse environments and could be responsible for the dramatic adaptation to a wide range of environments.

Self-incompatibility (SI) systems in flowering plants distinguish self- and non-self pollen to prevent inbreeding. While other SI systems rely on the self-recognition between specific male- and female-determinants, the Solanaceae family has a non-self recognition system resulting in the detoxification of female-determinants of S-ribonucleases (S-RNases), expressed in pistils, by multiple male-determinants of S-locus F-box proteins (SLFs), expressed in pollen. It is not known how many SLF components of this non-self recognition system there are in Solanaceae species, or how they evolved. We identified 16?20 SLFs in each S-haplotype in SI Petunia, from a total of 168 SLF sequences using large-scale next-generation sequencing and genomic polymerase chain reaction (PCR) techniques. We predicted the target S-RNases of SLFs by assuming that a particular S-allele must not have a conserved SLF that recognizes its own S-RNase, and validated these predictions by transformation experiments. A simple mathematical model confirmed that 16?20 SLF sequences would be adequate to recognize the vast majority of target S-RNases. We found evidence of gene conversion events, which we suggest are essential to the constitution of a non-self recognition system and also contribute to self-compatible mutations., [プレスリリース]植物が自分以外の花粉を選ぶための遺伝子セットを特定！ ～生物の非自己認識システムの進化の謎解明へ～, [Author's version title]Gene duplication and genetic exchange drive the evolution of S-RNase based self-incompatibility

To evaluate the effect of low-molecular-weight glutenin subunit (LMW-GS) genes on dough strength, locus-specific primers of LMW-GS genes and gliadin bands tightly linked to LMW-GS genes were analyzed in common wheat. Segregation analysis of the F2 progeny from a cross between ‘Haruhikari’, a good bread-making quality cultivar, and ‘Asakaze-komugi’, a poor bread-making quality cultivar, showed that dough strength significantly correlated with one amplified LMW-GS gene located at the Glu-B3 locus from ‘Haruhikari’. There was no specific reference to the gliadin bands identified as promising markers in the cross under study. The LMW-GS gene of ‘Haruhikari’ had a seven amino-acid deletion in the repetitive domain relative to ‘Asakaze-komugi’, as well as six amino-acid substitutions, three of which would be expected to cause changes in hydrophilicity. The presence of the LMW-GS gene and other LMW-GS genes tightly linked to it may affect the dough strength of wheat.

Branchio-oto (BO) syndrome is complicated with congenital preauricular fistulae, branchial fistulae (cysts), and hearing loss (sensorineural, conductive or mixed). As well as branchio-oto-renal (BOR) syndrome, it is known to be an autosomal dominant hereditary disorder. Since mutations in the EYA1 gene have been identified in both BO and BOR syndromes, mutation screening of this gene has been drawing attention as a genetic test to diagnose BOR/BO syndromes. In this study, we genetically investigated the presence of EYA1 mutations in a BO syndrome family in which we observed congenital preauricular fistulae, branchial fistulae (cysts) and hearing loss in four generations. Whereas there was a variety of phenotype expressions in this family, all subjects tested had a nonsense mutation (R264X) in exon 8 of the EYA1 gene. The present report adds further examples to support the usefulness of molecular genetic testing for the diagnosis of patients with BO syndrome.

Finding causal relationships between genotypic and phenotypic variation is a key focus of evolutionary biology, human genetics and plant breeding. To identify genome-wide patterns underlying trait diversity, we assembled a high-quality reference genome of Cardamine hirsuta, a close relative of the model plant Arabidopsis thaliana. We combined comparative genome and transcriptome analyses with the experimental tools available in C. hirsuta to investigate gene function and phenotypic diversification. Our findings highlight the prevalent role of transcription factors and tandem gene duplications in morphological evolution. We identified a specific role for the transcriptional regulators PLETHORA5/7 in shaping leaf diversity and link tandem gene duplication with differential gene expression in the explosive seed pod of C. hirsuta. Our work highlights the value of comparative approaches in genetically tractable species to understand the genetic basis for evolutionary change.

Ever since Darwin's pioneering research, the evolution of self-fertilisation (selfing) has been regarded as one of the most prevalent evolutionary transitions in flowering plants. A major mechanism to prevent selfing is the self-incompatibility (SI) recognition system, which consists of male and female specificity genes at the S-locus and SI modifier genes. Under conditions that favour selfing, mutations disabling the male recognition component are predicted to enjoy a relative advantage over those disabling the female component, because male mutations would increase through both pollen and seeds whereas female mutations would increase only through seeds. Despite many studies on the genetic basis of loss of SI in the predominantly selfing plant Arabidopsis thaliana, it remains unknown whether selfing arose through mutations in the female specificity gene (S-receptor kinase, SRK), male specificity gene (S-locus cysteine-rich protein, SCR; also known as S-locus protein 11, SP11) or modifier genes, and whether any of them rose to high frequency across large geographic regions. Here we report that a disruptive 213-base-pair (bp) inversion in the SCR gene (or its derivative haplotypes with deletions encompassing the entire SCR-A and a large portion of SRK-A) is found in 95% of European accessions, which contrasts with the genome-wide pattern of polymorphism in European A. thaliana. Importantly, interspecific crossings using Arabidopsis halleri as a pollen donor reveal that some A. thaliana accessions, including Wei-1, retain the female SI reaction, suggesting that all female components including SRK are still functional. Moreover, when the 213-bp inversion in SCR was inverted and expressed in transgenic Wei-1 plants, the functional SCR restored the SI reaction. The inversion within SCR is the first mutation disrupting SI shown to be nearly fixed in geographically wide samples, and its prevalence is consistent with theoretical predictions regarding the evolutionary advantage of mutations in male components.

Genome duplication with hybridization, or allopolyploidization, occurs commonly in plants, and is considered to be a strong force for generating new species. However, genome-wide quantification of homeolog expression ratios was technically hindered because of the high homology between homeologous gene pairs. To quantify the homeolog expression ratio using RNA-seq obtained from polyploids, a new method named HomeoRoq was developed, in which the genomic origin of sequencing reads was estimated using mismatches between the read and each parental genome. To verify this method, we first assembled the two diploid parental genomes of Arabidopsis halleri subsp. gemmifera and Arabidopsis lyrata subsp. petraea (Arabidopsis petraea subsp. umbrosa), then generated a synthetic allotetraploid, mimicking the natural allopolyploid Arabidopsis kamchatica. The quantified ratios corresponded well to those obtained by Pyrosequencing. We found that the ratios of homeologs before and after cold stress treatment were highly correlated (r = 0.870). This highlights the presence of nonstochastic polyploid gene regulation despite previous research identifying stochastic variation in expression. Moreover, our new statistical test incorporating overdispersion identified 226 homeologs (1.11% of 20 369 expressed homeologs) with significant ratio changes, many of which were related to stress responses. HomeoRoq would contribute to the study of the genes responsible for polyploid-specific environmental responses.