Characterization,Polymorphism, and Evolution of MHC Class II B Genes in Birds of Prey

During the last decade, the major histocompatibility complex (MHC) has received much attention in the fields of evolutionary and conservation biology because of its potential implications in many biological processes. New insights into the gene structure and evolution of MHC genes can be gained through study of additional lineages of birds not yet investigated at the genomic level. In this study, we characterized MHC class II B genes in five families of birds of prey (Accipitridae, Pandionidae, Strigidae, Tytonidae, and Falconidae). Using PCR approaches, we isolated genomic MHC sequences up to 1300 bp spanning exons 1 to 3 in 26 representatives of each raptor lineage, finding no stop codons or frameshift mutations in any coding region. A survey of diversity across the entirety of exon 2 in the lesser kestrel Falco naumanni reported 26 alleles in 21 individuals. Bayesian analysis revealed 21 positively selected amino acid sites, which suggests that the MHC genes described here are functional and probably expressed. Finally, through interlocus comparisons and phylogenetic analysis, we also discuss genetic evidence for concerted and transspecies evolution in the raptor MHC.     

Characterization and Phylogenetic Relationship of Prosimian MHC Class I Genes

MHC class I cDNA sequences from the most divergent primate group of extant primates compared to human, the suborder Strepsirrhini (prosimians), are described. The sequences are derived from the gray mouse lemur (Microcebus murinus) and the ring-tailed lemur (Lemur catta), which are members of the malagasy Lemuriformes, as well as from the pygmy slow loris (Nycticebus pygmaeus), a prosimian from East Asia. The M. murinus sequences have been analyzed in detail. Analysis of the expression level, G/C content, and synonymous vs. nonsynonymous substitution rates in the peptide-binding region codons suggests that these cDNA clones represent classical class I (class Ia) genes. According to Southern blot analysis, the genome of the gray mouse lemur might contain about 10 class I genes. In gene tree analysis, the strepsirrhine class Ia genes described here cluster significantly separately from the known class I genes of Catarrhini (humans, apes, Old World monkeys) and Platyrrhini (New World monkeys) species, suggesting that the class I loci of Simiiformes arose by gene duplications which occurred after the divergence of prosimians.     

Development of MHC Class I and II B Primers in Common Carp and its Molecular Characterization

The major histocompatibility complex (MHC) has an important role in immune response and is known as the most polymorphic locus in vertebrates. We developed three pairs of polymerase chain reaction primers of the alpha-2 domain (exon 3) of MHC class I and the beta-2 (exon 3) and beta-3 domains (exon 4) of MHC class II B gene in the German mirror common carp (Cyprinus carpio L.). We analyzed the three loci in a population of 65 individuals that had suffered the serious disease of gill rot. Five to six variable nucleotide sites and two to six variable amino acid sites (71.43%) were detected in the exon sequence of the sampled populations, indicating that many of them corresponded to amino acids involved in antigen recognition. Deviation from Hardy–Weinberg equilibrium and linkage disequilibrium were differentially found in some loci, which will be important for further study of disease resistance/susceptibility and population evolution.     

Class II MHC self-antigen presentation in human B and T lymphocytes

Human CD4(+) T cells process and present functional class II MHC-peptide complexes, but the endogenous peptide repertoire of these non-classical antigen presenting cells remains unknown. We eluted and sequenced HLA-DR-bound self-peptides presented by CD4(+) T cells in order to compare the T cell-derived peptide repertoire to sequences derived from genetically identical B cells. We identified several novel epitopes derived from the T cell-specific proteome, including fragments of CD4 and IL-2. While these data confirm that T cells can present peptides derived from the T-cell specific proteome, the vast majority of peptides sequenced after elution from MHC were derived from the common proteome. From this pool, we identified several identical peptide epitopes in the T and B cell repertoire derived from common endogenous proteins as well as novel endogenous epitopes with promiscuous binding. These findings indicate that the endogenous HLA-DR-bound peptide repertoire, regardless of APC type and across MHC isotype, is largely derived from the same pool of self-protein.     

Trans-Species Polymorphism and Selection in the MHC Class II DRA Genes of Domestic Sheep

Highly polymorphic genes with central roles in lymphocyte mediated immune surveillance are grouped together in the major histocompatibility complex (MHC) in higher vertebrates. Generally, across vertebrate species the class II MHC DRA gene is highly conserved with only limited allelic variation. Here however, we provide evidence of trans-species polymorphism at the DRA locus in domestic sheep (Ovis aries). We describe variation at the Ovar-DRA locus that is far in excess of anything described in other vertebrate species. The divergent DRA allele (Ovar-DRA*0201) differs from the sheep reference sequences by 20 nucleotides, 12 of which appear non-synonymous. Furthermore, DRA*0201 is paired with an equally divergent DRB1 allele (Ovar-DRB1*0901), which is consistent with an independent evolutionary history for the DR sub-region within this MHC haplotype. No recombination was observed between the divergent DRA and B genes in a range of breeds and typical levels of MHC class II DR protein expression were detected at the surface of leukocyte populations obtained from animals homozygous for the DRA*0201, DRB1*0901 haplotype. Bayesian phylogenetic analysis groups Ovar-DRA*0201 with DRA sequences derived from species within the Oryx and Alcelaphus genera rather than clustering with other ovine and caprine DRA alleles. Tests for Darwinian selection identified 10 positively selected sites on the branch leading to Ovar-DRA*0201, three of which are predicted to be associated with the binding of peptide antigen. As the Ovis, Oryx and Alcelaphus genera have not shared a common ancestor for over 30 million years, the DRA*0201 and DRB1*0901 allelic pair is likely to be of ancient origin and present in the founding population from which all contemporary domestic sheep breeds are derived. The conservation of the integrity of this unusual DR allelic pair suggests some selective advantage which is likely to be associated with the presentation of pathogen antigen to T-cells and the induction of protective immunity.     

Genetic Contribution of MHC Class II Genes in Susceptibility to West Nile Virus Infection

WNV is a zoonotic neurotropic flavivirus that has recently emerged globally as a significant cause of viral encephalitis. The last five years, 624 incidents of WNV infection have been reported in Greece. The risk for severe WNV disease increases among immunosuppressed individuals implying thus the contribution of the MHC locus to the control of WNV infection. In order to investigate a possible association of MHC class II genes, especially HLA-DPA1, HLA-DQA1, HLA-DRB1, we examined 105 WNV patients, including 68 cases with neuroinvasive disease and 37 cases with mild clinical phenotype, collected during the period from 2010 to2013, and 100 control individuals selected form the Greek population. Typing was performed for exon 2 for all three genes. DQA1*01:01 was considered to be "protective" against WNV infection (25.4% vs 40.1%, P = 0.004) while DQA1*01:02 was associated with increased susceptibility (48.0% vs 32.1%, P = 0.003). Protection against neuroinvasion was associated with the presence of DRB1*11:02 (4.99% vs 0.0%, P = 0.018). DRB1*16:02 was also absent from the control cohort (P = 0.016). Three additional population control groups were used in order to validate our results. No statistically significant association with the disease was found for HLA-DPA alleles. The results of the present study provide some evidence that MHC class II is involved in the response to WNV infection, outlining infection "susceptibility" and "CNS-high-risk" candidates. Furthermore, three new alleles were identified while the frequency of all alleles in the study was compared with worldwide data. The characterization of the MHC locus could help to estimate the risk for severe WNV cases in a country.     

Coevolution between MHC Class I and Antigen-Processing Genes in Salamanders

Proteins encoded by antigen-processing genes (APGs) provide major histocompatibility complex (MHC) class I (MHC-I) with antigenic peptides. In mammals, polymorphic multigenic MHC-I family is served by monomorphic APGs, whereas in certain nonmammalian species both MHC-I and APGs are polymorphic and coevolve within stable haplotypes. Coevolution was suggested as an ancestral gnathostome feature, presumably enabling only a single highly expressed classical MHC-I gene. In this view coevolution, while optimizing some aspects of adaptive immunity, would also limit its flexibility by preventing the expansion of classical MHC-I into a multigene family. However, some nonmammalian taxa, such as salamanders, have multiple highly expressed MHC-I genes, suggesting either that coevolution is relaxed or that it does not prevent the establishment of multigene MHC-I. To distinguish between these two alternatives, we use salamanders (30 species from 16 genera representing six families) to test, within a comparative framework, a major prediction of the coevolution hypothesis: the positive correlation between MHC-I and APG diversity. We found that MHC-I diversity explained both within-individual and species-wide diversity of two APGs, TAP1 and TAP2, supporting their coevolution with MHC-I, whereas no consistent effect was detected for the other three APGs (PSMB8, PSMB9, and TAPBP). Our results imply that although coevolution occurs in salamanders, it does not preclude the expansion of the MHC-I gene family. Contrary to the previous suggestions, nonmammalian vertebrates thus may be able to accommodate diverse selection pressures with flexibility granted by rapid expansion or contraction of the MHC-I family, while retaining the benefits of coevolution between MHC-I and TAPs.     

Sequence Polymorphism and Evolution of Three Cetacean MHC Genes

Sequence variability at three major histocompatibility complex (MHC) genes (DQB, DRA, and MHC-I) of cetaceans was investigated in order to get an overall understanding of cetacean MHC evolution. Little sequence variation was detected at the DRA locus, while extensive and considerable variability were found at the MHC-I and DQB loci. Phylogenetic reconstruction and sequence comparison revealed extensive sharing of identical MHC alleles among different species at the three MHC loci examined. Comparisons of phylogenetic trees for these MHC loci with the trees reconstructed only based on non-PBR sites revealed that allelic similarity/identity possibly reflected common ancestry and were not due to adaptive convergence. At the same time, trans-species evolution was also evidenced that the allelic diversity of the three MHC loci clearly pre-dated species divergence events according to the relaxed molecular clock. It may be the forces of balancing selection acting to maintain the high sequence variability and identical alleles in trans-specific manner at the MHC-I and DQB loci.     

MHC Class II haplotypes of Colombian Amerindian tribes

We analyzed 1041 individuals belonging to 17 Amerindian tribes of Colombia, Chimila, Bari and Tunebo (Chibcha linguistic family), Embera, Waunana (Choco linguistic family), Puinave and Nukak (Maku-Puinave linguistic families), Cubeo, Guanano, Tucano, Desano and Piratapuyo (Tukano linguistic family), Guahibo and Guayabero (Guayabero Linguistic Family), Curripaco and Piapoco (Arawak linguistic family) and Yucpa (Karib linguistic family). for MHC class II haplotypes (HLA-DRB1, DQA1, DQB1). Approximately 90% of the MHC class II haplotypes found among these tribes are haplotypes frequently encountered in other Amerindian tribes. Nonetheless, striking differences were observed among Chibcha and non-Chibcha speaking tribes. The DRB1*04:04, DRB1*04:11, DRB1*09:01 carrying haplotypes were frequently found among non-Chibcha speaking tribes, while the DRB1*04:07 haplotype showed significant frequencies among Chibcha speaking tribes, and only marginal frequencies among non-Chibcha speaking tribes. Our results suggest that the differences in MHC class II haplotype frequency found among Chibcha and non-Chibcha speaking tribes could be due to genetic differentiation in Mesoamerica of the ancestral Amerindian population into Chibcha and non-Chibcha speaking populations before they entered into South America.     

Evolutionary History of MHC Class I Genes in the Mammalian Order Perissodactyla

We carried out an analysis of partial sequences from expressed major histocompatibility complex (MHC) class I genes isolated from a range of equid species and more distantly related members of the mammalian order Perissodactyla. Phylogenetic analysis revealed a minimum of six groups, five of which contained genes and alleles that are found in equid species and one group specific to the rhinoceros. Four of the groups contained only one, or very few sequences, indicating the presence of relatively nonpolymorphic loci, while another group contained the majority of the equid sequences identified. These data suggest that a diversification of MHC genes took place after the split between the Equidae and the Rhinocerotidae yet before the speciation events within the genus Equus. Received: 17 November 1998 / Accepted: 7 April 1999     

Presentation of Phagocytosed Antigens by MHC Class I and II

Phagocytosis provides innate immune cells with a mechanism to take up and destroy pathogenic bacteria, apoptotic cells and other large particles. In some cases, however, peptide antigens from these particles are preserved for presentation in association with major histocompatibility complex (MHC) class I or class II molecules in order to stimulate antigen‐specific T cells. Processing and presentation of antigens from phagosomes presents a number of distinct challenges relative to antigens internalized by other means; while bacterial antigens were among the first discovered to be presented to T cells, analyses of the cellular mechanisms by which peptides from phagocytosed antigens assemble with MHC molecules and by which these complexes are then expressed at the plasma membrane have lagged behind those of conventional model soluble antigens. In this review, we cover recent advances in our understanding of these processes, including the unique cross‐presentation of phagocytosed antigens by MHC class I molecules, and in their control by signaling modalities in phagocytic cells.     

Class II MHC molecules and antigen enter the same vesicles during internalization by resting B lymphocytes

Efficient presentation of Ag by a B cell to a T cell requires that Ag bind to the Ag receptor (Ig) on the B cell, after which it is internalized into an acid compartment where it is modified and returned to the cell surface in the context of class II MHC molecules. It remains uncertain whether processed Ag binds to class II which has been internalized and recycled with Ag, or to nascent class II inside the cell. To determine if cell surface class II enters the same vesicles as Ag, or is excluded during internalization of Ag which is bound to the B cell receptor, 5- and 16-nm gold particles were labeled with anti-class II and anti-Ig, respectively. Cells were incubated at 37 degrees C and internalization of these particles was observed using electron microscopy. By 10 min, 60-75% of the B cell sections contained vesicles with gold particles inside them. Between 40 and 64% of these vesicles had both 5- and 16-nm particles. Maximum internalization occurred by 30-60 min, and by 2 hr the number of small and large particles on the B cell surface became constant or increased, respectively. Both kinds of particles moved from electron-lucent to electron-dense vesicles as the incubation time increased, although a portion of the anti-class II particles remained in electron-lucent vesicles. These data clearly show that labeled, cell surface class II is not selectively excluded from Ag-containing vesicles during Ag internalization. Thus, cointernalization of Ag and class II may represent a mechanism by which processed Ag meets class II.     

Class II MHC cDNAs in 15I5 B-congenic chickens

 cDNA was obtained from the bursae of Fabricius of chickens from six B-congenic lines developed at this laboratory and studied for expression of class II B-LB genes. Following cDNA amplification, cloning and sequencing, genes were assigned to B-LB loci based on characteristic DNA sequences, amino acid relatedness to characterized genes, and level of expression. Genes from the B-LBI, B-LBII, and B-LBVI loci were differentially expressed in chickens with the B ² , B ⁵ , B ¹³ , B ¹⁵ , or B ²¹ haplotypes. Chickens of all haplotypes expressed a B-LBII gene. Additional B-LB genes expressed included: B-LBI genes in the B ⁵ and B ¹⁹ haplotypes; a B-LBI/VI recombinant gene in the B ² haplotype; and a B-LBVI gene in the B ¹³ haplotype. The B-congenic lines have demonstrable differences in resistance to Marek’s disease (MD), and in responses to MD viral vaccines. This variability in disease resistance may be correlated with polymorphisms in the expressed B-LB genes, or with differential expression of genes at different loci. Received: 18 April 1997 / Revised: 29 September 1997     

Transport of Streptococcus pneumoniae capsular polysaccharide in MHC Class II tubules

Bacterial capsular polysaccharides are virulence factors and are considered T cell-independent antigens. However, the capsular polysaccharide Sp1 from Streptococcus pneumoniae serotype 1 has been shown to activate CD4(+) T cells in a major histocompatibility complex (MHC) class II-dependent manner. The mechanism of carbohydrate presentation to CD4(+) T cells is unknown. We show in live murine dendritic cells (DCs) that Sp1 translocates from lysosomal compartments to the plasma membrane in MHCII-positive tubules. Sp1 cell surface presentation results in reduction of self-peptide presentation without alteration of the MHCII self peptide repertoire. In DM-deficient mice, retrograde transport of Sp1/MHCII complexes resulting in T cell-dependent immune responses to the polysaccharide in vitro and in vivo is significantly reduced. The results demonstrate the capacity of a bacterial capsular polysaccharide antigen to use DC tubules as a vehicle for its transport as an MHCII/saccharide complex to the cell surface for the induction of T cell activation. Furthermore, retrograde transport requires the functional role of DM in self peptide-carbohydrate exchange. These observations open new opportunities for the design of vaccines against microbial encapsulated pathogens.     

Structural basis for the neutralization and specificity of Staphylococcal enterotoxin B against its MHC Class II binding site

Staphylococcal enterotoxin (SE) B is among the potent toxins produced by Staphylococcus aureus that cause toxic shock syndrome (TSS), which can result in multi-organ failure and death. Currently, neutralizing antibodies have been shown to be effective immunotherapeutic agents against this toxin, but the structural basis of the neutralizing mechanism is still unknown. In this study, we generated a neutralizing monoclonal antibody, 3E2, against SEB, and analyzed the crystal structure of the SEB-3E2 Fab complex. Crystallographic analysis suggested that the neutralizing epitope overlapped with the MHC II molecule binding site on SEB, and thus 3E2 could inhibit SEB function by preventing interaction with the MHC II molecule. Mutagenesis studies were done on SEB, as well as the related Staphylococcus aureus toxins SEA and SEC. These studies revealed that tyrosine (Y)46 and lysine (K)71 residues of SEB are essential to specific antibody–antigen recognition and neutralization. Substitution of Y at SEA glutamine (Q)49, which corresponds to SEB Y46, increased both 3E2’s binding to SEA in vitro and the neutralization of SEA in vivo. These results suggested that SEB Y46 is responsible for distinguishing SEB from SEA. These findings may be helpful for the development of antibody-based therapy for SEB-induced TSS.     

Structural basis for the neutralization and specificity of Staphylococcal enterotoxin B against its MHC Class II binding site

Staphylococcal enterotoxin (SE) B is among the potent toxins produced by Staphylococcus aureus that cause toxic shock syndrome (TSS), which can result in multi-organ failure and death. Currently, neutralizing antibodies have been shown to be effective immunotherapeutic agents against this toxin, but the structural basis of the neutralizing mechanism is still unknown. In this study, we generated a neutralizing monoclonal antibody, 3E2, against SEB, and analyzed the crystal structure of the SEB-3E2 Fab complex. Crystallographic analysis suggested that the neutralizing epitope overlapped with the MHC II molecule binding site on SEB, and thus 3E2 could inhibit SEB function by preventing interaction with the MHC II molecule. Mutagenesis studies were done on SEB, as well as the related Staphylococcus aureus toxins SEA and SEC. These studies revealed that tyrosine (Y)46 and lysine (K)71 residues of SEB are essential to specific antibody–antigen recognition and neutralization. Substitution of Y at SEA glutamine (Q)49, which corresponds to SEB Y46, increased both 3E2’s binding to SEA in vitro and the neutralization of SEA in vivo. These results suggested that SEB Y46 is responsible for distinguishing SEB from SEA. These findings may be helpful for the development of antibody-based therapy for SEB-induced TSS.     

Single-Molecule Motions of MHC Class II Rely on Bound Peptides

The major histocompatibility complex (MHC) class II protein can bind peptides of different lengths in the region outside the peptide-binding groove. Peptide-flanking residues (PFRs) contribute to the binding affinity of the peptide for MHC and change the immunogenicity of the peptide/MHC complex with regard to T cell receptor (TCR). The mechanisms underlying these phenomena are currently unknown. The molecular flexibility of the peptide/MHC complex may be an important determinant of the structures recognized by certain T cells. We used single-molecule x-ray analysis (diffracted x-ray tracking (DXT)) and fluorescence anisotropy to investigate these mechanisms. DXT enabled us to monitor the real-time Brownian motion of the peptide/MHC complex and revealed that peptides without PFRs undergo larger rotational motions than peptides with PFRs. Fluorescence anisotropy further revealed that peptides without PFRs exhibit slightly larger motions on the nanosecond timescale. These results demonstrate that peptides without PFRs undergo dynamic motions in the groove of MHC and consequently are able to assume diverse structures that can be recognized by T cells.     

Characterization of the Possible Roles for B Class MADS Box Genes in Regulation of Perianth Formation in Orchid

To investigate sepal/petal/lip formation in Oncidium Gower Ramsey, three paleoAPETALA3 genes, O. Gower Ramsey MADS box gene5 (OMADS5; clade 1), OMADS3 (clade 2), and OMADS9 (clade 3), and one PISTILLATA gene, OMADS8, were characterized. The OMADS8 and OMADS3 mRNAs were expressed in all four floral organs as well as in vegetative leaves. The OMADS9 mRNA was only strongly detected in petals and lips. The mRNA for OMADS5 was only strongly detected in sepals and petals and was significantly down-regulated in lip-like petals and lip-like sepals of peloric mutant flowers. This result revealed a possible negative role for OMADS5 in regulating lip formation. Yeast two-hybrid analysis indicated that OMADS5 formed homodimers and heterodimers with OMADS3 and OMADS9. OMADS8 only formed heterodimers with OMADS3, whereas OMADS3 and OMADS9 formed homodimers and heterodimers with each other. We proposed that sepal/petal/lip formation needs the presence of OMADS3/8 and/or OMADS9. The determination of the final organ identity for the sepal/petal/lip likely depended on the presence or absence of OMADS5. The presence of OMADS5 caused short sepal/petal formation. When OMADS5 was absent, cells could proliferate, resulting in the possible formation of large lips and the conversion of the sepal/petal into lips in peloric mutants. Further analysis indicated that only ectopic expression of OMADS8 but not OMADS5/9 caused the conversion of the sepal into an expanded petal-like structure in transgenic Arabidopsis (Arabidopsis thaliana) plants.The ABCDE model predicts the formation of any flower organ by the interaction of five classes of homeotic genes in plants (Yanofsky et al., 1990; Jack et al., 1992; Mandel et al., 1992; Goto and Meyerowitz, 1994; Jofuku et al., 1994; Pelaz et al., 2000, 2001; Theißen and Saedler, 2001; Pinyopich et al., 2003; Ditta et al., 2004; Jack, 2004). The A class genes control sepal formation. The A, B, and E class genes work together to regulate petal formation. The B, C, and E class genes control stamen formation. The C and E class genes work to regulate carpel formation, whereas the D class gene is involved in ovule development. MADS box genes seem to have a central role in flower development, because most ABCDE genes encode MADS box proteins (Coen and Meyerowitz, 1991; Weigel and Meyerowitz, 1994; Purugganan et al., 1995; Rounsley et al., 1995; Theißen and Saedler, 1995; Theißen et al., 2000; Theißen, 2001).The function of B group genes, such as APETALA3 (AP3) and PISTILLATA (PI), has been thought to have a major role in specifying petal and stamen development (Jack et al., 1992; Goto and Meyerowitz, 1994; Krizek and Meyerowitz, 1996; Kramer et al., 1998; Hernandez-Hernandez et al., 2007; Kanno et al., 2007; Whipple et al., 2007; Irish, 2009). In Arabidopsis (Arabidopsis thaliana), mutation in AP3 or PI caused identical phenotypes of second whorl petal conversion into a sepal structure and third flower whorl stamen into a carpel structure (Bowman et al., 1989; Jack et al., 1992; Goto and Meyerowitz, 1994). Similar homeotic conversions for petal and stamen were observed in the mutants of the AP3 and PI orthologs from a number of core eudicots such as Antirrhinum majus, Petunia hybrida, Gerbera hybrida, Solanum lycopersicum, and Nicotiana benthamiana (Sommer et al., 1990; Tröbner et al., 1992; Angenent et al., 1993; van der Krol et al., 1993; Yu et al., 1999; Liu et al., 2004; Vandenbussche et al., 2004; de Martino et al., 2006), from basal eudicot species such as Papaver somniferum and Aquilegia vulgaris (Drea et al., 2007; Kramer et al., 2007), as well as from monocot species such as Zea mays and Oryza sativa (Ambrose et al., 2000; Nagasawa et al., 2003; Prasad and Vijayraghavan, 2003; Yadav et al., 2007; Yao et al., 2008). This indicated that the function of the B class genes AP3 and PI is highly conserved during evolution.It has been thought that B group genes may have arisen from an ancestral gene through multiple gene duplication events (Doyle, 1994; Theißen et al., 1996, 2000; Purugganan, 1997; Kramer et al., 1998; Kramer and Irish, 1999; Lamb and Irish, 2003; Kim et al., 2004; Stellari et al., 2004; Zahn et al., 2005; Hernandez-Hernandez et al., 2007). In the gymnosperms, there was a single putative B class lineage that duplicated to generate the paleoAP3 and PI lineages in angiosperms (Kramer et al., 1998; Theißen et al., 2000; Irish, 2009). The paleoAP3 lineage is composed of AP3 orthologs identified in lower eudicots, magnolid dicots, and monocots (Kramer et al., 1998). Genes in this lineage contain the conserved paleoAP3- and PI-derived motifs in the C-terminal end of the proteins, which have been thought to be characteristics of the B class ancestral gene (Kramer et al., 1998; Tzeng and Yang, 2001; Hsu and Yang, 2002). The PI lineage is composed of PI orthologs that contain a highly conserved PI motif identified in most plant species (Kramer et al., 1998). Subsequently, there was a second duplication at the base of the core eudicots that produced the euAP3 and TM6 lineages, which have been subject to substantial sequence changes in eudicots during evolution (Kramer et al., 1998; Kramer and Irish, 1999). The paleoAP3 motif in the C-terminal end of the proteins was retained in the TM6 lineage and replaced by a conserved euAP3 motif in the euAP3 lineage of most eudicot species (Kramer et al., 1998). In addition, many lineage-specific duplications for paleoAP3 lineage have occurred in plants such as orchids (Hsu and Yang, 2002; Tsai et al., 2004; Kim et al., 2007; Mondragón-Palomino and Theißen, 2008, 2009; Mondragón-Palomino et al., 2009), Ranunculaceae, and Ranunculales (Kramer et al., 2003; Di Stilio et al., 2005; Shan et al., 2006; Kramer, 2009).Unlike the A or C class MADS box proteins, which form homodimers that regulate flower development, the ability of B class proteins to form homodimers has only been reported in gymnosperms and in the paleoAP3 and PI lineages of some monocots. For example, LMADS1 of the lily Lilium longiflorum (Tzeng and Yang, 2001), OMADS3 of the orchid Oncidium Gower Ramsey (Hsu and Yang, 2002), and PeMADS4 of the orchid Phalaenopsis equestris (Tsai et al., 2004) in the paleoAP3 lineage, LRGLOA and LRGLOB of the lily Lilium regale (Winter et al., 2002), TGGLO of the tulip Tulipa gesneriana (Kanno et al., 2003), and PeMADS6 of the orchid P. equestris (Tsai et al., 2005) in the PI lineage, and GGM2 of the gymnosperm Gnetum gnemon (Winter et al., 1999) were able to form homodimers that regulate flower development. Proteins in the euAP3 lineage and in most paleoAP3 lineages were not able to form homodimers and had to interact with PI to form heterodimers in order to regulate petal and stamen development in various plant species (Schwarz-Sommer et al., 1992; Tröbner et al., 1992; Riechmann et al., 1996; Moon et al., 1999; Winter et al., 2002; Kanno et al., 2003; Vandenbussche et al., 2004; Yao et al., 2008). In addition to forming dimers, AP3 and PI were able to interact with other MADS box proteins, such as SEPALLATA1 (SEP1), SEP2, and SEP3, to regulate petal and stamen development (Pelaz et al., 2000; Honma and Goto, 2001; Theißen and Saedler, 2001; Castillejo et al., 2005).Orchids are among the most important plants in the flower market around the world, and research on MADS box genes has been reported for several species of orchids during the past few years (Lu et al., 1993, 2007; Yu and Goh, 2000; Hsu and Yang, 2002; Yu et al., 2002; Hsu et al., 2003; Tsai et al., 2004, 2008; Xu et al., 2006; Guo et al., 2007; Kim et al., 2007; Chang et al., 2009). Unlike the flowers in eudicots, the nearly identical shape of the sepals and petals as well as the production of a unique lip in orchid flowers make them a very special plant species for the study of flower development. Four clades (1–4) of genes in the paleoAP3 lineage have been identified in several orchids (Hsu and Yang, 2002; Tsai et al., 2004; Kim et al., 2007; Mondragón-Palomino and Theißen, 2008, 2009; Mondragón-Palomino et al., 2009). Several works have described the possible interactions among these four clades of paleoAP3 genes and one PI gene that are involved in regulating the differentiation and formation of the sepal/petal/lip of orchids (Tsai et al., 2004; Kim et al., 2007; Mondragón-Palomino and Theißen, 2008, 2009). However, the exact mechanism that involves the orchid B class genes remains unclear and needs to be clarified by more experimental investigations.O. Gower Ramsey is a popular orchid with important economic value in cut flower markets. Only a few studies have been reported on the role of MADS box genes in regulating flower formation in this plant species (Hsu and Yang, 2002; Hsu et al., 2003; Chang et al., 2009). An AP3-like MADS gene that regulates both floral formation and initiation in transgenic Arabidopsis has been reported (Hsu and Yang, 2002). In addition, four AP1/AGAMOUS-LIKE9 (AGL9)-like MADS box genes have been characterized that show novel expression patterns and cause different effects on floral transition and formation in Arabidopsis (Hsu et al., 2003; Chang et al., 2009). Compared with other orchids, the production of a large and well-expanded lip and five small identical sepals/petals makes O. Gower Ramsey a special case for the study of the diverse functions of B class MADS box genes during evolution. Therefore, the isolation of more B class MADS box genes and further study of their roles in the regulation of perianth (sepal/petal/lip) formation during O. Gower Ramsey flower development are necessary. In addition to the clade 2 paleoAP3 gene OMADS3, which was previously characterized in our laboratory (Hsu and Yang, 2002), three more B class MADS box genes, OMADS5, OMADS8, and OMADS9, were characterized from O. Gower Ramsey in this study. Based on the different expression patterns and the protein interactions among these four orchid B class genes, we propose that the presence of OMADS3/8 and/or OMADS9 is required for sepal/petal/lip formation. Further sepal and petal formation at least requires the additional presence of OMADS5, whereas large lip formation was seen when OMADS5 expression was absent. Our results provide a new finding and information pertaining to the roles for orchid B class MADS box genes in the regulation of sepal/petal/lip formation.