Upendra Kumar Devisetty

Comparative transcriptomics reveals patterns of selection in domesticated and wild tomato

Significance One of the most important technological advances by humans is the domestication of plant species for the production of food. We have used high-throughput sequencing to identify changes in DNA sequence and gene expression that differentiate cultivated tomato and its wild relatives. We also identify hundreds of candidate genes that have evolved new protein sequences or have changed expression levels in response to natural selection in wild tomato relatives. Taken together, our analyses provide a snapshot of genome evolution under artificial and natural conditions. Although applied over extremely short timescales, artificial selection has dramatically altered the form, physiology, and life history of cultivated plants. We have used RNAseq to define both gene sequence and expression divergence between cultivated tomato and five related wild species. Based on sequence differences, we detect footprints of positive selection in over 50 genes. We also document thousands of shifts in gene-expression level, many of which resulted from changes in selection pressure. These rapidly evolving genes are commonly associated with environmental response and stress tolerance. The importance of environmental inputs during evolution of gene expression is further highlighted by large-scale alteration of the light response coexpression network between wild and cultivated accessions. Human manipulation of the genome has heavily impacted the tomato transcriptome through directed admixture and by indirectly favoring nonsynonymous over synonymous substitutions. Taken together, our results shed light on the pervasive effects artificial and natural selection have had on the transcriptomes of tomato and its wild relatives.

Accurate timekeeping is controlled by a cycling activator in Arabidopsis

Transcriptional feedback loops are key to circadian clock function in many organisms. Current models of the Arabidopsis circadian network consist of several coupled feedback loops composed almost exclusively of transcriptional repressors. Indeed, a central regulatory mechanism is the repression of evening-phased clock genes via the binding of morning-phased Myb-like repressors to evening element (EE) promoter motifs. We now demonstrate that a related Myb-like protein, REVEILLE8 (RVE8), is a direct transcriptional activator of EE-containing clock and output genes. Loss of RVE8 and its close homologs causes a delay and reduction in levels of evening-phased clock gene transcripts and significant lengthening of clock pace. Our data suggest a substantially revised model of the circadian oscillator, with a clock-regulated activator essential both for clock progression and control of clock outputs. Further, our work suggests that the plant clock consists of a highly interconnected, complex regulatory network rather than of coupled morning and evening feedback loops. DOI: http://dx.doi.org/10.7554/eLife.00473.001

Shade Avoidance Components and Pathways in Adult Plants Revealed by Phenotypic Profiling

Shade from neighboring plants limits light for photosynthesis; as a consequence, plants have a variety of strategies to avoid canopy shade and compete with their neighbors for light. Collectively the response to foliar shade is called the shade avoidance syndrome (SAS). The SAS includes elongation of a variety of organs, acceleration of flowering time, and additional physiological responses, which are seen throughout the plant life cycle. However, current mechanistic knowledge is mainly limited to shade-induced elongation of seedlings. Here we use phenotypic profiling of seedling, leaf, and flowering time traits to untangle complex SAS networks. We used over-representation analysis (ORA) of shade-responsive genes, combined with previous annotation, to logically select 59 known and candidate novel mutants for phenotyping. Our analysis reveals shared and separate pathways for each shade avoidance response. In particular, auxin pathway components were required for shade avoidance responses in hypocotyl, petiole, and flowering time, whereas jasmonic acid pathway components were only required for petiole and flowering time responses. Our phenotypic profiling allowed discovery of seventeen novel shade avoidance mutants. Our results demonstrate that logical selection of mutants increased success of phenotypic profiling to dissect complex traits and discover novel components.

Polymorphism identification and improved genome annotation of Brassica rapa through Deep RNA sequencing.

The mapping and functional analysis of quantitative traits in Brassica rapa can be greatly improved with the availability of physically positioned, gene-based genetic markers and accurate genome annotation. In this study, deep transcriptome RNA sequencing (RNA-Seq) of Brassica rapa was undertaken with two objectives: SNP detection and improved transcriptome annotation. We performed SNP detection on two varieties that are parents of a mapping population to aid in development of a marker system for this population and subsequent development of high-resolution genetic map. An improved Brassica rapa transcriptome was constructed to detect novel transcripts and to improve the current genome annotation. This is useful for accurate mRNA abundance and detection of expression QTL (eQTLs) in mapping populations. Deep RNA-Seq of two Brassica rapa genotypes—R500 (var. trilocularis, Yellow Sarson) and IMB211 (a rapid cycling variety)—using eight different tissues (root, internode, leaf, petiole, apical meristem, floral meristem, silique, and seedling) grown across three different environments (growth chamber, greenhouse and field) and under two different treatments (simulated sun and simulated shade) generated 2.3 billion high-quality Illumina reads. A total of 330,995 SNPs were identified in transcribed regions between the two genotypes with an average frequency of one SNP in every 200 bases. The deep RNA-Seq reassembled Brassica rapa transcriptome identified 44,239 protein-coding genes. Compared with current gene models of B. rapa, we detected 3537 novel transcripts, 23,754 gene models had structural modifications, and 3655 annotated proteins changed. Gaps in the current genome assembly of B. rapa are highlighted by our identification of 780 unmapped transcripts. All the SNPs, annotations, and predicted transcripts can be viewed at http://phytonetworks.ucdavis.edu/.

Modeling development and quantitative trait mapping reveal independent genetic modules for leaf size and shape.

Improved predictions of fitness and yield may be obtained by characterizing the genetic controls and environmental dependencies of organismal ontogeny. Elucidating the shape of growth curves may reveal novel genetic controls that single-time-point (STP) analyses do not because, in theory, infinite numbers of growth curves can result in the same final measurement. We measured leaf lengths and widths in Brassica rapa recombinant inbred lines (RILs) throughout ontogeny. We modeled leaf growth and allometry as function valued traits (FVT), and examined genetic correlations between these traits and aspects of phenology, physiology, circadian rhythms and fitness. We used RNA-seq to construct a SNP linkage map and mapped trait quantitative trait loci (QTL). We found genetic trade-offs between leaf size and growth rate FVT and uncovered differences in genotypic and QTL correlations involving FVT vs STPs. We identified leaf shape (allometry) as a genetic module independent of length and width and identified selection on FVT parameters of development. Leaf shape is associated with venation features that affect desiccation resistance. The genetic independence of leaf shape from other leaf traits may therefore enable crop optimization in leaf shape without negative effects on traits such as size, growth rate, duration or gas exchange.

The RAD51 and DMC1 homoeologous genes of bread wheat: cloning, molecular characterization and expression analysis

BackgroundMeiotic recombination in eukaryotes requires two homologues of the E. coli RecA proteins: Rad51 and Dmc1. Both proteins play important roles in the binding of single stranded DNA, homology search, strand invasion and strand exchange. Meiotic recombination has been well studied in Arabidopsis, rice, maize and the orthologues of RAD51 and DMC1 have been characterized. However genetic analysis of the RAD51 and DMC1 genes in bread wheat has been hampered due to the absence of complete sequence information and because of the existence of multiple copies of each gene in the hexaploid wheat genome.FindingsIn this study we have identified that TaRAD51 and TaDMC1 homoeologues are located on group 7 and group 5 chromosomes of hexaploid wheat, respectively. Comparative sequence analysis of cDNA derived from the TaRAD51 and TaDMC1 homoeologues revealed limited sequence divergence at both the nucleotide and the amino acid level. Indeed, comparisons between the predicted amino acid sequences of TaRAD51 and TaDMC1 and those of other eukaryotes reveal a high degree of evolutionary conservation. Despite the high degree of sequence conservation at the nucleotide level, genome-specific primers for cDNAs of TaRAD51 and TaDMC1 were developed to evaluate expression patterns of individual homoeologues during meiosis. QRT-PCR analysis showed that expression of the TaRAD51 and TaDMC1 cDNA homoeologues was largely restricted to meiotic tissue, with elevated levels observed during the stages of prophase I when meiotic recombination occurs. All three homoeologues of both strand-exchange proteins (TaRAD51 and TaDMC1) are expressed in wheat.ConclusionsBread wheat contains three expressed copies of each of the TaRAD51 and TaDMC1 homoeologues. While differences were detected between the three cDNA homoeologues of TaRAD51 as well as the three homoeologues of TaDMC1, it is unlikely that the predicted amino acid substitutions would have an effect on the protein structure, based on our three-dimensional structure prediction analyses. There are differences in the levels of expression of the three homoeologues of TaRAD51 and TaDMC1 as determined by QRT-PCR and if these differences are reflected at the protein level, bread wheat may be more dependent upon a particular homoeologue to achieve full fertility than all three equally.

Genetic architecture, biochemical underpinnings and ecological impact of floral UV patterning

Floral attraction traits can significantly affect pollinator visitation patterns, but adaptive evolution of these traits may be constrained by correlations with other traits. In some cases, molecular pathways contributing to floral attraction are well characterized, offering the opportunity to explore loci potentially underlying variation among individuals. Here, we quantify the range of variation in floral UV patterning (i.e. UV ‘bulls‐eye nectar guides) among crop and wild accessions of Brassica rapa. We then use experimental crosses to examine the genetic architecture, candidate loci and biochemical underpinnings of this patterning as well as phenotypic manipulations to test the ecological impact. We find qualitative variation in UV patterning between wild (commonly lacking UV patterns) and crop (commonly exhibiting UV patterns) accessions. Similar to the majority of crops, recombinant inbred lines (RILs) derived from an oilseed crop × WI fast‐plant® cross exhibit UV patterns, the size of which varies extensively among genotypes. In RILs, we further observe strong statistical‐genetic and QTL correlations within petal morphological traits and within measurements of petal UV patterning; however, correlations between morphology and UV patterning are weak or nonsignificant, suggesting that UV patterning is regulated and may evolve independently of overall petal size. HPLC analyses reveal a high concentration of sinapoyl glucose in UV‐absorbing petal regions, which, in concert with physical locations of UV‐trait QTLs, suggest a regulatory and structural gene as candidates underlying observed quantitative variation. Finally, insects prefer flowers with UV bulls‐eye patterns over those that lack patterns, validating the importance of UV patterning in pollen‐limited populations of B. rapa.

Using RNA-Seq for Genomic Scaffold Placement, Correcting Assemblies, and Genetic Map Creation in a Common Brassica rapa Mapping Population

Brassica rapa is a model species for agronomic, ecological, evolutionary, and translational studies. Here, we describe high-density SNP discovery and genetic map construction for a B. rapa recombinant inbred line (RIL) population derived from field collected RNA sequencing (RNA-Seq) data. This high-density genotype data enables the detection and correction of putative genome misassemblies and accurate assignment of scaffold sequences to their likely genomic locations. These assembly improvements represent 7.1–8.0% of the annotated B. rapa genome. We demonstrate how using this new resource leads to a significant improvement for QTL analysis over the current low-density genetic map. Improvements are achieved by the increased mapping resolution and by having known genomic coordinates to anchor the markers for candidate gene discovery. These new molecular resources and improvements in the genome annotation will benefit the Brassicaceae genomics community and may help guide other communities in fine-tuning genome annotations.

A genomic analysis of factors driving lincRNA diversification: lessons from plants

Transcriptomic analyses from across eukaryotes indicate that most of the genome is transcribed at some point in the developmental trajectory of an organism. One class of these transcripts is termed long intergenic noncoding RNAs (lincRNAs). Recently, attention has focused on understanding the evolutionary dynamics of lincRNAs, particularly their conservation within genomes. Here, we take a comparative genomic and phylogenetic approach to uncover factors influencing lincRNA emergence and persistence in the plant family Brassicaceae, to which Arabidopsis thaliana belongs. We searched 10 genomes across the family for evidence of > 5000 lincRNA loci from A. thaliana. From loci conserved in the genomes of multiple species, we built alignments and inferred phylogeny. We then used gene tree/species tree reconciliation to examine the duplication history and timing of emergence of these loci. Emergence of lincRNA loci appears to be linked to local duplication events, but, surprisingly, not whole genome duplication events (WGD), or transposable elements. Interestingly, WGD events are associated with the loss of loci for species having undergone relatively recent polyploidy. Lastly, we identify 1180 loci of the 6480 previously annotated A. thaliana lincRNAs (18%) with elevated levels of conservation. These conserved lincRNAs show higher expression, and are enriched for stress-responsiveness and cis-regulatory motifs known as conserved noncoding sequences (CNSs). These data highlight potential functional pathways and suggest that CNSs may regulate neighboring genes at both the genomic and transcriptomic level. In sum, we provide insight into processes that may influence lincRNA diversification by providing an evolutionary context for previously annotated lincRNAs.

Evolinc: A Tool for the Identification and Evolutionary Comparison of Long Intergenic Non-coding RNAs

Long intergenic non-coding RNAs (lincRNAs) are an abundant and functionally diverse class of eukaryotic transcripts. Reported lincRNA repertoires in mammals vary, but are commonly in the thousands to tens of thousands of transcripts, covering ~90% of the genome. In addition to elucidating function, there is particular interest in understanding the origin and evolution of lincRNAs. Aside from mammals, lincRNA populations have been sparsely sampled, precluding evolutionary analyses focused on their emergence and persistence. Here we present Evolinc, a two-module pipeline designed to facilitate lincRNA discovery and characterize aspects of lincRNA evolution. The first module (Evolinc-I) is a lincRNA identification workflow that also facilitates downstream differential expression analysis and genome browser visualization of identified lincRNAs. The second module (Evolinc-II) is a genomic and transcriptomic comparative analysis workflow that determines the phylogenetic depth to which a lincRNA locus is conserved within a user-defined group of related species. Here we validate lincRNA catalogs generated with Evolinc-I against previously annotated Arabidopsis and human lincRNA data. Evolinc-I recapitulated earlier findings and uncovered an additional 70 Arabidopsis and 43 human lincRNAs. We demonstrate the usefulness of Evolinc-II by examining the evolutionary histories of a public dataset of 5,361 Arabidopsis lincRNAs. We used Evolinc-II to winnow this dataset to 40 lincRNAs conserved across species in Brassicaceae. Finally, we show how Evolinc-II can be used to recover the evolutionary history of a known lincRNA, the human telomerase RNA (TERC). These latter analyses revealed unexpected duplication events as well as the loss and subsequent acquisition of a novel TERC locus in the lineage leading to mice and rats. The Evolinc pipeline is currently integrated in CyVerses Discovery Environment and is free for use by researchers.