OsLPR5 Encoding Ferroxidase Positively Regulates the Tolerance to Salt Stress in Rice.

Salinity is a major abiotic stress that harms rice growth and productivity. Low phosphate roots (LPRs) play a central role in Pi deficiency-mediated inhibition of primary root growth and have ferroxidase activity. However, the function of LPRs in salt stress response and tolerance in plants remains largely unknown. Here, we reported that the OsLPR5 was induced by NaCl stress and positively regulates the tolerance to salt stress in rice. Under NaCl stress, overexpression of OsLPR5 led to increased ferroxidase activity, more green leaves, higher levels of chlorophyll and lower MDA contents compared with the WT. In addition, OsLPR5 could promote the accumulation of cell osmotic adjustment substances and promote ROS-scavenging enzyme activities. Conversely, the mutant lpr5 had a lower ferroxidase activity and suffered severe damage under salt stress. Moreover, knock out of OsLPR5 caused excessive Na+ levels and Na+/K+ ratios. Taken together, our results exemplify a new molecular link between ferroxidase and salt stress tolerance in rice.

Deleterious Variants in Asian Rice and the Potential Cost of Domestication.

Many SNPs are predicted to encode deleterious amino acid variants. These slightly deleterious mutations can provide unique insights into population history, the dynamics of selection, and the genetic bases of phenotypes. This is especially true for domesticated species, where a history of bottlenecks and selection may affect the frequency of deleterious variants and signal a "cost of domestication". Here, we investigated the numbers and frequencies of deleterious variants in Asian rice (Oryza sativa), focusing on two varieties (japonica and indica) and their wild relative (O. rufipogon). We investigated three signals of a potential cost of domestication in Asian rice relative to O. rufipogon: an increase in the frequency of deleterious SNPs (dSNPs), an enrichment of dSNPs compared with synonymous SNPs (sSNPs), and an increased number of deleterious variants. We found evidence for all three signals, and domesticated individuals contained ∼3-4% more deleterious alleles than wild individuals. Deleterious variants were enriched within low recombination regions of the genome and experienced frequency increases similar to sSNPs within regions of putative selective sweeps. A characteristic feature of rice domestication was a shift in mating system from outcrossing to predominantly selfing. Forward simulations suggest that this shift in mating system may have been the dominant factor in shaping both deleterious and neutral diversity in rice.

Mutational bias is the driving force for shaping the synonymous codon usage pattern of alternatively spliced genes in rice (Oryza sativa L.).

Alternative splicing plays important roles in diverse aspects of plant development, metabolism, and stress responses. However, the regulatory mechanisms of alternative splicing of genes still remain incompletely elucidated, especially in plants. In this study, the synonymous codon usage pattern of alternatively spliced (AS) genes in rice was firstly explored using the combination of correspondence analysis (CA), internal CA, correlation and ANOVA analyses. The results show that alternatively and non-alternatively spliced (non-AS) genes have similar tendency for overall codon usage, but exhibit significant difference in 58 out of 64 codons. AS and non-AS genes are both under strong purifying selection, but the former ones have significant lower mutation rate and are prone to be enriched towards the chromosomal ends. In the group of AS genes, the variability in synonymous codon usage between genes is mainly due to the variations in GC content, CDS length, as well as gene functions. Mutational bias that accounts for 25.85 % of the total codon usage variability plays a major role in shaping the codon usage pattern of AS genes. In contrast, no obvious evidence is found for the contributions of translational selection, AS types, the conservation of AS events, and numbers of AS variants to the codon usage divergence between AS genes. These findings may be useful for further understanding the mechanisms of origination, differentiation and regulation of alternatively spliced genes in plants.

Genome-wide identification and evolutionary analysis of positively selected miRNA genes in domesticated rice.

The next-generation sequencing of tens to hundreds of plant genotypes made the uncovering of miRNA genes evolution available at the genome-wide level. Using the combinations of population genetics and evolutionary biology approaches, we have identified 21 miRNA loci having significant negative Tajima's D and Fu and Li's D* and F* values, of which 14 miRNAs (ps-miRNAs) showing clear signatures of positive selection in domesticated rice. The average sequence diversity (π) of the 21 miRNAs in cultivated rice is only 13.8 % of that in their wild progenitors. Interestingly, protein-coding genes immediately flanking these ps-miRNAs are apparently under weaker selective constraints. Totally, the 21 miRNAs are predicted to target 68 mRNA genes, of which 12 targets are estimated to have endured positive selection during rice evolution. In addition, the expression pattern and potential biological functions of ps-miRNAs targets are further investigated by searching published micro-array data and different mutant databases, respectively. We conclude that miRNAs, like protein-coding genes, should be crucial for driving rice evolution. These analyses may deepen our understanding on the miRNA genes evolution and functions during rice domestication.

Suppression of OsSAUR2 gene expression immobilizes soil arsenic bioavailability by modulating root exudation and rhizosphere microbial assembly in rice.

One of the factors influencing the behavior of arsenic (As) in environment is microbial-mediated As transformation. However, the detailed regulatory role of gene expression on the changes of root exudation, rhizosphere microorganisms, and soil As occurrence forms remains unclear. In this study, we evidence that loss-of-function of OsSAUR2 gene, a member of the SMALL AUXIN-UP RNA family in rice, results in significantly higher As uptake in roots but greatly lower As accumulation in grains via affecting the expression of OsLsi1, OsLsi2 in roots and OsABCC1 in stems. Further, the alteration of OsSAUR2 expression extensively affects the metabolomic of root exudation, and thereby leading to the variations in the composition of rhizosphere microbial communities in rice. The microbial community in the rhizosphere of Ossaur2 plants strongly immobilizes the occurrence forms of As in soil. Interestingly, Homovanillic acid (HA) and 3-Coumaric acid (CA), two differential metabolites screened from root exudation, can facilitate soil iron reduction, enhance As bioavailability, and stimulate As uptake and accumulation in rice. These findings add our further understanding in the relationship of OsSAUR2 expression with the release of root exudation and rhizosphere microbial assembly under As stress in rice, and provide potential rice genetic resources and root exudation in phytoremediation of As-contaminated paddy soil.

Mutational Bias and Natural Selection Driving the Synonymous Codon Usage of Single-Exon Genes in Rice (Oryza sativa L.).

The relative abundance of single-exon genes (SEGs) in higher plants is perplexing. Uncovering the synonymous codon usage pattern of SEGs will benefit for further understanding their underlying evolutionary mechanism in plants. Using internal correspondence analysis (ICA), we reveal a significant difference in synonymous codon usage between SEGs and multiple-exon genes (MEGs) in rice. But the effect is weak, accounting for only 2.61% of the total codon usage variability. SEGs and MEGs contain remarkably different base compositions, and are under clearly differential selective constraints, with the former having higher GC content, and evolving relatively faster during evolution. In the group of SEGs, the variability in synonymous codon usage among genes is partially due to the variations in GC content, gene function, and gene expression level, which accounts for 22.03%, 5.99%, and 3.32% of the total codon usage variability, respectively. Therefore, mutational bias and natural selection should work on affecting the synonymous codon usage of SEGs in rice. These findings may deepen our knowledge for the mechanisms of origination, differentiation and regulation of SEGs in plants.

De novo characterization of the Chinese fir (Cunninghamia lanceolata) transcriptome and analysis of candidate genes involved in cellulose and lignin biosynthesis.

BACKGROUND: Chinese fir (Cunninghamia lanceolata) is an important timber species that accounts for 20-30% of the total commercial timber production in China. However, the available genomic information of Chinese fir is limited, and this severely encumbers functional genomic analysis and molecular breeding in Chinese fir. Recently, major advances in transcriptome sequencing have provided fast and cost-effective approaches to generate large expression datasets that have proven to be powerful tools to profile the transcriptomes of non-model organisms with undetermined genomes. RESULTS: In this study, the transcriptomes of nine tissues from Chinese fir were analyzed using the Illumina HiSeq™ 2000 sequencing platform. Approximately 40 million paired-end reads were obtained, generating 3.62 gigabase pairs of sequencing data. These reads were assembled into 83,248 unique sequences (i.e. Unigenes) with an average length of 449 bp, amounting to 37.40 Mb. A total of 73,779 Unigenes were supported by more than 5 reads, 42,663 (57.83%) had homologs in the NCBI non-redundant and Swiss-Prot protein databases, corresponding to 27,224 unique protein entries. Of these Unigenes, 16,750 were assigned to Gene Ontology classes, and 14,877 were clustered into orthologous groups. A total of 21,689 (29.40%) were mapped to 119 pathways by BLAST comparison against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The majority of the genes encoding the enzymes in the biosynthetic pathways of cellulose and lignin were identified in the Unigene dataset by targeted searches of their annotations. And a number of candidate Chinese fir genes in the two metabolic pathways were discovered firstly. Eighteen genes related to cellulose and lignin biosynthesis were cloned for experimental validating of transcriptome data. Overall 49 Unigenes, covering different regions of these selected genes, were found by alignment. Their expression patterns in different tissues were analyzed by qRT-PCR to explore their putative functions. CONCLUSIONS: A substantial fraction of transcript sequences was obtained from the deep sequencing of Chinese fir. The assembled Unigene dataset was used to discover candidate genes of cellulose and lignin biosynthesis. This transcriptome dataset will provide a comprehensive sequence resource for molecular genetics research of C. lanceolata.

Novel miRNAs in the control of arsenite levels in rice.

In a species, the miRNA repertoire comprises many lineage- or species-specific miRNAs. Using deep sequencing, the whole transcriptome of seedling roots was studied and a total of 18 novel miRNAs in indica rice Minghui 86 were identified. Among these novel miRNAs six were up-regulated and two were down-regulated under arsenite stress. Quantitative real-time RT-PCR analysis revealed that miR6254 was predominantly expressed in roots relative to other tissues, whereas miR6250 and miR169i-3p were constitutively expressed in all tissues examined, with high abundance in roots (miR6250) or leaves (miR169i-3p). The miR6250 and miR169i-3p miRNAs also exhibited distinct expression patterns in rice cultivars Minghui 86 and Nipponbare at different time points after arsenite treatment. The predicted targets for these miRNAs included some protein kinases, DNA, or ATP-binding proteins. Besides arsenite, the expression of targets of miR6250 and miR6254 was also up- or down-regulated in response to abiotic environmental stresses, indicative of their involvement in regulation of plant adaptation. Three types of cis-elements involved in hormone, light, and stress response were found to occur frequently in the promoter regions. Interestingly, miR6254 was originally characterized to be an exonic miRNA located in the exon of AK101391, supporting the notion that miRNAs may also originate from some exons in plants.

Analyses of the sucrose synthase gene family in cotton: structure, phylogeny and expression patterns.

BACKGROUND: In plants, sucrose synthase (Sus) is widely considered as a key enzyme involved in sucrose metabolism. Several paralogous genes encoding different isozymes of Sus have been identified and characterized in multiple plant genomes, while limited information of Sus genes is available to date for cotton. RESULTS: Here, we report the molecular cloning, structural organization, phylogenetic evolution and expression profiles of seven Sus genes (GaSus1 to 7) identified from diploid fiber cotton (Gossypium arboreum). Comparisons between cDNA and genomic sequences revealed that the cotton GaSus genes were interrupted by multiple introns. Comparative screening of introns in homologous genes demonstrated that the number and position of Sus introns are highly conserved among Sus genes in cotton and other more distantly related plant species. Phylogenetic analysis showed that GaSus1, GaSus2, GaSus3, GaSus4 and GaSus5 could be clustered together into a dicot Sus group, while GaSus6 and GaSus7 were separated evenly into other two groups, with members from both dicot and monocot species. Expression profiles analyses of the seven Sus genes indicated that except GaSus2, of which the transcripts was undetectable in all tissues examined, and GaSus7, which was only expressed in stem and petal, the other five paralogues were differentially expressed in a wide ranges of tissues, and showed development-dependent expression profiles in cotton fiber cells. CONCLUSIONS: This is a comprehensive study of the Sus gene family in cotton plant. The results presented in this work provide new insights into the evolutionary conservation and sub-functional divergence of the cotton Sus gene family in response to cotton fiber growth and development.

The miR528-D3 Module Regulates Plant Height in Rice by Modulating the Gibberellin and Abscisic Acid Metabolisms.

Plant height, as one of the important agronomic traits of rice, is closely related to yield. In recent years, plant height-related genes have been characterized and identified, among which the DWARF3 (D3) gene is one of the target genes of miR528, and regulates rice plant height and tillering mainly by affecting strigolactone (SL) signal transduction. However, it remains unknown whether the miR528 and D3 interaction functions in controlling plant height, and the underlying regulatory mechanism in rice. In this study, we found that the plant height, internode length, and cell length of internodes of d3 mutants and miR528-overexpressing (OE-miR528) lines were greatly shorter than WT, D3-overexpressing (OE-D3), and miR528 target mimicry (OE-MIM528) transgenic plants. Knockout of D3 gene (d3 mutants) or miR528-overexpressing (OE-miR528) triggers a substantial reduction of gibberellin (GA) content, but a significant increase of abscisic acid (ABA) accumulation than in WT. The d3 and OE-miR528 transgenic plants were much more sensitive to GA, but less sensitive to ABA than WT. Moreover, the expression level of GA biosynthesis-related key genes, including OsCPS1, OsCPS2, OsKO2 and OsKAO was remarkably higher in OE-D3 plants, while the NECD2 expression, a key gene involved in ABA biosynthesis, was significantly higher in d3 mutants than in WT and OE-D3 plants. The results indicate that the miR528-D3 module negatively regulates plant height in rice by modulating the GA and ABA homeostasis, thereby further affecting the elongation of internodes, and resulting in lower plant height, which adds a new regulatory role to the D3-mediated plant height controlling in rice.

The miR528-AO Module Confers Enhanced Salt Tolerance in Rice by Modulating the Ascorbic Acid and Abscisic Acid Metabolism and ROS Scavenging.

The monocot lineage-specific miR528 was previously established as a multistress regulator. However, it remains largely unclear how miR528 participates in response to salinity stress in rice. Here, we show that miR528 positively regulates rice salt tolerance by down-regulating a gene encoding l-ascorbate oxidase (AO), thereby bolstering up the AO-mediated abscisic acid (ABA) synthesis and ROS scavenging. Overexpression of miR528 caused a substantial increase in ascorbic acid (AsA) and ABA contents but a significant reduction in ROS accumulation, resulting in the enhanced salt tolerance of rice plants. Conversely, knockdown of miR528 or overexpression of AO stimulated the expression of the AO gene, hence lowering the level of AsA, a critical antioxidant that promotes the ABA content but reduces the ROS level, and then compromising rice tolerance to salinity. Together, the findings reveal a novel mechanism of the miR528-AO module-mediated salt tolerance by modulating the processes of AsA and ABA metabolism as well as ROS detoxification, which adds a new regulatory role to the miR528-AO stress defense pathway in rice.

Functional divergence of the NIP III subgroup proteins involved altered selective constraints and positive selection.

BACKGROUND: Nod26-like intrinsic proteins (NIPs) that belong to the aquaporin superfamily are unique to plants. According to homology modeling and phylogenetic analysis, the NIP subfamily can be further divided into three subgroups with distinct biological functions (NIP I, NIP II, and NIP III). In some grasses, the NIP III subgroup proteins (NIP2s) were demonstrated to be permeable to solutes with larger diameter, such as silicic acid and arsenous acids. However, to date there is no data-mining or direct experimental evidences for the permeability of such larger solutes for dicot NIP2s, although they exhibit similar three-dimensional structures as those in grasses. It is therefore intriguing to investigate the molecular mechanisms that drive the evolution of plant NIP2s. RESULTS: The NIP III subgroup is more ancient with a divergence time that predates the monocot-dicot split. The proliferation of NIP2 genes in modern grass species is primarily attributed to whole genome and segmental chromosomal duplication events. The structure of NIP2 genes is relatively conserved, possessing five exons and four introns. All NIP2s possess an ar/R filter consisting of G, S, G, and R, except for the cucumber CsNIP2;2, where a small G in the H2 is substituted with the bulkier C residue. Our maximum likelihood analysis revealed that NIP2s, especially the loop A (LA) region, have undergone strong selective pressure for adaptive evolution. The analysis at the amino acid level provided strong statistical evidences for the functional divergence between monocot and dicot NIP III subgroup proteins. In addition, several SDPs (Specificity Determining Positions) responsible for functional specificity were predicted. CONCLUSIONS: The present study provides the first evidences of functional divergence between dicot and monocot NIP2s, and suggests that positive selection, as well as a radical shift of evolutionary rate at some critical amino acid sites is the primary driver. These findings will expand our understanding to evolutionary mechanisms driving the functional diversification of monocot and dicot NIP III subgroup proteins.

Evolutionary divergence of function and expression of laccase genes in plants.

Laccases (LACs) are versatile enzymes that catalyze oxidation of a wide range of substrates, thereby functioning in regulation of plant developmental processes and stress responses. However, with a few exceptions, the function of most LACs remains unclear in plants. In this study, we newly identified 4, 12, 22, 26, 27, 28 and 49 LAC genes for Physcomitrella patens, Amborella trichopoda, Zeamays, Ricinus communis, Vitis vinifera, Triticum aestivum and Glycine max, on the basis of exhaustive homologous sequence searches. In these plants, LACs differ greatly in sequence length and physical properties, such as molecular weight and theoretical isoelectric point (pI), but majority of them contain a signal peptide at their N-terminus. The originality of LACs could be traced back to as early as the emergence of moss. Plant LACs are clearly divided into seven distinct classes, where six ancient LACs should be present prior to the divergence of gymnosperms and angiosperms. Functional divergence analysis reveal that functional differentiation should occur among different groups of LACs because of altered selective constraints working on some critical amino acid sites (CAASs) within conserved laccase domains during evolution. Soybean and maize LACs have significantly different exon frequency (6.08 vs 4.82), and they are unevenly distributed and tend to form gene clusters on some chromosomes. Further analysis shows that the expansion of LAC gene family would be due toextensive tandem and chromosomal segmental duplications in the two plant species. Interestingly, *81.6% and 36.4% of soybean and maize LACs are potential targets of miRNAs, such as miR397a/b, miR408d, or miR528a/b etc. Both soybean and maize LACs are tissue specifically and developmental-specifically expressed, and are in response to different external abiotic and biotic stressors. These results suggest a diversity of functions of plant LAC genes, which will broaden our understanding and lay solid foundation for further investigating their biological functions in plants.

Divergence in function and expression of the NOD26-like intrinsic proteins in plants.

BACKGROUND: NOD26-like intrinsic proteins (NIPs) that belong to the aquaporin superfamily are plant-specific and exhibit a similar three-dimensional structure. Experimental evidences however revealed that functional divergence should have extensively occurred among NIP genes. It is therefore intriguing to further investigate the evolutionary mechanisms being responsible for the functional diversification of the NIP genes. To better understand this process, a comprehensive analysis including the phylogenetic, positive selection, functional divergence, and transcriptional analysis was carried out. RESULTS: The origination of NIPs could be dated back to the primitive land plants, and their diversification would be no younger than the emergence time of the moss P. patens. The rapid proliferation of NIPs in plants may be primarily attributed to the segmental chromosome duplication produced by polyploidy and tandem duplications. The maximum likelihood analysis revealed that NIPs should have experienced strong selective pressure for adaptive evolution after gene duplication and/or speciation, prompting the formation of distinct NIP groups. Functional divergence analysis at the amino acid level has provided strong statistical evidence for shifted evolutionary rate and/or radical change of the physiochemical properties of amino acids after gene duplication, and DIVERGE2 has identified the critical amino acid sites that are thought to be responsible for the divergence for further investigation. The expression of plant NIPs displays a distinct tissue-, cell-type-, and developmental specific pattern, and their responses to various stress treatments are quite different also. The differences in organization of cis-acting regulatory elements in the promoter regions may partially explain their distinction in expression. CONCLUSION: A number of analyses both at the DNA and amino acid sequence levels have provided strong evidences that plant NIPs have suffered a high divergence in function and expression during evolution, which is primarily attributed to the strong positive selection or a rapid change of evolutionary rate and/or physiochemical properties of some critical amino acid sites.

Dicer-like (DCL) proteins in plants.

Dicer and Dicer-like (DCL) proteins are key components in small RNA biogenesis. DCLs form a small protein family in plants whose diversification time dates to the emergence of mosses (Physcomitrella patens). DCLs are ubiquitously but not evenly expressed in tissues, at different developmental stages, and in response to environmental stresses. In Arabidopsis, AtDCL1, AtDCL2, and AtDCL4 exhibit similar expression pattern during the leaf or stem development, which is distinguished from AtDCL3. However, distinct expression profiles for all DCLs are found during the development of reproductive organs flower and seed. The grape VvDCL1 and VvDCL3 may act sequentially to face the fungi challenge. Overall, the responses of DCLs to drought, cold, and salt are quite different, indicating that plants might have specialized regulatory mechanism in response to different abiotic stresses. Further analysis of the promoter regions reveals a few of cis-elements that are hormone- and stress-responsive and developmental-related. However, gain and loss of cis-elements are frequent during evolution, and not only paralogous but also orthologous DCLs have dissimilar cis-element organization. In addition to cis-elements, AtDCL1 is probably regulated by both ath-miR162 and ath-miR414. Posterior analysis has identified some critical amino acid sites that are responsible for functional divergence between DCL family members. These findings provide new insights into understanding DCL protein functions.

Gene conversion and positive selection driving the evolution of the Caenorhabditis ssp. ZIM/HIM-8 protein family.

In C. elegans, four C2H2 zinc-finger proteins (ZIM-1, ZIM-2, ZIM-3, and HIM-8), which are arranged in tandem, mediate chromosome-specific pairing and synapsis during meiosis. The zim/him-8 genes from three Caenorhabditis species were contrasted in an effort to investigate the mechanisms driving their evolution. Here it is shown that the preservation of higher degree of sequence similarity in the N-terminal portion, particularly in several regions within the second exon between paralogous zim genes (especially between zim-1 and zim-3), is due to independent interparalogue gene conversions. However, the evolutionary force is not uniformly strong across species. The present data reveal that more frequent gene conversion events have occurred in C. elegans, whereas only gene conversions between zim-1 and zim-3 are detected in C. remanei. Although gene conversions are predicted to be present among zim-1, zim-2, and zim-3 in C. briggsae, the conversion tracts between zim-1/zim-2 and zim-2/zim-3 are very short. Moreover, positive selection analysis was performed on the basis of the significantly discordant phylogenies reconstructed using the N- and C-terminal sequences, respectively. Several codon sites located in the regions that are supposed not to have experienced gene conversions are predicted to be under the influence of positive selection. In comparison, stronger positive selection has acted on the C-terminal region relative to the N-terminal region. Thus, the zim/him-8 genes that evolve concertedly have also been shown to undergo adaptive diversifying selection.

The genome-wide dynamics of purging during selfing in maize.

Self-fertilization (also known as selfing) is an important reproductive strategy in plants and a widely applied tool for plant genetics and plant breeding. Selfing can lead to inbreeding depression by uncovering recessive deleterious variants, unless these variants are purged by selection. Here we investigated the dynamics of purging in a set of eleven maize lines that were selfed for six generations. We show that heterozygous, putatively deleterious single nucleotide polymorphisms are preferentially lost from the genome during selfing. Deleterious single nucleotide polymorphisms were lost more rapidly in regions of high recombination, presumably because recombination increases the efficacy of selection by uncoupling linked variants. Overall, heterozygosity decreased more slowly than expected, by an estimated 35% to 40% per generation instead of the expected 50%, perhaps reflecting pervasive associative overdominance. Finally, three lines exhibited marked decreases in genome size due to the purging of transposable elements. Genome loss was more likely to occur for lineages that began with larger genomes with more transposable elements and chromosomal knobs. These three lines purged an average of 398 Mb from their genomes, an amount equivalent to three Arabidopsis thaliana genomes per lineage, in only a few generations.

Molecular evolution of the MLO gene family in Oryza sativa and their functional divergence.

The present study identified 12 MLO genes in rice that were located on chromosomes 1, 2, 3, 4, 5, 6, 10, and 11 respectively without any obvious clustering. On a genome scale we showed that the expansion of rice MLO gene family was primarily attributed to segmental duplication produced by polyploidy, rather than through tandem amplification. Gene conversion events should also play important roles in the evolution of MLO genes. The results of relative rate ratio test and maximum likelihood analysis suggested that positive selection should have occurred after gene duplication and/or speciation, prompting the formation of distinct MLO subfamilies. Functional divergence analysis provided statistical evidence for shifted evolutionary rate after gene duplication. Compared to extracellular loop 3 and Ca(2+)-binding domain, much stronger functional constraints should impose on intracellular loop 2, although all of the three regions might be under purifying selection. The sliding window analysis of d(N)/d(S) ratio values identified one sequence region where strong functional constraints must impose on, and consequently should be crucial for functionality of MLO genes.

Evolution and functional divergence of monocarboxylate transporter genes in vertebrates.

Monocarboxylate transporters (MCTs) form a gene family with an ancient past. The identification of MCTs (MCHs) from bacteria, protozoa, fungi, invertebrates, as well as vertebrates, but not from plants and virus, allowed illuminating the phylogenetic and evolutionary history of this gene family. The significant expansion of vertebrate MCT genes should have primarily occurred after the divergence of vertebrates and invertebrates, but before the divergence time between ray-finned fish and mammals. The divergence of insect MCTs should have at least occurred in the common ancestor of fruit fly, beetle, and honeybee. Fungi monocarboxylate transporter homologues (MCHs) might evolve independently from an ancient ancestor. The results of functional divergence analysis provided statistical evidences for shifted evolutionary rate and/or changes of amino acid property after gene duplication. The sliding window analysis of the d(N)/d(S) ratio values showed that strong functional constraints must impose on the N- and C-terminal domains of vertebrate MCTs. These corresponding regions may play crucial roles for functionality of MCT proteins.

Identification of rice TUBBY-like genes and their evolution.

The identification of TUBBY-like genes in organisms ranging from single-celled to multicellular eukaryotes has allowed the phylogenetic history of this gene family to be traced back to the early evolutionary stages of eukaryote development. Rice TUBBY-like genes were located on chromosomes 1, 2, 3, 4, 5, 7, 8, 11 and 12 without any obvious clustering. On a genomic scale, it was revealed that the rice TUBBY-like gene family probably evolved mainly through segmental duplication produced by polyploidy. The altered selective constraints (or site-specific rate changes), related to functional divergence during protein evolution between plant and animal TUBBY-like genes, were statistically significant. Based on posterior probability analysis, five amino acid sites (103, 312, 315, 317 and 319) are thought to be responsible for functional divergence.