Extending the Toolkit for Beauty: Differential Co-Expression of DROOPING LEAF-Like and Class B MADS-Box Genes during Phalaenopsis Flower Development.

The molecular basis of orchid flower development is accomplished through a specific regulatory program in which the class B MADS-box AP3/DEF genes play a central role. In particular, the differential expression of four class B AP3/DEF genes is responsible for specification of organ identities in the orchid perianth. Other MADS-box genes (AGL6 and SEP-like) enrich the molecular program underpinning the orchid perianth development, resulting in the expansion of the original "orchid code" in an even more complex gene regulatory network. To identify candidates that could interact with the AP3/DEF genes in orchids, we conducted an in silico differential expression analysis in wild-type and peloric Phalaenopsis. The results suggest that a YABBY DL-like gene could be involved in the molecular program leading to the development of the orchid perianth, particularly the labellum. Two YABBY DL/CRC homologs are present in the genome of Phalaenopsis equestris, PeDL1 and PeDL2, and both express two alternative isoforms. Quantitative real-time PCR analyses revealed that both genes are expressed in column and ovary. In addition, PeDL2 is more strongly expressed the labellum than in the other tepals of wild-type flowers. This pattern is similar to that of the AP3/DEF genes PeMADS3/4 and opposite to that of PeMADS2/5. In peloric mutant Phalaenopsis, where labellum-like structures substitute the lateral inner tepals, PeDL2 is expressed at similar levels of the PeMADS2-5 genes, suggesting the involvement of PeDL2 in the development of the labellum, together with the PeMADS2-PeMADS5 genes. Although the yeast two-hybrid analysis did not reveal the ability of PeDL2 to bind the PeMADS2-PeMADS5 proteins directly, the existence of regulatory interactions is suggested by the presence of CArG-boxes and other MADS-box transcription factor binding sites within the putative promoter of the orchid DL2 gene.

Diversification of defensins and NLRs in Arabidopsis species by different evolutionary mechanisms.

BACKGROUND: Genes encoding proteins underlying host-pathogen co-evolution and which are selected for new resistance specificities frequently are under positive selection, a process that maintains diversity. Here, we tested the contribution of natural selection, recombination and transcriptional divergence to the evolutionary diversification of the plant defensins superfamily in three Arabidopsis species. The intracellular NOD-like receptor (NLR) family was used for comparison because positive selection has been well documented in its members. Similar to defensins, NLRs are encoded by a large and polymorphic gene family and many of their members are involved in the immune response. RESULTS: Gene trees of Arabidopsis defensins (DEFLs) show a high prevalence of clades containing orthologs. This indicates that their diversity dates back to a common ancestor and species-specific duplications did not significantly contribute to gene family expansion. DEFLs are characterized by a pervasive pattern of neutral evolution with infrequent positive and negative selection as well as recombination. In comparison, most NLR alignment groups are characterized by frequent occurrence of positive selection and recombination in their leucine-rich repeat (LRR) domain as well negative selection in their nucleotide-binding (NB-ARC) domain. While major NLR subgroups are expressed in pistils and leaves both in presence or absence of pathogen infection, the members of DEFL alignment groups are predominantly transcribed in pistils. Furthermore, conserved groups of NLRs and DEFLs are differentially expressed in response to Fusarium graminearum regardless of whether these genes are under positive selection or not. CONCLUSIONS: The present analyses of NLRs expands previous studies in Arabidopsis thaliana and highlights contrasting patterns of purifying and diversifying selection affecting different gene regions. DEFL genes show a different evolutionary trend, with fewer recombination events and significantly fewer instances of natural selection. Their heterogeneous expression pattern suggests that transcriptional divergence probably made the major contribution to functional diversification. In comparison to smaller families encoding pathogenesis-related (PR) proteins under positive selection, DEFLs are involved in a wide variety of processes that altogether might pose structural and functional trade-offs to their family-wide pattern of evolution.

Similarities between Reproductive and Immune Pistil Transcriptomes of Arabidopsis Species.

Independent lines of evidence suggest that members from ancient and polymorphic gene families such as defensins and receptor-like kinases mediate intercellular communication during both the immune response and reproduction. Here, we report a large-scale analysis to investigate the extent of overlap between these processes by comparing differentially expressed genes (DEGs) in the pistil transcriptomes of Arabidopsis thaliana and Arabidopsis halleri during self-pollination and interspecific pollination and during infection with Fusarium graminearum In both Arabidopsis species, the largest number of DEGs was identified in infected pistils, where genes encoding regulators of cell division and development were most frequently down-regulated. Comparison of DEGs between infection and various pollination conditions showed that up to 79% of down-regulated genes are shared between conditions and include especially defensin-like genes. Interspecific pollination of A.thaliana significantly up-regulated thionins and defensins. The significant overrepresentation of similar groups of DEGs in the transcriptomes of reproductive and immune responses of the pistil makes it a prime system in which to study the consequences of plant-pathogen interactions on fertility and the evolution of intercellular communication in pollination.

Expression of paralogous SEP-, FUL-, AG- and STK-like MADS-box genes in wild-type and peloric Phalaenopsis flowers.

The diverse flowers of Orchidaceae are the result of several major morphological transitions, among them the most studied is the differentiation of the inner median tepal into the labellum, a perianth organ key in pollinator attraction. Type A peloria lacking stamens and with ectopic labella in place of inner lateral tepals are useful for testing models on the genes specifying these organs by comparing their patterns of expression between wild-type and peloric flowers. Previous studies focused on DEFICIENS- and GLOBOSA-like MADS-box genes because of their conserved role in perianth and stamen development. The "orchid code" model summarizes this work and shows in Orchidaceae there are four paralogous lineages of DEFICIENS/AP3-like genes differentially expressed in each floral whorl. Experimental tests of this model showed the conserved, higher expression of genes from two specific DEF-like gene lineages is associated with labellum development. The present study tests whether eight MADS-box candidate SEP-, FUL-, AG-, and STK-like genes have been specifically duplicated in the Orchidaceae and are also differentially expressed in association with the distinct flower organs of Phalaenopsis hyb. "Athens." The gene trees indicate orchid-specific duplications. In a way analogous to what is observed in labellum-specific DEF-like genes, a two-fold increase in the expression of SEP3-like gene PhaMADS7 was measured in the labellum-like inner lateral tepals of peloric flowers. The overlap between SEP3-like and DEF-like genes suggests both are associated with labellum specification and similar positional cues determine their domains of expression. In contrast, the uniform messenger levels of FUL-like genes suggest they are involved in the development of all organs and their expression in the ovary suggests cell differentiation starts before pollination. As previously reported AG-like and STK-like genes are exclusively expressed in gynostemium and ovary, however no evidence for transcriptional divergence was found in the stage investigated. Gene expression suggests a developmental regulatory system based on the combined activity of duplicate MADS-box genes. We discuss its feasibility based on documented protein interactions and patterns of expression.

Perspectives on MADS-box expression during orchid flower evolution and development.

The diverse morphology of orchid flowers and their complex, often deceptive strategies to become pollinated have fascinated researchers for a long time. However, it was not until the 20th century that the ontogeny of orchid flowers, the genetic basis of their morphology and the complex phylogeny of Orchidaceae were investigated. In parallel, the improvement of techniques for in vitro seed germination and tissue culture, together with studies on biochemistry, physiology, and cytology supported the progress of what is now a highly productive industry of orchid breeding and propagation. In the present century both basic research in orchid flower evo-devo and the interest for generating novel horticultural varieties have driven the characterization of many members of the MADS-box family encoding key regulators of flower development. This perspective summarizes the picture emerging from these studies and discusses the advantages and limitations of the comparative strategy employed so far. I address the growing role of natural and horticultural mutants in these studies and the emergence of several model species in orchid evo-devo and genomics. In this context, I make a plea for an increasingly integrative approach.

High time for a roll call: gene duplication and phylogenetic relationships of TCP-like genes in monocots.

BACKGROUND AND AIMS: The TCP family is an ancient group of plant developmental transcription factors that regulate cell division in vegetative and reproductive structures and are essential in the establishment of flower zygomorphy. In-depth research on eudicot TCPs has documented their evolutionary and developmental role. This has not happened to the same extent in monocots, although zygomorphy has been critical for the diversification of Orchidaceae and Poaceae, the largest families of this group. Investigating the evolution and function of TCP-like genes in a wider group of monocots requires a detailed phylogenetic analysis of all available sequence information and a system that facilitates comparing genetic and functional information. METHODS: The phylogenetic relationships of TCP-like genes in monocots were investigated by analysing sequences from the genomes of Zea mays, Brachypodium distachyon, Oryza sativa and Sorghum bicolor, as well as EST data from several other monocot species. KEY RESULTS: All available monocot TCP-like sequences are associated in 20 major groups with an average identity ≥64 % and most correspond to well-supported clades of the phylogeny. Their sequence motifs and relationships of orthology were documented and it was found that 67 % of the TCP-like genes of Sorghum, Oryza, Zea and Brachypodium are in microsyntenic regions. This analysis suggests that two rounds of whole genome duplication drove the expansion of TCP-like genes in these species. CONCLUSIONS: A system of classification is proposed where putative or recognized monocot TCP-like genes are assigned to a specific clade of PCF-, CIN- or CYC/tb1-like genes. Specific biases in sequence data of this family that must be tackled when studying its molecular evolution and phylogeny are documented. Finally, the significant retention of duplicated TCP genes from Zea mays is considered in the context of balanced gene drive.

Conserved differential expression of paralogous DEFICIENS- and GLOBOSA-like MADS-box genes in the flowers of Orchidaceae: refining the 'orchid code'.

In flowering plants, class-B floral homeotic genes encode MADS-domain transcription factors, which are key in the specification of petal and stamen identity, and have two ancient clades: DEF-like and GLO-like genes. Many species have one gene of each clade, but orchids have typically four DEF-like genes, representing ancient gene clades 1, 2, 3 and 4. We tested the 'orchid code', a combinatorial genetic model suggesting that differences between the organs of the orchid perianth (outer tepals, inner lateral tepals and labellum) are generated by the combinatorial differential expression of four DEF-like genes. Our experimental test involves highly sensitive and specific measurements, with qRT-PCR of the expression of DEF- and GLO-like genes from the distantly related Vanilla planifolia and Phragmipedium longifolium, as well as from wild-type and peloric Phalaenopsis hybrid flowers. Our findings support the first 'orchid code' hypothesis, in that absence of clade-3 and -4 gene expression distinguishes the outer tepals from the inner tepals. In contrast to the original hypothesis, however, mRNA of both clade-3 and -4 genes accumulates in wild-type inner lateral tepals and the labellum, and in labellum-like inner lateral tepals of peloric flowers, albeit in different quantities. Our data suggest a revised hypothesis where high levels of clade-1 and -2, and low levels of clade-3 and -4, gene expression specify inner lateral tepals, whereas labellum development requires low levels of clade-1 and -2 expression and high levels of clade-3 and -4 expression.

Positive selection and ancient duplications in the evolution of class B floral homeotic genes of orchids and grasses.

BACKGROUND: Positive selection is recognized as the prevalence of nonsynonymous over synonymous substitutions in a gene. Models of the functional evolution of duplicated genes consider neofunctionalization as key to the retention of paralogues. For instance, duplicate transcription factors are specifically retained in plant and animal genomes and both positive selection and transcriptional divergence appear to have played a role in their diversification. However, the relative impact of these two factors has not been systematically evaluated. Class B MADS-box genes, comprising DEF-like and GLO-like genes, encode developmental transcription factors essential for establishment of perianth and male organ identity in the flowers of angiosperms. Here, we contrast the role of positive selection and the known divergence in expression patterns of genes encoding class B-like MADS-box transcription factors from monocots, with emphasis on the family Orchidaceae and the order Poales. Although in the monocots these two groups are highly diverse and have a strongly canalized floral morphology, there is no information on the role of positive selection in the evolution of their distinctive flower morphologies. Published research shows that in Poales, class B-like genes are expressed in stamens and in lodicules, the perianth organs whose identity might also be specified by class B-like genes, like the identity of the inner tepals of their lily-like relatives. In orchids, however, the number and pattern of expression of class B-like genes have greatly diverged. RESULTS: The DEF-like genes from Orchidaceae form four well-supported, ancient clades of orthologues. In contrast, orchid GLO-like genes form a single clade of ancient orthologues and recent paralogues. DEF-like genes from orchid clade 2 (OMADS3-like genes) are under less stringent purifying selection than the other orchid DEF-like and GLO-like genes. In comparison with orchids, purifying selection was less stringent in DEF-like and GLO-like genes from Poales. Most importantly, positive selection took place before the major organ reduction and losses in the floral axis that eventually yielded the zygomorphic grass floret. CONCLUSION: In DEF-like genes of Poales, positive selection on the region mediating interactions with other proteins or DNA could have triggered the evolution of the regulatory mechanisms behind the development of grass-specific reproductive structures. Orchidaceae show a different trend, where gene duplication and transcriptional divergence appear to have played a major role in the canalization and modularization of perianth development.

Why are orchid flowers so diverse? Reduction of evolutionary constraints by paralogues of class B floral homeotic genes.

BACKGROUND: The nearly 30 000 species of orchids produce flowers of unprecedented diversity. However, whether specific genetic mechanisms contributed to this diversity is a neglected topic and remains speculative. We recently published a theory, the 'orchid code', maintaining that the identity of the different perianth organs is specified by the combinatorial interaction of four DEF-like MADS-box genes with other floral homeotic genes. SCOPE: Here the developmental and evolutionary implications of our theory are explored. Specifically, it is shown that all frequent floral terata, including all peloric types, can be explained by monogenic gain- or-loss-of-function mutants, changing either expression of a DEF-like or CYC-like gene. Supposed dominance or recessiveness of mutant alleles is correlated with the frequency of terata in both cultivation and nature. Our findings suggest that changes in DEF- and CYC-like genes not only underlie terata but also the natural diversity of orchid species. We argue, however, that true changes in organ identity are rare events in the evolution of orchid flowers, even though we review some likely cases. CONCLUSIONS: The four DEF paralogues shaped floral diversity in orchids in a dramatic way by modularizing the floral perianth based on a complex series of sub- and neo-functionalization events. These genes may have eliminated constraints, so that different kinds of perianth organs could then evolve individually and thus often in dramatically different ways in response to selection by pollinators or by genetic drift. We therefore argue that floral diversity in orchids may be the result of an unprecedented developmental genetic predisposition that originated early in orchid evolution.

MADS about the evolution of orchid flowers.

Orchids have unique flowers involving three types of perianth organs: outer tepals, lateral inner tepals, and a lip. Expression studies indicate that the identity of these organs is specified by the combinatorial interaction of four different DEFICIENS-like MADS-box genes. We suggest that clarifying the evolution of these genes provides a rational framework for reconstructing the enigmatic origin and unique diversification of the orchid flower. For example, two rounds of gene duplications during early orchid evolution might have generated the genes that were probably recruited to distinguish the different types of perianth organs. This hypothesis suggests intriguing, experimentally testable mechanisms by which gene duplications followed by sub- and neo-functionalization events might have contributed to the evolutionary origin of morphological novelties in orchids - and well beyond.

Adaptive selection of mitochondrial complex I subunits during primate radiation.

Mammalian oxidative phosphorylation (OXPHOS) complexes I, III, IV and V are assembled from both mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) encoded subunits, with complex I encompassing 39 nDNA and seven mtDNA subunits. Yet the sequence variation of the mtDNA genes is more than ten fold greater than that of the nDNA encoded genes of the OXPHOS complexes and the mtDNA proteins have been found to be influenced by positive (adaptive) selection. To maintain a functional complex I, nDNA and mtDNA subunits must interact, implying that certain nDNA complex I genes may also have been influenced by positive selection. To determine if positive selection has influenced nDNA complex I genes, we analyzed the DNA sequences of all of the nDNA and mtDNA encoded complex I subunits from orangutan, gorilla, chimpanzee, human and all available vertebrate sequences. This revealed that three nDNA complex I genes (NDUFC2, NDUFA1, and NDUFA4) had significantly increased amino acid substitution rates by both PAML and Z-test, suggesting that they have been subjected to adaptive selection during primate radiation. Since all three of these subunits reside in the membrane domain of complex I along with the mtDNA subunits, we compared amino acid changes in these three nDNA genes with those of the mtDNA genes across species. Changes in the nDNA NDUFC2 cysteine 39 were found to correlate with those in the mtDNA ND5 cysteine 330. Therefore, adaptive selection has influenced some nDNA complex I genes and nDNA and mtDNA complex I genes may have co-evolved.

Gene conversion and the evolution of three leucine-rich repeat gene families in Arabidopsis thaliana.

The high number of duplicated genes in plant genomes provides a potential template for gene conversion and unequal crossing-over. Within a gene family these two processes can render all members homogeneous or generate diversity by reassorting variants among paralogs. The latter is especially feasible in families where gene diversity confers a selective advantage and thus conversion events are likely to be retained. Consequently, the most complete record of gene conversion is expected to be most evident in gene families commonly subjected to positive selection. Here, we describe the extent and characteristics of gene conversion and unequal crossing-over in the coding and noncoding regions of nucleotide-binding site leucine-rich repeat (NBS-LRR), receptor-like kinases (RLK), and receptor-like proteins (RLP) in the plant Arabidopsis thaliana. Members of these three gene families are associated with disease resistance and their pathogen-recognition domain is a documented target of positive selection. Our bioinformatic approach to study the major family features that may influence gene conversion revealed that in these families there is a significant association between the occurrence of gene conversion and high levels of sequence similarity, close physical clustering, gene orientation, and recombination rate. We discuss these results in the context of the overlap between gene conversion and positive selection during the evolutionary expansion of the NBS-LRR, RLK, and RLP gene families.

Patterns of positive selection in the complete NBS-LRR gene family of Arabidopsis thaliana.

Plant disease resistance genes have been shown to be subject to positive selection, particularly in the leucine rich repeat (LRR) region that may determine resistance specificity. We performed a genome-wide analysis of positive selection in members of the nucleotide binding site (NBS)-LRR gene family of Arabidopsis thaliana. Analyses were possible for 103 of 163 NBS-LRR nucleotide sequences in the genome, and the analyses uncovered substantial evidence of positive selection. Sites under positive selection were detected and identified for 10 sequence groups representing 53 NBS-LRR sequences. Functionally characterized Arabidopsis resistance genes were in these 10 groups, but several groups with extensive evidence of positive selection contained no previously characterized resistance genes. Amino acid residues under positive selection were identified, and these residues were mapped onto protein secondary structure. Positively selected positions were disproportionately located in the LRR domain (P < 0.001), particularly a nine-amino acid beta-strand submotif that is likely to be solvent exposed. However, a substantial proportion (30%) of positively selected sites were located outside LRRs, suggesting that regions other than the LRR may function in determining resistance specificity. Because of the unusual sequence variability in the LRRs of this class of proteins, secondary-structure analysis identifies LRRs that are not identified by similarity analyses alone. LRRs also contain substantial indel variation, suggesting elasticity in LRR length could also influence resistance specificity.

The PE-PGRS glycine-rich proteins of Mycobacterium tuberculosis: a new family of fibronectin-binding proteins?

A clone was isolated by screening of a cosmid library of Mycobacterium tuberculosis with an oligonucleotide designed from the N-terminal sequence of a previously reported proline-rich protein. Characterization of the 4481 bp insert showed the presence of polymorphic CG-repetitive sequences (PGRSs) with an ORF of 2.7 kb, encoding a 81.3 kDa protein (PE-PGRS81). Southern blot analysis and BLAST-p searches revealed several homologous sequences in the genome of M. tuberculosis. The deduced amino acid sequence was highly similar to a stretch of about 98 residues in the N-terminus present in several members of the PE-PGRS family available in the GenBank database, including 100% identity with the partial amino acid sequence of the potential protein encoded by orf3' as well as with the Rv0278c sequence. A neighbour-joining analysis of the 99 PE-PGRS sequences available in the database indicated that PE-PGRS81 is included in a group where its closest relatives are the sequences orf3', Rv0278c, Rv0279c, Rv1759c, Rv3652 and Rv0747. Probing with the complete coding regions of PE-PGRS81 and Rv1759c in Southern blot assays, on samples of genomic DNA from M. tuberculosis H37Rv, Mycobacterium bovis BCG and M. tuberculosis clinical isolates, showed a complex hybridization pattern for all strains. This shows the existence of intrastrain PGRS variability as reported for other PGRS members. In contrast, probing with the short conserved N-terminal region of Rv1759c reduced the hybridization to a single band. This marker allowed identification of M. tuberculosis clinical strains that lack Rv1759c. A recombinant C-terminal fragment of Rv1759c showed fibronectin-binding properties and was recognized by sera from patients infected with M. tuberculosis, suggesting that at least this member of the PE-PGRS is expressed in tuberculosis infection.