Duplicate and diverge: the evolution of plant genome microstructure

Genome-wide characterization of Alfin-like (AL) genes in apple and functional identification of MdAL4 in response to drought stress

Eleven Alfin-like (AL) genes were obtained from apple and MdAL4 was selected for improving drought stress tolerance of transgenic apple callus and Arabidopsis. Drought is an important environmental factor affecting plant growth all over the world. Alfin-like (AL) have well-documented functions in abiotic stress response, but their drought stress tolerance in apple (Malus domestica) are poorly understood. According to the transcriptome data, 11 MdAL genes containing conserved Alfin and PHD-finger domain were identified in apple and divided into three subgroups with a total of 35 members from different species. Subsequently, gene structures, conserved amino acid sequences, promoter cis-acting elements, and gene evolution events were analyzed. Based on differential expression of MdALs in response to abiotic stresses, MdAL4, which was highly expressed under drought, was further cloned and investigated. MdAL4 encoding nuclear-localized protein conferred enhanced drought tolerance in overexpressing transgenic calli of apple 'Orin'. Moreover, the ectopic expression of MdAL4 improved the drought tolerance of transgenic Arabidopsis, as judged from remarkably decreased malonaldehyde (MDA) content and electrolyte leakage in MdAL4 overexpressing plants relative to WT. Furthermore, MdAL4 possibly could bind to promoter regions of ROS-scavenging and stress-related genes to improve drought tolerance. Additionally, we found in silico evidence that three proteins containing the WD40 domain that interact with MdAL4. Based on these results, MdAL4 was identified as a positive regulator for improving drought stress of apple.

Genome-Wide Characterization of the Aquaporin Gene Family in Radish and Functional Analysis of RsPIP2-6 Involved in Salt Stress

Aquaporins (AQPs) constitute a highly diverse family of channel proteins that transport water and neutral solutes. AQPs play crucial roles in plant development and stress responses. However, the characterization and biological functions of RsAQPs in radish (Raphanus sativus L.) remain elusive. In this study, 61 non-redundant members of AQP-encoding genes were identified from the radish genome database and located on nine chromosomes. Radish AQPs (RsAQPs) were divided into four subfamilies, including 21 plasma membrane intrinsic proteins (PIPs), 19 tonoplast intrinsic proteins (TIPs), 16 NOD-like intrinsic proteins (NIPs), and 5 small basic intrinsic proteins (SIPs), through phylogenetic analysis. All RsAQPs contained highly conserved motifs (motifs 1 and 4) and transmembrane regions, indicating the potential transmembrane transport function of RsAQPs. Tissue- and stage-specific expression patterns of AQP gene analysis based on RNA-seq data revealed that the expression levels of PIPs were generally higher than TIPs, NIPs, and SIPs in radish. In addition, quantitative real-time polymerase chain reaction (qRT-PCR) revealed that seven selected RsPIPs, according to our previous transcriptome data (e.g., RsPIP1-3, 1-6, 2-1, 2-6, 2-10, 2-13, and 2-14), exhibited significant upregulation in roots of salt-tolerant radish genotype. In particular, the transcriptional levels of RsPIP2-6 dramatically increased after 6 h of 150 mM NaCl treatment during the taproot thickening stage. Additionally, overexpression of RsPIP2-6 could enhance salt tolerance by Agrobacterium rhizogenes-mediated transgenic radish hairy roots, which exhibited the mitigatory effects of plant growth reduction, leaf relative water content (RWC) reduction and alleviation of O2- in cells, as shown by nitro blue tetrazolium (NBT) staining, under salt stress. These findings are helpful for deeply dissecting the biological function of RsAQPs on the salt stress response, facilitating practical application and genetic improvement of abiotic stress resistance in radish.

Genome-Wide Identification and Analysis of the Polycomb Group Family in Medicago truncatula

Polycomb group (PcG) proteins, which are important epigenetic regulators, play essential roles in the regulatory networks involved in plant growth, development, and environmental stress responses. Currently, as far as we know, no comprehensive and systematic study has been carried out on the PcG family in Medicago truncatula. In the present study, we identified 64 PcG genes with distinct gene structures from the M. truncatula genome. All of the PcG genes were distributed unevenly over eight chromosomes, of which 26 genes underwent gene duplication. The prediction of protein interaction network indicated that 34 M. truncatula PcG proteins exhibited protein-protein interactions, and MtMSI1;4 and MtVRN2 had the largest number of protein-protein interactions. Based on phylogenetic analysis, we divided 375 PcG proteins from 27 species into three groups and nine subgroups. Group I and Group III were composed of five components from the PRC1 complex, and Group II was composed of four components from the PRC2 complex. Additionally, we found that seven PcG proteins in M. truncatula were closely related to the corresponding proteins of Cicer arietinum. Syntenic analysis revealed that PcG proteins had evolved more conservatively in dicots than in monocots. M. truncatula had the most collinearity relationships with Glycine max (36 genes), while collinearity with three monocots was rare (eight genes). The analysis of various types of expression data suggested that PcG genes were involved in the regulation and response process of M. truncatula in multiple developmental stages, in different tissues, and for various environmental stimuli. Meanwhile, many differentially expressed genes (DEGs) were identified in the RNA-seq data, which had potential research value in further studies on gene function verification. These findings provide novel and detailed information on the M. truncatula PcG family, and in the future it would be helpful to carry out related research on the PcG family in other legumes.

Genome-wide characterization and expression profiling of PDI family gene reveals function as abiotic and biotic stress tolerance in Chinese cabbage (Brassica rapa ssp. pekinensis)

Background

Protein disulfide isomerase (PDI) and PDI-like proteins contain thioredoxin domains that catalyze protein disulfide bond, inhibit aggregation of misfolded proteins, and function in isomerization during protein folding in endoplasmic reticulum and responses during abiotic stresses.Chinese cabbage is widely recognized as an economically important, nutritious vegetable, but its yield is severely hampered by various biotic and abiotic stresses. Because of, it is prime need to identify those genes whose are responsible for biotic and abiotic stress tolerance. PDI family genes are among of them.

Results

We have identified 32 PDI genes from the Br135K microarray dataset, NCBI and BRAD database, and in silico characterized their sequences. Expression profiling of those genes was performed using cDNA of plant samples imposed to abiotic stresses; cold, salt, drought and ABA (Abscisic Acid) and biotic stress; Fusarium oxysporum f. sp. conglutinans infection. The Chinese cabbage PDI genes were clustered in eleven groups in phylogeny. Among them, 15 PDI genes were ubiquitously expressed in various organs, while 24 PDI genes were up-regulated under salt and drought stress. By contrast, cold and ABA stress responsive gene number were ten and nine, respectively. In case of F. oxysporum f. sp. conglutinans infection 14 BrPDI genes were highly up-regulated. Interestingly, BrPDI1–1 gene was identified as putative candidate against abiotic (salt and drought) and biotic stresses, BrPDI5–2 gene for ABA stress, and BrPDI1–4, 6–1 and 9–2 were putative candidate genes for both cold and chilling injury stresses.

Conclusions

Our findings help to elucidate the involvement of PDI genes in stress responses, and they lay the foundation for functional genomics in future studies and molecular breeding of Brassica rapa crops. The stress-responsive PDI genes could be potential resources for molecular breeding of Brassica crops resistant to biotic and abiotic stresses.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4277-2) contains supplementary material, which is available to authorized users.

Molecular characterisation and expression profiling of calcineurin B-like (CBL) genes in Chinese cabbage under abiotic stresses

Calcium signals act as a second messenger in plant responses to various abiotic stresses, which regulate a range of physiological processes. Calcium-binding proteins, like calcineurin B-like (CBL) proteins, belong to a unique group of calcium sensors that play a role in calcium signalling. However, their identities and functions are unknown in Chinese cabbage. In this study, 17 CBL genes were identified from the Brassica rapa L. (Chinese cabbage) database and Br135K microarray datasets. They were used to construct a phylogenetic tree with known CBL proteins of other species. Analysis of genomic distribution and evolution revealed different gene duplication in Chinese cabbage compared to Arabidopsis. The microarray expression analysis showed differential expression of BrCBL genes at various temperatures. Organ-specific expression was observed by RT-PCR, and qRT-PCR analyses revealed responsiveness of BrCBL genes to cold, drought and salt stresses. Our findings confirm that CBL genes are involved in calcium signalling and regulate responses to environmental stimuli, suggesting this family gene have crucial role to play in plant responses to abiotic stresses. The results facilitate selection of candidate genes for further functional characterisation. In addition, abiotic stress-responsive genes reported in this study might be exploited for marker-aided backcrossing of Chinese cabbage.

Genome-wide expression profiling of aquaporin genes confer responses to abiotic and biotic stresses in Brassica rapa

BACKGROUND

Plants contain a range of aquaporin (AQP) proteins, which act as transporter of water and nutrient molecules through living membranes. AQPs also participate in water uptake through the roots and contribute to water homeostasis in leaves.

RESULTS

In this study, we identified 59 AQP genes in the B. rapa database and Br135K microarray dataset. Phylogenetic analysis revealed four distinct subfamilies of AQP genes: plasma membrane intrinsic proteins (PIPs), tonoplast intrinsic proteins (TIPs), NOD26-like intrinsic proteins (NIPs) and small basic intrinsic proteins (SIPs). Microarray analysis showed that the majority of PIP subfamily genes had differential transcript abundance between two B. rapa inbred lines Chiifu and Kenshin that differ in their susceptibility to cold. In addition, all BrPIP genes showed organ-specific expression. Out of 22 genes, 12, 7 and 17 were up-regulated in response to cold, drought and salt stresses, respectively. In addition, 18 BrPIP genes were up-regulated under ABA treatment and 4 BrPIP genes were up-regulated upon F. oxysporum f. sp. conglutinans infection. Moreover, all BrPIP genes showed down-regulation under waterlogging stress, reflecting likely the inactivation of AQPs controlling symplastic water movement.

CONCLUSIONS

This study provides a comprehensive analysis of AQPs in B. rapa and details the expression of 22 members of the BrPIP subfamily. These results provide insight into stress-related biological functions of each PIP gene of the AQP family, which will promote B. rapa breeding programs.

Transcriptome analysis of the variations between autotetraploid Paulownia tomentosa and its diploid using high-throughput sequencing

Timber properties of autotetraploid Paulownia tomentosa are heritable with whole genome duplication, but the molecular mechanisms for the predominant characteristics remain unclear. To illuminate the genetic basis, high-throughput sequencing technology was used to identify the related unigenes. 2677 unigenes were found to be significantly differentially expressed in autotetraploid P. tomentosa. In total, 30 photosynthesis-related, 21 transcription factor-related, and 22 lignin-related differentially expressed unigenes were detected, and the roles of the peroxidase in lignin biosynthesis, MYB DNA-binding proteins, and WRKY proteins associated with the regulation of relevant hormones are extensively discussed. The results provide transcriptome data that may bring a new perspective to explain the polyploidy mechanism in the long growth cycle of plants and offer some help to the future Paulownia breeding.

Duplicate and conquer: multiple homologs of PHOSPHORUS-STARVATION TOLERANCE1 enhance phosphorus acquisition and sorghum performance on low-phosphorus soils

Low soil phosphorus (P) availability is a major constraint for crop production in tropical regions. The rice (Oryza sativa) protein kinase, PHOSPHORUS-STARVATION TOLERANCE1 (OsPSTOL1), was previously shown to enhance P acquisition and grain yield in rice under P deficiency. We investigated the role of homologs of OsPSTOL1 in sorghum (Sorghum bicolor) performance under low P. Association mapping was undertaken in two sorghum association panels phenotyped for P uptake, root system morphology and architecture in hydroponics and grain yield and biomass accumulation under low-P conditions, in Brazil and/or in Mali. Root length and root surface area were positively correlated with grain yield under low P in the soil, emphasizing the importance of P acquisition efficiency in sorghum adaptation to low-P availability. SbPSTOL1 alleles reducing root diameter were associated with enhanced P uptake under low P in hydroponics, whereas Sb03g006765 and Sb03g0031680 alleles increasing root surface area also increased grain yield in a low-P soil. SbPSTOL1 genes colocalized with quantitative trait loci for traits underlying root morphology and dry weight accumulation under low P via linkage mapping. Consistent allelic effects for enhanced sorghum performance under low P between association panels, including enhanced grain yield under low P in the soil in Brazil, point toward a relatively stable role for Sb03g006765 across genetic backgrounds and environmental conditions. This study indicates that multiple SbPSTOL1 genes have a more general role in the root system, not only enhancing root morphology traits but also changing root system architecture, which leads to grain yield gain under low-P availability in the soil.

Molecular evolutionary analysis of the Alfin-like protein family in Arabidopsis lyrata, Arabidopsis thaliana, and Thellungiella halophila

In previous studies, the Alfin1 gene, a transcription factor, enhanced salt tolerance in alfalfa, primarily through altering gene expression levels in the root. Here, we examined the molecular evolution of the Alfin-like (AL) proteins in two Arabidopsis species (A. lyrata and A. thaliana) and a salt-tolerant close relative Thellungiella halophila. These AL-like proteins could be divided into four groups and the two known DUF3594 and PHD-finger domains had co-evolved within each group of genes, irrespective of species, due to gene duplication events in the common ancestor of all three species while gene loss was observed only in T. halophila. To detect whether natural selection acted in the evolution of AL genes, we calculated synonymous substitution ratios (dn/ds) and codon usage statistics, finding positive selection operated on four branches and significant differences in biased codon usage in the AL family between T. halophila and A. lyrata or A. thaliana. Distinctively, only the AL7 branch was under positive selection on the PHD-finger domain and the three members on the branch showed the smallest difference when codon bias was evaluated among the seven clusters. Functional analysis based on transgenic overexpression lines and T-DNA insertion mutants indicated that salt-stress-induced AtAL7 could play a negative role in salt tolerance of A. thaliana, suggesting that adaptive evolution occurred in the members of AL gene family.

Use of digital gene expression to discriminate gene expression differences in early generations of resynthesized Brassica napus and its diploid progenitors

Background

Polyploidy is an important evolutionary mechanism in flowering plants that often induces immediate extensive changes in gene expression through genomic merging and doubling. Brassica napus L. is one of the most economically important polyploid oil crops and has been broadly studied as an example of polyploid crop. RNA-seq is a recently developed technique for transcriptome study, which could be in choice for profiling gene expression pattern in polyploids.

Results

We examined the global gene expression patterns of the first four generations of resynthesized B. napus (F₁–F₄), its diploid progenitors B. rapa and B. oleracea, and natural B. napus using digital gene expression analysis. Almost 42 million clean tags were generated using Illumina technology to produce the expression data for 25959 genes, which account for 63% of the annotated B. rapa genome. More than 56% of the genes were transcribed from both strands, which indicate the importance of RNA-mediated gene regulation in polyploidization. Tag mapping of the B. rapa genome generated 19023, 18547, 24383, 20659, 18881, 20692, and 19955 annotated genes for the B. rapa, B. oleracea, F₁–F₄ of synthesized B. napus, and natural B. napus libraries, respectively. The unambiguous tag-mapped genes in the libraries were functionally categorized via gene ontological analysis. Thousands of differentially expressed genes (DEGs) were identified and revealed the substantial changes in F₁–F₄. Among the 20 most DEGs are DNA binding/transcription factor, cyclin-dependent protein kinase, epoxycarotenoid dioxygenase, and glycine-rich protein. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the DEGs suggested approximately 120 biological pathways.

Conclusions

The systematic deep sequencing analysis provided a comprehensive understanding of the transcriptome complexity of early generations of synthesized B. napus. This information broadens our understanding of the mechanisms of B. napus polyploidization and contributes to molecular and genetic research by enriching the Brassica database.

[Bioinformatics analysis of candidate sequences between two flanking markers linked with resistance gene Ht1 to Exserohilum turcicum in maize]

In order to get Ht1 candidate sequence(s), bioinformatics method was employed to analyze sequences between two flanking markers umc22a and umc122a linked with Ht1. Sixty-three open reading frame (ORF) sequences were found, fourteen of which encoded protein domains. The amino acid sequences encoded by these ORF sequences were blasted with the 24 cloned resistance gene amino acid sequences, and a phylogenetic tree was constructed. Among the fourteen ORF sequences, gpm565a with nucleotide binding site (NBS), containing a lot of random coils, had both high identity and close relationship with Xal resistant to bacterial blight, and was predicted to be related with Ht1. The remaining thirteen ORF sequences had no enough evidence to show any relationships with Ht1 due to lack of conserved domains, low identity or distant relationships.

Transferability of wheat microsatellites to diploid Aegilops species and determination of chromosomal localizations of microsatellites in the S genome

Overall, 253 genomic wheat (Triticum aestivum) microsatellite markers were studied for their transferability to the diploid species Aegilops speltoides, Aegilops longissima, and Aegilops searsii, representing the S genome. In total, 88% of all the analyzed primer pairs of markers derived from the B genome of hexaploid wheat amplified DNA fragments in the genomes of the studied species. The transferability of simple sequence repeat (SSR) markers of the T. aestivum A and D genomes totaled 74%. Triticum aestivum-Ae. speltoides, T. aestivum-Ae. longissima, and T. aestivum-Ae. searsii chromosome addition lines allowed us to determine the chromosomal localizations of 103 microsatellite markers in the Aegilops genomes. The majority of them were localized to homoeologous chromosomes in the genome of Aegilops. Several instances of nonhomoeologous localization of T. aestivum SSR markers in the Aegilops genome were considered to be either amplification of other loci or putative translocations. The results of microsatellite analysis were used to study phylogenetic relationships among the 3 species of the Sitopsis section (Ae. speltoides, Ae. longissima, and Ae. searsii) and T. aestivum. The dendrogram obtained generally reflects the current views on phylogenetic relationships among these species.

Phylogenetic and functional analysis of Arabidopsis RCI2 genes

Six new Arabidopsis thaliana genes (AtRCI2C-H) have been identified that show high homology to AtRCI2A and AtRCI2B. Sequence comparisons revealed that AtRCI2-related genes are widely spread among very different organisms, including other plant species, prokaryotes, fungi, and simply organized animals, and are also organized in gene families. Most RCI2 genes show a similar exon-intron organization, which indicates that they have been structurally conserved during evolution, and encode small, highly hydrophobic proteins containing two putative transmembrane domains. Consistently, the majority of AtRCI2 proteins localize in the plasma membrane. RCI2 proteins exhibit an elevated level of sequence similarity and seem to have evolved from a common ancestor. In spite of their high similarity, conserved subcellular localization, and common origin, experimental evidence is presented suggesting that different RCI2 proteins may have distinct functional roles. Thus, as previously demonstrated for AtRCI2A and AtRCI2B, the newly identified AtRCI2 genes (AtRCI2C-H) are differentially regulated in Arabidopsis organs and in response to abiotic stresses and ABA treatment. Furthermore, only the AtRCI2 proteins that do not contain the C-terminal hydrophilic tail (i.e. AtRCI2A-C and AtRCI2H) are able to complement for the loss of the yeast AtRCI2-related gene PMP3. On the basis of these results, different aspects on the evolution and roles of RCI2 genes are discussed.

Differential lineage-specific amplification of transposable elements is responsible for genome size variation in Gossypium

The DNA content of eukaryotic nuclei (C-value) varies approximately 200,000-fold, but there is only a approximately 20-fold variation in the number of protein-coding genes. Hence, most C-value variation is ascribed to the repetitive fraction, although little is known about the evolutionary dynamics of the specific components that lead to genome size variation. To understand the modes and mechanisms that underlie variation in genome composition, we generated sequence data from whole genome shotgun (WGS) libraries for three representative diploid (n = 13) members of Gossypium that vary in genome size from 880 to 2460 Mb (1C) and from a phylogenetic outgroup, Gossypioides kirkii, with an estimated genome size of 588 Mb. Copy number estimates including all dispersed repetitive sequences indicate that 40%-65% of each genome is composed of transposable elements. Inspection of individual sequence types revealed differential, lineage-specific expansion of various families of transposable elements among the different plant lineages. Copia-like retrotransposable element sequences have differentially accumulated in the Gossypium species with the smallest genome, G. raimondii, while gypsy-like sequences have proliferated in the lineages with larger genomes. Phylogenetic analyses demonstrated a pattern of lineage-specific amplification of particular subfamilies of retrotransposons within each species studied. One particular group of gypsy-like retrotransposon sequences, Gorge3 (Gossypium retrotransposable gypsy-like element), appears to have undergone a massive proliferation in two plant lineages, accounting for a major fraction of genome-size change. Like maize, Gossypium has undergone a threefold increase in genome size due to the accumulation of LTR retrotransposons over the 5-10 Myr since its origin.

Comparative genomics of Brassica oleracea and Arabidopsis thaliana reveal gene loss, fragmentation, and dispersal after polyploidy

We sequenced 2.2 Mb representing triplicated genome segments of Brassica oleracea, which are each paralogous with one another and homologous with a segmentally duplicated region of the Arabidopsis thaliana genome. Sequence annotation identified 177 conserved collinear genes in the B. oleracea genome segments. Analysis of synonymous base substitution rates indicated that the triplicated Brassica genome segments diverged from a common ancestor soon after divergence of the Arabidopsis and Brassica lineages. This conclusion was corroborated by phylogenetic analysis of protein families. Using A. thaliana as an outgroup, 35% of the genes inferred to be present when genome triplication occurred in the Brassica lineage have been lost, most likely via a deletion mechanism, in an interspersed pattern. Genes encoding proteins involved in signal transduction or transcription were not found to be significantly more extensively retained than those encoding proteins classified with other functions, but putative proteins predicted in the A. thaliana genome were underrepresented in B. oleracea. We identified one example of gene loss from the Arabidopsis lineage. We found evidence for the frequent insertion of gene fragments of nuclear genomic origin and identified four apparently intact genes in noncollinear positions in the B. oleracea and A. thaliana genomes.

Comparative analysis of methylthioalkylmalate synthase (MAM) gene family and flanking DNA sequences in Brassica oleracea and Arabidopsis thaliana

Gene BoGSL-PRO is associated with presence of 3-carbon side-chain glucosinolates (GSL). This gene is a member of the methylthioalkylmalate synthase (MAM) gene family. A BAC clone of Brassica oleracea, B21F5, containing this gene, was sequenced, annotated and compared to its corresponding region in Arabidopsis thaliana. Twelve protein-coding genes and 10 transposable elements were found in this clone. The corresponding region in A. thaliana chromosome I has 14 genes and no transposable elements. Analysis of MAM gene family in both species, which also include genes controlling 4-carbon side-chain GSL, separated the genes in two groups based on exon numbers and function. Phylogenetic analysis of the amino acid sequences encoded by these genes suggest that these two groups were produced by a duplication that must have occurred before the divergence of the Rosid and Asterid lineages of angiosperms. Comparison with putative orthologs from several prokaryotes further suggest that the members of the gene family with 10 exons, which encode proteins involved in 4-carbon side-chain GSL biosynthesis, were derived via truncation of the 3' end from ancestral genes more similar in length to those with 12 exons, which encode proteins involved in 3-carbon side-chain GSL biosynthesis. Lower gene density in B. oleracea compared to A. thaliana is due in part to presence of transposable elements (TE) mostly in inter-genic regions.

Molecular evolution of the rice miR395 gene family

MicroRNAs (miRNAs) are 20-22 nucleotide non-coding RNAs that play important roles in plant and animal development. They are usually processed from larger precursors that can form stem-loop structures. Among 20 miRNA families that are conserved between Arabidopsis and rice, the rice miR395 gene family was unique because it was organized into compact clusters that could be transcribed as one single transcript. We show here that in fact this family had four clusters of total 24 genes. Three of these clusters were segmental duplications. They contained miR395 genes of both 120 bp and 66 bp long. However, only the latter was repeatedly duplicated. The fourth cluster contained miR395 genes of two different sizes that could be the consequences of intergenic recombination of genes from the first three clusters. On each cluster, both 1-duplication and 2-duplication histories were observed based on the sequence similarity between miR395 genes, some of which were nearly identical suggesting a recent origin. This was supported by a miR395 locus survey among several species of the genus Oryza, where two clusters were only found in species with an AA genome, the genome of the cultivated rice. A comparative study of the genomic organization of Medicago truncatula miR395 gene family showed significant expansion of intergenic spaces indicating that the originally clustered genes were drifting away from each other. The diverse genomic organizations of a conserved microRNA gene family in different plant genomes indicated that this important negative gene regulation system has undergone dramatic tune-ups in plant genomes.

Comparative analysis of a transposon-rich Brassica oleracea BAC clone with its corresponding sequence in A. thaliana

We compared the sequence of a 96.7 Kb-long BAC clone (B 19 N 3) from Brassica oleracea (broccoli) with its corresponding regions in Arabidopsis thaliana. B 19 N 3 contains eight genes and 15 transposable elements (TEs). The first two genes in this clone, Bo 1 and Bo 2, have its corresponding region at the end of chromosome V of Arabidopsis (24 Mb). The third gene, Bo 3, corresponds to an ortholog at the opposite end (2.6 Mb) of the same chromosome. The other five genes, Bo 4 to Bo 8 also have a corresponding region on the same chromosome but at 7.7 Mb . These five genes are colinear with those found in the corresponding region of Arabidopsis, which contains, however, 15 genes. Therefore, a cluster of 10 genes is missing in B. oleracea clone (B 19 N 3). All five genes in common have the same order and orientation in the genomes of both species. Their 36 exons constituting the eight homologous genes have high conservation in size and sequence identity in both species. Among these, there is a major gene involved in aliphatic glucosinolate biosynthesis, Bo GSL-ELONG (Bo 4). Similar to A. thaliana, this gene, has a tandem duplicate, Bo 5. A contig for this region was constructed by primer walking and BAC-end-sequencing, revealing general gene colinearity between both species. During the 20 million years separating A. thaliana from B. oleracea from a common ancestor both genomes have diverged by chromosomal rearrangements and differential TE activity. These events, in addition to changes in chromosome number are responsible for the evolution of the genomes of both species. In spite of these changes, both species conserve general colinearity for their corresponding genes.

Phenotypic and genotypic expression of self-incompatibility haplotypes in Arabidopsis lyrata suggests unique origin of alleles in different dominance classes

The highly divergent alleles of the SRK gene in outcrossing Arabidopsis lyrata have provided important insights into the evolutionary history of self-incompatibility (SI) alleles and serve as an ideal model for studies of the evolutionary and molecular interactions between alleles in cell-cell recognition systems in general. One tantalizing question is how new specificities arise in systems that require coordination between male and female components. Allelic recruitment via gene conversion has been proposed as one possibility, based on the division of DNA sequences at the SRK locus into two distinctive groups: (1) sequences whose relationships are not well resolved and display the long branch lengths expected for a gene under balancing selection (Class A); and (2) sequences falling into a well-supported group with shorter branch lengths (Class B) that are closely related to an unlinked paralogous locus. The purpose of this study was to determine if differences in phenotype (site of expression assayed using allele-specific reverse transcription-polymerase chain reaction) or function (dominance relationships assayed through controlled pollinations) accompany the sequence-based classification. Expression of Class A alleles was restricted to floral tissues, as predicted for genes involved in the SI response. In contrast, Class B alleles, despite being tightly linked to the SI phenotype, were unexpectedly expressed in both leaves and floral tissues; the same pattern found for a related unlinked paralogous sequence. Whereas Class A included haplotypes in three different dominance classes, all Class B haplotypes were found to be recessive to all except one Class A haplotype. In addition, mapping of expression and dominance patterns onto an S-domain-based genealogy suggested that allelic dominance may be determined more by evolutionary history than by frequency-dependent selection for lowered dominance as some theories suggest. The possibility that interlocus gene conversion might have contributed to allelic diversity is discussed.

Comparative analysis of a Brassica BAC clone containing several major aliphatic glucosinolate genes with its corresponding Arabidopsis sequence

We compared the sequence of a 101-kb-long bacterial artificial chromosome (BAC) clone (B21H13) from Brassica oleracea with its homologous region in Arabidopsis thaliana. This clone contains a gene family involved in the synthesis of aliphatic glucosinolates. The A. thaliana homologs for this gene family are located on chromosome IV and correspond to three 2-oxoglutarate-dependent dioxygenase (AOP) genes. We found that B21H13 harbors 23 genes, whereas the equivalent region in Arabidopsis contains 37 genes. All 23 common genes have the same order and orientation in both Brassica and Arabidopsis. The 16 missing genes in the broccoli BAC clone were arranged in two major blocks of 5 and 7 contiguous genes, two singletons, and a twosome. The 118 exons comprising these 23 genes have high conservation between the two species. The arrangement of the AOP gene family in A. thaliana is as follows: AOP3 (GS-OHP) - AOP2 (GS-ALK) - pseudogene - AOP1. In contrast, in B. oleracea (broccoli and collard), two of the genes are duplicated and the third, AOP3, is missing. The remaining genes are arranged as follows: Bo-AOP2.1 (BoGSL-ALKa) - pseudogene - AOP2.2 (BoGSL-ALKb) - AOP1.1 - AOP1.2. When the survey was expanded to other Brassica accessions, we found variation in copy number and sequence for the Brassica AOP2 homologs. This study confirms that extensive rearrangements have taken place during the evolution of the Brassicacea at both gene and chromosomal levels.

Conservation of the microstructure of genome segments in Brassica napus and its diploid relatives

The cultivated Brassica species are the group of crops most closely related to Arabidopsis thaliana (Arabidopsis). They represent models for the application in crops of genomic information gained in Arabidopsis and provide an opportunity for the investigation of polyploid genome formation and evolution. The scientific literature contains contradictory evidence for the dynamics of the evolution of polyploid genomes. We aimed at overcoming the inherent complexity of Brassica genomes and clarify the effects of polyploidy on the evolution of genome microstructure in specific segments of the genome. To do this, we have constructed bacterial artificial chromosome (BAC) libraries from genomic DNA of B. rapa subspecies trilocularis (JBr) and B. napus var Tapidor (JBnB) to supplement an existing BAC library from B. oleracea. These allowed us to analyse both recent polyploidization (under 10,000 years in B. napus) and more ancient polyploidization events (ca. 20 Myr for B. rapa and B. oleracea relative to Arabidopsis), with an analysis of the events occurring on an intermediate time scale (over the ca. 4 Myr since the divergence of the B. rapa and B. oleracea lineages). Using the Arabidopsis genome sequence and clones from the JBr library, we have analysed aspects of gene conservation and microsynteny between six regions of the genome of B. rapa with the homoeologous regions of the genomes of B. oleracea and Arabidopsis. Extensive divergence of gene content was observed between the B. rapa paralogous segments and their homoeologous segments within the genome of Arabidopsis. A pattern of interspersed gene loss was identified that is similar, but not identical, to that observed in B. oleracea. The conserved genes show highly conserved collinearity with their orthologues across genomes, but a small number of species-specific rearrangements were identified. Thus the evolution of genome microstructure is an ongoing process. Brassica napus is a recently formed polyploid resulting from the hybridization of B. rapa (containing the Brassica A genome) and B. oleracea (containing the Brassica C genome). Using clones from the JBnB library, we have analysed the microstructure of the corresponding segments of the B. napus genome. The results show that there has been little or no change to the microstructure of the analysed segments of the Brassica A and C genomes as a consequence of the hybridization event forming natural B. napus. The observations indicate that, upon polyploid formation, these segments of the genome did not undergo a burst of evolution discernible at the scale of microstructure.

Analyses of LTR-retrotransposon structures reveal recent and rapid genomic DNA loss in rice

We initially analyzed 11 families of low- and middle-copy-number long terminal repeat (LTR) retrotransposons in rice to determine how their structures have diverged from their predicted ancestral forms. These elements, many highly fragmented, were identified on the basis of sequence homology and structural characteristics. The 11 families, totaling 1000 elements, have copy numbers ranging from 1 to 278. Less than one-quarter of these elements are intact, whereas the remaining are solo LTRs and variously truncated fragments. We also analyzed two highly repetitive families (Osr8 and Osr30) of LTR retrotransposons and observed the same results. Our data indicate that unequal homologous recombination and illegitimate recombination are primarily responsible for LTR-retrotransposon removal. Further analysis suggests that most of the detectable LTR retrotransposons in rice inserted less than 8 million years ago, and have now lost over two-thirds of their encoded sequences. Hence, we predict that the half-life of LTR-retrotransposon sequences in rice is less than 6 million years. Moreover, our data demonstrate that at least 22% (97 Mb) of the current rice genome is comprised of LTR-retrotransposon sequences, and that more than 190 Mb of LTR-retrotransposon sequences have been deleted from the rice genome in the last 8 million years.

Characterisation of an Arabidopsis-Leptosphaeria maculans pathosystem: resistance partially requires camalexin biosynthesis and is independent of salicylic acid, ethylene and jasmonic acid signalling

Out of 168 Arabidopsis accessions screened with isolates of Leptosphaeria maculans, one (An-1) showed clear disease symptoms. In order to identify additional components involved in containment of L. maculans in Arabidopsis, a screen for L. maculans-susceptible (lms) mutants was performed. Eleven lms mutants were isolated, which displayed differential susceptibility responses to L. maculans. lms1 was crossed with Columbia (Col-0) and Ws-0, and mapping data for both populations showed the highest linkage to a region on chromosome 2. Reduced levels of PR-1 and PDF1.2 expression were found in lms1 compared to wild-type plants 48 h after pathogen inoculation. In contrast, the lms1 mutant displayed upregulation of either marker gene upon chemical treatment, possibly as an effect of an altered ethylene (ET) response. To assess the contribution of different defence pathways, genotypes implicated in salicylic acid (SA) signalling plants expressing the bacterial salicylate hydroxylase (nahG) gene, non-expressor of PR1 (npr1)-1 and phytoalexin-deficient (pad4-1), jasmonic acid (JA) signalling (coronatine insensitive (coi)1-16, enhanced disease susceptibility (eds)8-1 and jasmonic acid resistant (jar)1-1) and ET signalling (eds4-1, ethylene insensitive (ein)2, ein3-1 and ethylene resistant (etr)1-1) were screened. All the genotypes screened were as resistant as wild-type plants, demonstrating the dispensability of the pathways in L. maculans resistance. When mutants implicated in cell death responses were assayed, responsive to antagonist 1 (ran1)-1 exhibited a weak susceptible phenotype, whereas accelerated cell death (acd)1-20 showed a rapid lesion development. Camalexin is only partially responsible for L. maculans containment in Arabidopsis, as pad3-1 and enhanced susceptibility to Alternaria (esa)1 clearly showed a susceptible response while wild-type levels of camalexin were present in An-1 and lms1. The data presented point to the existence of multiple defence mechanisms controlling the containment of L. maculans in Arabidopsis.

Sequence-based alignment of sorghum chromosome 3 and rice chromosome 1 reveals extensive conservation of gene order and one major chromosomal rearrangement

The completed rice genome sequence will accelerate progress on the identification and functional classification of biologically important genes and serve as an invaluable resource for the comparative analysis of grass genomes. In this study, methods were developed for sequence-based alignment of sorghum and rice chromosomes and for refining the sorghum genetic/physical map based on the rice genome sequence. A framework of 135 BAC contigs spanning approximately 33 Mbp was anchored to sorghum chromosome 3. A limited number of sequences were collected from 118 of the BACs and subjected to BLASTX analysis to identify putative genes and BLASTN analysis to identify sequence matches to the rice genome. Extensive conservation of gene content and order between sorghum chromosome 3 and the homeologous rice chromosome 1 was observed. One large-scale rearrangement was detected involving the inversion of an approximately 59 cM block of the short arm of sorghum chromosome 3. Several small-scale changes in gene collinearity were detected, indicating that single genes and/or small clusters of genes have moved since the divergence of sorghum and rice. Additionally, the alignment of the sorghum physical map to the rice genome sequence allowed sequence-assisted assembly of an approximately 1.6 Mbp sorghum BAC contig. This streamlined approach to high-resolution genome alignment and map building will yield important information about the relationships between rice and sorghum genes and genomic segments and ultimately enhance our understanding of cereal genome structure and evolution.

Comparative mapping between potato (Solanum tuberosum) and Arabidopsis thaliana reveals structurally conserved domains and ancient duplications in the potato genome

A genetic map of potato (Solanum tuberosum) was constructed based on 293 restriction fragment length polymorphism (RFLP) markers including 31 EST markers of Arabidopsis. The in silico comparison of all marker sequences with the Arabidopsis genomic sequence resulted in 189 markers that detected in Arabidopsis 787 loci with sequence conservation. Based on conserved linkage between groups of at least three different markers on the genetic map of potato and the physical map of Arabidopsis, 90 putative syntenic blocks were identified covering 41% of the potato genetic map and 50% of the Arabidopsis physical map. The existence and distribution of syntenic blocks suggested a higher degree of structural conservation in some parts of the potato genome when compared to others. Syntenic blocks were redundant: most potato syntenic blocks were related to several Arabidopsis genome segments and vice versa. Some duplicated potato syntenic blocks correlated well with ancient segmental duplications in Arabidopsis. Syntenic relationships between different genomic segments of potato and the same segment of the Arabidopsis genome indicated that potato genome evolution included ancient intra- and interchromosomal duplications. The partial genome coveridge and the redundancy of syntenic blocks limits the use of synteny for functional comparisons between the crop species potato and the model plant Arabidopsis.

Syntenic relationships between Medicago truncatula and Arabidopsis reveal extensive divergence of genome organization

Arabidopsis and Medicago truncatula represent sister clades within the dicot subclass Rosidae. We used genetic map-based and bacterial artificial chromosome sequence-based approaches to estimate the level of synteny between the genomes of these model plant species. Mapping of 82 tentative orthologous gene pairs reveals a lack of extended macrosynteny between the two genomes, although marker collinearity is frequently observed over small genetic intervals. Divergence estimates based on non-synonymous nucleotide substitutions suggest that a majority of the genes under analysis have experienced duplication in Arabidopsis subsequent to divergence of the two genomes, potentially confounding synteny analysis. Moreover, in cases of localized synteny, genetically linked loci in M. truncatula often share multiple points of synteny with Arabidopsis; this latter observation is consistent with the large number of segmental duplications that compose the Arabidopsis genome. More detailed analysis, based on complete sequencing and annotation of three M. truncatula bacterial artificial chromosome contigs suggests that the two genomes are related by networks of microsynteny that are often highly degenerate. In some cases, the erosion of microsynteny could be ascribed to the selective gene loss from duplicated loci, whereas in other cases, it is due to the absence of close homologs of M. truncatula genes in Arabidopsis.

Structural divergence of chromosomal segments that arose from successive duplication events in the Arabidopsis genome

Using the extensive segmental duplications of the Arabidopsis thaliana genome, a comparative study of homoeologous segments occurring in chromosomes 1, 2, 4 and 5 was performed. The gene-by-gene BLASTP approach was applied to identify duplicated genes in homoeologues. The levels of synonymous substitutions between duplicated coding sequences suggest that these regions were formed by at least two rounds of duplications. Moreover, remnants of even more ancient duplication events were recognised by a whole-genome study. We describe a subchromosomal organisation of genes, including the tandemly repeated genes, and the distribution of transposable elements (TEs). In certain cases, evidence of the possible mechanisms of structural rearrangements within the segments could be found. We provide a probable scenario of the rearrangements that took place during the evolution of the homoeologous regions. Furthermore, on the basis of the comparative analysis of the chromosomal segments in the Columbia and Landsberg erecta accessions, an additional structural variation in the A.thaliana genome is described. Analysis of the segments, spanning 7 Mb or 5.6% of the genome, permitted us to propose a model of evolution at the subchromosomal level.

Comparative genome and proteome analysis of Anopheles gambiae and Drosophila melanogaster

Comparison of the genomes and proteomes of the two diptera Anopheles gambiae and Drosophila melanogaster, which diverged about 250 million years ago, reveals considerable similarities. However, numerous differences are also observed; some of these must reflect the selection and subsequent adaptation associated with different ecologies and life strategies. Almost half of the genes in both genomes are interpreted as orthologs and show an average sequence identity of about 56%, which is slightly lower than that observed between the orthologs of the pufferfish and human (diverged about 450 million years ago). This indicates that these two insects diverged considerably faster than vertebrates. Aligned sequences reveal that orthologous genes have retained only half of their intron/exon structure, indicating that intron gains or losses have occurred at a rate of about one per gene per 125 million years. Chromosomal arms exhibit significant remnants of homology between the two species, although only 34% of the genes colocalize in small "microsyntenic" clusters, and major interarm transfers as well as intra-arm shuffling of gene order are detected.

Insights into cereal genomes from two draft genome sequences of rice

Draft genome sequences have been reported for two subspecies of rice. The drafts include the sequences of an estimated 99% of all rice genes and provide major advances in our understanding of the content and complexity of cereal genomes in general and the rice genome in particular.

Fine mapping in tomato using microsynteny with the Arabidopsis genome: the Diageotropica (Dgt) locus

BACKGROUND

The Arabidopsis thaliana genome sequence provides a catalog of reference genes applicable to comparative microsynteny analysis of other species, facilitating map-based cloning in economically important crops. We have applied such an analysis to the tomato expressed sequence tag (EST) database to expedite high-resolution mapping of the Diageotropica (Dgt) gene within the distal end of chromosome 1 in tomato (Lycopersicon esculentum).

RESULTS

A BLAST search of the Arabidopsis database with nucleotide sequences of markers that flank the tomato dgt locus revealed regions of microsynteny between the distal end of chromosome 1 in tomato, two regions of Arabidopsis chromosome 4, and one on chromosome 2. Tomato ESTs homeologous to Arabidopsis gene sequences within those regions were converted into co-dominant molecular markers via cleaved amplified polymorphic sequence (CAPS) analysis and scored against an informative backcross mapping population. Six new microsyntenic EST (MEST) markers were rapidly identified in the dgt region, two of which further defined the placement of the Dgt gene and permitted the selection of a candidate tomato bacterial artificial chromosome clone for sequence analysis.

CONCLUSIONS

Microsynteny-based comparative mapping combined with CAPS analysis of recombinant plants rapidly and economically narrowed the dgt mapping region from 0.8 to 0.15 cM. This approach should contribute to developing high-density maps of molecular markers to target-specific regions for positional cloning and marker-assisted selection in a variety of plants.

The probability of preservation of a newly arisen gene duplicate

Newly emerging data from genome sequencing projects suggest that gene duplication, often accompanied by genetic map changes, is a common and ongoing feature of all genomes. This raises the possibility that differential expansion/contraction of various genomic sequences may be just as important a mechanism of phenotypic evolution as changes at the nucleotide level. However, the population-genetic mechanisms responsible for the success vs. failure of newly arisen gene duplicates are poorly understood. We examine the influence of various aspects of gene structure, mutation rates, degree of linkage, and population size (N) on the joint fate of a newly arisen duplicate gene and its ancestral locus. Unless there is active selection against duplicate genes, the probability of permanent establishment of such genes is usually no less than 1/(4N) (half of the neutral expectation), and it can be orders of magnitude greater if neofunctionalizing mutations are common. The probability of a map change (reassignment of a key function of an ancestral locus to a new chromosomal location) induced by a newly arisen duplicate is also generally >1/(4N) for unlinked duplicates, suggesting that recurrent gene duplication and alternative silencing may be a common mechanism for generating microchromosomal rearrangements responsible for postreproductive isolating barriers among species. Relative to subfunctionalization, neofunctionalization is expected to become a progressively more important mechanism of duplicate-gene preservation in populations with increasing size. However, even in large populations, the probability of neofunctionalization scales only with the square of the selective advantage. Tight linkage also influences the probability of duplicate-gene preservation, increasing the probability of subfunctionalization but decreasing the probability of neofunctionalization.

Conservation of Microstructure between a Sequenced Region of the Genome of Rice and Multiple Segments of the Genome of Arabidopsis thaliana

The nucleotide sequence was determined for a 340-kb segment of rice chromosome 2, revealing 56 putative protein-coding genes. This represents a density of one gene per 6.1 kb, which is higher than was reported for a previously sequenced segment of the rice genome. Sixteen of the putative genes were supported by matches to ESTs. The predicted products of 29 of the putative genes showed similarity to known proteins, and a further 17 genes showed similarity only to predicted or hypothetical proteins identified in genome sequence data. The region contains a few transposable elements: one retrotransposon, and one transposon. The segment of the rice genome studied had previously been identified as representing a part of rice chromosome 2 that may be homologous to a segment of Arabidopsis chromosome 4. We confirmed the conservation of gene content and order between the two genome segments. In addition, we identified a further four segments of the Arabidopsis genome that contain conserved gene content and order. In total, 22 of the 56 genes identified in the rice genome segment were represented in this set of Arabidopsis genome segments, with at least five genes present, in conserved order, in each segment. These data are consistent with the hypothesis that theArabidopsis genome has undergone multiple duplication events. Our results demonstrate that conservation of the genome microstructure can be identified even between monocot and dicot species. However, the frequent occurrence of duplication, and subsequent microstructure divergence, within plant genomes may necessitate the integration of subsets of genes present in multiple redundant segments to deduce evolutionary relationships and identify orthologous genes.

Conservation of microstructure between a sequenced region of the genome of rice and multiple segments of the genome of Arabidopsis thaliana

The nucleotide sequence was determined for a 340-kb segment of rice chromosome 2, revealing 56 putative protein-coding genes. This represents a density of one gene per 6.1 kb, which is higher than was reported for a previously sequenced segment of the rice genome. Sixteen of the putative genes were supported by matches to ESTs. The predicted products of 29 of the putative genes showed similarity to known proteins, and a further 17 genes showed similarity only to predicted or hypothetical proteins identified in genome sequence data. The region contains a few transposable elements: one retrotransposon, and one transposon. The segment of the rice genome studied had previously been identified as representing a part of rice chromosome 2 that may be homologous to a segment of Arabidopsis chromosome 4. We confirmed the conservation of gene content and order between the two genome segments. In addition, we identified a further four segments of the Arabidopsis genome that contain conserved gene content and order. In total, 22 of the 56 genes identified in the rice genome segment were represented in this set of Arabidopsis genome segments, with at least five genes present, in conserved order, in each segment. These data are consistent with the hypothesis that the Arabidopsis genome has undergone multiple duplication events. Our results demonstrate that conservation of the genome microstructure can be identified even between monocot and dicot species. However, the frequent occurrence of duplication, and subsequent microstructure divergence, within plant genomes may necessitate the integration of subsets of genes present in multiple redundant segments to deduce evolutionary relationships and identify orthologous genes.