Development of a mobile, high-throughput, and low-cost image-based plant growth phenotyping system.

Nondestructive plant phenotyping forms a key technique for unraveling molecular processes underlying plant development and response to the environment. While the emergence of high-throughput phenotyping facilities can further our understanding of plant development and stress responses, their high costs greatly hinder scientific progress. To democratize high-throughput plant phenotyping, we developed sets of low-cost image- and weight-based devices to monitor plant shoot growth and evapotranspiration. We paired these devices to a suite of computational pipelines for integrated and straightforward data analysis. The developed tools were validated for their suitability for large genetic screens by evaluating a cowpea (Vigna unguiculata) diversity panel for responses to drought stress. The observed natural variation was used as an input for a genome-wide association study, from which we identified nine genetic loci that might contribute to cowpea drought resilience during early vegetative development. The homologs of the candidate genes were identified in Arabidopsis (Arabidopsis thaliana) and subsequently evaluated for their involvement in drought stress by using available T-DNA insertion mutant lines. These results demonstrate the varied applicability of this low-cost phenotyping system. In the future, we foresee these setups facilitating the identification of genetic components of growth, plant architecture, and stress tolerance across a wide variety of plant species.

Syntrichia ruralis: emerging model moss genome reveals a conserved and previously unknown regulator of desiccation in flowering plants.

Water scarcity, resulting from climate change, poses a significant threat to ecosystems. Syntrichia ruralis, a dryland desiccation-tolerant moss, provides valuable insights into survival of water-limited conditions. We sequenced the genome of S. ruralis, conducted transcriptomic analyses, and performed comparative genomic and transcriptomic analyses with existing genomes and transcriptomes, including with the close relative S. caninervis. We took a genetic approach to characterize the role of an S. ruralis transcription factor, identified in transcriptomic analyses, in Arabidopsis thaliana. The genome was assembled into 12 chromosomes encompassing 21 169 protein-coding genes. Comparative analysis revealed copy number and transcript abundance differences in known desiccation-associated gene families, and highlighted genome-level variation among species that may reflect adaptation to different habitats. A significant number of abscisic acid (ABA)-responsive genes were found to be negatively regulated by a MYB transcription factor (MYB55) that was upstream of the S. ruralis ortholog of ABA-insensitive 3 (ABI3). We determined that this conserved MYB transcription factor, uncharacterized in Arabidopsis, acts as a negative regulator of an ABA-dependent stress response in Arabidopsis. The new genomic resources from this emerging model moss offer novel insights into how plants regulate their responses to water deprivation.

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs.

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

An identical-by-descent segment harbors a 12-bp insertion determining fruit softening during domestication and speciation in Pyrus.

BACKGROUND: Although the wild relatives of pear originated in southwest China, this fruit crop was independently domesticated and improved in Asia and Europe, and there are major phenotypic differences (e.g., maturity and fruit firmness) between Asian and European pears.  RESULTS: In this study, we examined the genomes of 113 diverse pear accessions using an identity-by-descent (IBD) approach to investigate how historical gene flow has shaped fruit firmness traits in Asian and European pears. We found a 3-Mbp IBD-enriched region (IBD-ER) that has undergone "convergent domestication" in both the Asian and European pear lineages, and a genome-wide association study (GWAS) of fruit firmness phenotypes strongly implicated the TRANSLOCON AT THE INNER CHLOROPLAST ENVELOPE55 (TIC55) locus within this 3-Mbp IBD-ER. Furthermore, we identified a tandem duplication that includes a 12-bp insertion located in the first exon of TIC55 that is uniquely present in Asian pears, and expression analysis showed that the pear TIC55 gene is highly expressed in Asian pear, while it is weakly or not expressed in European pear; this could contribute to the differences in fruit firmness between Asian and European pear fruits. CONCLUSIONS: Our findings provide insights into how pear fruit softening has been impacted during domestication, and we identified candidate genes associated with fruit softening that can contribute to the breeding and improvement of pear and other fruit crops.

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae.

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

The spinach YY genome reveals sex chromosome evolution, domestication, and introgression history of the species.

BACKGROUND: Spinach (Spinacia oleracea L.) is a dioecious species with an XY sex chromosome system, but its Y chromosome has not been fully characterized. Our knowledge about the history of its domestication and improvement remains limited. RESULTS: A high-quality YY genome of spinach is assembled into 952 Mb in six pseudo-chromosomes. By a combination of genetic mapping, Genome-Wide Association Studies, and genomic analysis, we characterize a 17.42-Mb sex determination region (SDR) on chromosome 1. The sex chromosomes of spinach evolved when an insertion containing sex determination genes occurred, followed by a large genomic inversion about 1.98 Mya. A subsequent burst of SDR-specific repeats (0.1-0.15 Mya) explains the large size of this SDR. We identify a Y-specific gene, NRT1/PTR 6.4 which resides in this insertion, as a strong candidate for the sex determination or differentiation factor. Resequencing of 112 spinach genomes reveals a severe domestication bottleneck approximately 10.87 Kya, which dates the domestication of spinach 7000 years earlier than the archeological record. We demonstrate that a strong selection signal associated with internode elongation and leaf area expansion is associated with domestication of edibility traits in spinach. We find that several strong genomic introgressions from the wild species Spinacia turkestanica and Spinacia tetrandra harbor desirable alleles of genes related to downy mildew resistance, frost resistance, leaf morphology, and flowering-time shift, which likely contribute to spinach improvement. CONCLUSIONS: Analysis of the YY genome uncovers evolutionary forces shaping nascent sex chromosome evolution in spinach. Our findings provide novel insights about the domestication and improvement of spinach.

Identification of structural variation and polymorphisms of a sex co-segregating scaffold in spinach.

Spinach is a common vegetable, and dioecy is maintained by a pair of XY sex chromosomes. Due to limited genomic resources and its highly repetitive genome, limited studies were conducted to investigate the genomic landscape of the region near sex-determining loci. In this study, we screened the structure variations (SVs) between Y-linked contigs and a 1.78-Mb X scaffold (Super_scaffold 66), which enabled the development of 12 sex co-segregating DNA markers. These markers were tested in one F1 mapping population and 40 spinach accessions, which comprised 692 individual plants with the strong sex linkage pattern. In addition, we found that Super_scaffold 66 was highly repetitive along with the enriched LTR-RTs insertions and decreased microsatellite distribution compared with the rest genome, which matches extremely low gene density featured by only nine annotated genes. Synteny analysis between Y contigs and Superscaffold_66 revealed a 340-Kb accumulative Y contig (non-continuous) and a 500-Kb X counterpart along with SVs and wide-spread tandem duplications. Among the nine genes, one ABC transporter gene revealed noticeable SVs between Y contig and X counterpart, as an approximate 5-Kb recent Gypsy LTR-RT insertion in the Y-linked allele, but not the X allele. The gene paucity, SVs, and sex-linked polymorphisms attributed to the recombination suppression. We proposed that Super_scaffold 66 is part of the non-recombining region containing the sex determination genes. The spread of 12 sex co-segregating markers from this 1.78 Mb genomic region indicated the existence and expansion of sex determination region during progression of the Y chromosome.

Cloning and physical localization of male-biased repetitive DNA sequences in Spinacia oleracea (Amaranthaceae).

Spinach (Spinacia oleracea Linnaeus, 1753) is an ideal material for studying molecular mechanisms of early-stage sex chromosome evolution in dioecious plants. Degenerate oligonucleotide-primed polymerase chain reaction (DOP-PCR) technique facilitates the retrotransposon-relevant studies by enriching specific repetitive DNA sequences from a micro-dissected single chromosome. We conducted genomic subtractive hybridization to screen sex-biased DNA sequences by using the DOP-PCR amplification products of micro-dissected spinach Y chromosome. The screening yielded 55 male-biased DNA sequences with 30 576 bp in length, of which, 32 DNA sequences (12 049 bp) contained repeat DNA sequences, including LTR/Copia, LTR/Gypsy, simple repeats, and DNA/CMC-EnSpm. Among these repetitive DNA sequences, four DNA sequences that contained a fragment of Ty3-gypsy retrotransposons (SP73, SP75, SP76, and SP77) were selected as fluorescence probes to hybridization on male and female spinach karyotypes. Fluorescence in situ hybridization (FISH) signals of SP73 and SP75 were captured mostly on the centromeres and their surrounding area for each homolog. Hybridization signals primarily appeared near the putative centromeres for each homologous chromosome pair by using SP76 and SP77 probes for FISH, and sporadic signals existed on the long arms. Results can be served as a basis to study the function of repetitive DNA sequences in sex chromosome evolution in spinach.

Construction of high-density genetic maps defined sex determination region of the Y chromosome in spinach.

Spinach (Spinacia olracea L.) is a dioecious leafy vegetable with a highly repetitive genome of around 990 Mb, which is challenging for de-novo genome assembly. In our study, a segregating F1 (double pseudo-testcross) population from 'Viroflay' × 'Cornell-NO. 9' was used for genetic mapping by resequencing genotyping. In the paternal 'Cornell-NO. 9' map, 212,414 SNPs were mapped, and the total linkage distance was 476.83 cM; the maternal 'Viroflay' map included 29,282 SNPs with 401.28 cM total genetic distance. Both paternal and maternal maps have the expected number of six linkage groups (LGs). A non-recombining region with 5678 SNPs (39 bin markers) co-segregates with sex type which located at 45.2 cM of LG1 in the 'Cornell-NO. 9' map while indicates the sex determination region (SDR). Integration of two maps into a consensus map guided us to anchor additional 1242 contigs to six pseudomolecules from the published reference genome, which improved additional 233 Mb (23.4%) assembly based on spinach estimated genome size. Particularly, the X counterpart of SDR in our assembly is estimated around 18.4 Mb which locates at the largest chromosome, as consensus with sex-biased FISH signals from previous cytogenetics studies. The region is featured by reduced gene density, higher percentage of repetitive sequences, and no recombination. Our linkage maps provide the resource for improving spinach genome de-novo assembly and identification of sex-determining genes in spinach.

Identification of Candidate Auxin Response Factors Involved in Pomegranate Seed Coat Development.

Auxin response factors (ARFs) are transcription factors, regulating the auxin signaling pathways involved in plant development and related processes. In this study, we performed the genome-wide identification and characterization of ARFs in pomegranate and compared them with ARFs from three other species. Seventeen PgrARFs were identified and clustered into four groups, according to their phylogenetic relationship with the remaining 59 ARFs. A recent whole-genome duplication event in pomegranate may have contributed to the expansion and diversification of PgrARFs. Genomic truncation and variant splicing mechanisms contributed to the divergence of PgrARFs, a conclusion that was supported by different exon-intron structures of genes and incomplete conserved domains of PgrARFs in a specific phylogenetic group (group III). Interestingly, the absence of motifs from certain PgrARF genes corresponded to their low transcription levels, which contrasted to the highly expressed PgrARFs with intact motifs. Specifically, PgrARF1 and PgrARF2 highly expressed in both inner and outer seed coat, and phylogenetically related to Arabidopsis orthologs which mediates cell divisions in seed coat. We infer these two PgrARFs might involve in seed coat development through cell divisions in response to auxin regulation. These findings provided information on the characteristics and evolutionary relationships of PgrARFs, but also shed lights on their potential roles during seed coat development in pomegranate.

Genome-Wide Analysis of Transposable Elements and Satellite DNAs in Spinacia Species to Shed Light on Their Roles in Sex Chromosome Evolution.

Sex chromosome evolution has mostly been studied in species with heteromorphic sex chromosomes. The Spinacia genus serves as an ideal model for investigating evolutionary mechanisms underlying the transition from homomorphic to heteromorphic sex chromosomes. Among evolutionary factors, repetitive sequences play multiple roles in sex chromosome evolution while their forces have not been fully explored in Spinacia species. Here, we identified major repetitive sequence classes in male and female genomes of Spinacia species and their ancestral relative sugar beet to elucidate the evolutionary processes of sex chromosome evolution using next-generation sequencing (NGS) data. Comparative analysis revealed that the repeat elements of Spinacia species are considerably higher than of sugar beet, especially the Ty3/Gypsy and Ty1/Copia retrotransposons. The long terminal repeat retroelements (LTR) Angela, Athila, and Ogre may be accounted for the higher proportion of repeats in the spinach genome. Comparison of the repeats proportion between female and male genomes of three Spinacia species indicated the different representation in Spinacia tetrandra samples but not in the S. oleracea or S. turkestanica samples. From these results, we speculated that emergence of repetitive DNA sequences may correlate the formation of sex chromosome and the transition from homomorphic sex chromosomes to heteromorphic sex chromosomes as heteromorphic sex chromosomes exclusively existed in Spinacia tetrandra. Three novel sugar beet-specific satellites were identified and confirmed by fluorescence in situ hybridization (FISH); six out of eight new spinach-specific satellites were mapped to the short arm of sex chromosomes. A total of 141 copies of SolSat01-171-s were found in the sex determination region (SDR). Thus, the accumulation of satellite DNA on the short arm of chromosome 1 may be involved in the sex chromosome evolution in Spinacia species. Our study provides a fundamental resource for understanding repeat sequences in Spinacia species and their roles in sex chromosome evolution.

Differential gene expression among three sex types reveals a MALE STERILITY 1 (CpMS1) for sex differentiation in papaya.

BACKGROUND: Carica papaya is a trioecious plant species with a genetic sex-determination system defined by sex chromosomes. Under unfavorable environmental conditions male and hermaphrodite exhibit sex-reversal. Previous genomic research revealed few candidate genes for sex differentiation in this species. Nevertheless, more analysis is still needed to identify the mechanism responsible for sex flower organ development in papaya. RESULTS: The aim of this study was to identify differentially expressed genes among male, female and hermaphrodite flowers in papaya during early (pre-meiosis) and later (post-meiosis) stages of flower development. RNA-seq was used to evaluate the expression of differentially expressed genes and RT-qPCR was used to verify the results. Putative functions of these genes were analyzed based on their homology with orthologs in other plant species and their expression patterns. We identified a Male Sterility 1 gene (CpMS1) highly up-regulated in male and hermaphrodite flower buds compared to female flower buds, which expresses in small male flower buds (3-8 mm), and that might be playing an important role in male flower organ development due to its homology to MS1 genes previously identified in other plants. This is the first study in which the sex-biased expression of genes related to tapetum development in the anther developmental pathway is being reported in papaya. Besides important transcription factors related to flower organ development and flowering time regulation, we identified differential expression of genes that are known to participate in ABA, ROS and auxin signaling pathways (ABA-8-hydroxylases, AIL5, UPBEAT 1, VAN3-binding protein). CONCLUSIONS: CpMS1 was expressed in papaya male and hermaphrodite flowers at early stages, suggesting that this gene might participate in male flower organ development processes, nevertheless, this gene cannot be considered a sex-determination gene. Due to its homology with other plant MS1 proteins and its expression pattern, we hypothesize that this gene participates in anther development processes, like tapetum and pollen development, downstream gender specification. Further gene functional characterization studies in papaya are required to confirm this hypothesis. The role of ABA and ROS signaling pathways in papaya flower development needs to be further explored as well.

Genome-Wide Identification, Evolution and Functional Divergence of MYB Transcription Factors in Chinese White Pear (Pyrus bretschneideri).

The MYB superfamily is large and functionally diverse in plants. To date, MYB family genes have not yet been identified in Chinese white pear (Pyrus bretschneideri), and their functions remain unclear. In this study, we identified 231 genes as candidate MYB genes and divided them into four subfamilies. The R2R3-MYB (PbrMYB) family shared an R2R3 domain with 104 amino acid residues, including five conserved tryptophan residues. The Pbr MYB family was divided into 37 functional subgroups including 33 subgroups which contained both MYB genes of Rosaceae plants and AtMYB genes, and four subgroups which included only Rosaceae MYB genes or AtMYB genes. PbrMYB genes with similar functions clustered into the same subgroup, indicating functional conservation. We also found that whole-genome duplication (WGD) and dispersed duplications played critical roles in the expansion of the MYB family. The 87 Pbr MYB duplicated gene pairs dated back to the two WGD events. Purifying selection was the primary force driving Pbr MYB gene evolution. The 15 gene pairs presented 1-7 codon sites under positive selection. A total of 147 expressed genes were identified from RNA-sequencing data of fruit, and six Pbr MYB members in subgroup C1 were identified as important candidate genes in the regulation of lignin synthesis by quantitative real-time PCR analysis. Further correlation analysis revealed that six PbrMYBs were significantly correlated with five structural gene families (F5H, HCT, CCR, POD and C3'H) in the lignin pathway. The phylogenetic, evolution and expression analyses of the MYB gene family in Chinese white pear establish a solid foundation for future comprehensive functional analysis of Pbr MYB genes.