Comparative transcriptomics and genomic patterns of discordance in Capsiceae (Solanaceae).

The integration of genomics and phylogenetics allows new insight into the structure of gene tree discordance, the relationships among gene position, gene history, and rate of evolution, as well as the correspondence of gene function, positive selection, and gene ontology enrichment across lineages. We explore these issues using the tribe Capsiceae (Solanaceae), which is comprised of the genera Lycianthes and Capsicum (peppers). In combining the annotated genomes of Capsicum with newly sequenced transcriptomes of four species of Lycianthes and Capsicum, we develop phylogenies for 6747 genes, and construct a backbone species tree using both concordance and explicit phylogenetic network approaches. We quantify phylogenetic discordance among individual gene trees, measure their rates of synonymous and nonsynonymous substitution, and test whether they were positively selected along any branch of the phylogeny. We then map these genes onto the annotated Capsicum genome and test whether rates of evolution, gene history, and gene ontology vary significantly with gene position. We observed substantial discordance among gene trees. A bifurcating species tree placing Capsicum within a paraphyletic Lycianthes was supported over all phylogenetic networks. Rates of synonymous and nonsynonymous substitution varied 41-fold and 130-fold among genes, respectively, and were significantly lower in pericentromeric regions. We found that results of concordance tree analyses vary depending on the subset of genes used, and that genes within the pericentromeric regions only capture a portion of the observed discordance. We identified 787 genes that have been positively selected throughout the diversification history of Capsiceae, and discuss the importance of these genes as targets for investigation of economically important traits in the domesticated peppers.

Identification of QTLs for capsaicinoids, fruit quality, and plant architecture-related traits in an interspecific Capsicum RIL population.

Quantitative trait loci (QTL) analyses in pepper are common for horticultural, disease resistance, and fruit quality traits; although none of the studies to date have used sequence-based markers associated with genes. In this study we measured plant architectural, phenological, and fruit quality traits in a pepper mapping population consisting of 92 recombinant inbred lines derived from a cross between Capsicum frutescens acc. 2814-6 and C. annuum var. NuMexRNAKY. Phenotypic measurements were correlated to loci in a high-density EST-based genetic map. In total, 96 QTL were identified for 38 traits, including 12 QTL associated with capsaicinoid levels. Twenty-one loci showed correlation among seemingly unrelated phenotypic categories, highlighting tight linkage or shared genetics between previously unassociated traits in pepper.

Reference quality assembly of the 3.5-Gb genome of Capsicum annuum from a single linked-read library.

Linked-Read sequencing technology has recently been employed successfully for de novo assembly of human genomes, however, the utility of this technology for complex plant genomes is unproven. We evaluated the technology for this purpose by sequencing the 3.5-gigabase (Gb) diploid pepper (Capsicum annuum) genome with a single Linked-Read library. Plant genomes, including pepper, are characterized by long, highly similar repetitive sequences. Accordingly, significant effort is used to ensure that the sequenced plant is highly homozygous and the resulting assembly is a haploid consensus. With a phased assembly approach, we targeted a heterozygous F1 derived from a wide cross to assess the ability to derive both haplotypes and characterize a pungency gene with a large insertion/deletion. The Supernova software generated a highly ordered, more contiguous sequence assembly than all currently available C. annuum reference genomes. Over 83% of the final assembly was anchored and oriented using four publicly available de novo linkage maps. A comparison of the annotation of conserved eukaryotic genes indicated the completeness of assembly. The validity of the phased assembly is further demonstrated with the complete recovery of both 2.5-Kb insertion/deletion haplotypes of the PUN1 locus in the F1 sample that represents pungent and nonpungent peppers, as well as nearly full recovery of the BUSCO2 gene set within each of the two haplotypes. The most contiguous pepper genome assembly to date has been generated which demonstrates that Linked-Read library technology provides a tool to de novo assemble complex highly repetitive heterozygous plant genomes. This technology can provide an opportunity to cost-effectively develop high-quality genome assemblies for other complex plants and compare structural and gene differences through accurate haplotype reconstruction.

Arabidopsis SHORT INTEGUMENTS 2 is a mitochondrial DAR GTPase.

The Arabidopsis short integuments 2-1 (sin2-1) mutant produces ovules with short integuments due to early cessation of cell division in these structures. SIN2 was isolated and encodes a putative GTPase sharing features found in the novel DAR GTPase family. DAR proteins share a signature DAR motif and a unique arrangement of the four conserved GTPase G motifs. We found that DAR GTPases are present in all examined prokaryotes and eukaryotes and that they have diversified into four paralogous lineages in higher eukaryotes. Eukaryotic members of the SIN2 clade of DAR GTPases have been found to localize to mitochondria and are related to eubacterial proteins that facilitate essential steps in biogenesis of the large ribosomal subunit. We propose a similar role for SIN2 in mitochondria. A sin2 insertional allele has ovule effects similar to sin2-1, but more pronounced pleiotropic effects on vegetative and floral development. The diverse developmental effects of the mitochondrial SIN2 GTPase support a mitochondrial role in the regulation of multiple developmental pathways.

Regions Underlying Population Structure and the Genomics of Organ Size Determination in Capsicum annuum.

Fruits, as an important part of the human diet, have been under strong selection during domestication. In general, continued directed selection has led to varieties having larger fruit with greater shape variation and tremendous increases in fruit mass. Common cultivated peppers ( L.) are found in a wide range of sizes and shapes. Analysis of genetic relatedness and population structure has shown that the large-fruited, nonpungent types have reduced diversity and comprise a highly structured group. To explore this population structure, a statistical method for detecting fixation within subpopulations was applied to a set of 21 pungent and 19 nonpungent lines that represent the pepper breeding germplasm. We have identified 17 blocks within the pepper genome that are conserved among nonpungent large-fruited varieties. To determine if these regions were fixed by selection on fruit size or pungency, quantitative trait loci (QTLs) from seven studies along with capsaicin biosynthesis genes and homologs of organ size regulatory genes were mapped onto the current pepper genome assembly. Of the 17 fixed regions, 14 overlapped with fruit size or shape QTLs. There were seven putative organ size regulators and seven capsaicin biosynthetic genes within these regions. This work defines genomic regions that underly structure within the nonpungent pepper germplasm and QTLs or genes that may have been selected for during the development of large-fruited nonpungent pepper varieties.

Characterization of Capsicum annuum genetic diversity and population structure based on parallel polymorphism discovery with a 30K unigene Pepper GeneChip.

The widely cultivated pepper, Capsicum spp., important as a vegetable and spice crop world-wide, is one of the most diverse crops. To enhance breeding programs, a detailed characterization of Capsicum diversity including morphological, geographical and molecular data is required. Currently, molecular data characterizing Capsicum genetic diversity is limited. The development and application of high-throughput genome-wide markers in Capsicum will facilitate more detailed molecular characterization of germplasm collections, genetic relationships, and the generation of ultra-high density maps. We have developed the Pepper GeneChip® array from Affymetrix for polymorphism detection and expression analysis in Capsicum. Probes on the array were designed from 30,815 unigenes assembled from expressed sequence tags (ESTs). Our array design provides a maximum redundancy of 13 probes per base pair position allowing integration of multiple hybridization values per position to detect single position polymorphism (SPP). Hybridization of genomic DNA from 40 diverse C. annuum lines, used in breeding and research programs, and a representative from three additional cultivated species (C. frutescens, C. chinense and C. pubescens) detected 33,401 SPP markers within 13,323 unigenes. Among the C. annuum lines, 6,426 SPPs covering 3,818 unigenes were identified. An estimated three-fold reduction in diversity was detected in non-pungent compared with pungent lines, however, we were able to detect 251 highly informative markers across these C. annuum lines. In addition, an 8.7 cM region without polymorphism was detected around Pun1 in non-pungent C. annuum. An analysis of genetic relatedness and diversity using the software Structure revealed clustering of the germplasm which was confirmed with statistical support by principle components analysis (PCA) and phylogenetic analysis. This research demonstrates the effectiveness of parallel high-throughput discovery and application of genome-wide transcript-based markers to assess genetic and genomic features among Capsicum annuum.

An Arabidopsis F-box protein acts as a transcriptional co-factor to regulate floral development.

Plants flower in response to both environmental and endogenous signals. The Arabidopsis LEAFY (LFY) transcription factor is crucial in integrating these signals, and acts in part by activating the expression of multiple floral homeotic genes. LFY-dependent activation of the homeotic APETALA3 (AP3) gene requires the activity of UNUSUAL FLORAL ORGANS (UFO), an F-box component of an SCF ubiquitin ligase, yet how this regulation is effected has remained unclear. Here, we show that UFO physically interacts with LFY both in vitro and in vivo, and this interaction is necessary to recruit UFO to the AP3 promoter. Furthermore, a transcriptional repressor domain fused to UFO reduces endogenous LFY activity in plants, supporting the idea that UFO acts as part of a transcriptional complex at the AP3 promoter. Moreover, chemical or genetic disruption of proteasome activity compromises LFY-dependent AP3 activation, indicating that protein degradation is required to promote LFY activity. These results define an unexpected role for an F-box protein in functioning as a DNA-associated transcriptional co-factor in regulating floral homeotic gene expression. These results suggest a novel mechanism for promoting flower development via protein degradation and concomitant activation of the LFY transcription factor. This mechanism may be widely conserved, as homologs of UFO and LFY have been identified in a wide array of plant species.

The nuclear genome of Brachypodium distachyon: analysis of BAC end sequences.

Due in part to its small genome (approximately 350 Mb), Brachypodium distachyon is emerging as a model system for temperate grasses, including important crops like wheat and barley. We present the analysis of 10.9% of the Brachypodium genome based on 64,696 bacterial artificial chromosome (BAC) end sequences (BES). Analysis of repeat DNA content in BES revealed that approximately 11.0% of the genome consists of known repetitive DNA. The vast majority of the Brachypodium repetitive elements are LTR retrotransposons. While Bare-1 retrotransposons are common to wheat and barley, Brachypodium repetitive element sequence-1 (BRES-1), closely related to Bare-1, is also abundant in Brachypodium. Moreover, unique Brachypodium repetitive element sequences identified constitute approximately 7.4% of its genome. Simple sequence repeats from BES were analyzed, and flanking primer sequences for SSR detection potentially useful for genetic mapping are available at http://brachypodium.pw.usda.gov . Sequence analyses of BES indicated that approximately 21.2% of the Brachypodium genome represents coding sequence. Furthermore, Brachypodium BES have more significant matches to ESTs from wheat than rice or maize, although these species have similar sizes of EST collections. A phylogenetic analysis based on 335 sequences shared among seven grass species further revealed a closer relationship between Brachypodium and Triticeae than Brachypodium and rice or maize.

ABERRANT TESTA SHAPE encodes a KANADI family member, linking polarity determination to separation and growth of Arabidopsis ovule integuments.

The Arabidopsis aberrant testa shape (ats) mutant produces a single integument instead of the two integuments seen in wild-type ovules. Cellular anatomy and patterns of marker gene expression indicate that the single integument results from congenital fusion of the two integuments of the wild type. Isolation of the ATS locus showed it to encode a member of the KANADI (KAN) family of putative transcription factors, previously referred to as KAN4. ATS was expressed at the border between the two integuments at the time of their initiation, with expression later confined to the abaxial layer of the inner integument. In an inner no outer (ino) mutant background, where an outer integument does not form, the ats mutation led to amorphous inner integument growth. The kan1kan2 double mutant exhibits a similar amorphous growth of the outer integument without affecting inner integument growth. We hypothesize that ATS and KAN1/KAN2 play similar roles in the specification of polarity in the inner and outer integuments, respectively, that parallel the known roles of KAN proteins in promoting abaxial identity during leaf development. INO and other members of the YABBY gene family have been hypothesized to have similar parallel roles in outer integument and leaf development. Together, these two hypotheses lead us to propose a model for normal integument growth that also explains the described mutant phenotypes.

Regulation of APETALA3 floral homeotic gene expression by meristem identity genes.

The Arabidopsis APETALA3 (AP3) floral homeotic gene is required for specifying petal and stamen identities, and is expressed in a spatially limited domain of cells in the floral meristem that will give rise to these organs. Here we show that the floral meristem identity genes LEAFY (LFY) and APETALA1 (AP1) are required for the activation of AP3. The LFY transcription factor binds to a sequence, with dyad symmetry, that lies within a region of the AP3 promoter required for early expression of AP3. Mutation of this region abolishes LFY binding in vitro and in yeast one hybrid assays, but has no obvious effect on AP3 expression in planta. Experiments using a steroid-inducible form of LFY show that, in contrast to its direct transcriptional activation of other floral homeotic genes, LFY acts in both a direct and an indirect manner to regulate AP3 expression. This LFY-induced expression of AP3 depends in part on the function of the APETALA1 (AP1) floral homeotic gene, since mutations in AP1 reduce LFY-dependent induction of AP3 expression. LFY therefore appears to act through several pathways, one of which is dependent on AP1 activity, to regulate AP3 expression.