Functional analysis and tissue-differential expression of four FAD2 genes in amphidiploid Brassica napus derived from Brassica rapa and Brassica oleracea

Fatty acid desaturase 2 (FAD2), which resides in the endoplasmic reticulum (ER), plays a crucial role in producing linoleic acid (18:2) through catalyzing the desaturation of oleic acid (18:1) by double bond formation at the delta 12 position. FAD2 catalyzes the first step needed for the production of polyunsaturated fatty acids found in the glycerolipids of cell membranes and the triacylglycerols in seeds. In this study, four FAD2 genes from amphidiploid Brassica napus genome were isolated by PCR amplification, with their enzymatic functions predicted by sequence analysis of the cDNAs. Fatty acid analysis of budding yeast transformed with each of the FAD2 genes showed that whereas BnFAD2-1, BnFAD2-2, and BnFAD2-4 are functional enzymes, and BnFAD2-3 is nonfunctional. The four FAD2 genes of B. napus originated from synthetic hybridization of its diploid progenitors Brassica rapa and Brassica oleracea, each of which has two FAD2 genes identical to those of B. napus. The BnFAD2-3 gene of B. napus, a nonfunctional pseudogene mutated by multiple nucleotide deletions and insertions, was inherited from B. rapa. All BnFAD2 isozymes except BnFAD2-3 localized to the ER. Nonfunctional BnFAD2-3 localized to the nucleus and chloroplasts. Four BnFAD2 genes can be classified on the basis of their expression patterns.     

Hemocyanin gene family evolution in spiders (Araneae), with implications for phylogenetic relationships and divergence times in the infraorder Mygalomorphae

Hemocyanins are multimeric copper-containing hemolymph proteins involved in oxygen binding and transport in all major arthropod lineages. Most arachnids have seven primary subunits (encoded by paralogous genes a–g), which combine to form a 24-mer (4 × 6) quaternary structure. Within some spider lineages, however, hemocyanin evolution has been a dynamic process with extensive paralog duplication and loss. We have obtained hemocyanin gene sequences from numerous representatives of the spider infraorders Mygalomorphae and Araneomorphae in order to infer the evolution of the hemocyanin gene family and estimate spider relationships using these conserved loci. Our hemocyanin gene tree is largely consistent with the previous hypotheses of paralog relationships based on immunological studies, but reveals some discrepancies in which paralog types have been lost or duplicated in specific spider lineages. Analyses of concatenated hemocyanin sequences resolved deep nodes in the spider phylogeny and recovered a number of clades that are supported by other molecular studies, particularly for mygalomorph taxa. The concatenated data set is also used to estimate dates of higher-level spider divergences and suggests that the diversification of extant mygalomorphs preceded that of extant araneomorphs. Spiders are diverse in behavior and respiratory morphology, and our results are beneficial for comparative analyses of spider respiration. Lastly, the conserved hemocyanin sequences allow for the inference of spider relationships and ancient divergence dates.     

Small inverted repeats drive mitochondrial genome evolution in Lake Baikal sponges

Demosponges, the largest and most diverse class in the phylum Porifera, possess mitochondrial DNA (mtDNA) markedly different from that in other animals. Although several studies investigated evolution of demosponge mtDNA among major lineages of the group, the changes within these groups remain largely unexplored. Recently we determined mitochondrial genomic sequence of the Lake Baikal sponge Lubomirskia baicalensis and described proliferation of small inverted repeats (hairpins) that occurred in it since the divergence between L. baicalensis and the most closely related cosmopolitan freshwater sponge Ephydatia muelleri. Here we report mitochondrial genomes of three additional species of Lake Baikal sponges: Swartschewskia papyracea, Rezinkovia echinata and Baikalospongia intermedia morpha profundalis (Demospongiae, Haplosclerida, Lubomirskiidae) and from a more distantly related freshwater sponge Corvomeyenia sp. (Demospongiae, Haplosclerida, Metaniidae). We use these additional sequences to explore mtDNA evolution in Baikalian sponges, paying particular attention to the variation in the rates of nucleotide substitutions and the distribution of hairpins, abundant in these genomes. We show that most of the changes in Lubomirskiidae mitochondrial genomes are due to insertion/deletion/duplication of these elements rather than single nucleotide substitutions. Thus inverted repeats can act as an important force in evolution of mitochondrial genome architecture and be a valuable marker for population- and species-level studies in this group. In addition, we infer (((Rezinkovia+Lubomirskia)+Swartschewskia)+Baikalospongia) phylogeny for the family Lubomirskiidae based on the analysis of mitochondrial coding sequences from freshwater sponges.     

Molecular evolution of two duplicated CDPK genes CPK7 and CPK12 in grass species: a case study in wheat (Triticum aestivum L.)

Gene duplication contributes to the expansion of gene families and subsequent functional diversification. Calcium-dependent protein kinases (CDPKs) are members of an important calcium sensor family involved in abiotic and biotic stress signaling in plants. We report here the molecular evolution and expression analysis of a pair of duplicated CDPK genes CPK7 and CPK12 that arose in the common ancestor of grass species. With higher nonsynonymous/synonymous ratios (dN/dS, or ω), CPK12 genes appear to diverge more rapidly than CPK7s, suggesting relaxed selection constraints on CPK12s. Sliding window analysis revealed increased dN and ω values at N-terminal regions and the calcium-binding EF hand loops. Likelihood analyses using various models in PAML 4.0 showed purifying selection on both CPK7 and CPK12 lineages. In addition to the divergence in cis-element combinations on their promoters, functional divergence of CPK7 and CPK12 genes was also observed in wheat where TaCPK7 was found to respond to drought (PEG), salt (NaCl), cold, and hydrogen peroxide (H(2)O(2)) while TaCPK12 responded only to the treatment of ABA, a feature that may complement or expand TaCPK7-mediated stress signaling networks of wheat. The contrasting expression patterns of CPK7 and CPK12 genes under stress conditions were also observed in rice, suggesting conservative functional evolution of these genes. Since no positive selection was detected between the two lineages, the divergence of CPK7 and CPK12 genes should be ascribed to subfunctionalization, rather than neofunctionalization. Thus, our work demonstrates another case of evolutionary employment of duplicated genes via subfunctionalization for better adaptation.     

Evolutionary analysis of three gibberellin oxidase genes in rice, Arabidopsis, and soybean

GAs are plant hormones that play fundamental roles in plant growth and development. GA2ox, GA3ox, and GA20ox are three key enzymes in GA biosynthesis. These enzymes belong to the 2OG-Fe (II) oxygenase superfamily and are independently encoded by different gene families. To date, genome-wide comparative analyses of GA oxidases in plant species have not been thoroughly carried out. In the present work, 61 GA oxidase family genes from rice (Oryza sativa), Arabidopsis, and soybean (Glycine max) were identified and a full study of these genes including phylogenetic tree construction, gene structure, gene family expansion and analysis of functional motifs was performed. Based on phylogeny, most of the GA oxidases were divided into four subgroups that reflected functional classifications. Intron/intron average length of GA oxidase genes in rice analysis revealed that GA oxidase genes in rice experienced substantial evolutionary divergence. Segmental duplication events were mainly found in soybean genome. However, in rice and Arabidopsis, no single expansion pattern exhibited dominance, indicating that GA oxidase genes from these species might have been subjected to a more complex evolutionary mechanism. In addition, special functional motifs were discovered in GA20ox, GA3ox, and GA2ox, which suggested that different functional motifs are associated with differences in protein function. Taken together our results suggest that GA oxidase family genes have undergone divergent evolutionary routes, especially at the monocot-dicot split, with dynamic evolution occurring in Arabidopsis thaliana and soybean.     

SINEs, evolution and genome structure in the opossum

Short INterspersed Elements (SINEs) are non-autonomous retrotransposons, usually between 100 and 500 base pairs (bp) in length, which are ubiquitous components of eukaryotic genomes. Their activity, distribution, and evolution can be highly informative on genomic structure and evolutionary processes. To determine recent activity, we amplified more than one hundred SINE1 loci in a panel of 43 M. domestica individuals derived from five diverse geographic locations. The SINE1 family has expanded recently enough that many loci were polymorphic, and the SINE1 insertion-based genetic distances among populations reflected geographic distance. Genome-wide comparisons of SINE1 densities and GC content revealed that high SINE1 density is associated with high GC content in a few long and many short spans. Young SINE1s, whether fixed or polymorphic, showed an unbiased GC content preference for insertion, indicating that the GC preference accumulates over long time periods, possibly in periodic bursts. SINE1 evolution is thus broadly similar to human Alu evolution, although it has an independent origin. High GC content adjacent to SINE1s is strongly correlated with bias towards higher AT to GC substitutions and lower GC to AT substitutions. This is consistent with biased gene conversion, and also indicates that like chickens, but unlike eutherian mammals, GC content heterogeneity (isochore structure) is reinforced by substitution processes in the M. domestica genome. Nevertheless, both high and low GC content regions are apparently headed towards lower GC content equilibria, possibly due to a relative shift to lower recombination rates in the recent Monodelphis ancestral lineage. Like eutherians, metatherian (marsupial) mammals have evolved high CpG substitution rates, but this is apparently a convergence in process rather than a shared ancestral state.     

The evolution and functional divergence of the beta-carotene oxygenase gene family in teleost fish—Exemplified by Atlantic salmon

In mammals, two carotenoid cleaving oxygenases are known; beta-carotene 15,15′-monooxygenase (BCMO1) and beta-carotene 9′,10′-oxygenase (BCO2). BCMO1 is a key enzyme in vitamin A synthesis by symmetrically cleaving beta-carotene into 2 molecules of all-trans-retinal, while BCO2 is responsible for asymmetric cleavage of a broader range of carotenoids. Here, we show that the Atlantic salmon beta-carotene oxygenase (bco) gene family contains 5 members, three bco2 and two bcmo1 paralogs. Using public sequence databases, multiple bco genes were also found in several additional teleost species. Phylogenetic analysis indicates that bco2a and bco2b originate from the teleost fish specific genome duplication (FSGD or 3R), while the third and more distant paralog, bco2 like, might stem from a prior duplication event in the teleost lineage. The two bcmo1 paralogs (bcmo1 and bcmo1 like) appear to be the result of an ancient duplication event that took place before the divergence of ray-finned (Actinopterygii) and lobe-finned fish (Sarcopterygii), with subsequent nonfunctionalization and loss of one Sarcopterygii paralog. Gene expression analysis of the bcmo1 and bco2 paralogs in Atlantic salmon reveals regulatory divergence with tissue specific expression profiles, suggesting that the beta-carotene oxygenase subtypes have evolved functional divergences. We suggest that teleost fish have evolved and maintained an extended repertoire of beta-carotene oxygenases compared to the investigated Sarcopterygii species, and hypothesize that the main driver behind this functional divergence is the exposure to a diverse set of carotenoids in the aquatic environment.     

Extremely slow rate of evolution in the HOX cluster revealed by comparison between Tanzanian and Indonesian coelacanths

Coelacanths are known as "living fossils" because their morphology has changed very little from that in the fossil record. To elucidate why coelacanths have evolved so slowly is thus of primary importance in evolutionary biology. In the present study, we determined the entire sequence of the HOX cluster of the Tanzanian coelacanth (Latimeria chalumnae) and compared it with that of the Indonesian coelacanth (L. menadoensis), which was available in the literature. The most intriguing result was the extremely small genetic divergence between the two coelacanths. The synonymous divergence of the HOX coding region between the two coelacanths was estimated to be 0.07%, which is ~11-fold smaller than that of human-chimp. When we applied the estimated divergence time of the two coelacanths of 6 million years ago (MYA) and 30 MYA, which were proposed in independent mitochondrial DNA analyses, the synonymous substitution rate of the coelacanth HOX cluster was estimated to be ~11-fold and 56-fold smaller than that of human-chimp, respectively. Thus, the present study implies that the reduction of the nucleotide substitution rate in coelacanth HOX genes may account for the conservation of coelacanth morphology during evolution.