Genome-wide quantification of homeolog expression ratio revealed nonstochastic gene regulation in synthetic allopolyploid Arabidopsis

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 response

Phenotypic variation of a new synthetic allotetraploid Arabidopsis

The phenotypic variation of vegetative organs and reproductive organs of newly synthesized and natural Arabidopsis kamchatica genotypes was investigated in both a controlled environment and a natural environment in an experimental garden. When we compared the variation of their leaf shape as a vegetative organ, the synthetic A. kamchatica individuals grown in the garden showed larger variation compared with the individuals incubated in a growth chamber, suggesting enhanced phenotypic variation in a natural fluctuating environment. In contrast, the natural A. kamchatica genotypes did not show significant change in variation by growth condition. The phenotypic variation of floral organs by growth condition was much smaller in both synthetic and natural A. kamchatica genotypes, and the difference in variation width between the growth chamber and the garden was not significant in each genotype as well as among genotypes. The higher phenotypic variation in synthetic leaf may imply flexible transcriptomic regulation of a newly synthesized polyploid compared with a natural polyploid

Plant trichomes and a single gene GLABRA1 contribute to insect community composition on field-grown Arabidopsis thaliana

BACKGROUND: Genetic variation in plants alters insect abundance and community structure in the field; however, little is known about the importance of a single gene among diverse plant genotypes. In this context, Arabidopsis trichomes provide an excellent system to discern the roles of natural variation and a key gene, GLABRA1, in shaping insect communities. In this study, we transplanted two independent glabrous mutants (gl1-1 and gl1-2) and 17 natural accessions of Arabidopsis thaliana to two localities in Switzerland and Japan. RESULTS: Fifteen insect species inhabited the plant accessions, with the insect community composition significantly attributed to variations among plant accessions. The total abundance of leaf-chewing herbivores was negatively correlated with trichome density at both field sites, while glucosinolates had variable effects on leaf chewers between the sites. Interestingly, there was a parallel tendency for the abundance of leaf chewers to be higher on gl1-1 and gl1-2 than on their different parental accessions, Ler-1 and Col-0, respectively. Furthermore, the loss of function in the GLABRA1 gene significantly decreased the resistance of plants to the two predominant chewers; flea beetles and turnip sawflies. CONCLUSIONS: Overall, our results indicate that insect community composition significantly varies among A. thaliana accessions across two distant field sites, with GLABRA1 playing a key role in altering the abundance of leaf-chewing herbivores. Given that such a trichome variation is widely observed in Brassicaceae plants, the present study exemplifies the community-wide effect of a single plant gene on crucifer-feeding insects in the field

A low-coverage 3′ RNA-seq to detect homeolog expression in polyploid wheat

Although allopolyploid species are common among natural and crop species, it is not easy to distinguish duplicated genes, known as homeologs, during their genomic analysis. Yet, cost-efficient RNA sequencing (RNA-seq) is to be developed for large-scale transcriptomic studies such as time-series analysis and genome-wide association studies in allopolyploids. In this study, we employed a 3′ RNA-seq utilizing 3′ untranslated regions (UTRs) containing frequent mutations among homeologous genes, compared to coding sequence. Among the 3′ RNA-seq protocols, we examined a low-cost method Lasy-Seq using an allohexaploid bread wheat, Triticum aestivum. HISAT2 showed the best performance for 3′ RNA-seq with the least mapping errors and quick computational time. The number of detected homeologs was further improved by extending 1 kb of the 3′ UTR annotation. Differentially expressed genes in response to mild cold treatment detected by the 3′ RNA-seq were verified with high-coverage conventional RNA-seq, although the latter detected more differentially expressed genes. Finally, downsampling showed that even a 2 million sequencing depth can still detect more than half of expressed homeologs identifiable by the conventional 32 million reads. These data demonstrate that this low-cost 3′ RNA-seq facilitates large-scale transcriptomic studies of allohexaploid wheat and indicate the potential application to other allopolyploid species

Drop on a Bent Fibre

Inspired by the huge droplets attached on cypress tree leaf tips after rain, we find that a bent fibre can hold significantly more water in the corner than a horizontally placed fibre (typically up to three times or more). The maximum volume of the liquid that can be trapped is remarkably affected by the bending angle of the fibre and surface tension of the liquid. We experimentally find the optimal included angle ($\sim {36}{^\circ}$) that holds the most water. Analytical and semi-empirical models are developed to explain these counter-intuitive experimental observations and predict the optimal angle. The data and models could be useful for designing microfluidic and fog harvesting devices

Female sterility associated with increased clonal propagation suggests a unique combination of androdioecy and asexual reproduction in populations of Cardamine amara (Brassicaceae)

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 Brassicacea

Evolution of self-compatibility in Arabidopsis by a mutation in the male specificity gene

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 assembly and annotation ofArabidopsis halleri, a model for heavy metal hyperaccumulation and evolutionary ecology

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 makes 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 genome-wide studies of A. halleri as well as the allopolyploid Arabidopsis kamchatica derived from A. lyrata and A. halleri

Dominance in self-compatibility between subgenomes of allopolyploid Arabidopsis kamchatica shown by transgenic restoration of self-incompatibility

The evolutionary transition to self-compatibility facilitates polyploid speciation. In Arabidopsis relatives, the self-incompatibility system is characterized by epigenetic dominance modifiers, among which small RNAs suppress the expression of a recessive SCR/SP11 haplogroup. Although the contribution of dominance to polyploid self-compatibility is speculated, little functional evidence has been reported. Here we employ transgenic techniques to the allotetraploid plant A. kamchatica. We find that when the dominant SCR-B is repaired by removing a transposable element insertion, self-incompatibility is restored. This suggests that SCR was responsible for the evolution of self-compatibility. By contrast, the reconstruction of recessive SCR-D cannot restore self-incompatibility. These data indicate that the insertion in SCR-B conferred dominant self-compatibility to A. kamchatica. Dominant self-compatibility supports the prediction that dominant mutations increasing selfing rate can pass through Haldane’s sieve against recessive mutations. The dominance regulation between subgenomes inherited from progenitors contrasts with previous studies on novel epigenetic mutations at polyploidization termed genome shock

Geographic Variation in Genomic Signals of Admixture Between Two Closely Related European Sepsid Fly Species.

UNLABELLED The extent of interspecific gene flow and its consequences for the initiation, maintenance, and breakdown of species barriers in natural systems remain poorly understood. Interspecific gene flow by hybridization may weaken adaptive divergence, but can be overcome by selection against hybrids, which may ultimately promote reinforcement. An informative step towards understanding the role of gene flow during speciation is to describe patterns of past gene flow among extant species. We investigate signals of admixture between allopatric and sympatric populations of the two closely related European dung fly species Sepsis cynipsea and S. neocynipsea (Diptera: Sepsidae). Based on microsatellite genotypes, we first inferred a baseline demographic history using Approximate Bayesian Computation. We then used genomic data from pooled DNA of natural and laboratory populations to test for past interspecific gene flow based on allelic configurations discordant with the inferred population tree (ABBA-BABA test with D-statistic). Comparing the detected signals of gene flow with the contemporary geographic relationship among interspecific pairs of populations (sympatric vs. allopatric), we made two contrasting observations. At one site in the French Cevennes, we detected an excess of past interspecific gene flow, while at two sites in Switzerland we observed lower signals of past microsatellite genotypes gene flow among populations in sympatry compared to allopatric populations. These results suggest that the species boundaries between these two species depend on the past and/or present eco-geographic context in Europe, which indicates that there is no uniform link between contemporary geographic proximity and past interspecific gene flow in natural populations. SUPPLEMENTARY INFORMATION The online version contains supplementary material available at 10.1007/s11692-023-09612-5