Molecular evolution of the clustered <Emphasis Type="Italic">MIC</Emphasis>-<Emphasis Type="Italic">3</Emphasis> multigene family of <Emphasis Type="Italic">Gossypium</Emphasis> species

The Gossypium MIC-3 (Meloidogyne Induced Cotton-3) gene family is of great interest for molecular evolutionary studies because of its uniqueness to Gossypium species, multi-gene content, clustered localization, and root-knot nematode resistance-associated features. Molecular evolution of the MIC-3 gene family was studied in 15 tetraploid and diploid Gossypium genotypes that collectively represent seven phylogenetically distinct genomes. Synonymous (d_S) and non-synonymous (d_N) nucleotide substitution rates suggest that the second of the two exons of the MIC-3 genes has been under strong positive selection pressure, while the first exon has been under strong purifying selection to preserve function. Based on nucleotide substitution rates, we conclude that MIC-3 genes are evolving by a birth-and-death process and that a ‘gene amplification’ mechanism has helped to retain all duplicate copies, which best fits with the “bait and switch” model of R-gene evolution. The data indicate MIC-3 gene duplication events occurred at various rates, once per 1 million years (MY) in the allotetraploids, once per ~2 MY in the A/F genome clade, and once per ~8 MY in the D-genome clade. Variations in the MIC-3 gene family seem to reflect evolutionary selection for increased functional stability, while also expanding the capacity to develop novel “switch” pockets for responding to diverse pests and pathogens. Such evolutionary roles are congruent with the hypothesis that members of this unique resistance gene family provide fitness advantages in Gossypium.     

Molecular evolution of paclitaxel biosynthetic genes <Emphasis Type="Italic">TS</Emphasis> and <Emphasis Type="Italic">DBAT</Emphasis> of <Emphasis Type="Italic">Taxus</Emphasis> species

Evolutionary patterns of sequence divergence were analyzed in genes from the conifer genus Taxus (yew), encoding paclitaxel biosynthetic enzymes taxadiene synthase (TS) and 10-deacetylbaccatin III-10β-O-acetyltransferase (DBAT). N-terminal fragments of TS, full-length DBAT and internal transcribed spacer (ITS) were amplified from 15 closely related Taxus species and sequenced. Premature stop codons were not found in TS and DBAT sequences. Codon usage bias was not found, suggesting that synonymous mutations are selectively neutral. TS and DBAT gene trees are not consistent with the ITS tree, where species formed monophyletic clades. In fact, for both genes, alleles were sometimes shared across species and parallel amino acid substitutions were identified. While both TS and DBAT are, overall, under purifying selection, we identified a number of amino acids of TS under positive selection based on inference using maximum likelihood models. Positively selected amino acids in the N-terminal region of TS suggest that this region might be more important for enzyme function than previously thought. Moreover, we identify lineages with significantly elevated rates of amino acid substitution using a genetic algorithm. These findings demonstrate that the pattern of adaptive paclitaxel biosynthetic enzyme evolution can be documented between closely related Taxus species, where species-specific taxane metabolism has evolved recently.     

Genomewide analysis of <Emphasis Type="Italic">ABCB</Emphasis>s with a focus on <Emphasis Type="Italic">ABCB1</Emphasis> and <Emphasis Type="Italic">ABCB19</Emphasis> in <Emphasis Type="Italic">Malus domestica</Emphasis>

The B subfamily of ATP-binding cassette (ABC) proteins (ABCB) plays a vital role in auxin efflux. However, no systematic study has been done in apple. In this study, we performed genomewide identification and expression analyses of the ABCB family in Malus domestica for the first time. We identified a total of 25 apple ABCBs that were divided into three clusters based on the phylogenetic analysis. Most ABCBs within the same cluster demonstrated a similar exon–intron organization. Additionally, the digital expression profiles of ABCB genes shed light on their functional divergence. ABCB1 and ABCB19 are two well-studied auxin efflux carrier genes, and we found that their expression levels are higher in young shoots of M106 than in young shoots of M9. Since young shoots are the main source of auxin synthesis and auxin efflux involves in tree height control. This suggests that ABCB1 and ABCB19 may also take a part in the auxin efflux and tree height control in apple.     

Identification and expression of amphioxus <Emphasis Type="Italic">AmphiSmad1/5/8</Emphasis> and <Emphasis Type="Italic">AmphiSmad4</Emphasis>

Smad family proteins are identified as intracellular signal mediators of the TGF-β superfamily. In this study, we identified two novel members of the Smad family, termed as AmphiSmad1/5/8 and AmphiSmad4, from Chinese amphioxus. Both AmphiSmad1/5/8 and AmphiSmad4 showed a typical domain structure of Smad proteins consisting of conserved MH1 and MH2 domains. Phylogenetic analysis placed AmphiSmad1/5/8 in the Smad1, 5 and 8 subgroup of the R-Smad subfamily, and AmphiSmad4 in the Co-Smad subfamily. The spatial and temporal gene expression patterns of AmphiSmad1/5/8 and AmphiSmad4 showed that they may be involved in the embryonic development of notochord, myotome and alimentary canal, and may help to establish the specification of dorsal-ventral axis of amphioxus. Moreover, AmphiSmad1/5/8 and AmphiSmad4 showed extensive distribution in all adult tissues examined, suggesting that these two genes may play important roles in the morphogenesis of a variety of tissues especially notochord and gonad.     

Molecular dissection of <Emphasis Type="Italic">Bombyx mori</Emphasis> nucleopolyhedrovirus <Emphasis Type="Italic">orf8</Emphasis> gene

Viruses including baculoviruses are obligatory parasites, as their genomes do not encode all the proteins required for replication. Therefore, viruses have evolved to exploit the behavior and the physiology of their hosts and often coevolved with their hosts over millions of years. Recent comparative analyses of complete genome sequences of baculoviruses revealed the patterns of gene acquisitions and losses that have occurred during baculovirus evolution. In addition, knowledge of virus genes has also provided understanding of the mechanism of baculovirus infection including replication, species-specific virulence and host range. The Bm8 gene of Bombyx mori nucleopolyhedrovirus (NPV) and its homologues are found only in group I NPV genomes. The Autographa californica NPV Ac16 gene is a homologue of Bm8 and, encodes a viral structural protein. It has been shown that Bm8/Ac16 interacts with baculoviral and cellular proteins. Bm8/Ac16 interacts with baculoviral IE1 that is facilitated by coiled coil domains, and the interaction with IE1 is important for Bm8 function. Ac16 also forms a complex with viral FP25 and cellular actin and associates with membranes via palmitoylation. These data suggested that this gene family encodes a multifunctional protein that accomplishes specific needs of group I NPVs.      

WRKY gene family evolution in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

The Arabidopsis thaliana WRKY proteins are characterized by a sequence of 60 amino acids including WRKY domain. It is well established that these proteins are involved in the regulation of various physiological programs unique to plants including pathogen defense, senescence and response to environmental stresses, which attracts attention of the scientific community as to how this family might have evolved. We tried to satisfy this curiosity and analyze reasons for duplications of these gene sequences leading to their diversified gene actions. The WRKY sequences available in Arabidopsis thaliana were used to evaluate selection pressure following duplication events. A phylogenetic tree was constructed and the WRKY family was divided into five sub-families. After that, tests were conducted to decide whether positive or purified selection played key role in these events. Our results suggest that purifying selection played major role during the evolution of this family. Some amino acid changes were also detected in specific branches of phylogeny suggesting that relaxed constraints might also have contributed to functional divergence among sub-families. Sites relaxed from purifying selection were identified and mapped onto the structural and functional regions of the WRKY1 protein. These analyses will enhance our understanding of the precise role played by natural selection to create functional diversity in WRKY family.     

A novel family of mitochondrial proteins is represented by the <Emphasis Type="Italic">Drosophila</Emphasis> genes <Emphasis Type="Italic">slmo</Emphasis>, <Emphasis Type="Italic">preli-like</Emphasis> and <Emphasis Type="Italic">real-time</Emphasis>

Mitochondria play essential roles in development and disease. The characterisation of mitochondrial proteins is therefore of particular importance. The slowmo (slmo) gene of Drosophila melanogaster has been shown to encode a novel type of mitochondrial protein, and is essential in the developing central nervous system. The Slmo protein contains a conserved PRELI/MSF1p domain, found in proteins from a wide variety of eukaryotic organisms. However, the function of the proteins of this family is currently unknown. In this study, the evolutionary relationships between members of the PRELI/MSF1p family are described, and we present the first analysis of two novel Drosophila genes predicted to encode proteins of this type. The first of these, preli-like (prel), is expressed ubiquitously during embryonic development, whilst the second, real-time (retm), is expressed dynamically in the developing gut and central nervous system. retm encodes a member of a novel conserved subclass of larger PRELI/MSF1p domain proteins, which also contain the CRAL-TRIO motif thought to mediate the transport of small hydrophobic ligands. Here we provide evidence that, like Slmo, both the Prel and Retm proteins are localised to the mitochondria, indicating that the function of the PRELI/MSF1p domain is specific to this organelle.Edited by P. Simpson     

Organization of phenylalanine ammonia lyase (<Emphasis Type="Italic">PAL</Emphasis>), acidic <Emphasis Type="Italic">PR-5</Emphasis> and osmotin-like (<Emphasis Type="Italic">OSM</Emphasis>) defence-response gene families in the potato genome

Defence-response (DR) genes are candidates for the genetic functions underlying quantitative resistance to plant pathogens. The organization of three DR gene families encoding phenylalanine ammonia lyase (PAL), acidic PR-(pathogenesis-related) protein 5, and basic PR-5, or osmotin-like (OSM), proteins was studied in the potato genome. A bacterial artificial chromosome (BAC) library containing ~50,000 clones was constructed from high-molecular weight genomic DNA of the diploid potato clone PD59, a hybrid between Solanum tuberosum and S. phureja. BAC clones carrying one or more copies of the DR genes were identified and characterized by Southern hybridization, sequence analysis and genetic mapping. PAL, acidic PR-5 and OSM (basic PR-5) genes were all organized into gene families of varying complexity. The PAL gene family consisted of at least 16 members, several of which were physically linked. Four acidic PR-5 homologous were localized to a 45-kb segment on potato chromosome XII. One of these, PR-5/319, codes for the acidic thaumatin-like protein C found in intercellular fluids of potato. Nine OSM genes were organized at two loci: eight form a 90-kb cluster on chromosome VIII, and a single gene was found on chromosome XI. The topology of a phylogenetic tree based on PR-5 and OSM protein sequences from Solanaceae suggests a mode of evolution for these gene families. The results will form the basis for further studies on the potential role of these defence-related loci in quantitative resistance to pathogens.     

Functional analysis of <Emphasis Type="Italic">Xa3/Xa26</Emphasis> family members in rice resistance to <Emphasis Type="Italic">Xanthomonas</Emphasis><Emphasis Type="Italic">oryzae</Emphasis> pv. <Emphasis Type="Italic">oryzae</Emphasis>

Plant disease resistant (R) genes are frequently clustered in the genome. The diversity of members in a complex R-gene family may provide variation in resistance specificity. Rice Xa3/Xa26, conferring resistance to Xanthomonas oryzae pv. oryzae (Xoo) encodes a leucine-rich repeat (LRR) receptor kinase-type protein and belongs to a multigene family, consisting of Xa3/Xa26, MRKa, MRKc and MRKd in rice cultivar Minghui 63. MRKa and MRKc are intact genes, while MRKd is a pseudogene. Complementary analyses showed that MRKa and MRKc could not mediate resistance to Xoo when regulated by their native promoters, but MRKa not MRKc conferred partial resistance to Xoo when regulated by a strong constitutive promoter. Plants carrying truncated XA3/XA26, which lacked the kinase domain, were compromised in their resistance to Xoo. However, the kinase domain of MRKa could partially restore the function of the truncated XA3/XA26 in resistance. MRKa and MRKc showed similar expression pattern as Xa3/Xa26, which expressed only in the vascular systems of different tissues. The expressional characteristic of MRKa and MRKc perfectly fits the function of genes conferring resistance to Xoo, a vascular pathogen. These results suggest that although MRKa and MRKc cannot mediate bacterial blight resistance nowadays, they may be once effective genes for Xoo resistance. Their expressional characteristic and sequence similarity to Xa3/Xa26 will provide templates for generating novel recognition specificity to face the evolution of Xoo. In addition, both LRR and kinase domains encoded by Xa3/Xa26 and MRKa are the functional determinants and MRKa-mediated resistance is dosage-dependent. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comparative genomic analysis of CIPK gene family in <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Populus</Emphasis>

Calcium serves as a second messenger in various signal transduction pathways in plants. CBL-interacting protein kinases (CIPKs), which have a variety of functions, are involved in calcium signal transduction. Previous, the studies on CIPK family members focused on Arabidopsis and rice. Here, we present a comparative genomic analysis of the CIPK gene family in Arabidopsis and poplar, a model tree species. Twenty-seven potential CIPKs were identified from poplar using genome-wide analysis. Like the CIPK gene family from Arabidopsis, CIPK genes from poplar were also divided into intron-free and intron-harboring groups. In the intron-harboring group, the intron distribution of CIPKs is rather conserved during the genome evolutionary process. Many homologous gene pairs were found in the CIPK gene family, indicating duplication events might contribute to the amplification of this gene family. The phylogenetic comparison of CIPKs in combination with intron distribution analysis revealed that CIPK genes from both Arabidopsis and poplar might have an ancient origin, which formed earlier than the separation of these two eudicot species. Our genomic and bioinformatic analysis will provide an important foundation for further functional dissection of the CBL-CIPK signaling network in poplars. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evolution of the <Emphasis Type="Italic">Groucho</Emphasis>/<Emphasis Type="Italic">Tle</Emphasis> gene family: gene organization and duplication events

The Groucho/Tle family of corepressor proteins has important roles in development and in adult tissue in both Protostomes and Deuterostomes. In Drosophila, a single member of this family has been identified. Unlike in Protostomes, most Deutrostomes contain more than two full-length Tle genes. In this study, I analyse the genomic organization and phylogenetic relationship between the long and short forms of the Groucho/Tle family members in Chordata. The genomic location and sequence similarities suggest that Aes/Grg5 and Tle6/Grg6 arose from duplication of the Tle2 gene; each evolved independently and acquired new functions as negative regulators of the other Tle proteins. Based on these data, a model for Groucho/Tle gene evolution is proposed. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evolution of <Emphasis Type="Italic">Hox3</Emphasis> and <Emphasis Type="Italic">ftz</Emphasis> in arthropods: insights from the crustacean <Emphasis Type="Italic">Daphnia pulex</Emphasis>

The Drosophila melanogaster genes zerknüllt (zen) and fushi tarazu (ftz) are members of the Hox gene family whose roles have changed significantly in the insect lineage and thus provide an opportunity to study the mechanisms underlying the functional evolution of Hox proteins. We have studied the expression of orthologs of zen (DpuHox3) and ftz (Dpuftz) in the crustacean Daphnia pulex (Branchiopoda), both of which show a dynamic expression pattern. DpuHox3 is expressed in a complex pattern in early embryogenesis, with the most anterior boundary of expression lying at the anterior limit of the second antennal segment as well as a ring of expression around the embryo. In later embryos, DpuHox3 expression is restricted to the mesoderm of mandibular limb buds. Dpuftz is first expressed in a ring around the embryo following the posterior limit of the mandibular segment. Later, Dpuftz is restricted to the posterior part of the mandibular segment. This is the first report of expression of a Hox3 ortholog in a crustacean, and together with Dpuftz data, the results presented here show that Hox3 and ftz have retained a Hox-like expression pattern in crustaceans. This is in accordance with the proposed model of Hox3 and ftz evolution in arthropods and allows a more precise pinpointing of the loss of ftz “Hox-like behaviour”: in the lineage between the Branchiopoda and the basal insect Thysanura.     

The evolutionary history of PDR in <Emphasis Type="Italic">Brachypodium</Emphasis><Emphasis Type="Italic">distachyon</Emphasis> polyploids

The ATP-binding cassette transporter genes include the pleiotropic drug resistance (PDR) family found only in fungi and plants. These transporters transport toxic compounds across biological membranes. Here, we investigated the evolution of the PDR1 gene in Brachypodium distachyon, a widely distributed temperate grass species that belongs to the Poaceae (Gramineae) family, which also contains the domesticated cereal crops. Because this species has multiple ploidy levels, investigating PDR1 evolution in B. distachyon will offer insights into the formation and evolution of polyploidy. From 23 B. distachyon ecotypes, 39 PDR1 homologs were identified. All ecotypes had either one or two PDR1 copies. Based on restriction site analysis, the PDR1 homologs were classified as E or H type. All but one diploid and tetraploid ecotypes had only a single H type PDR1. All but one hexaploid ecotypes had both an E and a H type PDR1. Phylogenetic analysis revealed that each type formed a well-supported cluster. The two PDR1 types appeared to evolve differently. These different evolutionary patterns could indicate a difference in age between the two types or might indicate different mutation rates or selection pressures on the two types. The phylogenetic analysis also revealed that the hexaploid ecotypes shared a genomic origin for their E type PDR1, but there were multiple origins for hexaploid H type PDR1 homologs. Overall, the results suggest that tetraploid and hexaploid might be misnomers in B. distachyon and suggest a complex polyploidization history during B. distachyon evolution.     

Characterization of multiple <Emphasis Type="Italic">CYP9A</Emphasis> genes in the silkworm, <Emphasis Type="Italic">Bombyx mori</Emphasis>

Based on the advances in the silkworm genome project, a new genome-wide analysis of cytochrome P450 genes was performed, focusing mainly on gene duplication. All four CYP9A subfamily members from the silkworm, Bombyx mori, were cloned by RT-PCR and designated CYP9A19–CYP9A22 by the P450 Nomenclature Committee. They each contain an open reading frame of 1,593 bp in length and encode a putative polypeptide of 531 amino acids. Both nucleic acid and amino acid sequences share very high identities with one another. The typical motifs of insect cytochrome P450, including the heme-binding region, helix-C, helix-I, helix-K, and PERF, show high sequence conservation among the multiple proteins. Alignment with their cDNA sequences revealed that these paralogues share identical gene structures, each comprising ten exons and nine introns of variable sizes. The locations of their introns (all nine introns follow the GT–AG rule) are absolutely conserved. CYP9A19, CYP9A20, and CYP9A21 form a tandem cluster on chromosome 17, whereas CYP9A22 is separated from the cluster by four tandem alcohol-dehydrogenase-like genes. Their phylogenetic relationships and structural comparisons indicated that these paralogues arose as the results of gene duplication events. RT-PCR detected their mRNAs in different “first line of defense” tissues, as well as in several other organs, suggesting diverse functions. Tissue-selective expression also indicates their functional divergence. The identified CYP9A genes have not yet been found outside the Lepidoptera, and are probably unique to the Lepidoptera. They show high sequence and structural similarities to each other, indicating that the Lepidoptera-specific P450s may be of functional importance. This analysis constitutes the first report of the clustering, spatial organization, and functional divergence of P450 in the silkworm.     

<Emphasis Type="Italic">SnRK2</Emphasis> Homologs in <Emphasis Type="Italic">Gossypium</Emphasis> and <Emphasis Type="Italic">GhSnRK2.6</Emphasis> Improved Salt Tolerance in Transgenic Upland Cotton and <Emphasis Type="Italic">Arabidopsis</Emphasis>

SnRK2s are a large family of plant-specific protein kinases, which play important roles in multiple abiotic stress responses in various plant species. But the family in Gossypium has not been well studied. Here, we identified 13, 10, and 13 members of the SnRK2 family from Gossypium raimondii, Gossypium arboreum, and Gossypium hirsutum, respectively, and analyzed the locations of SnRK2 homologs in chromosomes based on genome data of cotton species. Phylogenetic tree analysis of SnRK2 proteins showed that these families were classified into three groups. All SnRK2 genes were comprised of nine exons and eight introns, and the exon distributions and the intron phase of homolog genes among different cotton species were analogous. Moreover, GhSnRK2.6 was overexpressed in Arabidopsis and upland cotton, respectively. Under salt treatment, overexpressed Arabidopsis could maintain higher biomass accumulation than wild-type plants, and GhSnRK2.6 overexpression in cotton exhibited higher germination rate than the control. So, the gene GhSnRK2.6 could be utilized in cotton breeding for salt tolerance.     

Functional constraint and divergence in the G protein family in <Emphasis Type="Italic">Caenorhabditis elegans</Emphasis> and <Emphasis Type="Italic">Caenorhabditis briggsae</Emphasis>

Part of the challenge of the post-genomic world is to identify functional elements within the wide array of information generated by genome sequencing. Although cross-species comparisons and investigation of rates of sequence divergence are an efficient approach, the relationship between sequence divergence and functional conservation is not clear. Here, we use a comparative approach to examine questions of evolutionary rates and conserved function within the guanine nucleotide-binding protein (G protein) gene family in nematodes of the genus Caenorhabditis. In particular, we show that, in cases where the Caenorhabditis elegans ortholog shows a loss-of-function phenotype, G protein genes of C. elegans and Caenorhabditis briggsae diverge on average three times more slowly than G protein genes that do not exhibit any phenotype when mutated in C. elegans, suggesting that genes with loss of function phenotypes are subject to stronger selective constraints in relation to their function in both species. Our results also indicate that selection is as strong on G proteins involved in environmental perception as it is on those controlling other important processes. Finally, using phylogenetic footprinting, we identify a conserved non-coding motif present in multiple copies in the genomes of four species of Caenorhabditis. The presence of this motif in the same intron in the gpa-1 genes of C. elegans, C. briggsae and Caenorhabditis remanei suggests that it plays a role in the regulation of gpa-1, as well as other loci.Electronic Supplementary Material Supplementary material is available for this article at