Natural variation in GL1 and its effects on trichome density in Arabidopsis thaliana

Transcriptome Analysis Reveals Key Genes Involved in Trichome Formation in Pepper (Capsicum annuum L.)

Trichomes are specialized organs located in the plant epidermis that play important defense roles against biotic and abiotic stresses. However, the mechanisms regulating the development of pepper epidermal trichomes and the related regulatory genes at the molecular level are not clear. Therefore, we performed transcriptome analyses of A114 (less trichome) and A115 (more trichome) to dig deeper into the genes involved in the regulatory mechanisms of epidermal trichome development in peppers. In this study, the epidermal trichome density of A115 was found to be higher by phenotypic observation and was highest in the leaves at the flowering stage. A total of 39,261 genes were quantified by RNA-Seq, including 11,939 genes not annotated in the previous genome analysis and 18,833 differentially expressed genes. Based on KEGG functional enrichment, it was found that DEGs were mainly concentrated in three pathways: plant-pathogen interaction, MAPK signaling pathway-plant, and plant hormone signal transduction. We further screened the DEGs associated with the development of epidermal trichomes in peppers, and the expression of the plant signaling genes GID1B-like (Capana03g003488) and PR-6 (Capana09g001847), the transcription factors MYB108 (Capana05g002225) and ABR1-like (Capana04g001261), and the plant resistance genes PGIP-like (Capana09g002077) and At5g49770 (Capana08g001721) in the DEGs were higher at A115 compared to A114, and were highly expressed in leaves at the flowering stage. In addition, based on the WGCNA results and the establishment of co-expression networks showed that the above genes were highly positively correlated with each other. The transcriptomic data and analysis of this study provide a basis for the study of the regulatory mechanisms of pepper epidermal trichomes.

The GaKAN2, a KANADI transcription factor, modulates stem trichomes in Gossypium arboreum

KEY MESSAGE

GaKAN2, a member of the KANADI family, was found to be widely expressed in the cotton tissues and regulates trichome development through complex pathways. Cotton trichomes are believed to be the defense barrier against insect pests. Cotton fiber and trichomes are single-cell epidermal extensions with shared regulatory mechanisms. Despite several studies underlying mechanism of trichome development remains elusive. The KANADI is one of the key transcription factors (TFs) family, regulating Arabidopsis trichomes growth. However, the function of KANADI genes in cotton remains unknown. In the current study genome-wide scanning, transcriptomic analysis, gene silencing, subcellular localization, and yeast two-hybrid techniques were employed to decipher the function of KANADI TFs family genes in cotton crop. A total of 7 GaKAN genes were found in the Gossypium arboreum. Transcriptomic data revealed that these genes were significantly expressed in stem and root. Moreover, GaKAN2 was widely expressed in other tissues also. Subsequently, we selected GaKAN2 to validate the function of KANADI genes. Silencing of GaKAN2 resulted in a 24.99% decrease in single-cell trichomes and an 11.33% reduction in internodal distance, indicating its potential role in regulating trichomes and plant growth. RNA-Seq analysis elucidated that GaSuS and GaERS were the downstream genes of GaKAN2. The transcriptional activation and similarity in silencing phenotype between GaKAN2 and GaERS suggested that GaKAN2 regulates trichomes development through GaERS. Moreover, KEGG analysis revealed that a significant number of genes were enriched in the biosynthesis of secondary metabolites and plant hormone signal transduction pathways, thereby suggesting that GaKAN2 regulates the stem trichomes and plant growth. The GFP subcellular localization and yeast transcriptional activation analysis elucidated that GaKAN2 was located in the nucleus and capable of regulating the transcription of downstream genes. This study elucidated the function and characteristics of the KANADI gene family in cotton, providing a fundamental basis for further research on GaKAN2 gene in cotton plant trichomes and plant developmental processes.

Transcriptome Analysis Reveals Coexpression Networks and Hub Genes Involved in Papillae Development in Lilium auratum

Lilium is a genus of important ornamental plants with many colouring pattern variations. Lilium auratum is the parent of Oriental hybrid lilies. A typical feature of L. auratum is the presence of red-orange special raised spots named papillae on the interior tepals. Unlike the usual raised spots, the papillae are slightly rounded or connected into sheets and usually have hairy tips. To elucidate the potential genes regulating papillae development in L. auratum, we performed high-throughput sequencing of its tepals at different stages. Genes involved in the flavonoid biosynthesis pathway were significantly enriched during the colouration of the papillae, and CHS, F3H, F3'H, FLS, DFR, ANS, and UFGT were significantly upregulated. To identify the key genes involved in the papillae development of L. auratum, we performed weighted gene coexpression network analysis (WGCNA) and further analysed four modules. In total, 51, 24, 1, and 6 hub genes were identified in four WGCNA modules, MEbrown, MEyellow, MEpurple, and MEred, respectively. Then, the coexpression networks were constructed, and important genes involved in trichome development and coexpressed with anthocyanin biosynthesis genes, such as TT8, TTG1, and GEM, were identified. These results indicated that the papillae are essentially trichomes that accumulate anthocyanins. Finally, we randomly selected 12 hub genes for qRT-PCR analysis to verify the accuracy of our RNA-Seq analysis. Our results provide new insights into the papillae development in L. auratum flowers.

Mechanodetection of neighbor plants elicits adaptive leaf movements through calcium dynamics

Plants detect their neighbors via various cues, including reflected light and touching of leaf tips, which elicit upward leaf movement (hyponasty). It is currently unknown how touch is sensed and how the signal is transferred from the leaf tip to the petiole base that drives hyponasty. Here, we show that touch-induced hyponasty involves a signal transduction pathway that is distinct from light-mediated hyponasty. We found that mechanostimulation of the leaf tip upon touching causes cytosolic calcium ([Ca2+]cyt induction in leaf tip trichomes that spreads towards the petiole. Both perturbation of the calcium response and the absence of trichomes reduce touch-induced hyponasty. Finally, using plant competition assays, we show that touch-induced hyponasty is adaptive in dense stands of Arabidopsis. We thus establish a novel, adaptive mechanism regulating hyponastic leaf movement in response to mechanostimulation by neighbors in dense vegetation.

Glandular trichome development, morphology, and maturation are influenced by plant age and genotype in high THC-containing cannabis (Cannabis sativa L.) inflorescences

BACKGROUND

Glandular capitate trichomes which form on bract tissues of female inflorescences of high THC-containing Cannabis sativa L. plants are important sources of terpenes and cannabinoids. The influence of plant age and cannabis genotype on capitate trichome development, morphology, and maturation has not been extensively studied. Knowledge of the various developmental changes that occur in trichomes over time and the influence of genotype and plant age on distribution, numbers, and morphological features should lead to a better understanding of cannabis quality and consistency.

METHODS

Bract tissues of two genotypes-"Moby Dick" and "Space Queen"-were examined from 3 weeks to 8 weeks of flower development using light and scanning electron microscopy. Numbers of capitate trichomes on upper and lower bract surfaces were recorded at different positions within the inflorescence. Observations on distribution, extent of stalk formation, glandular head diameter, production of resin, and extent of dehiscence and senescence were made at various time points. The effects of post-harvesting handling and drying on trichome morphology were examined in an additional five genotypes.

RESULTS

Two glandular trichome types-bulbous and capitate (sessile or stalked)-were observed. Capitate trichome numbers and stalk length were significantly (P = 0.05) greater in "Space Queen" compared to "Moby Dick" at 3 and 6 weeks of flower development. Significantly more stalked-capitate trichomes were present on lower compared to upper bract surfaces at 6 weeks in both genotypes, while sessile-capitate trichomes predominated at 3 weeks. Epidermal and hypodermal cells elongated to different extents during stalk formation, producing significant variation in length (from 20 to 1100 μm). Glandular heads ranged from 40 to 110 μm in diameter. Maturation of stalked-capitate glandular heads was accompanied by a brown color development, reduced UV autofluorescence, and head senescence and dehiscence. Secreted resinous material from glandular heads appeared as droplets on the cuticular surface that caused many heads to stick together or collapse. Trichome morphology was affected by the drying process.

CONCLUSION

Capitate trichome numbers, development, and degree of maturation were influenced by cannabis genotype and plant age. The observations of trichome development indicate that asynchronous formation leads to different stages of trichome maturity on bracts. Trichome stalk lengths also varied between the two genotypes selected for study as well as over time. The variability in developmental stage and maturation between genotypes can potentially lead to variation in total cannabinoid levels in final product. Post-harvest handling and drying were shown to affect trichome morphology.

CRISPR/Cas9-mediated genome editing of RsGL1a and RsGL1b in radish (Raphanus sativus L.)

The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) is a powerful tool widely used for genome editing in various organisms, including plants. It introduces and facilitates the study of rare genetic mutations in a short time and is a potent tool to assist in plant molecular breeding. Radish (Raphanus sativus L.) is an important Brassicaceae vegetable cultivated and consumed worldwide. However, the application of the CRISPR/Cas9 system is limited by the absence of an efficient transformation system in radish. This study aimed to establish a CRISPR/Cas9 system in radish employing the Agrobacterium-mediated genetic transformation system reported recently. For this purpose, we performed genome editing using the CRISPR/Cas9 system targeting the GLABRA1 (GL1) orthologs, RsGL1a and RsGL1b, that induces leaf trichome formation in radish. A Cas9/single guide RNA (sgRNA) vector with a common sgRNA corresponding to RsGL1a and RsGL1b was transferred. A total of eight T₀ plants were analyzed, of which six (editing efficiency 75%) had a mutated RsGL1a, five (62.5%) had a mutated RsGL1b, and five showed mutations in both RsGL1a and RsGL1b. Most mutations in T₀ plants were short (<3 bp) deletions or insertions, causing frameshift mutations that might produce non-functional proteins. Chimeric mutations were detected in several T₀ generation plants. In the T₁ generation, the hairless phenotype was observed only in plants with knockout mutations in both RsGL1a and RsGL1b. The majority of mutant alleles in T₀ plants, with the exception of the chimeric mutant plants detected, were stably inherited in the T₁ generation. In conclusion, we successfully knocked out RsGL1a and RsGL1b using the CRISPR/Cas9 system and demonstrated that both RsGL1a and RsGL1b independently contribute to the induction of leaf trichome formation in radish. In this study, genome-edited plants without T-DNA, which are useful as breeding material, were obtained. The findings prove the feasibility of genome editing in radish using a CRISPR/Cas9 system that could accelerate its molecular breeding to improve agronomically desirable traits.

Differential environmental and genomic architectures shape the natural diversity for trichome patterning and morphology in different Arabidopsis organs

Despite the adaptive and taxonomic relevance of the natural diversity for trichome patterning and morphology, the molecular and evolutionary mechanisms underlying these traits remain mostly unknown, particularly in organs other than leaves. In this study, we address the ecological, genetic and molecular bases of the natural variation for trichome patterning and branching in multiple organs of Arabidopsis (Arabidopsis thaliana). To this end, we characterized a collection of 191 accessions and carried out environmental and genome-wide association (GWA) analyses. Trichome amount in different organs correlated negatively with precipitation in distinct seasons, thus suggesting a precise fit between trichome patterning and climate throughout the Arabidopsis life cycle. In addition, GWA analyses showed small overlapping between the genes associated with different organs, indicating partly independent genetic bases for vegetative and reproductive phases. These analyses identified a complex locus on chromosome 2, where two adjacent MYB genes (ETC2 and TCL1) displayed differential effects on trichome patterning in several organs. Furthermore, analyses of transgenic lines carrying different natural alleles demonstrated that TCL1 accounts for the variation for trichome patterning in all organs, and for stem trichome branching. By contrast, two other MYB genes (TRY and GL1), mainly showed effects on trichome patterning or branching, respectively.

Transcriptome profiling reveals key genes in regulation of the tepal trichome development in Lilium pumilum D.C

A number of potential genes and pathways involved in tepal trichome development were identified in a natural lily mutant by transcriptome analysis and were confirmed with trichome and trichomeless species. Trichome is a specialized structure found on the surface of the plant with an important function in survival against abiotic and biotic stress. It is also an important economic trait in crop breeding. Extensive research has investigated the foliar trichome in model plants (Arabidopsis and tomato). However, the developmental mechanism of tepal trichome remains elusive. Lilium pumilum is an edible ornamental bulb and a good breeding parent possessing cold and salt-alkali resistance. Here, we found a natural mutant of Lilium pumilum grown on a highland whose tepals are covered by trichomes. Our data indicate that trichomes of the mutant are multicellular and branchless. Notably, stomata are also developed on the tepal of the mutant as well, suggesting there may be a correlation between trichome and stomata regulation. Furthermore, we isolated 27 differentially expressed genes (DEGs) by comparing the transcriptome profiling between the natural mutant and the wild type. These 27 genes belong to 4 groups: epidermal cell cycle and division, trichome morphogenesis, stress response, and transcription factors. Quantitative real-time PCR in Lilium pumilum (natural mutant and the wild type) and other lily species (Lilium leichtlinii var. maximowiczii/trichome; Lilium davidii var. willmottiae/, trichomeless) confirmed the validation of RNA-seq data and identified several trichome-related genes.

A novel gene LbHLH from the halophyte Limonium bicolor enhances salt tolerance via reducing root hair development and enhancing osmotic resistance

BACKGROUND

Identifying genes involved in salt tolerance in the recretohalophyte Limonium bicolor could facilitate the breeding of crops with enhanced salt tolerance. Here we cloned the previously uncharacterized gene LbHLH and explored its role in salt tolerance.

RESULTS

The 2,067-bp open reading frame of LbHLH encodes a 688-amino-acid protein with a typical helix-loop-helix (HLH) domain. In situ hybridization showed that LbHLH is expressed in salt glands of L. bicolor. LbHLH localizes to the nucleus, and LbHLH is highly expressed during salt gland development and in response to NaCl treatment. To further explore its function, we heterologously expressed LbHLH in Arabidopsis thaliana under the 35S promoter. The overexpression lines showed significantly increased trichome number and reduced root hair number. LbHLH might interact with GLABRA1 to influence trichome and root hair development, as revealed by yeast two-hybrid analysis. The transgenic lines showed higher germination percentages and longer roots than the wild type under NaCl treatment. Analysis of seedlings grown on medium containing sorbitol with the same osmotic pressure as 100 mM NaCl demonstrated that overexpressing LbHLH enhanced osmotic resistance.

CONCLUSION

These results indicate that LbHLH enhances salt tolerance by reducing root hair development and enhancing osmotic resistance under NaCl stress.

MYB transcription factors drive evolutionary innovations in Arabidopsis fruit trichome patterning

Both inter- and intra-specific diversity has been described for trichome patterning in fruits, which is presumably involved in plant adaptation. However, the mechanisms underlying this developmental trait have been hardly addressed. Here we examined natural populations of Arabidopsis (Arabidopsis thaliana) that develop trichomes in fruits and pedicels, phenotypes previously not reported in the Arabidopsis genus. Genetic analyses identified five loci, MALAMBRUNO 1-5 (MAU1-5), with MAU2, MAU3, and MAU5 showing strong epistatic interactions that are necessary and sufficient to display these traits. Functional characterization of these three loci revealed cis-regulatory mutations in TRICHOMELESS1 and TRIPTYCHON, as well as a structural mutation in GLABRA1. Therefore, the multiple mechanisms controlled by three MYB transcription factors of the core regulatory network for trichome patterning have jointly been modulated to trigger trichome development in fruits. Furthermore, analyses of worldwide accessions showed that these traits and mutations only occur in a highly differentiated relict lineage from the Iberian Peninsula. In addition, these traits and alleles were associated with low spring precipitation, which suggests that trichome development in fruits and pedicels might be involved in climatic adaptation. Thus, we show that the combination of synergistic mutations in a gene regulatory circuit has driven evolutionary innovations in fruit trichome patterning in Arabidopsis.

Construction of a Full-Length cDNA Over-Expressing Library to Identify Valuable Genes from Populus tomentosa

Poplar wood is the main source of renewable biomass energy worldwide, and is also considered to be a model system for studying woody plants. The Full-length cDNA Over-eXpressing (FOX) gene hunting system is an effective method for generating gain-of-function mutants. Large numbers of novel genes have successfully been identified from many herbaceous plants according to the phenotype of gain-of-function mutants under normal or abiotic stress conditions using this system. However, the system has not been used for functional gene identification with high-throughput mutant screening in woody plants. In this study, we constructed a FOX library from the Chinese white poplar, Populus tomentosa. The poplar cDNA library was constructed into the plant expression vector pEarleyGate101 and further transformed into Arabidopsis thaliana (thale cress). We collected 1749 T1 transgenic plants identified by PCR. Of these, 593 single PCR bands from different transgenic lines were randomly selected for sequencing, and 402 diverse sequences of poplar genes were isolated. Most of these genes were involved in photosynthesis, environmental adaptation, and ribosome biogenesis based on Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation. We characterized in detail two mutant lines carrying PtoCPCa or PtoWRKY13 cDNA insertions. Phenotypic characterization showed that overexpression of these genes in A. thaliana affected trichome development or secondary cell wall (SCW) deposition, respectively. Together, the Populus-FOX-Arabidopsis library generated in our experiments will be helpful for efficient discovery of novel genes in poplar.

Differentiation in the genetic basis of stem trichome development between cultivated tetraploid cotton species

BACKGROUND

Cotton stem trichomes and seed fibers are each single celled structures formed by protrusions of epidermal cells, and were found sharing the overlapping molecular mechanism. Compared with fibers, cotton stem trichomes are more easily observed, but the molecular mechanisms underlying their development are still poorly understood.

RESULTS

In this study, Gossypium hirsutum (Gh) and G. barbadense (Gb) were found to differ greatly in percentages of varieties/accessions with glabrous stems and in trichome density, length, and number per trichopore. Gh varieties normally had long singular and clustered trichomes, while Gb varieties had short clustered trichomes. Genetic mapping using five F₂ populations from crosses between glabrous varieties and those with different types of stem trichomes revealed that much variation among stem trichome phenotypes could be accounted for by different combinations of genes/alleles on Chr. 06 and Chr. 24. The twenty- six F₁ generations from crosses between varieties with different types of trichomes had varied phenotypes, further suggesting that the trichomes of tetraploid cotton were controlled by different genes/alleles. Compared to modern varieties, a greater proportion of Gh wild accessions were glabrous or had shorter and denser trichomes; whereas a smaller proportion of Gb primitive accessions had glabrous stems. A close correlation between fuzz fiber number and stem trichome density was observed in both Gh and Gb primitive accessions and modern varieties.

CONCLUSION

Based on these findings, we hypothesize that stem trichomes evolved in parallel with seed fibers during the domestication of cultivated tetraploid cotton. In addition, the current results illustrated that stem trichome can be used as a morphological index of fiber quality in cotton conventional breeding.

Analysis and review of trichomes in plants

BACKGROUND

Trichomes play a key role in the development of plants and exist in a wide variety of species.

RESULTS

In this paper, it was reviewed that the structure and morphology characteristics of trichomes, alongside the biological functions and classical regulatory mechanisms of trichome development in plants. The environment factors, hormones, transcription factor, non-coding RNA, etc., play important roles in regulating the initialization, branching, growth, and development of trichomes. In addition, it was further investigated the atypical regulation mechanism in a non-model plant, found that regulating the growth and development of tea (Camellia sinensis) trichome is mainly affected by hormones and the novel regulation factors.

CONCLUSIONS

This review further displayed the complex and differential regulatory networks in trichome initiation and development, provided a reference for basic and applied research on trichomes in plants.

Functional identification of PsMYB57 involved in anthocyanin regulation of tree peony

BACKGROUND

R2R3 myeloblastosis (MYB) genes are widely distributed in plants and comprise one of the largest transcription factor gene families. They play important roles in the regulatory networks controlling development, metabolism, and stress responses. Researches on functional genes in tree peony are still in its infancy. To date, few MYB genes have thus far been reported.

RESULTS

In this study, we constructed a comprehensive reference gene set by transcriptome sequencing to obtain R2R3 MYB genes. The transcriptomes of eight different tissues were sequenced, and 92,837 unigenes were obtained with an N50 of 1662 nt. A total of 48,435 unigenes (77.98%) were functionally annotated in public databases. Based on the assembly, we identified 57 R2R3 MYB genes containing full-length open reading frames, which clustered into 35 clades by phylogenetic analysis. PsMYB57 clustered with anthocyanin regulation genes in Arabidopsis and was mainly transcribed in the buds and young leaves. The overexpression of PsMYB57 induced anthocyanin accumulation in tobacco, and four detected anthocyanin structural genes, including NtCHS, NtF3'H, NtDFR, and NtANS, were upregulated. The two endogenous bHLH genes NtAn1a and NtAn1b were also upregulated and may work in combination with PsMYB57 in regulating anthocyanin structural genes.

CONCLUSIONS

Our study offers a useful reference to the selection of candidate MYB genes for further functional studies in tree peony. Function analysis of PsMYB57 is helpful to understand the color accumulation in vegetative organs of tree peony. PsMYB57 is also a promising resource to improve plant color in molecular breeding.

AthCNV: A Map of DNA Copy Number Variations in the Arabidopsis Genome

Copy number variations (CNVs) greatly contribute to intraspecies genetic polymorphism and phenotypic diversity. Recent analyses of sequencing data for >1000 Arabidopsis (Arabidopsis thaliana) accessions focused on small variations and did not include CNVs. Here, we performed genome-wide analysis and identified large indels (50 to 499 bp) and CNVs (500 bp and larger) in these accessions. The CNVs fully overlap with 18.3% of protein-coding genes, with enrichment for evolutionarily young genes and genes involved in stress and defense. By combining analysis of both genes and transposable elements (TEs) affected by CNVs, we revealed that the variation statuses of genes and TEs are tightly linked and jointly contribute to the unequal distribution of these elements in the genome. We also determined the gene copy numbers in a set of 1060 accessions and experimentally validated the accuracy of our predictions by multiplex ligation-dependent probe amplification assays. We then successfully used the CNVs as markers to analyze population structure and migration patterns. Finally, we examined the impact of gene dosage variation triggered by a CNV spanning the SEC10 gene on SEC10 expression at both the transcript and protein levels. The catalog of CNVs, CNV-overlapping genes, and their genotypes in a top model dicot will stimulate the exploration of the genetic basis of phenotypic variation.

Genome-wide association study reveals new genes involved in leaf trichome formation in polyploid oilseed rape (Brassica napus L.)

Leaf trichomes protect against various biotic and abiotic stresses in plants. However, there is little knowledge about this trait in oilseed rape (Brassica napus). Here, we demonstrated that hairy leaves were less attractive to Plutella xylostella larvae than glabrous leaves. We established a core germplasm collection with 290 accessions for a genome-wide association study (GWAS) of the leaf trichome trait in oilseed rape. We compared the transcriptomes of the shoot apical meristem (SAM) between hairy- and glabrous-leaf genotypes to narrow down the candidate genes identified by GWAS. The single nucleotide polymorphisms and the different transcript levels of BnaA.GL1.a, BnaC.SWEET4.a, BnaC.WAT1.a and BnaC.WAT1.b corresponded to the divergence of the hairy- and glabrous-leaf phenotypes, indicating the role of sugar and/or auxin signalling in leaf trichome initiation. The hairy-leaf SAMs had lower glucose and sucrose contents but higher expression of putative auxin responsive factors than the glabrous-leaf SAMs. Spraying of exogenous auxin (8 μm) increased leaf trichome number in certain genotypes, whereas spraying of sucrose (1%) plus glucose (6%) slightly repressed leaf trichome initiation. These data contribute to the existing knowledge about the genetic control of leaf trichomes and would assist breeding towards the desired leaf surface type in oilseed rape.

Genetic and physiological mechanisms of freezing tolerance in locally adapted populations of a winter annual

PREMISE

Despite myriad examples of local adaptation, the phenotypes and genetic variants underlying such adaptive differentiation are seldom known. Recent work on freezing tolerance and local adaptation in ecotypes of Arabidopsis thaliana from Italy and Sweden provides an essential foundation for uncovering the genotype-phenotype-fitness map for an adaptive response to a key environmental stress.

METHODS

We examined the consequences of a naturally occurring loss-of-function (LOF) mutation in an Italian allele of the gene that encodes the transcription factor CBF2, which underlies a major freezing-tolerance locus. We used four lines with a Swedish genetic background, each containing a LOF CBF2 allele. Two lines had introgression segments containing the Italian CBF2 allele, and two contained deletions created using CRISPR-Cas9. We used a growth chamber experiment to quantify freezing tolerance and gene expression before and after cold acclimation.

RESULTS

Freezing tolerance was lower in the Italian (11%) compared to the Swedish (72%) ecotype, and all four experimental CBF2 LOF lines had reduced freezing tolerance compared to the Swedish ecotype. Differential expression analyses identified 10 genes for which all CBF2 LOF lines, and the IT ecotype had similar patterns of reduced cold responsive expression compared to the SW ecotype.

CONCLUSIONS

We identified 10 genes that are at least partially regulated by CBF2 that may contribute to the differences in cold-acclimated freezing tolerance between the Italian and Swedish ecotypes. These results provide novel insight into the molecular and physiological mechanisms connecting a naturally occurring sequence polymorphism to an adaptive response to freezing conditions.

Transcriptional Variation in Glucosinolate Biosynthetic Genes and Inducible Responses to Aphid Herbivory on Field-Grown Arabidopsis thaliana

Recently, increasing attempts have been made to understand how plant genes function in natura. In this context, transcriptional profiles represent plant physiological status in response to environmental stimuli. Herein, we combined high-throughput RNA-Seq with insect survey data on 19 accessions of Arabidopsis thaliana grown at a field site in Switzerland. We found that genes with the gene ontology (GO) annotations of "glucosinolate biosynthetic process" and "response to insects" were most significantly enriched, and the expression of these genes was highly variable among plant accessions. Nearly half of the total expression variation in the glucosinolate biosynthetic genes (AOPs, ESM1, ESP, and TGG1) was explained by among-accession variation. Of these genes, the expression level of AOP3 differed among Col-0 accession individuals depending on the abundance of the mustard aphid (Lipaphis erysimi). We also found that the expression of the major cis-jasmone activated gene CYP81D11 was positively correlated with the number of flea beetles (Phyllotreta striolata and Phyllotreta atra). Combined with the field RNA-Seq data, bioassays confirmed that AOP3 was up-regulated in response to attack by mustard aphids. The combined results from RNA-Seq and our ecological survey illustrate the feasibility of using field transcriptomics to detect an inducible defense, providing a first step towards an in natura understanding of biotic interactions involving phenotypic plasticity.

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 spatiotemporally regulated transcriptional complex underlies heteroblastic development of leaf hairs in Arabidopsis thaliana

Heteroblasty refers to a phenomenon that a plant produces morphologically or functionally different lateral organs in an age-dependent manner. In the model plant Arabidopsis thaliana, the production of trichomes (epidermal leaf hairs) on the abaxial (lower) side of leaves is a heteroblastic mark for the juvenile-to-adult transition. Here, we show that the heteroblastic development of abaxial trichomes is regulated by a spatiotemporally regulated complex comprising the leaf abaxial fate determinant (KAN1) and the developmental timer (miR172-targeted AP2-like proteins). We provide evidence that a short-distance chromatin loop brings the downstream enhancer element into close association with the promoter elements of GL1, which encodes a MYB transcription factor essential for trichome initiation. During juvenile phase, the KAN1-AP2 repressive complex binds to the downstream sequence of GL1 and represses its expression through chromatin looping. As plants age, the gradual reduction in AP2-like protein levels leads to decreased amount of the KAN1-AP2 complex, thereby licensing GL1 expression and the abaxial trichome initiation. Our results thus reveal a novel molecular mechanism by which a heteroblastic trait is governed by integrating age and leaf polarity cue in plants.

Drought adaptation in Arabidopsis thaliana by extensive genetic loss-of-function

Interdisciplinary syntheses are needed to scale up discovery of the environmental drivers and molecular basis of adaptation in nature. Here we integrated novel approaches using whole genome sequences, satellite remote sensing, and transgenic experiments to study natural loss-of-function alleles associated with drought histories in wild Arabidopsis thaliana. The genes we identified exhibit population genetic signatures of parallel molecular evolution, selection for loss-of-function, and shared associations with flowering time phenotypes in directions consistent with longstanding adaptive hypotheses seven times more often than expected by chance. We then confirmed predicted phenotypes experimentally in transgenic knockout lines. These findings reveal the importance of drought timing to explain the evolution of alternative drought tolerance strategies and further challenge popular assumptions about the adaptive value of genetic loss-of-function in nature. These results also motivate improved species-wide sequencing efforts to better identify loss-of-function variants and inspire new opportunities for engineering climate resilience in crops.

Water shortages caused by droughts lead to crop losses that affect billions of people around the world each year. By discovering how wild plants adapt to drought, it may be possible to identify traits and genes that help to improve the growth of crop plants when water is scarce. It has been suggested that plants have adapted to droughts by flowering at times of the year when droughts are less likely to occur. For example, if droughts are more likely to happen in spring, the plants may delay flowering until the summer.

Arabidopsis thaliana is a small plant that is found across Eurasia, Africa and North America, including in areas that are prone to drought at different times of the year. Individual plants of the same species may carry different versions of the same gene (known as alleles). Some of these alleles may not work properly and are referred to as loss-of-function alleles. Monroe et al. investigated whether A. thaliana plants carry any loss-of-function alleles that are associated with droughts happening in the spring or summer, and whether they are linked to when those plants will flower.

Monroe et al. analyzed satellite images collected over the last 30 years to measure when droughts have occurred. Next, they searched genome sequences of Arabidopsis thaliana for alleles that might help the plants to adapt to droughts in the spring or summer. Combining the two approaches revealed that loss-of-function alleles associated with spring droughts were strongly predicted to be associated with the plants flowering later in the year. Similarly, loss-of-function alleles associated with summer droughts were predicted to be associated with the plants flowering earlier in the year.

These findings support the idea that plants can adapt to drought by changing when they produce flowers, and suggest that loss-of-function alleles play a major role in this process. New techniques for editing genes mean it is easier than ever to generate new loss-of-function alleles in specific genes. Therefore, the results presented by Monroe et al. may help researchers to develop new varieties of crop plants that are better adapted to droughts.

Distributions of Mutational Effects and the Estimation of Directional Selection in Divergent Lineages of Arabidopsis thaliana

Mutations are crucial to evolution, providing the ultimate source of variation on which natural selection acts. Due to their key role, the distribution of mutational effects on quantitative traits is a key component to any inference regarding historical selection on phenotypic traits. In this paper, we expand on a previously developed test for selection that could be conducted assuming a Gaussian mutation effect distribution by developing approaches to also incorporate any of a family of heavy-tailed Laplace distributions of mutational effects. We apply the test to detect directional natural selection on five traits along the divergence of Columbia and Landsberg lineages of Arabidopsis thaliana, constituting the first test for natural selection in any organism using quantitative trait locus and mutation accumulation data to quantify the intensity of directional selection on a phenotypic trait. We demonstrate that the results of the test for selection can depend on the mutation effect distribution specified. Using the distributions exhibiting the best fit to mutation accumulation data, we infer that natural directional selection caused divergence in the rosette diameter and trichome density traits of the Columbia and Landsberg lineages.

Fine-scale frequency differentiation along a herbivory gradient in the trichome dimorphism of a wild Arabidopsis

Geographic variation is commonly observed in plant resistance traits, where plant species might experience different selection pressure across a heterogeneous landscape. Arabidopsis halleri subsp. gemmifera is dimorphic for trichome production, generating two morphs, trichome-producing (hairy) and trichomeless (glabrous) plants. Trichomes of A. halleri are known to confer resistance against the white butterfly, cabbage sawfly, and brassica leaf beetle, but not against flea beetles. We combined leaf damage, microclimate, and microsatellite loci data of 26 A. halleri populations in central Japan, to explore factors responsible for fine-scale geographic variation in the morph frequency. We found that hairy plants were less damaged than glabrous plants within populations, but the among-site variation was the most significant source of variation in the individual-level damage. Fixation index (Gst″) of a putative trichome locus exhibited a significant divergence along population-level damage with an exception of an outlier population, inferring the local adaptation to herbivory. Notably, this outlier was a population wherein our previous study reported a balancing role of the brassica leaf beetle Phaedon brassicae on the morph frequency. This differentiation of the trichome locus was unrelated to neutral genetic differentiation (evaluated by Gst″ of microsatellite loci) and meteorological factors (including temperature and solar radiation). The present findings, combined with those of our previous work, provide suggestive evidence that herbivore-driven divergence and occasional outbreak of a specific herbivore have jointly contributed to the ecogeographic pattern in the frequency of two morphs.

The genome sequence of Barbarea vulgaris facilitates the study of ecological biochemistry

The genus Barbarea has emerged as a model for evolution and ecology of plant defense compounds, due to its unusual glucosinolate profile and production of saponins, unique to the Brassicaceae. One species, B. vulgaris, includes two ‘types’, G-type and P-type that differ in trichome density, and their glucosinolate and saponin profiles. A key difference is the stereochemistry of hydroxylation of their common phenethylglucosinolate backbone, leading to epimeric glucobarbarins. Here we report a draft genome sequence of the G-type, and re-sequencing of the P-type for comparison. This enables us to identify candidate genes underlying glucosinolate diversity, trichome density, and study the genetics of biochemical variation for glucosinolate and saponins. B. vulgaris is resistant to the diamondback moth, and may be exploited for “dead-end” trap cropping where glucosinolates stimulate oviposition and saponins deter larvae to the extent that they die. The B. vulgaris genome will promote the study of mechanisms in ecological biochemistry to benefit crop resistance breeding.

Genome-Wide Survey Reveals Transcriptional Differences Underlying the Contrasting Trichome Phenotypes of Two Sister Desert Poplars

Trichomes, which are widely used as an important diagnostic characteristic in plant species delimitation, play important roles in plant defense and adaptation to adverse environments. In this study, we used two sister poplar species, Populus pruinosa and Populus euphratica—which have, respectively, dense and sparse trichomes—to examine the genetic differences associated with these contrasting phenotypes. The results showed that 42 and 45 genes could be identified as candidate genes related to trichomes in P. pruinosa and P. euphratica, respectively; most of these genes possessed high degrees of diversification in their coding sequences, but they were similar in intron/exon structure in the two species. We also found that most of the candidate trichome genes were expressed at higher levels in P. pruinosa, which has dense trichomes, than in P. euphratica, where there are few trichomes. Based on analyses of transcriptional profiles, a total of 195 genes, including many transcription factors, were found to show distinct differences in expression. The results of gene function annotation suggested that the genes identified as having contrasting levels of expression level are mainly associated with trichome elongation, ATPase activity, and hormone transduction. Changes in the expression of these and other related genes with high sequence diversification may have contributed to the contrast in the pattern of trichome phenotypes between the two species.

The genetics of extreme microgeographic adaptation: an integrated approach identifies a major gene underlying leaf trichome divergence in Yellowstone Mimulus guttatus

Microgeographic adaptation provides a particularly interesting context for understanding the genetic basis of phenotypic divergence and may also present unique empirical challenges. In particular, plant adaptation to extreme soil mosaics may generate barriers to gene flow or shifts in mating system that confound simple genomic scans for adaptive loci. Here, we combine three approaches - quantitative trait locus (QTL) mapping of candidate intervals in controlled crosses, population resequencing (PoolSeq) and analyses of wild recombinant individuals - to investigate one trait associated with Mimulus guttatus (yellow monkeyflower) adaptation to geothermal soils in Yellowstone National Park. We mapped a major QTL causing dense leaf trichomes in thermally adapted plants to a <50-kb region of linkage Group 14 (Tr14) previously implicated in trichome divergence between independent M. guttatus populations. A PoolSeq scan of Tr14 region revealed a cluster of six genes, coincident with the inferred QTL peak, with high allele frequency differences sufficient to explain observed phenotypic differentiation. One of these, the R2R3 MYB transcription factor Migut.N02661, is a plausible functional candidate and was also strongly associated (r²  = 0.27) with trichome phenotype in analyses of wild-collected admixed individuals. Although functional analyses will be necessary to definitively link molecular variants in Tr14 with trichome divergence, our analyses are a major step in that direction. They point to a simple, and parallel, genetic basis for one axis of Mimulus guttatus adaptation to an extreme habitat, suggest a broadly conserved genetic basis for trichome variation across flowering plants and pave the way for further investigations of this challenging case of microgeographic incipient speciation.

Exploring Genetic Diversity in Plants Using High-Throughput Sequencing Techniques

Food security has emerged as an urgent concern because of the rising world population. To meet the food demands of the near future, it is required to improve the productivity of various crops, not just of staple food crops. The genetic diversity among plant populations in a given species allows the plants to adapt to various environmental conditions. Such diversity could therefore yield valuable traits that could overcome the food-security challenges. To explore genetic diversity comprehensively and to rapidly identify useful genes and/or allele, advanced high-throughput sequencing techniques, also called next-generation sequencing (NGS) technologies, have been developed. These provide practical solutions to the challenges in crop genomics. Here, we review various sources of genetic diversity in plants, newly developed genetic diversity-mining tools synergized with NGS techniques, and related genetic approaches such as quantitative trait locus analysis and genome-wide association study.

A fragment substitution in the promoter of CsHDZIV11/CsGL3 is responsible for fruit spine density in cucumber (Cucumis sativus L.)

The indel in the promoter of CsHDZIV11 co-segregates with fruit spine density and could be used for molecular breeding in cucumber. Fruit spine density is an important quality trait for marketing in cucumber (Cucumis sativus L.). However, the molecular basis of fruit spine density in cucumber remains unclear. In this study, we isolated a mutant, few spines 1 (fs1), from CNS2 (wild type, WT), a North China-type cucumber with a high density of fruit spines. Genetic analysis showed that fs1 was controlled by a single recessive Mendelian factor. Bulked segregant analysis combined with genome resequencing were used for mapping fs1 in the F2 population derived from a cross between the fs1 mutant and WT, and it was located on chromosome 6 through association analysis. To develop more polymorphic markers to locate fs1, another F2 population was constructed from the cross between fs1 and 'Chinese long' 9930. Then, fs1 was narrowed down to a 110.4-kb genomic region containing 25 annotated genes. A fragment substitution was identified in the promoter region of Csa6M514870 between fs1 and WT. This fragment in fs1 was also present in wild cucumber. Csa6M514870 encodes a PDF2-related protein, a homeodomain-leucine zipper IV transcription factor (CsHDZIV11/CsGL3) sharing high identity and similarity with proteins related to trichome formation or epidermal cell differentiation. Quantitative reverse-transcription PCR revealed a higher expression level of CsHDZIV11 in young fruits from fs1 compared to WT. A molecular marker based on this indel co-segregated with the spine density. This work provides a solid foundation not only for understanding the molecular mechanism of fruit spine density, but also for molecular breeding in cucumber.

The Genetics Underlying Natural Variation in the Biotic Interactions of Arabidopsis thaliana: The Challenges of Linking Evolutionary Genetics and Community Ecology

In the context of global change, predicting the responses of plant communities in an ever-changing biotic environment calls for a multipronged approach at the interface of evolutionary genetics and community ecology. However, our understanding of the genetic basis of natural variation involved in mediating biotic interactions, and associated adaptive dynamics of focal plants in their natural communities, is still in its infancy. Here, we review the genetic and molecular bases of natural variation in the response to biotic interactions (viruses, bacteria, fungi, oomycetes, herbivores, and plants) in the model plant Arabidopsis thaliana as well as the adaptive value of these bases. Among the 60 identified genes are a number that encode nucleotide-binding site leucine-rich repeat (NBS-LRR)-type proteins, consistent with early examples of plant defense genes. However, recent studies have revealed an extensive diversity in the molecular mechanisms of defense. Many types of genetic variants associate with phenotypic variation in biotic interactions, even among the genes of large effect that tend to be identified. In general, we found that (i) balancing selection rather than directional selection explains the observed patterns of genetic diversity within A. thaliana and (ii) the cost/benefit tradeoffs of adaptive alleles can be strongly dependent on both genomic and environmental contexts. Finally, because A. thaliana rarely interacts with only one biotic partner in nature, we highlight the benefit of exploring diffuse biotic interactions rather than tightly associated host-enemy pairs. This challenge would help to improve our understanding of coevolutionary quantitative genetics within the context of realistic community complexity.

Genome-wide identification and characterization of R2R3MYB family in Solanum lycopersicum

The R2R3MYB proteins comprise one of the largest families of transcription factors and play regulatory roles in developmental processes and defense responses in plants. However, there has been relatively little effort to systematically carry out comprehensive genomic and functional analyses of these genes in tomato (Solanum lycopersicum L.), a reference species for Solanaceae plants, and the model plant for fruit development. In this study, a total of 121 R2R3MYB genes were identified in the tomato genome released recently and further classified into 29 subgroups based on the phylogenetic analysis of the complete protein sequences. Phylogenetic comparison of the members of this superfamily among tomato, Arabidopsis, grape, rice, poplar, soybean, cucumber and apple revealed that the putative functions of some tomato R2R3MYB proteins were clustered into the Arabidopsis functional clades. The chromosome distribution pattern revealed that tomato R2R3MYB genes were enriched on several chromosomes and 52 % of the family members were tandemly duplicated genes. Tissue specificity or different expression levels of SlR2R3MYBs in different tissues suggested differential regulation of tissue development as well as metabolic regulation. The transcript abundance level analysis during abiotic conditions identified a group of R2R3MYB genes that responded to one or more treatments suggesting that the SlR2R3MYBs played major roles in the plant response to abiotic conditions and involved in signal transduction pathways. This study not only provides a solid foundation for further functional dissection of tomato R2R3MYB family genes, but may also be profitable for, in the future, the improvement of tomato stress tolerance and fruit quality.

Molecular basis of natural variation and environmental control of trichome patterning

Trichomes are differentiated epidermal cells on above ground organs of nearly all land plants. They play important protective roles as structural defenses upon biotic attacks such as herbivory, oviposition and fungal infections, and against abiotic stressors such as drought, heat, freezing, excess of light, and UV radiation. The pattern and density of trichomes is highly variable within natural population suggesting tradeoffs between traits positively affecting fitness such as resistance and the costs of trichome production. The spatial distribution of trichomes is regulated through a combination of endogenous developmental programs and external signals. This review summarizes the current understanding on the molecular basis of the natural variation and the role of phytohormones and environmental stimuli on trichome patterning.

The genetic architecture of constitutive and induced trichome density in two new recombinant inbred line populations of Arabidopsis thaliana: phenotypic plasticity, epistasis, and bidirectional leaf damage response

BACKGROUND

Herbivory imposes an important selective pressure on plants. In Arabidopsis thaliana leaf trichomes provide a key defense against insect herbivory; however, trichome production incurs a fitness cost in the absence of herbivory. Previous work on A. thaliana has shown an increase in trichome density in response to leaf damage, suggesting a mechanism by which the cost associated with constitutively high trichome density might be mitigated; however, the genetic basis of trichome density induction has not been studied.

RESULTS

Here, we describe the mapping of quantitative trait loci (QTL) for constitutive and damage induced trichome density in two new recombinant inbred line populations of A. thaliana; mapping for constitutive and induced trichome density also allowed for the investigation of damage response (plasticity) QTL. Both novel and previously identified QTL for constitutive trichome density and the first QTL for induced trichome density and response are identified. Interestingly, two of the four parental accessions and multiple RILs in each population exhibited lower trichome density following leaf damage, a response not previously described in A. thaliana. Importantly, a single QTL was mapped for the response phenotype and allelic variation at this locus appears to determine response trajectory in RILs. The data also show that epistatic interactions are a significant component of the genetic architecture of trichome density.

CONCLUSIONS

Together, our results provide further insights into the genetic architecture of constitutive trichome density and new insights into induced trichome density in A. thaliana specifically and to our understanding of the genetic underpinnings of natural variation generally.

Comparison of five major trichome regulatory genes in Brassica villosa with orthologues within the Brassicaceae

Coding sequences for major trichome regulatory genes, including the positive regulators GLABRA 1(GL1), GLABRA 2 (GL2), ENHANCER OF GLABRA 3 (EGL3), and TRANSPARENT TESTA GLABRA 1 (TTG1) and the negative regulator TRIPTYCHON (TRY), were cloned from wild Brassica villosa, which is characterized by dense trichome coverage over most of the plant. Transcript (FPKM) levels from RNA sequencing indicated much higher expression of the GL2 and TTG1 regulatory genes in B. villosa leaves compared with expression levels of GL1 and EGL3 genes in either B. villosa or the reference genome species, glabrous B. oleracea; however, cotyledon TTG1 expression was high in both species. RNA sequencing and Q-PCR also revealed an unusual expression pattern for the negative regulators TRY and CPC, which were much more highly expressed in trichome-rich B. villosa leaves than in glabrous B. oleracea leaves and in glabrous cotyledons from both species. The B. villosa TRY expression pattern also contrasted with TRY expression patterns in two diploid Brassica species, and with the Arabidopsis model for expression of negative regulators of trichome development. Further unique sequence polymorphisms, protein characteristics, and gene evolution studies highlighted specific amino acids in GL1 and GL2 coding sequences that distinguished glabrous species from hairy species and several variants that were specific for each B. villosa gene. Positive selection was observed for GL1 between hairy and non-hairy plants, and as expected the origin of the four expressed positive trichome regulatory genes in B. villosa was predicted to be from B. oleracea. In particular the unpredicted expression patterns for TRY and CPC in B. villosa suggest additional characterization is needed to determine the function of the expanded families of trichome regulatory genes in more complex polyploid species within the Brassicaceae.