The Arabidopsis Plant Intracellular Ras-group LRR (PIRL) Family and the Value of Reverse Genetic Analysis for Identifying Genes that Function in Gametophyte Development

Expression analysis of genes encoding extracellular leucine-rich repeat proteins in Arabidopsis thaliana

Leucine-rich repeat (LRR)-containing proteins have been identified in diverse species, including plants. The diverse intracellular and extracellular LRR variants are responsible for numerous biological processes. We analyzed the expression patterns of Arabidopsis thaliana extracellular LRR (AtExLRR) genes, 10 receptor-like proteins, and 4 additional genes expressing the LRR-containing protein by a promoter: β-glucuronidase (GUS) study. According to in silico expression studies, several AtExLRR genes were expressed in a tissue- or stage-specific and abiotic/hormone stress-responsive manner, indicating their potential participation in specific biological processes. Based on the promoter: GUS assay, AtExLRRs were expressed in different cells and organs. A quantitative real-time PCR investigation revealed that the expressions of AtExLRR3 and AtExLRR9 were distinct under various abiotic stress conditions. This study investigated the potential roles of extracellular LRR proteins in plant growth, development, and response to various abiotic stresses.

Integrated "-omics" analysis highlights the role of brassinosteroid signaling and antioxidant machinery underlying improved rice seed longevity during artificial aging treatment

Seed longevity is a critical characteristic in agriculture, yet the specific genes/proteins responsible for this trait and the molecular mechanisms underlying reduced longevity during seed aging remain largely elusive. Here we report the comparative proteome and metabolome profiling of three rice cultivars exhibiting varying degrees of aging tolerance: Dharial, an aging-tolerant cultivar; Ilmi, an aging-sensitive cultivar; and A2, a moderately aging-tolerant cultivar developed from the crossbreeding of Dharial and Ilmi. Artificial aging treatment (AAT) markedly reduced the germination percentage and enhanced the activities of antioxidant enzymes in all the cultivars. Further, proteomics results showed a key role of the ubiquitin (Ub)-proteasome pathway in the degradation of damaged proteins during AAT while other proteases were majorly reduced. In addition, proteins associated with energy production and protein synthesis were strongly reduced in Ilmi while these were majorly increased in A2 and Dharial. These, along with metabolomics results, suggest that Ub-proteasome mediated protein degradation during AAT results in the accumulation of free amino acids in Ilmi while tolerant cultivars potentially utilize those for energy production and synthesis of stress-related proteins, especially hsp20/alpha-crystallin family protein. Additionally, both Dharial and A2 seem to activate brassinosteroid signaling and suppress jasmonate signaling which initiates a signaling cascade that allows accumulation of enzymatic and non-enzymatic antioxidants for efficient detoxification of aging-induced ROS. Taken together, these results provide an in-depth understanding of the aging-induced changes in rice seeds and highlight key pathways responsible for maintaining seed longevity during AAT.

Fine mapping and candidate gene analysis of qFL-A12-5: a fiber length-related QTL introgressed from Gossypium barbadense into Gossypium hirsutum

The fiber length-related qFL-A12-5 identified in CSSLs introgressed from Gossypium barbadense into Gossypium hirsutum was fine-mapped to an 18.8 kb region on chromosome A12, leading to the identification of the GhTPR gene as a potential regulator of cotton fiber length. Fiber length is a key determinant of fiber quality in cotton, and it is a key target of artificial selection for breeding and domestication. Although many fiber length-related quantitative trait loci have been identified, there are few reports on their fine mapping or candidate gene validation, thus hampering efforts to understand the mechanistic basis of cotton fiber development. Our previous study identified the qFL-A12-5 associated with superior fiber quality on chromosome A12 in the chromosome segment substitution line (CSSL) MBI7747 (BC₄F_3:5). A single segment substitution line (CSSL-106) screened from BC₆F₂ was backcrossed to construct a larger segregation population with its recurrent parent CCRI45, thus enabling the fine mapping of 2852 BC₇F₂ individuals using denser simple sequence repeat markers to narrow the qFL-A12-5 to an 18.8 kb region of the genome, in which six annotated genes were identified in Gossypium hirsutum. Quantitative real-time PCR and comparative analyses led to the identification of GH_A12G2192 (GhTPR) encoding a tetratricopeptide repeat-like superfamily protein as a promising candidate gene for qFL-A12-5. A comparative analysis of the protein-coding regions of GhTPR among Hai1, MBI7747, and CCRI45 revealed two non-synonymous mutations. The overexpression of GhTPR resulted in longer roots in Arabidopsis, suggesting that GhTPR may regulate cotton fiber development. These results provide a foundation for future efforts to improve cotton fiber length.

Genetic Adaptation of Siberian Larch (Larix sibirica Ledeb.) to High Altitudes

Forest trees growing in high altitude conditions offer a convenient model for studying adaptation processes. They are subject to a whole range of adverse factors that are likely to cause local adaptation and related genetic changes. Siberian larch (Larix sibirica Ledeb.), whose distribution covers different altitudes, makes it possible to directly compare lowland with highland populations. This paper presents for the first time the results of studying the genetic differentiation of Siberian larch populations, presumably associated with adaptation to the altitudinal gradient of climatic conditions, based on a joint analysis of altitude and six other bioclimatic variables, together with a large number of genetic markers, single nucleotide polymorphisms (SNPs), obtained from double digest restriction-site-associated DNA sequencing (ddRADseq). In total, 25,143 SNPs were genotyped in 231 trees. In addition, a dataset of 761 supposedly selectively neutral SNPs was assembled by selecting SNPs located outside coding regions in the Siberian larch genome and mapped to different contigs. The analysis using four different methods (PCAdapt, LFMM, BayeScEnv and RDA) revealed 550 outlier SNPs, including 207 SNPs whose variation was significantly correlated with the variation of some of environmental factors and presumably associated with local adaptation, including 67 SNPs that correlated with altitude based on either LFMM or BayeScEnv and 23 SNPs based on both of them. Twenty SNPs were found in the coding regions of genes, and 16 of them represented non-synonymous nucleotide substitutions. They are located in genes involved in the processes of macromolecular cell metabolism and organic biosynthesis associated with reproduction and development, as well as organismal response to stress. Among these 20 SNPs, nine were possibly associated with altitude, but only one of them was identified as associated with altitude by all four methods used in the study, a nonsynonymous SNP in scaffold_31130 in position 28092, a gene encoding a cell membrane protein with uncertain function. Among the studied populations, at least two main groups (clusters), the Altai populations and all others, were significantly genetically different according to the admixture analysis based on any of the three SNP datasets as follows: 761 supposedly selectively neutral SNPs, all 25,143 SNPs and 550 adaptive SNPs. In general, according to the AMOVA results, genetic differentiation between transects or regions or between population samples was relatively low, although statistically significant, based on 761 neutral SNPs (F_ST = 0.036) and all 25,143 SNPs (F_ST = 0.017). Meanwhile, the differentiation based on 550 adaptive SNPs was much higher (F_ST = 0.218). The data showed a relatively weak but highly significant linear correlation between genetic and geographic distances (r = 0.206, p = 0.001).

Expression analysis of plant intracellular Ras-group related leucine-rich repeat proteins (PIRLs) in Arabidopsis thaliana

Arabidopsis thaliana contains a family of nine genes known as plant intracellular Ras-group related leucine-rich repeat (LRR) proteins (PIRLs). These are structurally similar to animals and fungal LRR proteins and play important roles in developmental pathways. However, to date, no detailed tissue-specific expression analysis of these PIRLs has been performed. Therefore, in this study, we generated promoter:GUS transgenic plants for the nine A. thaliana PIRL genes and identified their expression patterns in seedlings and floral organs at different developmental stages. Most PIRL members showed expression in the root apical region and in the vascular tissue of primary and lateral roots. Shoot apex-specific expression was recorded for PIRL1 and PIRL8. Furthermore, PIRL1, PIRL3, PIRL5, PIRL6, and PIRL7 showed distinct expression patterns in flowers, especially in pollen and anthers. In addition, co-expression network analysis identified cases where PIRLs were co-expressed with other genes known to have specific functions related to growth and development. Taken together, the tissue-specific expression patterns of PIRL genes improve our understanding of the functions of this gene family in plant growth and development.

•

PIRL constituting gene family in A. thaliana is widely distributed among plants.

•

PIRL1, 2, 3, and 6 were strongly expressed in anther and/or pollen, consistent with their function in pollen development.

•

PIRL7 was distinctively expressed in pollen and pollen tube, suggesting its role in pollen function.

The molecular signatures of compatible and incompatible pollination in Arabidopsis

Background

Fertilization in flowering plants depends on the early contact and acceptance of pollen grains by the receptive papilla cells of the stigma. Deciphering the specific transcriptomic response of both pollen and stigmatic cells during their interaction constitutes an important challenge to better our understanding of this cell recognition event.

Results

Here we describe a transcriptomic analysis based on single nucleotide polymorphisms (SNPs) present in two Arabidopsis thaliana accessions, one used as female and the other as male. This strategy allowed us to distinguish 80% of transcripts according to their parental origins. We also developed a tool which predicts male/female specific expression for genes without SNP. We report an unanticipated transcriptional activity triggered in stigma upon incompatible pollination and show that following compatible interaction, components of the pattern-triggered immunity (PTI) pathway are induced on the female side.

Conclusions

Our work unveils the molecular signatures of compatible and incompatible pollinations both at the male and female side. We provide invaluable resource and tools to identify potential new molecular players involved in pollen-stigma interaction.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07503-7.

Construction of the First SNP-Based Linkage Map Using Genotyping-by-Sequencing and Mapping of the Male-Sterility Gene in Leaf Chicory

We report the first high-density linkage map construction through genotyping-by-sequencing (GBS) in leaf chicory (Cichorium intybus subsp. intybus var. foliosum, 2n = 2x = 18) and the SNP-based fine mapping of the linkage group region carrying a recessive gene responsible for male-sterility (ms1). An experimental BC₁ population, segregating for the male sterility trait, was specifically generated and 198 progeny plants were preliminary screened through a multiplexed SSR genotyping analysis for the identification of microsatellite markers linked to the ms1 locus. Two backbone SSR markers belonging to linkage group 4 of the available Cichorium consensus map were found genetically associated to the ms1 gene at 5.8 and 12.1 cM apart. A GBS strategy was then used to produce a high-density SNP-based linkage map, containing 727 genomic loci organized into 9 linkage groups and spanning a total length of 1,413 cM. 13 SNPs proved to be tightly linked to the ms1 locus based on a subset of 44 progeny plants analyzed. The map position of these markers was further validated by sequence-specific PCR experiments using an additional set of 64 progeny plants, enabling to verify that four of them fully co-segregated with male-sterility. A mesosynteny analysis revealed that 10 genomic DNA sequences encompassing the 13 selected SNPs of chicory mapped in a peripheral region of chromosome 5 of lettuce (Lactuca sativa L.) spanning about 18 Mbp. Since a MYB103-like gene, encoding for a transcription factor involved in callose dissolution of tetrads and exine development of microspores, was found located in the same chromosomal region, this orthologous was chosen as candidate for male-sterility. The amplification and sequencing of its CDS using accessions with contrasting phenotypes/genotypes (i.e., 4 male sterile mutants, ms1ms1, and 4 male fertile inbreds, Ms1Ms1) enabled to detect an INDEL of 4 nucleotides in its second exon, responsible for an anticipated stop codon in the male sterile mutants. This polymorphism was subsequently validated through allele-specific PCR assays and found to fully co-segregate with male-sterility, using 64 progeny plants of the same mapping BC₁ population. Overall, our molecular data could be practically exploited for genotyping plant materials and for marker-assisted breeding schemes in leaf chicory.

Arabidopsis PIRL6 Is Essential for Male and Female Gametogenesis and Is Regulated by Alternative Splicing

Plant intracellular Ras-group leucine-rich repeat (LRR) proteins (PIRLs) are related to Ras-interacting animal LRR proteins that participate in developmental cell signaling. Systematic knockout analysis has implicated some members of the Arabidopsis (Arabidopsis thaliana) PIRL family in pollen development. However, for PIRL6, no bona fide knockout alleles have been recovered, suggesting that it may have an essential function in both male and female gametophytes. To test this hypothesis, we investigated PIRL6 expression and induced knockdown by RNA interference. Knockdown triggered defects in gametogenesis, resulting in abnormal pollen and early developmental arrest in the embryo sac. Consistent with this, PIRL6 was expressed in gametophytes: functional transcripts were detected in wild-type flowers but not in sporocyteless (spl) mutant flowers, which do not produce gametophytes. A genomic PIRL6-GFP fusion construct confirmed expression in both pollen and the embryo sac. Interestingly, PIRL6 is part of a convergent overlapping gene pair, a scenario associated with an increased likelihood of alternative splicing. We detected multiple alternative PIRL6 mRNAs in vegetative organs and spl mutant flowers, tissues that lacked the functionally spliced transcript. cDNA sequencing revealed that all contained intron sequences and premature termination codons. These alternative mRNAs accumulated in the nonsense-mediated decay mutant upf3, indicating that they are normally subjected to degradation. Together, these results demonstrate that PIRL6 is required in both male and female gametogenesis and suggest that sporophytic expression is negatively regulated by unproductive alternative splicing. This posttranscriptional mechanism may function to minimize PIRL6 protein expression in sporophyte tissues while allowing the overlapping adjacent gene to remain widely transcribed.

Large-scale identification of wheat genes resistant to cereal cyst nematode Heterodera avenae using comparative transcriptomic analysis

Background

Cereal cyst nematode Heterodera avenae, an important soil-borne pathogen in wheat, causes numerous annual yield losses worldwide, and use of resistant cultivars is the best strategy for control. However, target genes are not readily available for breeding resistant cultivars. Therefore, comparative transcriptomic analyses were performed to identify more applicable resistance genes for cultivar breeding.

Methods

The developing nematodes within roots were stained with acid fuchsin solution. Transcriptome assemblies and redundancy filteration were obtained by Trinity, TGI Clustering Tool and BLASTN, respectively. Gene Ontology annotation was yielded by Blast2GO program, and metabolic pathways of transcripts were analyzed by Path_finder. The ROS levels were determined by luminol-chemiluminescence assay. The transcriptional gene expression profiles were obtained by quantitative RT-PCR.

Results

The RNA-sequencing was performed using an incompatible wheat cultivar VP1620 and a compatible control cultivar WEN19 infected with H. avenae at 24 h, 3 d and 8 d. Infection assays showed that VP1620 failed to block penetration of H. avenae but disturbed the transition of developmental stages, leading to a significant reduction in cyst formation. Two types of expression profiles were established to predict candidate resistance genes after developing a novel strategy to generate clean RNA-seq data by removing the transcripts of H. avenae within the raw data before assembly. Using the uncoordinated expression profiles with transcript abundance as a standard, 424 candidate resistance genes were identified, including 302 overlapping genes and 122 VP1620-specific genes. Genes with similar expression patterns were further classified according to the scales of changed transcript abundances, and 182 genes were rescued as supplementary candidate resistance genes. Functional characterizations revealed that diverse defense-related pathways were responsible for wheat resistance against H. avenae. Moreover, phospholipase was involved in many defense-related pathways and localized in the connection position. Furthermore, strong bursts of reactive oxygen species (ROS) within VP1620 roots infected with H. avenae were induced at 24 h and 3 d, and eight ROS-producing genes were significantly upregulated, including three class III peroxidase and five lipoxygenase genes.

Conclusions

Large-scale identification of wheat resistance genes were processed by comparative transcriptomic analysis. Functional characterization showed that phospholipases associated with ROS production played vital roles in early defense responses to H. avenae via involvement in diverse defense-related pathways as a hub switch. This study is the first to investigate the early defense responses of wheat against H. avenae, not only provides applicable candidate resistance genes for breeding novel wheat cultivars, but also enables a better understanding of the defense mechanisms of wheat against H. avenae.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2037-8) contains supplementary material, which is available to authorized users.