Extensive gene tree discordance and hemiplasy shaped the genomes of North American columnar cacti.

Few clades of plants have proven as difficult to classify as cacti. One explanation may be an unusually high level of convergent and parallel evolution (homoplasy). To evaluate support for this phylogenetic hypothesis at the molecular level, we sequenced the genomes of four cacti in the especially problematic tribe Pachycereeae, which contains most of the large columnar cacti of Mexico and adjacent areas, including the iconic saguaro cactus (Carnegiea gigantea) of the Sonoran Desert. We assembled a high-coverage draft genome for saguaro and lower coverage genomes for three other genera of tribe Pachycereeae (Pachycereus, Lophocereus, and Stenocereus) and a more distant outgroup cactus, Pereskia We used these to construct 4,436 orthologous gene alignments. Species tree inference consistently returned the same phylogeny, but gene tree discordance was high: 37% of gene trees having at least 90% bootstrap support conflicted with the species tree. Evidently, discordance is a product of long generation times and moderately large effective population sizes, leading to extensive incomplete lineage sorting (ILS). In the best supported gene trees, 58% of apparent homoplasy at amino sites in the species tree is due to gene tree-species tree discordance rather than parallel substitutions in the gene trees themselves, a phenomenon termed "hemiplasy." The high rate of genomic hemiplasy may contribute to apparent parallelisms in phenotypic traits, which could confound understanding of species relationships and character evolution in cacti.

A modified GC-specific MAKER gene annotation method reveals improved and novel gene predictions of high and low GC content in Oryza sativa.

BACKGROUND: Accurate structural annotation depends on well-trained gene prediction programs. Training data for gene prediction programs are often chosen randomly from a subset of high-quality genes that ideally represent the variation found within a genome. One aspect of gene variation is GC content, which differs across species and is bimodal in grass genomes. When gene prediction programs are trained on a subset of grass genes with random GC content, they are effectively being trained on two classes of genes at once, and this can be expected to result in poor results when genes are predicted in new genome sequences. RESULTS: We find that gene prediction programs trained on grass genes with random GC content do not completely predict all grass genes with extreme GC content. We show that gene prediction programs that are trained with grass genes with high or low GC content can make both better and unique gene predictions compared to gene prediction programs that are trained on genes with random GC content. By separately training gene prediction programs with genes from multiple GC ranges and using the programs within the MAKER genome annotation pipeline, we were able to improve the annotation of the Oryza sativa genome compared to using the standard MAKER annotation protocol. Gene structure was improved in over 13% of genes, and 651 novel genes were predicted by the GC-specific MAKER protocol. CONCLUSIONS: We present a new GC-specific MAKER annotation protocol to predict new and improved gene models and assess the biological significance of this method in Oryza sativa. We expect that this protocol will also be beneficial for gene prediction in any organism with bimodal or other unusual gene GC content.

A G-Box-Like Motif Is Necessary for Transcriptional Regulation by Circadian Pseudo-Response Regulators in Arabidopsis.

PSEUDO-RESPONSE REGULATORs (PRRs) play overlapping and distinct roles in maintaining circadian rhythms and regulating diverse biological processes, including the photoperiodic control of flowering, growth, and abiotic stress responses. PRRs act as transcriptional repressors and associate with chromatin via their conserved C-terminal CCT (CONSTANS, CONSTANS-like, and TIMING OF CAB EXPRESSION 1 [TOC1/PRR1]) domains by a still-poorly understood mechanism. Here, we identified genome-wide targets of PRR9 using chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) and compared them with PRR7, PRR5, and TOC1/PRR1 ChIP-seq data. We found that PRR binding sites are located within genomic regions of low nucleosome occupancy and high DNase I hypersensitivity. Moreover, conserved noncoding regions among Brassicaceae species are enriched around PRR binding sites, indicating that PRRs associate with functionally relevant cis-regulatory regions. The PRRs shared a significant number of binding regions, and our results indicate that they coordinately restrict the expression of target genes to around dawn. A G-box-like motif was overrepresented at PRR binding regions, and we showed that this motif is necessary for mediating transcriptional regulation of CIRCADIAN CLOCK ASSOCIATED 1 and PRR9 by the PRRs. Our results further our understanding of how PRRs target specific promoters and provide an extensive resource for studying circadian regulatory networks in plants.

Fine mapping of the stem rust resistance gene SrTA10187.

KEY MESSAGE: SrTA10187 was fine-mapped to a 1.1 cM interval, candidate genes were identified in the region of interest, and molecular markers were developed for marker-assisted selection and Sr gene pyramiding. Stem rust (Puccinia graminis f. sp. tritici, Pgt) races belonging to the Ug99 (TTKSK) race group pose a serious threat to global wheat (Triticum aestivum L.) production. To improve Pgt host resistance, the Ug99-effective resistance gene SrTA10187 previously identified in Aegilops tauschii Coss. was introgressed into wheat, and mapped to the short arm of wheat chromosome 6D. In this study, high-resolution mapping of SrTA10187 was done using a population of 1,060 plants. Pgt resistance was screened using race QFCSC. PCR-based SNP and STS markers were developed from genotyping-by-sequencing tags and SNP sequences available in online databases. SrTA10187 segregated as expected in a 3:1 ratio of resistant to susceptible individuals in three out of six BC3F2 families, and was fine-mapped to a 1.1 cM region on wheat chromosome 6DS. Marker context sequence was aligned to the reference Ae. tauschii genome to identify the physical region encompassing SrTA10187. Due to the size of the corresponding region, candidate disease resistance genes could not be identified with confidence. Comparisons with the Ae. tauschii genetic map developed by Luo et al. (PNAS 110(19):7940-7945, 2013) enabled identification of a discrete genetic locus and a BAC minimum tiling path of the region spanning SrTA10187. Annotation of pooled BAC library sequences led to the identification of candidate genes in the region of interest-including a single NB-ARC-LRR gene. The shorter genetic interval and flanking KASP™ and STS markers developed in this study will facilitate marker-assisted selection, gene pyramiding, and positional cloning of SrTA10187.

Definition of early transcriptional circuitry involved in light-induced reversal of PIF-imposed repression of photomorphogenesis in young Arabidopsis seedlings.

Light signals perceived by the phytochromes induce the transition from skotomorphogenic to photomorphogenic development (deetiolation) in dark-germinated seedlings. Evidence that a quadruple mutant (pifq) lacking four phytochrome-interacting bHLH transcription factors (PIF1, 3, 4, and 5) is constitutively photomorphogenic in darkness establishes that these factors sustain the skotomorphogenic state. Moreover, photoactivated phytochromes bind to and induce rapid degradation of the PIFs, indicating that the photoreceptor reverses their constitutive activity upon light exposure, initiating photomorphogenesis. Here, to define the modes of transcriptional regulation and cellular development imposed by the PIFs, we performed expression profile and cytological analyses of pifq mutant and wild-type seedlings. Dark-grown mutant seedlings display cellular development that extensively phenocopies wild-type seedlings grown in light. Similarly, 80% of the gene expression changes elicited by the absence of the PIFs in dark-grown pifq seedlings are normally induced by prolonged light in wild-type seedlings. By comparing rapidly light-responsive genes in wild-type seedlings with those responding in darkness in the pifq mutant, we identified a subset, enriched in transcription factor-encoding genes, that are potential primary targets of PIF transcriptional regulation. Collectively, these data suggest that the transcriptional response elicited by light-induced PIF proteolysis is a major component of the mechanism by which the phytochromes pleiotropically regulate deetiolation and that at least some of the rapidly light-responsive genes may comprise a transcriptional network directly regulated by the PIF proteins.