Loss-of-function of SAWTOOTH 1 affects leaf dorsiventrality genes to promote leafy heads in lettuce.

The mechanisms underlying leafy heads in vegetables are poorly understood. Here, we cloned a quantitative trait locus (QTL) controlling leafy heads in lettuce (Lactuca sativa). The QTL encodes a transcription factor, SAWTOOTH 1 (LsSAW1), which has a BEL1-like homeodomain and is a homolog of Arabidopsis thaliana. A 1-bp deletion in Lssaw1 contributes to the development of leafy heads. Laser-capture microdissection and RNA-sequencing showed that LsSAW1 regulates leaf dorsiventrality and loss-of-function of Lssaw1 downregulates the expression of many adaxial genes but upregulates abaxial genes. LsSAW1 binds to the promoter region of the adaxial gene ASYMMETRIC LEAVES 1 (LsAS1) to upregulate its expression. Overexpression of LsAS1 compromised the effects of Lssaw1 on heading. LsSAW1 also binds to the promoter region of the abaxial gene YABBY 1 (LsYAB1), but downregulates its expression. Overexpression of LsYAB1 led to bending leaves in LsSAW1 genotypes. LsSAW1 directly interacts with KNOTTED 1 (LsKN1), which is necessary for leafy heads in lettuce. RNA-seq data showed that LsSAW1 and LsKN1 exert antagonistic effects on the expression of thousands of genes. LsSAW1 compromises the ability of LsKN1 to repress LsAS1. Our results suggest that downregulation or loss-of-function of adaxial genes and upregulation of abaxial genes allow for the development of leafy heads.

A WOX homolog disrupted by a transposon led to the loss of spines and contributed to the domestication of lettuce.

The loss of spines is one of the most important domestication traits for lettuce (Lactuca sativa). However, the genetics and regulation of spine development in lettuce remain unclear. We examined the genetics of spines in lettuce using a segregating population derived from a cross between cultivated and wild lettuce (Lactuca serriola). A gene encoding WUSCHEL-related homeobox transcription factor, named as WOX-SPINE1 (WS1), was identified as the candidate gene controlling the spine development in lettuce, and its function on spines was verified. A CACTA transposon was found to be inserted into the first exon of the ws1 allele, knocking out its function and leading to the lack of spines in cultivated lettuce. All lettuce cultivars investigated have the nonfunctional ws1 gene, and a selection sweep was found at the WS1 locus, suggesting its important role in lettuce domestication. The expression levels of WS1 were associated with the density of spines among different accessions of wild lettuce. At least two independent loss-of-function mutations in the ws1 gene caused the loss of spines in wild lettuce. These findings provide new insights into the development of spines and facilitate the exploitation of wild genetic resources in future lettuce breeding programs.

Alternative splicing triggered by the insertion of a CACTA transposon attenuates LsGLK and leads to the development of pale-green leaves in lettuce.

Lettuce (Lactuca sativa) is one of the most important vegetable crops in the world. As a leafy vegetable, the polymorphism of lettuce leaves from dark to pale green is an important trait. However, the genetic and molecular mechanisms underlying such variations remain poorly understood. In this study, one major locus controlling the polymorphism of dark- and pale-green leaves in lettuce was identified using genome-wide association studies (GWAS). This locus was then fine mapped to an interval of 5375 bp on chromosome 4 using a segregating population containing 2480 progeny. Only one gene, homologous to the GLK genes in Arabidopsis and other plants, is present in the candidate region. A complementation test confirmed that the candidate gene, LsGLK, contributes to the variation of dark- and pale-green leaves. Sequence analysis showed that a CACTA transposon of 7434 bp was inserted 10 bp downstream of the stop codon of LsGLK, followed by a duplication of a 1826-bp fragment covering exons 3-6 of the LsGLK gene. The transposon insertion did not change the expression level of the LsGLK gene. However, because of alternative splicing, only 6% of the transcripts produced from the transposon insertion were wild-type transcripts, which led to the production of pale-green leaves. An evolutionary analysis revealed that the insertion of the CACTA transposon occurred in cultivated lettuce and might have been selected in particular cultivars to satisfy the diverse demands of consumers. In this study, we demonstrated that a transposon insertion near a gene may affect its splicing and consequently generate phenotypic variations.

Seed color in lettuce is determined by the LsTT2, LsCHS, and Ls2OGD genes from the flavonoid biosynthesis pathway.

KEY MESSAGE: The mutated LsTT2 and Ls2OGD genes are responsible for white seeds and yellow seeds in lettuce, respectively. Three LsCHS genes are involved in the biosynthesis of flavonoid in seed coats. Lettuce seeds have several different colors, including black, yellow, and white. The genetic mechanisms underlying color variations of lettuce seeds remain unknown. We used genome-wide association studies (GWAS) and map-based cloning approaches to clone genes controlling the color of lettuce seeds. LsTT2, which encodes an R2R3-MYB transcription factor and is homologous to the TT2 gene in Arabidopsis, was shown to be the causal gene for the variation of black and white seeds in lettuce. A point mutation leads to the lack of stop codon in the LsTT2 transcript, resulting in white seeds. Knockout of the LsTT2 gene converted black seeds to white seeds. The locus controlling yellow seeds was mapped to Chromosome 2. Knockout of two 2-oxoglutarate-dependent dioxygenases (2OGD) genes from the candidate region converted black seeds to yellow seeds, suggesting that these two 2OGD proteins catalyze the conversion of yellow metabolites to black metabolites. We also showed that three LsCHS genes from the candidate region are associated with flavonoid biosynthesis in seeds. Knockout mutants of the three LsCHS genes decreased color intensity. This study provides new insights into the regulation of flavonoid biosynthesis in plants.

Upregulation of a KN1 homolog by transposon insertion promotes leafy head development in lettuce.

Leafy head is a unique type of plant architecture found in some vegetable crops, with leaves bending inward to form a compact head. The genetic and molecular mechanisms underlying leafy head in vegetables remain poorly understood. We genetically fine-mapped and cloned a major quantitative trait locus controlling heading in lettuce. The candidate gene (LsKN1) is a homolog of knotted 1 (KN1) from Zea mays Complementation and CRISPR/Cas9 knockout experiments confirmed the role of LsKN1 in heading. In heading lettuce, there is a CACTA-like transposon inserted into the first exon of LsKN1 (LsKN1▽). The transposon sequences act as a promoter rather than an enhancer and drive high expression of LsKN1▽. The enhanced expression of LsKN1▽ is necessary but not sufficient for heading in lettuce. Data from ChIP-sequencing, electrophoretic mobility shift assays, and dual luciferase assays indicate that the LsKN1▽ protein binds the promoter of LsAS1 and down-regulates its expression to alter leaf dorsoventrality. This study provides insight into plant leaf development and will be useful for studies on heading in other vegetable crops.

The upregulated LsKN1 gene transforms pinnately to palmately lobed leaves through auxin, gibberellin, and leaf dorsiventrality pathways in lettuce.

Leaf shape represents a vital agronomic trait for leafy vegetables such as lettuce. Some lettuce cultivars produce lobed leaves, varying from pinnately to palmately lobed, but the genetic mechanisms remain unclear. In this study, we cloned one major quantitative trait locus (QTL) controlling palmately lobed leaves. The candidate gene, LsKN1, encodes a homeobox transcription factor, and has been shown previously to be critical for the development of leafy heads in lettuce. The LsKN1 allele that is upregulated by the insertion of a transposon promotes the development of palmately lobed leaves. We demonstrated that LsKN1 upregulated LsCUC2 and LsCUC3 through different mechanisms, and their upregulation was critical for the development of palmately lobed leaves. LsKN1 binds the promoter of LsPID to promote auxin biosynthesis, which positively contributes to the development of palmately lobed leaves. In contrast, LsKN1 suppresses GA biosynthesis to promote palmately lobed leaves. LsKN1 also binds to the promoter of LsAS1, a dorsiventrality gene, to downregulate its expression. Overexpression of the LsAS1 gene compromised the effects of the LsKN1 gene changing palmately to pinnately lobed leaves. Our study illustrated that the upregulated LsKN1 gene led to palmately lobed leaves in lettuce by integrating several downstream pathways, including auxin, gibberellin, and leaf dorsiventrality pathways.

LsNRL4 enhances photosynthesis and decreases leaf angles in lettuce.

Lettuce (Lactuca sativa) is one of the most important vegetables worldwide and an ideal plant for producing protein drugs. Both well-functioning chloroplasts that perform robust photosynthesis and small leaf angles that enable dense planting are essential for high yields. In this study, we used an F2 population derived from a cross between a lettuce cultivar with pale-green leaves and large leaf angles to a cultivar with dark-green leaves and small leaf angles to clone LsNRL4, which encodes an NPH3/RPT2-Like (NRL) protein. Unlike other NRL proteins in lettuce, the LsNRL4 lacks the BTB domain. Knockout mutants engineered using CRISPR/Cas9 and transgenic lines overexpressing LsNRL4 verified that LsNRL4 contributes to chloroplast development, photosynthesis and leaf angle. The LsNRL4 gene was not present in the parent with pale-green leaves and enlarged leaf angles. Loss of LsNRL4 results in the enlargement of chloroplasts, decreases in the amount of cellular space allocated to chloroplasts and defects in secondary cell wall biosynthesis in lamina joints. Overexpressing LsNRL4 significantly improved photosynthesis and decreased leaf angles. Indeed, the plant architecture of the overexpressing lines is ideal for dense planting. In summary, we identified a novel NRL gene that enhances photosynthesis and influences plant architecture. Our study provides new approaches for the breeding of lettuce that can be grown in dense planting in the open field or in modern plant factories. LsNRL4 homologues may also be used in other crops to increase photosynthesis and improve plant architecture.

Two TOBAMOVIRUS MULTIPLICATION 2A homologs in tobacco control asymptomatic response to tobacco mosaic virus.

The most common response of a host to pathogens is arguably the asymptomatic response. However, the genetic and molecular mechanisms responsible for asymptomatic responses to pathogens are poorly understood. Here we report on the genetic cloning of two genes controlling the asymptomatic response to tobacco mosaic virus (TMV) in cultivated tobacco (Nicotiana tabacum). These two genes are homologous to tobamovirus multiplication 2A (TOM2A) from Arabidopsis, which was shown to be critical for the accumulation of TMV. Expression analysis indicates that the TOM2A genes might play fundamental roles in plant development or in responses to stresses. Consistent with this hypothesis, a null allele of the TOM2A ortholog in tomato (Solanum lycopersicum) led to the development of bent branches and a high tolerance to both TMV and tomato mosaic virus (ToMV). However, the TOM2A ortholog in Nicotiana glauca did not account for the asymptomatic response to TMV in N. glauca. We showed that TOM2A family is plant-specific and originated from Chlorophyte, and the biological functions of TOM2A orthologs to promote TMV accumulation are highly conserved in the plant kingdom-in both TMV host and nonhost species. In addition, we showed that the interaction between tobacco TOM1 and TOM2A orthologs in plant species is conserved, suggesting a conserved nature of TOM1-TOM2A module in promoting TMV multiplication in plants. The tradeoff between host development, the resistance of hosts to pathogens, and their influence on gene evolution are discussed. Our results shed light on mechanisms that contribute to asymptomatic responses to viruses in plants and provide approaches for developing TMV/ToMV-resistant crops.

Up-regulation of LsKN1 promotes cytokinin and suppresses gibberellin biosynthesis to generate wavy leaves in lettuce.

Lettuce (Lactuca sativa) is one of the most popular vegetables worldwide, and diverse leaf shapes, including wavy leaves, are important commercial traits. In this study, we examined the genetics of wavy leaves using an F2 segregating population, and identified a major QTL controlling wavy leaves. The candidate region contained LsKN1, which has previously been shown to be indispensable for leafy heads in lettuce. Complementation tests and knockout experiments verified the function of LsKN1 in producing wavy leaves. The LsKN1∇ allele, which has the insertion of a transposon and has previously been shown to control leafy heads, promoted wavy leaves in our population. Transposition of the CACTA transposon from LsKN1 compromised its function for wavy leaves. High expression of LsKN1 up-regulated several key genes associated with cytokinin (CK) to increase the content in the leaves, whereas it down-regulated the expression of genes in the gibberellin (GA) biosynthesis pathway to decrease the content. Application of CK to leaves enhanced the wavy phenotype, while application of GA dramatically flattened the leaves. We conclude that the changes in CK and GA contents that result from high expression of LsKN1 switch determinate cells to indeterminate, and consequently leads to the development of wavy leaves.

Frequent gain- and loss-of-function mutations of the BjMYB113 gene accounted for leaf color variation in Brassica juncea.

BACKGROUND: Mustard (Brassica juncea) is an important economic vegetable, and some cultivars have purple leaves and accumulate more anthocyanins than the green. The genetic and evolution of purple trait in mustard has not been well studied. RESULT: In this study, free-hand sections and metabolomics showed that the purple leaves of mustard accumulated more anthocyanins than green ones. The gene controlling purple leaves in mustard, Mustard Purple Leaves (MPL), was genetically mapped and a MYB113-like homolog was identified as the candidate gene. We identified three alleles of the MYB113-like gene, BjMYB113a from a purple cultivar, BjMYB113b and BjMYB113c from green cultivars. A total of 45 single nucleotide polymorphisms (SNPs) and 8 InDels were found between the promoter sequences of the purple allele BjMYB113a and the green allele BjMYB113b. On the other hand, the only sequence variation between the purple allele BjMYB113a and the green allele BjMYB113c is an insertion of 1,033-bp fragment in the 3'region of BjMYB113c. Transgenic assay and promoter activity studies showed that the polymorphism in the promoter region was responsible for the up-regulation of the purple allele BjMYB113a and high accumulation of anthocyanin in the purple cultivar. The up-regulation of BjMYB113a increased the expression of genes in the anthocyanin biosynthesis pathway including BjCHS, BjF3H, BjF3'H, BjDFR, BjANS and BjUGFT, and consequently led to high accumulation of anthocyanin. However, the up-regulation of BjMYB113 was compromised by the insertion of 1,033-bp in 3'region of the allele BjMYB113c. CONCLUSIONS: Our results contribute to a better understanding of the genetics and evolution of the BjMYB113 gene controlling purple leaves and provide useful information for further breeding programs of mustard.

Characterization of four polymorphic genes controlling red leaf colour in lettuce that have undergone disruptive selection since domestication.

Anthocyanins protect plants from biotic and abiotic stressors and provide great health benefits to consumers. In this study, we cloned four genes (Red Lettuce Leaves 1 to 4: RLL1 to RLL4) that contribute to colour variations in lettuce. The RLL1 gene encodes a bHLH transcription factor, and a 5-bp deletion in some cultivars abolishes its function to activate the anthocyanin biosynthesis pathway. The RLL2 gene encodes an R2R3-MYB transcription factor, which was derived from a duplication followed by mutations in its promoter region. The RLL3 gene encodes an R2-MYB transcription factor, which down-regulates anthocyanin biosynthesis through competing with RLL2 for interaction with RLL1; a mis-sense mutation compromises the capacity of RLL3 to bind RLL1. The RLL4 gene encodes a WD-40 transcription factor, homologous to the RUP genes suppressing the UV-B signal transduction pathway in Arabidopsis; a mis-sense mutation in rll4 attenuates its suppressing function, leading to a high concentration of anthocyanins. Sequence analysis of the RLL1-RLL4 genes from wild and cultivated lettuce showed that their function-changing mutations occurred after domestication. The mutations in rll1 disrupt anthocyanin biosynthesis, while the mutations in RLL2, rll3 and rll4 activate anthocyanin biosynthesis, showing disruptive selection for leaf colour during domestication of lettuce. The characterization of multiple polymorphic genes in this study provides the necessary molecular resources for the rational breeding of lettuce cultivars with distinct levels of red pigments and green cultivars with high levels of health-promoting flavonoids.

Independent activation of the BoMYB2 gene leading to purple traits in Brassica oleracea.

KEY MESSAGE: Transposon insertion and point mutation independently activated the BoMYB2 gene in three purple cultivars of Brassica oleracea including kale, kohlrabi, and cabbage. Several varieties of B. oleracea have both green and purple cultivars. In this study, the causal genes for the purple traits in kale, kohlrabi and cabbage were cloned using map-based cloning approach. The purple traits in all three varieties were mapped to the same locus as the BoMYB2 gene in cauliflower. Surprisingly, the insertion of Harbinger transposon of BoMYB2 in cauliflower was not found in purple kale, kohlrabi and cabbage. Sequencing of the BoMYB2 gene in purple kale and purple kohlrabi discovered a 7606 bp CACTA-like transposon in its promoter region. Transient assay and promoter activity study showed that the insertion upregulated the expression of the BoMYB2 gene. On the other hand, the activation of the BoMYB2 gene in purple cabbage was caused by point mutation and/or 1-bp insertion in its promoter region. Sequence analysis of the BoMYB2 gene in different varieties suggested that the activating events most likely occurred independently after the divergence of cabbage, cauliflower, and kale/kohlrabi. Our results not only contribute to a better understanding of anthocyanin inheritance in B. oleracea, but also provide useful information for future hybrid breeding of purple cultivars through combination of different functional alleles of the BoMYB2 gene.

P-MITE: a database for plant miniature inverted-repeat transposable elements.

Miniature inverted-repeat transposable elements (MITEs) are prevalent in eukaryotic species including plants. MITE families vary dramatically and usually cannot be identified based on homology. In this study, we de novo identified MITEs from 41 plant species, using computer programs MITE Digger, MITE-Hunter and/or Repetitive Sequence with Precise Boundaries (RSPB). MITEs were found in all, but one (Cyanidioschyzon merolae), species. Combined with the MITEs identified previously from the rice genome, >2.3 million sequences from 3527 MITE families were obtained from 41 plant species. In general, higher plants contain more MITEs than lower plants, with a few exceptions such as papaya, with only 538 elements. The largest number of MITEs is found in apple, with 237 302 MITE sequences. The number of MITE sequences in a genome is significantly correlated with genome size. A series of databases (plant MITE databases, P-MITE), available online at http://pmite.hzau.edu.cn/django/mite/, was constructed to host all MITE sequences from the 41 plant genomes. The databases are available for sequence similarity searches (BLASTN), and MITE sequences can be downloaded by family or by genome. The databases can be used to study the origin and amplification of MITEs, MITE-derived small RNAs and roles of MITEs on gene and genome evolution.

Molecular Modeling Studies of 11ß-Hydroxysteroid Dehydrogenase Type 1 Inhibitors through Receptor-Based 3D-QSAR and Molecular Dynamics Simulations.

11ß-Hydroxysteroid dehydrogenase type 1 (11ß-HSD1) is a potential target for the treatment of numerous human disorders, such as diabetes, obesity, and metabolic syndrome. In this work, molecular modeling studies combining molecular docking, 3D-QSAR, MESP, MD simulations and free energy calculations were performed on pyridine amides and 1,2,4-triazolopyridines as 11ß-HSD1 inhibitors to explore structure-activity relationships and structural requirement for the inhibitory activity. 3D-QSAR models, including CoMFA and CoMSIA, were developed from the conformations obtained by docking strategy. The derived pharmacophoric features were further supported by MESP and Mulliken charge analyses using density functional theory. In addition, MD simulations and free energy calculations were employed to determine the detailed binding process and to compare the binding modes of inhibitors with different bioactivities. The binding free energies calculated by MM/PBSA showed a good correlation with the experimental biological activities. Free energy analyses and per-residue energy decomposition indicated the van der Waals interaction would be the major driving force for the interactions between an inhibitor and 11ß-HSD1. These unified results may provide that hydrogen bond interactions with Ser170 and Tyr183 are favorable for enhancing activity. Thr124, Ser170, Tyr177, Tyr183, Val227, and Val231 are the key amino acid residues in the binding pocket. The obtained results are expected to be valuable for the rational design of novel potent 11ß-HSD1 inhibitors.

The I2 resistance gene homologues in Solanum have complex evolutionary patterns and are targeted by miRNAs.

BACKGROUND: Several resistance traits, including the I2 resistance against tomato fusarium wilt, were mapped to the long arm of chromosome 11 of Solanum. However, the structure and evolution of this locus remain poorly understood. RESULTS: Comparative analysis showed that the structure and evolutionary patterns of the I2 locus vary considerably between potato and tomato. The I2 homologues from different Solanaceae species usually do not have orthologous relationship, due to duplication, deletion and frequent sequence exchanges. At least 154 sequence exchanges were detected among 76 tomato I2 homologues, but sequence exchanges between I2 homologues in potato is less frequent. Previous study showed that I2 homologues in potato were targeted by miR482. However, our data showed that I2 homologues in tomato were targeted by miR6024 rather than miR482. Furthermore, miR6024 triggers phasiRNAs from I2 homologues in tomato. Sequence analysis showed that miR6024 was originated after the divergence of Solanaceae. We hypothesized that miR6024 and miR482 might have facilitated the expansion of the I2 family in Solanaceae species, since they can minimize their potential toxic effects by down-regulating their expression. CONCLUSIONS: The I2 locus represents a most divergent resistance gene cluster in Solanum. Its high divergence was partly due to frequent sequence exchanges between homologues. We propose that the successful expansion of I2 homologues in Solanum was at least partially attributed to miRNA mediated regulation.

Frequent loss of lineages and deficient duplications accounted for low copy number of disease resistance genes in Cucurbitaceae.

BACKGROUND: The sequenced genomes of cucumber, melon and watermelon have relatively few R-genes, with 70, 75 and 55 copies only, respectively. The mechanism for low copy number of R-genes in Cucurbitaceae genomes remains unknown. RESULTS: Manual annotation of R-genes in the sequenced genomes of Cucurbitaceae species showed that approximately half of them are pseudogenes. Comparative analysis of R-genes showed frequent loss of R-gene loci in different Cucurbitaceae species. Phylogenetic analysis, data mining and PCR cloning using degenerate primers indicated that Cucurbitaceae has limited number of R-gene lineages (subfamilies). Comparison between R-genes from Cucurbitaceae and those from poplar and soybean suggested frequent loss of R-gene lineages in Cucurbitaceae. Furthermore, the average number of R-genes per lineage in Cucurbitaceae species is approximately 1/3 that in soybean or poplar. Therefore, both loss of lineages and deficient duplications in extant lineages accounted for the low copy number of R-genes in Cucurbitaceae. No extensive chimeras of R-genes were found in any of the sequenced Cucurbitaceae genomes. Nevertheless, one lineage of R-genes from Trichosanthes kirilowii, a wild Cucurbitaceae species, exhibits chimeric structures caused by gene conversions, and may contain a large number of distinct R-genes in natural populations. CONCLUSIONS: Cucurbitaceae species have limited number of R-gene lineages and each genome harbors relatively few R-genes. The scarcity of R-genes in Cucurbitaceae species was due to frequent loss of R-gene lineages and infrequent duplications in extant lineages. The evolutionary mechanisms for large variation of copy number of R-genes in different plant species were discussed.

Miniature inverted-repeat transposable elements (MITEs) have been accumulated through amplification bursts and play important roles in gene expression and species diversity in Oryza sativa.

Miniature inverted-repeat transposable elements (MITEs) are predicted to play important roles on genome evolution. We developed a BLASTN-based approach for de novo identification of MITEs and systematically analyzed MITEs in rice genome. The genome of rice cultivar Nipponbare (Oryza sativa ssp. japonica) harbors 178,533 MITE-related sequences classified into 338 families. Pairwise nucleotide diversity and phylogenetic tree analysis indicated that individual MITE families were resulted from one or multiple rounds of amplification bursts. The timing of amplification burst varied considerably between different MITE families or subfamilies. MITEs are associated with 23,623 (58.2%) genes in rice genome. At least 7,887 MITEs are transcribed and more than 3,463 were transcribed with rice genes. The MITE sequences transcribed with rice coding genes form 1,130 pairs of potential natural sense/antisense transcripts. MITEs generate 23.5% (183,837 of 781,885) of all small RNAs identified from rice. Some MITE families generated small RNAs mainly from the terminals, while other families generated small RNAs predominantly from the central region. More than half (51.8%) of the MITE-derived small RNAs were generated exclusively by MITEs located away from genes. Genome-wide analysis showed that genes associated with MITEs have significantly lower expression than genes away from MITEs. Approximately 14.8% of loci with full-length MITEs have presence/absence polymorphism between rice cultivars 93-11 (O. sativa ssp. indica) and Nipponbare. Considering that different sets of genes may be regulated by MITE-derived small RNAs in different genotypes, MITEs provide considerable diversity for O. sativa.

Dynamic nucleotide-binding site and leucine-rich repeat-encoding genes in the grass family.

The proper use of resistance genes (R genes) requires a comprehensive understanding of their genomics and evolution. We analyzed genes encoding nucleotide-binding sites and leucine-rich repeats in the genomes of rice (Oryza sativa), maize (Zea mays), sorghum (Sorghum bicolor), and Brachypodium distachyon. Frequent deletions and translocations of R genes generated prevalent presence/absence polymorphism between different accessions/species. The deletions were caused by unequal crossover, homologous repair, nonhomologous repair, or other unknown mechanisms. R gene loci identified from different genomes were mapped onto the chromosomes of rice cv Nipponbare using comparative genomics, resulting in an integrated map of 495 R loci. Sequence analysis of R genes from the partially sequenced genomes of an African rice cultivar and 10 wild accessions suggested that there are many additional R gene lineages in the AA genome of Oryza. The R genes with chimeric structures (termed type I R genes) are diverse in different rice accessions but only account for 5.8% of all R genes in the Nipponbare genome. In contrast, the vast majority of R genes in the rice genome are type II R genes, which are highly conserved in different accessions. Surprisingly, pseudogene-causing mutations in some type II lineages are often conserved, indicating that their conservations were not due to their functions. Functional R genes cloned from rice so far have more type II R genes than type I R genes, but type I R genes are predicted to contribute considerable diversity in wild species. Type I R genes tend to reduce the microsynteny of their flanking regions significantly more than type II R genes, and their flanking regions have slightly but significantly lower G/C content than those of type II R genes.

The agonist binding mechanism of human CB2 receptor studied by molecular dynamics simulation, free energy calculation and 3D-QSAR studies.

CB2-selective agonists have drawn attention in drug discovery, since CB2 becomes a promising target for the treatment of neuropathic pain without psychoactive or other CNS-related side effects. However, the lack of experimental data of the 3D structures of human cannabinoid receptors hampers the understanding of the binding modes between ligands and CB2 by traditional methods. In the present work, combinational molecular modeling studies including flexible docking, MD simulations and free energy calculations were performed to investigate the interaction modes and mechanism of CB2-unselective agonist CP55940 and CB2-selective agonist GW842166X, separately binding with the homology model of CB2 in a DPPC/TIP3P simulated membrane environment. The binding free energies calculated by MM-PBSA method give an explanation for the activity differences of the studied ligands. Binding free energies decomposition by MM-GBSA method shows that the van der Waals interaction is the dominant driving force during the binding process. Our MD simulations demonstrate that Phe197 could be a critical residue for the binding of CB2-selective agonists. Furthermore, by using the MD simulated binding conformer as a template, the 3D-QSAR studies were performed with the comparative molecular field analysis (CoMFA) approach on a set of GW842166X analogues. A combinational exploration of both CoMFA steric and potential contour maps for CB2 affinities and the MD studied interaction modes sheds light on the structural requirements for CB2 agonists and serves as a basis for the design of novel CB2 agonists.