Comparative genome analyses of Arabidopsis spp.: inferring chromosomal rearrangement events in the evolutionary history of A. thaliana

Repeat turnover meets stable chromosomes: repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets

Sugar beet and its wild relatives share a base chromosome number of nine and similar chromosome morphologies. Yet, interspecific breeding is impeded by chromosome and sequence divergence that is still not fully understood. Since repetitive DNAs are among the fastest evolving parts of the genome, we investigated, if repeatome innovations and losses are linked to chromosomal differentiation and speciation. We traced genome and chromosome-wide evolution across 13 beet species comprising all sections of the genera Beta and Patellifolia. For this, we combined short and long read sequencing, flow cytometry, and cytogenetics to build a comprehensive framework that spans the complete scale from DNA to chromosome to genome. Genome sizes and repeat profiles reflect the separation into three gene pools with contrasting evolutionary patterns. Among all repeats, satellite DNAs harbor most genomic variability, leading to fundamentally different centromere architectures, ranging from chromosomal uniformity in Beta and Patellifolia to the formation of patchwork chromosomes in Corollinae/Nanae. We show that repetitive DNAs are causal for the genome expansions and contractions across the beet genera, providing insights into the genomic underpinnings of beet speciation. Satellite DNAs in particular vary considerably between beet genomes, leading to the evolution of distinct chromosomal setups in the three gene pools, likely contributing to the barriers in beet breeding. Thus, with their isokaryotypic chromosome sets, beet genomes present an ideal system for studying the link between repeats, genomic variability, and chromosomal differentiation and provide a theoretical fundament for understanding barriers in any crop breeding effort.

Synergid cell calcium oscillations refine understanding of FERONIA/LORELEI signaling during interspecific hybridization

KEY MESSAGE

Pollen tubes from closely related species and mutants lacking pollen tube MYB transcription factors are able to initiate FER/LRE-dependent synergid cell calcium oscillations. Reproductive isolation leads to the evolution of new species; however, the molecular mechanisms that maintain reproductive barriers between sympatric species are not well defined. In flowering plants, sperm cells are immotile and are delivered to female gametes by the pollen grain. After landing on the stigmatic surface, the pollen grain germinates a polarized extension, the pollen tube, into floral tissue. After growing via polar extension to the female gametes and shuttling its cargo of sperm cells through its cytoplasm, the pollen tube signals its arrival and identity to synergid cells that flank the egg. If signaling is successful, the pollen tube and receptive synergid cell burst, and sperm cells are released for fusion with female gametes. To better understand cell-cell recognition during reproduction and how reproductive barriers are maintained between closely related species, pollen tube-initiated synergid cell calcium ion dynamics were examined during interspecific crosses. It was observed that interspecific pollen tubes successfully trigger synergid cell calcium oscillations-a hallmark of reproductive success-but signaling fails downstream of key signaling genes and sperm are not released. This work further defines pollen tube-synergid cell signaling as a critical block to interspecific hybridization and suggests that the FERONIA/LORELEI signaling mechanism plays multiple parallel roles during pollen tube reception.

Cytological heterogeneity of heterochromatin among 10 sequenced Drosophila species

In Drosophila chromosomal rearrangements can be maintained and are associated with karyotypic variability among populations from different geographic localities. The abundance of variability in gene arrangements among chromosomal arms is even greater when comparing more distantly related species and the study of these chromosomal changes has provided insights into the evolutionary history of species in the genus. In addition, the sequencing of genomes of several Drosophila species has offered the opportunity to establish the global pattern of genomic evolution, at both genetic and chromosomal level. The combined approaches of comparative analysis of syntenic blocks and direct physical maps on polytene chromosomes have elucidated changes in the orientation of genomic sequences and the difference between heterochromatic and euchromatic regions. Unfortunately, the centromeric heterochromatic regions cannot be studied using the cytological maps of polytene chromosomes because they are underreplicated and therefore reside in the chromocenter. In Drosophila melanogaster, a cytological map of the heterochromatin has been elaborated using mitotic chromosomes from larval neuroblasts. In the current work, we have expanded on that mapping by producing cytological maps of the mitotic heterochromatin in an additional 10 sequenced Drosophila species. These maps highlight 2 apparently different paths, for the evolution of the pericentric heterochromatin between the subgenera Sophophora and Drosophila. One path leads toward a progressive complexity of the pericentric heterochromatin (Sophophora) and the other toward a progressive simplification (Drosophila). These maps are also useful for a better understanding how karyotypes have been altered by chromosome arm reshuffling during evolution.

Species-specific partial gene duplication in Arabidopsis thaliana evolved novel phenotypic effects on morphological traits under strong positive selection

Gene duplication is increasingly recognized as an important mechanism for the origination of new genes, as revealed by comparative genomic analysis. However, how new duplicate genes contribute to phenotypic evolution remains largely unknown, especially in plants. Here, we identified the new gene EXOV, derived from a partial gene duplication of its parental gene EXOVL in Arabidopsis thaliana. EXOV is a species-specific gene that originated within the last 3.5 million years and shows strong signals of positive selection. Unexpectedly, RNA-sequencing analyses revealed that, despite its young age, EXOV has acquired many novel direct and indirect interactions in which the parental gene does not engage. This observation is consistent with the high, selection-driven substitution rate of its encoded protein, in contrast to the slowly evolving EXOVL, suggesting an important role for EXOV in phenotypic evolution. We observed significant differentiation of morphological changes for all phenotypes assessed in genome-edited and T-DNA insertional single mutants and in double T-DNA insertion mutants in EXOV and EXOVL. We discovered a substantial divergence of phenotypic effects by principal component analyses, suggesting neofunctionalization of the new gene. These results reveal a young gene that plays critical roles in biological processes that underlie morphological evolution in A. thaliana.

Using carrot centromeric repeats to study karyotype relationships in the genus Daucus (Apiaceae)

Background

In the course of evolution, chromosomes undergo evolutionary changes; thus, karyotypes may differ considerably among groups of organisms, even within closely related taxa. The genus Daucus seems to be a promising model for exploring the dynamics of karyotype evolution. It comprises some 40 wild species and the cultivated carrot, a crop of great economic significance. However, Daucus species are very diverse morphologically and genetically, and despite extensive research, the taxonomic and phylogenetic relationships between them have still not been fully resolved. Although several molecular cytogenetic studies have been conducted to investigate the chromosomal structure and karyotype evolution of carrot and other Daucus species, detailed karyomorphological research has been limited to carrot and only a few wild species. Therefore, to better understand the karyotype relationships within Daucus, we (1) explored the chromosomal distribution of carrot centromeric repeats (CentDc) in 34 accessions of Daucus and related species by means of fluorescence in situ hybridization (FISH) and (2) performed detailed karyomorphological analysis in 16 of them.

Results

We determined the genomic organization of CentDc in 26 accessions of Daucus (belonging to both Daucus I and II subclades) and one accession of closely related species. The CentDc repeats were present in the centromeric regions of all chromosomes of 20 accessions (representing 11 taxa). In the other Daucus taxa, the number of chromosome pairs with CentDc signals varied depending on the species, yet their centromeric localization was conserved. In addition, precise chromosome measurements performed in 16 accessions showed the inter- and intraspecific karyological relationships among them.

Conclusions

The presence of the CentDc repeats in the genomes of taxa belonging to both Daucus subclades and one outgroup species indicated the ancestral status of the repeat. The results of our study provide useful information for further evolutionary, cytotaxonomic, and phylogenetic research on the genus Daucus and may contribute to a better understanding of the dynamic evolution of centromeric satellites in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07853-2.

Metal and Metalloid Toxicity in Plants: An Overview on Molecular Aspects

Worldwide, the effects of metal and metalloid toxicity are increasing, mainly due to anthropogenic causes. Soil contamination ranks among the most important factors, since it affects crop yield, and the metals/metalloids can enter the food chain and undergo biomagnification, having concomitant effects on human health and alterations to the environment. Plants have developed complex mechanisms to overcome these biotic and abiotic stresses during evolution. Metals and metalloids exert several effects on plants generated by elements such as Zn, Cu, Al, Pb, Cd, and As, among others. The main strategies involve hyperaccumulation, tolerance, exclusion, and chelation with organic molecules. Recent studies in the omics era have increased knowledge on the plant genome and transcriptome plasticity to defend against these stimuli. The aim of the present review is to summarize relevant findings on the mechanisms by which plants take up, accumulate, transport, tolerate, and respond to this metal/metalloid stress. We also address some of the potential applications of biotechnology to improve plant tolerance or increase accumulation.

Multiple and independent rearrangements revealed by comparative cytogenetic mapping in the dysploid Leptostachyus group (Phaseolus L., Leguminosae)

Polyploidy and dysploidy have been reported as the main events in karyotype evolution of plants. In the genus Phaseolus L. (2n = 22), a small monophyletic group of three species, the Leptostachyus group, presents a dysploid karyotype with 2n = 20. It was shown in Phaseolus leptostachyus that the dysploidy was caused by a nested chromosome fusion (NCF) accompanied by several translocations, suggesting a high rate of karyotype evolution in the group. To verify if this karyotype restructuring was a single event or occurred progressively during the evolution of this group, we analysed P. macvaughii, sister to Phaseolus micranthus + P. leptostachyus. Twenty-four genomic clones of P. vulgaris previously mapped on P. leptostachyus, in addition to the 5S and 35S rDNA probes, were used for fluorescence in situ hybridization. Only a single rearrangement was common to the two species: the nested chromosome fusion (NCF) involving chromosomes 10 and 11. The translocation of chromosome 2 is not the same found in P. leptostachyus, and pericentric inversions in chromosomed 3 and 4 were exclusive of P. macvaughii. The other rearrangements observed in P. leptostachyus were not shared with this species, suggesting that they occurred after the separation of these lineages. The presence of private rearrangements indicates a progressive accumulation of karyotype changes in the Leptostachyus group instead of an instant genome-wide repatterning.

Rolling Circle Amplification (RCA)-Mediated Genome-Wide ihpRNAi Mutant Library Construction in Brassica napus

With the successful completion of genomic sequencing for Brassica napus, identification of novel genes, determination of functions performed by genes, and exploring the molecular mechanisms underlying important agronomic traits were challenged. Mutagenesis-based functional genomics techniques including chemical, physical, and insertional mutagenesis have been used successfully in the functional characterization of genes. However, these techniques had their disadvantages and inherent limitations for allopolyploid Brassica napus, which contained a large number of homologous and redundant genes. Long intron-spliced hairpin RNA (ihpRNA) constructs which contained inverted repeats of the target gene separated by an intron, had been shown to be very effective in triggering RNAi in plants. In the present study, the genome-wide long ihpRNA library of B. napus was constructed with the rolling circle amplification (RCA)-mediated technology. Using the phytoene desaturase (PDS) gene as a target control, it was shown that the RCA-mediated long ihpRNA construct was significantly effective in triggering gene silence in B. napus. Subsequently, the resultant long ihpRNA library was transformed into B. napus to produce corresponding RNAi mutants. Among the obtained transgenic ihpRNA population of B. napus, five ihpRNA lines with observable mutant phenotypes were acquired including alterations in the floral model and the stamen development. The target genes could be quickly identified using specific primers. These results showed that the RCA-mediated ihpRNA construction method was effective for the genome-wide long ihpRNA library of B. napus, therefore providing a platform for study of functional genomics in allopolyploid B. napus.

High-Density Genetic Linkage Mapping of Lepidium Based on Genotyping-by-Sequencing SNPs and Segregating Contig Tag Haplotypes

Lepidium campestre has been targeted for domestication as future oilseed and catch crop. Three hundred eighty plants comprising genotypes of L. campestre, Lepidium heterophyllum, and their interspecific F2 mapping population were genotyped using genotyping by sequencing (GBS), and the generated polymorphic markers were used for the construction of high-density genetic linkage map. TASSEL-GBS, a reference genome-based pipeline, was used for this analysis using a draft L. campestre whole genome sequence. The analysis resulted in 120,438 biallelic single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) above 0.01. The construction of genetic linkage map was conducted using MSTMap based on phased SNPs segregating in 1:2:1 ratio for the F2 individuals, followed by genetic mapping of segregating contig tag haplotypes as dominant markers against the linkage map. The final linkage map consisted of eight linkage groups (LGs) containing 2,330 SNP markers and spanned 881 Kosambi cM. Contigs (10,302) were genetically mapped to the eight LGs, which were assembled into pseudomolecules that covered a total of ∼120.6 Mbp. The final size of the pseudomolecules ranged from 9.4 Mbp (LG-4) to 20.4 Mpb (LG-7). The following major correspondence between the eight Lepidium LGs (LG-1 to LG-8) and the five Arabidopsis thaliana (At) chromosomes (Atx-1-Atx-5) was revealed through comparative genomics analysis: LG-1&2_Atx-1, LG-3_Atx-2&3, LG-4_Atx-2, LG-5_Atx-2&Atx-3, LG-6_Atx-4&5, LG-7_Atx-4, and LG-8_Atx-5. This analysis revealed that at least 66% of the sequences of the LGs showed high collinearity with At chromosomes. The sequence identity between the corresponding regions of the LGs and At chromosomes ranged from 80.6% (LG-6) to 86.4% (LG-8) with overall mean of 82.9%. The map positions on Lepidium LGs of the homologs of 24 genes that regulate various traits in A. thaliana were also identified. The eight LGs revealed in this study confirm the previously reported (1) haploid chromosome number of eight in L. campestre and L. heterophyllum and (2) chromosomal fusion, translocation, and inversion events during the evolution of n = 8 karyotype in ancestral species shared by Lepidium and Arabidopsis to n = 5 karyotype in A. thaliana. This study generated highly useful genomic tools and resources for Lepidium that can be used to accelerate its domestication.

Ancestral Reconstruction of Karyotypes Reveals an Exceptional Rate of Nonrandom Chromosomal Evolution in Sunflower

Mapping the chromosomal rearrangements between species can inform our understanding of genome evolution, reproductive isolation, and speciation. Here, we present a novel algorithm for identifying regions of synteny in pairs of genetic maps, which is implemented in the accompanying R package syntR. The syntR algorithm performs as well as previous ad hoc methods while being systematic, repeatable, and applicable to mapping chromosomal rearrangements in any group of species. In addition, we present a systematic survey of chromosomal rearrangements in the annual sunflowers, which is a group known for extreme karyotypic diversity. We build high-density genetic maps for two subspecies of the prairie sunflower, Helianthus petiolaris ssp. petiolaris and H. petiolaris ssp. fallax Using syntR, we identify blocks of synteny between these two subspecies and previously published high-density genetic maps. We reconstruct ancestral karyotypes for annual sunflowers using those synteny blocks and conservatively estimate that there have been 7.9 chromosomal rearrangements per million years, a high rate of chromosomal evolution. Although the rate of inversion is even higher than the rate of translocation in this group, we further find that every extant karyotype is distinguished by between one and three translocations involving only 8 of the 17 chromosomes. This nonrandom exchange suggests that specific chromosomes are prone to translocation and may thus contribute disproportionately to widespread hybrid sterility in sunflowers. These data deepen our understanding of chromosome evolution and confirm that Helianthus has an exceptional rate of chromosomal rearrangement that may facilitate similarly rapid diversification.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Computational Sequence Analysis of Inversion Breakpoint Regions in the Cactophilic Drosophila mojavensis Lineage

We investigated rates of chromosomal evolution in Drosophila mojavensis using whole-genome sequence information from D. mojavensis, Drosophila buzzatii, and Drosophila virilis. Drosophila mojavensis is a cactophilic species of the repleta group living under extreme ecological conditions in the deserts of the Southwestern United States and Northwestern México. The genome of D. buzzatii, another member of the repleta group, was recently sequenced and the largest scaffolds anchored to all chromosomes using diverse procedures. Chromosome organization between D. mojavensis and D. buzzatii was compared using MUMmer and GRIMM software. Our results corroborate previous cytological analyses that indicated chromosome 2 differed between these 2 species by 10 inversions, chromosomes X and 5 differed by one inversion each, and chromosome 4 was homosequential. In contrast, we found that chromosome 3 differed by 5 inversions instead of the expected 2 that were previously inferred by cytological analyses. Thirteen of these inversions occurred in the D. mojavensis lineage: 12 are fixed and one of them is a polymorphic inversion previously described in populations from Sonora and Baja California, México. We previously investigated the breakpoints of chromosome 2 inversions fixed in D. mojavensis. Here we characterized the breakpoint regions of the 5 inversions found in chromosome 3 in order to infer the molecular mechanism that generated each inversion and its putative functional consequences. Overall, our results reveal a number of gene alterations at the inversion breakpoints with putative adaptive consequences that point to natural selection as the cause for fast chromosomal evolution in D. mojavensis.

Mutations in the CDS and promoter of BjuA07.CLV1 cause a multilocular trait in Brassica juncea

Multilocular trait has recently attracted considerable attention for its potential to increase yield. Our previous studies indicated that two genes (Bjln1 and Bjln2) are responsible for multilocular siliques in Brassica juncea and the Bjln1 gene has been delimited to a 208-kb region. In present study, the Bjln1 gene was successfully isolated using the map-based cloning method. Complementation test indicated that the BjuA07.CLV1 (equivalent to BjLn1) could rescue the multilocular phenotype and generate bilocular siliques. Two amino acids changes at positions 28 and 63 in BjuA07.clv1 as well as a 702-bp deletion in its promoter have been proved to affect the carpel numbers. Microscopic analyses suggested that BjuA07.CLV1 is involved in the maintenance of shoot and floral meristem size. The expression level of BjuA07.clv1 was significantly reduced in the SAM. Furthermore, WUS, CLV2, CLV3, RPK2 and POL, key genes in the CLV/WUS signal pathway, showed lower expression level in the multilocular plants. These data suggest that the mutations in the CDS and promoter of BjuA07.clv1 reduced its function and expression level, which disturbed CLV/WUS signal pathway, thereby leading to the enlargement of the shoot and floral meristem and resulting in the multilocular siliques.

A naturally occurring epiallele associates with leaf senescence and local climate adaptation in Arabidopsis accessions

Epigenetic variation has been proposed to facilitate adaptation to changing environments, but evidence that natural epialleles contribute to adaptive evolution has been lacking. Here we identify a retrotransposon, named “NMR19” (naturally occurring DNA methylation variation region 19), whose methylation and genomic location vary among Arabidopsis thaliana accessions. We classify NMR19 as NMR19-4 and NMR19-16 based on its location, and uncover NMR19-4 as an epiallele that controls leaf senescence by regulating the expression of PHEOPHYTIN PHEOPHORBIDE HYDROLASE (PPH). We find that the DNA methylation status of NMR19-4 is stably inherited and independent of genetic variation. In addition, further analysis indicates that DNA methylation of NMR19-4 correlates with local climates, implying that NMR19-4 is an environmentally associated epiallele. In summary, we discover a novel epiallele, and provide mechanistic insights into its origin and potential function in local climate adaptation.

Epigenetic variation underlies various aspects of phenotypic diversity of plants. Here, He et al show a naturally occurring epiallele controls Arabidopsis leaf senescence by regulating the expression of PHEOPHYTIN PHEOPHORBIDE HYDROLASE (PPH), and is associated with local climate adaptation.

Chromosome Evolution in Connection with Repetitive Sequences and Epigenetics in Plants

Chromosome evolution is a fundamental aspect of evolutionary biology. The evolution of chromosome size, structure and shape, number, and the change in DNA composition suggest the high plasticity of nuclear genomes at the chromosomal level. Repetitive DNA sequences, which represent a conspicuous fraction of every eukaryotic genome, particularly in plants, are found to be tightly linked with plant chromosome evolution. Different classes of repetitive sequences have distinct distribution patterns on the chromosomes. Mounting evidence shows that repetitive sequences may play multiple generative roles in shaping the chromosome karyotypes in plants. Furthermore, recent development in our understanding of the repetitive sequences and plant chromosome evolution has elucidated the involvement of a spectrum of epigenetic modification. In this review, we focused on the recent evidence relating to the distribution pattern of repetitive sequences in plant chromosomes and highlighted their potential relevance to chromosome evolution in plants. We also discussed the possible connections between evolution and epigenetic alterations in chromosome structure and repatterning, such as heterochromatin formation, centromere function, and epigenetic-associated transposable element inactivation.

Proteomic Changes Between Populus Allotriploids and Diploids Revealed Using an iTRAQ-based Quantitative Approach

Background

Polyploid breeding is a powerful approach for Populus genetic improve-ment because polyploid trees have valuable characteristics, including better timber quality and a higher degree of stress resistance compared with their full-sib diploids. However, the genetic mech-anism underlying this phenomenon remains unknown.

Objective

To better understand the proteomic changes between Populus allotriploids and diploids, we examined the proteomic profiles of allotriploid and diploid Populus by iTRAQ labeling coupled with two-dimensional liquid chromatography and MALDI-TOF/TOF mass spectrometry (MS).

Method

iTRAQ labeling coupled with two-dimensional liquid chromatography and MALDI-TOF/TOF mass spectrometry (MS).

Results

Between the Populus allotriploid and the full-sib diploid, 932 differentially expressed proteins (DEPs) were identified. These DEPs were primarily involved in stress, defense, transportation, transcriptional and/or translational modification, and energy production. The pathway analysis indi-cated that most of the DEPs were implicated in carbohydrate transport and metabolism, nitrogen me-tabolism and glycolysis, and the ribosome assembly pathway. These data suggest high protein di-vergence between Populus allotriploids and diploids, and rapid changes during hybridization.

Conclusion

The results provide new data for further understanding of the mechanisms of polyploid trees that generally display increased height growth compared with their full-sib diploids.

Re-sequencing transgenic plants revealed rearrangements at T-DNA inserts, and integration of a short T-DNA fragment, but no increase of small mutations elsewhere

KEY MESSAGE

Transformation resulted in deletions and translocations at T-DNA inserts, but not in genome-wide small mutations. A tiny T-DNA splinter was detected that probably would remain undetected by conventional techniques. We investigated to which extent Agrobacterium tumefaciens-mediated transformation is mutagenic, on top of inserting T-DNA. To prevent mutations due to in vitro propagation, we applied floral dip transformation of Arabidopsis thaliana. We re-sequenced the genomes of five primary transformants, and compared these to genomic sequences derived from a pool of four wild-type plants. By genome-wide comparisons, we identified ten small mutations in the genomes of the five transgenic plants, not correlated to the positions or number of T-DNA inserts. This mutation frequency is within the range of spontaneous mutations occurring during seed propagation in A. thaliana, as determined earlier. In addition, we detected small as well as large deletions specifically at the T-DNA insert sites. Furthermore, we detected partial T-DNA inserts, one of these a tiny 50-bp fragment originating from a central part of the T-DNA construct used, inserted into the plant genome without flanking other T-DNA. Because of its small size, we named this fragment a T-DNA splinter. As far as we know this is the first report of such a small T-DNA fragment insert in absence of any T-DNA border sequence. Finally, we found evidence for translocations from other chromosomes, flanking T-DNA inserts. In this study, we showed that next-generation sequencing (NGS) is a highly sensitive approach to detect T-DNA inserts in transgenic plants.

Draft Genome Sequence of Mentha longifolia and Development of Resources for Mint Cultivar Improvement

The genus Mentha encompasses mint species cultivated for their essential oils, which are formulated into a vast array of consumer products. Desirable oil characteristics and resistance to the fungal disease Verticillium wilt are top priorities for the mint industry. However, cultivated mints have complex polyploid genomes and are sterile. Breeding efforts, therefore, require the development of genomic resources for fertile mint species. Here, we present draft de novo genome and plastome assemblies for a wilt-resistant South African accession of Mentha longifolia (L.) Huds., a diploid species ancestral to cultivated peppermint and spearmint. The 353 Mb genome contains 35 597 predicted protein-coding genes, including 292 disease resistance gene homologs, and nine genes determining essential oil characteristics. A genetic linkage map ordered 1397 genome scaffolds on 12 pseudochromosomes. More than two million simple sequence repeats were identified, which will facilitate molecular marker development. The M. longifolia genome is a valuable resource for both metabolic engineering and molecular breeding. This is exemplified by employing the genome sequence to clone and functionally characterize the promoters in a peppermint cultivar, and demonstrating the utility of a glandular trichome-specific promoter to increase expression of a biosynthetic gene, thereby modulating essential oil composition.

Altered Transcription and Neofunctionalization of Duplicated Genes Rescue the Harmful Effects of a Chimeric Gene in Brassica napus

Chimeric genes contribute to the evolution of diverse functions in plants and animals. However, new chimeric genes also increase the risk of developmental defects. Here, we show that the chimeric gene Brassica napus male sterile 4 (Bnams4^b ) is responsible for genic male sterility in the widely used canola line 7365A (Bnams3 ms3ms4^bms4^b ). Bnams4^b originated via exon shuffling ∼4.6 million years ago. It causes defects in the normal functions of plastids and induces aborted anther formation and/or albino leaves and buds. Evidence of the age of the mutation, its tissue expression pattern, and its sublocalization indicated that it coevolved with BnaC9.Tic40 (BnaMs3). In Arabidopsis thaliana, Bnams4^b results in complete male sterility that can be rescued by BnaC9.Tic40, suggesting that BnaC9.Tic40 might restore fertility through effects on protein level. Another suppressor gene, Bnams4^a , rescues sterility by reducing the level of transcription of Bnams4^b Our results suggest that Brassica plants have coevolved altered transcription patterns and neofunctionalization of duplicated genes that can block developmental defects resulting from detrimental chimeric genes.

Microsynteny and phylogenetic analysis of tandemly organised miRNA families across five members of Brassicaceae reveals complex retention and loss history

Plant genomes are characterized by the presence of large miRNA gene families which are few in number. The expansion of miRNA families is thought to be driven by gene and genome duplication. Some members of these miRNA gene families are tandemly arranged and their analysis is of interest because such organisation may indicate origin through tandem duplication and also to investigate whether some such tandem clusters have similar expression patterns, and whether these are regulated through a common set of cis-regulatory elements (eg. promoters and enhancers). As a first step, we undertake a comprehensive study using micro-synteny analyses of tandemly organised miRNA families across the Brassicaceae spanning an evolutionary time scale of ca. 45 million years, among Arabidopsis, Capsella, Brassica and Thellungiella species, to address the following questions: Are most miRNA gene families present as tandem clusters? To what extent are these tandem patterns retained? To what extent can family sizes be ascribed to genome duplication? Our analysis of thirteen tandemly organised miRNA families revealed that synteny is largely conserved among Arabidopsis thaliana, A. lyrata and Capsella rubella, which form a clade spanning approximately between 6.2-9.8 my (Acarkan et al., 2000) [1]. On the other hand, comparison of sequences from these species with Brassica rapa, B. oleracea and Thellungiella halophila, which form a separate clade spanning 31 my (Franzke et al., 2011)[2] reveals many differences. The latter clade reveals several paralogous duplications that probably resulted from whole genome duplication, as well as disrupted synteny. Phylogenetic analyses of precursor sequences generally support the history inferred from synteny analysis. Synteny and phylogenetic analysis of six members of the tandemly organised miR169 family suggest that the Brassicaceae ancestral state consisted of a "dimer as a unit" which may have undergone direct local duplication to retain the transcriptional orientation followed by lineage specific changes. MiR169, to the best of our knowledge, is one of the largest tandemly organised miRNA gene family across plant kingdom and further analysis should reveal the generality of this pattern of evolution. The conserved organisation of miR395A-B-C and miR395 D-E-F as two clusters on same chromosome/scaffold across A. thaliana, B. rapa and salsuginea demonstrates retention of the large chromosomal segment across the two lineages. MiRNA family miR845 was detected only in Arabidopsis species and Thellungiella indicating a complex loss and retention history. MiR447A-B family was only found in A. thaliana indicating that it is a species-specific gene family of recent origin.

Chromosome number reduction in the sister clade of Carica papaya with concomitant genome size doubling

PREMISE OF THE STUDY

Caricaceae include six genera and 34 species, among them papaya, a model species in plant sex chromosome research. The family was held to have a conserved karyotype with 2n = 18 chromosomes, an assumption based on few counts. We examined the karyotypes and genome size of species from all genera to test for possible cytogenetic variation.

METHODS

We used fluorescent in situ hybridization using standard telomere, 5S, and 45S rDNA probes. New and published data were combined with a phylogeny, molecular clock dating, and C values (available for ∼50% of the species) to reconstruct genome evolution.

KEY RESULTS

The African genus Cylicomorpha, which is sister to the remaining Caricaceae (all neotropical), has 2n = 18, as do the species in two other genera. A Mexican clade of five species that includes papaya, however, has 2n = 18 (papaya), 2n = 16 (Horovitzia cnidoscoloides), and 2n = 14 (Jarilla caudata and J. heterophylla; third Jarilla not counted), with the phylogeny indicating that the dysploidy events occurred ∼16.6 and ∼5.5 million years ago and that Jarilla underwent genome size doubling (∼450 to 830-920 Mbp/haploid genome). Pericentromeric interstitial telomere repeats occur in both Jarilla adjacent to 5S rDNA sites, and the variability of 5S rDNA sites across all genera is high.

CONCLUSIONS

On the basis of outgroup comparison, 2n = 18 is the ancestral number, and repeated chromosomal fusions with simultaneous genome size increase as a result of repetitive elements accumulating near centromeres characterize the papaya clade. These results have implications for ongoing genome assemblies in Caricaceae.

Synteny and comparative analysis of miRNA retention, conservation, and structure across Brassicaceae reveals lineage- and sub-genome-specific changes

The recent availability of genome sequences together with syntenic block information for Brassicaceae offers an opportunity to study microRNA (miRNA) evolution across this family. We employed a synteny-based comparative genomics strategy to unambiguously identify miRNA homologs from the genome sequence of members of Brassicaceae. Such an analysis of miRNA across Brassicaceae allowed us to classify miRNAs as conserved, lineage-, karyotype- and sub-genome-specific. The differential loss of miRNA from sub-genomes in polyploid genomes of Brassica rapa and Brassica oleracea shows that miRNA also follows the rules of gene fractionation as observed in the case of protein-coding genes. The study of mature and miR* region of precursors revealed instances of in-dels and SNPs which reflect the evolutionary history of the genomes. High level of conservation in miR* regions in some cases points to their functional relevance which needs to be further investigated. We further show that sequence and length variability in precursor sequences can affect the free energy and foldback structure of miRNA which may ultimately affect their biogenesis and expression in the biological system.

Functional analysis of the three HMA4 copies of the metal hyperaccumulator Arabidopsis halleri

In Arabidopsis halleri, the HMA4 gene has an essential function in Zn/Cd hypertolerance and hyperaccumulation by mediating root-to-shoot translocation of metals. Constitutive high expression of AhHMA4 results from a tandem triplication and cis-activation of the promoter of all three copies. The three AhHMA4 copies possess divergent promoter sequences, but highly conserved coding sequences, and display identical expression profiles in the root and shoot vascular system. Here, an AhHMA4::GFP fusion was expressed under the control of each of the three A. halleri HMA4 promoters in a hma2hma4 double mutant of A. thaliana to individually examine the function of each AhHMA4 copy. The protein showed non-polar localization at the plasma membrane of the root pericycle cells of both A. thaliana and A. halleri. The expression of each AhHMA4::GFP copy complemented the severe Zn-deficiency phenotype of the hma2hma4 mutant by restoring root-to-shoot translocation of Zn. However, each copy had a different impact on metal homeostasis in the A. thaliana genetic background: AhHMA4 copies 2 and 3 were more highly expressed and provided higher Zn tolerance in roots and accumulation in shoots than copy 1, and AhHMA4 copy 3 also increased Cd tolerance in roots. These data suggest a certain extent of functional differentiation among the three A. halleri HMA4 copies, stemming from differences in expression levels rather than in expression profile. HMA4 is a key node of the Zn homeostasis network and small changes in expression level can have a major impact on Zn allocation to root or shoot tissues.

Zinc triggers a complex transcriptional and post-transcriptional regulation of the metal homeostasis gene FRD3 in Arabidopsis relatives

In Arabidopsis thaliana, FRD3 (FERRIC CHELATE REDUCTASE DEFECTIVE 3) plays a central role in metal homeostasis. FRD3 is among a set of metal homeostasis genes that are constitutively highly expressed in roots and shoots of Arabidopsis halleri, a zinc hyperaccumulating and hypertolerant species. Here, we examined the regulation of FRD3 by zinc in both species to shed light on the evolutionary processes underlying the evolution of hyperaccumulation in A. halleri. We combined gene expression studies with the use of β-glucuronidase and green fluorescent protein reporter constructs to compare the expression profile and transcriptional and post-transcriptional regulation of FRD3 in both species. The AtFRD3 and AhFRD3 genes displayed a conserved expression profile. In A. thaliana, alternative transcription initiation sites from two promoters determined transcript variants that were differentially regulated by zinc supply in roots and shoots to favour the most highly translated variant under zinc-excess conditions. In A. halleri, a single transcript variant with higher transcript stability and enhanced translation has been maintained. The FRD3 gene thus undergoes complex transcriptional and post-transcriptional regulation in Arabidopsis relatives. Our study reveals that a diverse set of mechanisms underlie increased gene dosage in the A. halleri lineage and illustrates how an environmental challenge can alter gene regulation.

Genomic variation in Arabidopsis: tools and insights from next-generation sequencing

The release of a reference genome for Arabidopsis thaliana in 2000 has been an enormous boon for the study of plant genetics. Less than a decade later, however, a revolution in sequencing technology had enabled rapid and inexpensive re-sequencing of whole A. thaliana genomes. Large-scale efforts to characterize natural genomic variation in A. thaliana have revealed remarkable intra-specific variation in this species, ranging from single-nucleotide differences to large structural rearrangements. The partitioning of this variation by geography and local adaptation has been described using powerful new methods and tools. Simultaneously, an ambitious research agenda has emerged to sequence 1001 A. thaliana lines from around the world, while sequencing of related species is enabling powerful evolutionary genomic analyses. In this review, I summarize recent progress in genomic analysis of natural variation in A. thaliana and its close relatives. This progress has set the stage for the emergence of Arabidopsis as a model genus for evolutionary and functional genomics.

Telomere-centric genome repatterning determines recurring chromosome number reductions during the evolution of eukaryotes

Whole-genome duplication (WGD) is central to the evolution of many eukaryotic genomes, in particular rendering angiosperm (flowering plant) genomes much less stable than those of animals. Following repeated duplication/triplication(s), angiosperm chromosome numbers have usually been restored to a narrow range, as one element in a 'diploidization' process that re-establishes diploid heredity. In several angiosperms affected by WGD, we show that chromosome number reduction (CNR) is best explained by intra- and/or inter-chromosomal crossovers to form new chromosomes that utilize the existing telomeres of 'invaded' and centromeres of 'invading' chromosomes, the alternative centromeres and telomeres being lost. Comparison with the banana (Musa acuminata) genome supports a 'fusion model' for the evolution of rice (Oryza sativa) chromosomes 2 and 3, implying that the grass common ancestor had seven chromosomes rather than the five implied by a 'fission model.' The 'invading' and 'invaded' chromosomes are frequently homoeologs, originating from duplication of a common ancestral chromosome and with greater-than-average DNA-level correspondence to one another. Telomere-centric CNR following recursive WGD in plants is also important in mammals and yeast, and may be a general mechanism of restoring small linear chromosome numbers in higher eukaryotes.

Evolution of DNA methylation patterns in the Brassicaceae is driven by differences in genome organization

DNA methylation is an ancient molecular modification found in most eukaryotes. In plants, DNA methylation is not only critical for transcriptionally silencing transposons, but can also affect phenotype by altering expression of protein coding genes. The extent of its contribution to phenotypic diversity over evolutionary time is, however, unclear, because of limited stability of epialleles that are not linked to DNA mutations. To dissect the relative contribution of DNA methylation to transposon surveillance and host gene regulation, we leveraged information from three species in the Brassicaceae that vary in genome architecture, Capsella rubella, Arabidopsis lyrata, and Arabidopsis thaliana. We found that the lineage-specific expansion and contraction of transposon and repeat sequences is the main driver of interspecific differences in DNA methylation. The most heavily methylated portions of the genome are thus not conserved at the sequence level. Outside of repeat-associated methylation, there is a surprising degree of conservation in methylation at single nucleotides located in gene bodies. Finally, dynamic DNA methylation is affected more by tissue type than by environmental differences in all species, but these responses are not conserved. The majority of DNA methylation variation between species resides in hypervariable genomic regions, and thus, in the context of macroevolution, is of limited phenotypic consequence.

DNA methylation is an epigenetic mark that has received a great deal of attention in plants because it can be stably transmitted across generations. However, the rate of DNA methylation change, or epimutation, is greater than that of DNA mutation. In addition, different from DNA sequence, DNA methylation can vary within an individual in response to developmental or environmental cues. Whether altered characters can be passed on to the next generation via directed modifications in DNA methylation is a question of great interest. We have compared how DNA methylation changes between species, tissues, and environments using three closely related crucifers as examples. We found that DNA methylation is different between roots and shoots and changes with temperatures, but that such changes are not conserved across species. Moreover, most of the methylated sites are not conserved between species. This suggests that DNA methylation may respond to immediate fluctuations in the environment, but this response is not retained over long evolutionary periods. Thus, in contrast to transcriptional responses, conserved epigenetic responses at the level of DNA methylation are not widespread. Instead, the patterns of DNA methylation are largely determined by the evolution of genome structure, and responsive loci are likely short-lived accidents of this process.

Disentangling relationships among the diploid members of the intricate genus Knautia (Caprifoliaceae, Dipsacoideae)

The genus Knautia (Caprifoliaceae, Dipsacoideae) encompasses 40-60 species mainly distributed in western Eurasia, with highest species diversity in the Alps and the Balkan Peninsula. It is traditionally regarded as one of the taxonomically most challenging European genera due to the widespread occurrence of polyploidy, the high incidence of hybridisation and the maintenance of morphologically intermediate forms. A prerequisite for assessing the complex spatiotemporal diversification of a polyploid group is a comprehensive hypothesis of the phylogenetic relationships among its diploid members. To this end, DNA sequence data (nrDNA ITS and plastid petN(ycf6)-psbM) combined with AFLP fingerprinting were performed on 148 diploid populations belonging to 35 taxa. Phylogenies obtained by maximum parsimony and Bayesian analyses were used to test the monophyly of the genus and its three sections Trichera, Tricheroides and Knautia, to provide insights into its evolutionary history and to test previous hypotheses of inter- and intrasectional classification. Both nuclear and chloroplast datasets support the monophyly of Knautia and its three sections, with ambiguous placement of K. cf. degenii. The majority of species belong to the nearly exclusively perennial section Trichera (x=10). Within section Trichera all markers revealed largely unresolved phylogenetic relationships suggesting rapid radiation and recent range expansion. In addition, extensive sharing of plastid haplotypes across taxa and wide geographic ranges of plastid haplotypes and ribotype groups were observed. The molecular data are partly at odds with the traditional informal grouping of taxa within section Trichera. Whereas the traditional groups of K. dinarica, K. drymeia and K. montana can be maintained, the new, smaller and well supported Midzorensis and Pancicii Groups as well as the SW European Group are separated from the heterogeneous traditional K. longifolia group. The former groups of K. arvensis, K. dalmatica, K. fleischmannii and K. velutina are clearly polyphyletic. Their diploid members have to be rearranged into the Xerophytic Group, the Carinthiaca Group, and the Northern and Southern Arvensis Groups. The annual sections Tricheroides (x=10) and Knautia (x=8) with only a few taxa are resolved in the ITS and plastid trees on long branches as early diverging lineages within the genus.

Two young MicroRNAs originating from target duplication mediate nitrogen starvation adaptation via regulation of glucosinolate synthesis in Arabidopsis thaliana

Nitrogen is an essential macronutrient required for plant growth and development. A number of genes respond to nitrogen starvation conditions. However, the functions of most of these nitrogen starvation-responsive genes are unclear. Our recent survey suggested that many microRNAs (miRNAs) are responsive to nitrogen starvation in Arabidopsis thaliana. Here, we identified a new miRNA (miR5090) from the complementary transcript of the MIR826 gene. Further investigation uncovered that both miRNA genes recently evolved from the inverse duplication of their common target gene, ALKENYL HYDROXALKYL PRODUCING2 (AOP2). Similar to miR826, miR5090 is induced by nitrogen starvation. By contrast, the AOP2 transcript level was negatively correlated with miR826 and miR5090 under nitrogen starvation. GUS-fused AOP2 expression suggested that AOP2 was posttranscriptionally suppressed by miR826 and miR5090. miRNA transgenic plants with significantly low AOP2 expression accumulated fewer Met-derived glucosinolates, phenocopying the aop2 mutants. Most glucosinolate synthesis-associated genes were repressed under nitrogen starvation conditions. Furthermore, miRNA transgenic plants with less glucosinolate displayed enhanced tolerance to nitrogen starvation, including high biomass, more lateral roots, increased chlorophyll, and decreased anthocyanin. Meanwhile, nitrogen starvation-responsive genes were up-regulated in transgenic plants, implying improved nitrogen uptake activity. Our study reveals a mechanism by which Arabidopsis thaliana regulates the synthesis of glucosinolates to adapt to environmental changes in nitrogen availability.

Distinct evolutionary strategies in the GGPPS family from plants

Multiple geranylgeranyl diphosphate synthases (GGPPS) for biosynthesis of geranylgeranyl diphosphate (GGPP) exist in plants. GGPP is produced in the isoprenoid pathway and is a central precursor for various primary and specialized plant metabolites. Therefore, its biosynthesis is an essential regulatory point in the isoprenoid pathway. We selected 119 GGPPSs from 48 species representing all major plant lineages, based on stringent homology criteria. After the diversification of land plants, the number of GGPPS paralogs per species increases. Already in the moss Physcomitrella patens, GGPPS appears to be encoded by multiple paralogous genes. In gymnosperms, neofunctionalization of GGPPS may have enabled optimized biosynthesis of primary and specialized metabolites. Notably, lineage-specific expansion of GGPPS occurred in land plants. As a representative species we focused here on Arabidopsis thaliana, which retained the highest number of GGPPS paralogs (twelve) among the 48 species we considered in this study. Our results show that the A. thaliana GGPPS gene family is an example of evolution involving neo- and subfunctionalization as well as pseudogenization. We propose subfunctionalization as one of the main mechanisms allowing the maintenance of multiple GGPPS paralogs in A. thaliana genome. Accordingly, the changes in the expression patterns of the GGPPS paralogs occurring after gene duplication led to developmental and/or condition specific functional evolution.

Next-generation sequencing, FISH mapping and synteny-based modeling reveal mechanisms of decreasing dysploidy in Cucumis

In the large Cucurbitaceae genus Cucumis, cucumber (C. sativus) is the only species with 2n = 2x = 14 chromosomes. The majority of the remaining species, including melon (C. melo) and the sister species of cucumber, C. hystrix, have 2n = 2x = 24 chromosomes, implying a reduction from n = 12 to n = 7. To understand the underlying mechanisms, we investigated chromosome synteny among cucumber, C. hystrix and melon using integrated and complementary approaches. We identified 14 inversions and a C. hystrix lineage-specific reciprocal inversion between C. hystrix and melon. The results reveal the location and orientation of 53 C. hystrix syntenic blocks on the seven cucumber chromosomes, and allow us to infer at least 59 chromosome rearrangement events that led to the seven cucumber chromosomes, including five fusions, four translocations, and 50 inversions. The 12 inferred chromosomes (AK1-AK12) of an ancestor similar to melon and C. hystrix had strikingly different evolutionary fates, with cucumber chromosome C1 apparently resulting from insertion of chromosome AK12 into the centromeric region of translocated AK2/AK8, cucumber chromosome C3 originating from a Robertsonian-like translocation between AK4 and AK6, and cucumber chromosome C5 originating from fusion of AK9 and AK10. Chromosomes C2, C4 and C6 were the result of complex reshuffling of syntenic blocks from three (AK3, AK5 and AK11), three (AK5, AK7 and AK8) and five (AK2, AK3, AK5, AK8 and AK11) ancestral chromosomes, respectively, through 33 fusion, translocation and inversion events. Previous results (Huang, S., Li, R., Zhang, Z. et al., , Nat. Genet. 41, 1275-1281; Li, D., Cuevas, H.E., Yang, L., Li, Y., Garcia-Mas, J., Zalapa, J., Staub, J.E., Luan, F., Reddy, U., He, X., Gong, Z., Weng, Y. 2011a, BMC Genomics, 12, 396) showing that cucumber C7 stayed largely intact during the entire evolution of Cucumis are supported. Results from this study allow a fine-scale understanding of the mechanisms of dysploid chromosome reduction that has not been achieved previously.

Exploiting a Reference Genome in Terms of Duplications: The Network of Paralogs and Single Copy Genes in Arabidopsis thaliana

Arabidopsis thaliana became the model organism for plant studies because of its small diploid genome, rapid lifecycle and short adult size. Its genome was the first among plants to be sequenced, becoming the reference in plant genomics. However, the Arabidopsis genome is characterized by an inherently complex organization, since it has undergone ancient whole genome duplications, followed by gene reduction, diploidization events and extended rearrangements, which relocated and split up the retained portions. These events, together with probable chromosome reductions, dramatically increased the genome complexity, limiting its role as a reference. The identification of paralogs and single copy genes within a highly duplicated genome is a prerequisite to understand its organization and evolution and to improve its exploitation in comparative genomics. This is still controversial, even in the widely studied Arabidopsis genome. This is also due to the lack of a reference bioinformatics pipeline that could exhaustively identify paralogs and singleton genes. We describe here a complete computational strategy to detect both duplicated and single copy genes in a genome, discussing all the methodological issues that may strongly affect the results, their quality and their reliability. This approach was used to analyze the organization of Arabidopsis nuclear protein coding genes, and besides classifying computationally defined paralogs into networks and single copy genes into different classes, it unraveled further intriguing aspects concerning the genome annotation and the gene relationships in this reference plant species. Since our results may be useful for comparative genomics and genome functional analyses, we organized a dedicated web interface to make them accessible to the scientific community

Development of a next-generation NIL library in Arabidopsis thaliana for dissecting complex traits

Background

The identification of the loci and specific alleles underlying variation in quantitative traits is an important goal for evolutionary biologists and breeders. Despite major advancements in genomics technology, moving from QTL to causal alleles remains a major challenge in genetics research. Near-isogenic lines are the ideal raw material for QTL validation, refinement of QTL location and, ultimately, gene discovery.

Results

In this study, a population of 75 Arabidopsis thaliana near-isogenic lines was developed from an existing recombinant inbred line (RIL) population derived from a cross between physiologically divergent accessions Kas-1 and Tsu-1. First, a novel algorithm was developed to utilize genome-wide marker data in selecting RILs fully isogenic to Kas-1 for a single chromosome. Seven such RILs were used in 2 generations of crossing to Tsu-1 to create BC1 seed. BC1 plants were genotyped with SSR markers so that lines could be selected that carried Kas-1 introgressions, resulting in a population carrying chromosomal introgressions spanning the genome. BC1 lines were genotyped with 48 genome-wide SSRs to identify lines with a targeted Kas-1 introgression and the fewest genomic introgressions elsewhere. 75 such lines were selected and genotyped at an additional 41 SNP loci and another 930 tags using 2b-RAD genotyping by sequencing. The final population carried an average of 1.35 homozygous and 2.49 heterozygous introgressions per line with average introgression sizes of 5.32 and 5.16 Mb, respectively. In a simple case study, we demonstrate the advantage of maintaining heterozygotes in our library whereby fine-mapping efforts are conducted simply by self-pollination. Crossovers in the heterozygous interval during this single selfing generation break the introgression into smaller, homozygous fragments (sub-NILs). Additionally, we utilize a homozygous NIL for validation of a QTL underlying stomatal conductance, a low heritability trait.

Conclusions

The present results introduce a new and valuable resource to the Brassicaceae research community that enables rapid fine-mapping of candidate loci in parallel with QTL validation. These attributes along with dense marker coverage and genome-wide chromosomal introgressions make this population an ideal starting point for discovery of genes underlying important complex traits of agricultural and ecological significance.

Identification, expression, and comparative genomic analysis of the IPT and CKX gene families in Chinese cabbage (Brassica rapa ssp. pekinensis)

BACKGROUND

Cytokinins (CKs) have significant roles in various aspects of plant growth and development, and they are also involved in plant stress adaptations. The fine-tuning of the controlled CK levels in individual tissues, cells, and organelles is properly maintained by isopentenyl transferases (IPTs) and cytokinin oxidase/dehydrogenases (CKXs). Chinese cabbage is one of the most economically important vegetable crops worldwide. The whole genome sequencing of Brassica rapa enables us to perform the genome-wide identification and functional analysis of the IPT and CKX gene families.

RESULTS

In this study, a total of 13 BrIPT genes and 12 BrCKX genes were identified. The gene structures, conserved domains and phylogenetic relationships were analyzed. The isoelectric point, subcellular localization and glycosylation sites of the proteins were predicted. Segmental duplicates were found in both BrIPT and BrCKX gene families. We also analyzed evolutionary patterns and divergence of the IPT and CKX genes in the Cruciferae family. The transcription levels of BrIPT and BrCKX genes were analyzed to obtain an initial picture of the functions of these genes. Abiotic stress elements related to adverse environmental stimuli were found in the promoter regions of BrIPT and BrCKX genes and they were confirmed to respond to drought and high salinity conditions. The effects of 6-BA and ABA on the expressions of BrIPT and BrCKX genes were also investigated.

CONCLUSIONS

The expansion of BrIPT and BrCKX genes after speciation from Arabidopsis thaliana is mainly attributed to segmental duplication events during the whole genome triplication (WGT) and substantial duplicated genes are lost during the long evolutionary history. Genes produced by segmental duplication events have changed their expression patterns or may adopted new functions and thus are obtained. BrIPT and BrCKX genes respond well to drought and high salinity stresses, and their transcripts are affected by exogenous hormones, such as 6-BA and ABA, suggesting their potential roles in abiotic stress conditions and regulatory mechanisms of plant hormone homeostasis. The appropriate modulation of endogenous CKs levels by IPT and CKX genes is a promising approach for developing economically important high-yielding and high-quality stress-tolerant crops in agriculture.

An Ultra-High-Density, Transcript-Based, Genetic Map of Lettuce

We have generated an ultra-high-density genetic map for lettuce, an economically important member of the Compositae, consisting of 12,842 unigenes (13,943 markers) mapped in 3696 genetic bins distributed over nine chromosomal linkage groups. Genomic DNA was hybridized to a custom Affymetrix oligonucleotide array containing 6.4 million features representing 35,628 unigenes of Lactuca spp. Segregation of single-position polymorphisms was analyzed using 213 F7:8 recombinant inbred lines that had been generated by crossing cultivated Lactuca sativa cv. Salinas and L. serriola acc. US96UC23, the wild progenitor species of L. sativa The high level of replication of each allele in the recombinant inbred lines was exploited to identify single-position polymorphisms that were assigned to parental haplotypes. Marker information has been made available using GBrowse to facilitate access to the map. This map has been anchored to the previously published integrated map of lettuce providing candidate genes for multiple phenotypes. The high density of markers achieved in this ultradense map allowed syntenic studies between lettuce and Vitis vinifera as well as other plant species.

Gene duplication within the Green Lineage: the case of TEL genes

Recent years have witnessed a breathtaking increase in the availability of genome sequence data, providing evidence of the highly duplicate nature of eukaryotic genomes. Plants are exceptional among eukaryotic organisms in that duplicate loci compose a large fraction of their genomes, partly because of the frequent occurrence of polyploidy (or whole-genome duplication) events. Tandem gene duplication and transposition have also contributed to the large number of duplicated genes in plant genomes. Evolutionary analyses allowed the dynamics of duplicate gene evolution to be studied and several models were proposed. It seems that, over time, many duplicated genes were lost and some of those that were retained gained new functions and/or expression patterns (neofunctionalization) or subdivided their functions and/or expression patterns between them (subfunctionalization). Recent studies have provided examples of genes that originated by duplication with successive diversification within plants. In this review, we focused on the TEL (TERMINAL EAR1-like) genes to illustrate such mechanisms. Emerged from the mei2 gene family, these TEL genes are likely to be land plant-specific. Phylogenetic analyses revealed one or two TEL copies per diploid genome. TEL gene degeneration and loss in several Angiosperm species such as in poplar and maize seem to have occurred. In Arabidopsis thaliana, whose genome experienced at least three polyploidy events followed by massive gene loss and genomic reorganization, two TEL genes were retained and two new shorter TEL-like (MCT) genes emerged. Molecular and expression analyses suggest for these genes sub- and neofunctionalization events, but confirmation will come from their functional characterization.

Signals of speciation within Arabidopsis thaliana in comparison with its relatives

The species within the now well-defined Arabidopsis genus provide biological materials suitable to investigate speciation and the development of reproductive isolation barriers between related species. Even within the model species A. thaliana, genetic differentiation between populations due to environmental adaptation or demographic history can lead to cases where hybrids between accessions are non-viable. Experimental evidence supports the importance of genome duplications and genetic epistatic interactions in the occurrence of reproductive isolation. Other examples of adaptation to specific environments can be found in Arabidopsis relatives where hybridization and chromosome doubling lead to new amphidiploid species. Molecular signals of speciation found in the Arabidopsis genus should provide a better understanding of speciation processes in plants from a genetic, molecular and evolutionary perspective.

Functional alleles of the flowering time regulator FRIGIDA in the Brassica oleracea genome

BACKGROUND

Plants adopt different reproductive strategies as an adaptation to growth in a range of climates. In Arabidopsis thaliana FRIGIDA (FRI) confers a vernalization requirement and thus winter annual habit by increasing the expression of the MADS box transcriptional repressor FLOWERING LOCUS C (FLC). Variation at FRI plays a major role in A. thaliana life history strategy, as independent loss-of-function alleles that result in a rapid-cycling habit in different accessions, appear to have evolved many times. The aim of this study was to identify and characterize orthologues of FRI in Brassica oleracea.

RESULTS

We describe the characterization of FRI from Brassica oleracea and identify the two B. oleracea FRI orthologues (BolC.FRI.a and BolC.FRI.b). These show extensive amino acid conservation in the central and C-terminal regions to FRI from other Brassicaceae, including A. thaliana, but have a diverged N-terminus. The genes map to two of the three regions of B. oleracea chromosomes syntenic to part of A. thaliana chromosome 5 suggesting that one of the FRI copies has been lost since the ancient triplication event that formed the B. oleracea genome. This genomic position is not syntenic with FRI in A. thaliana and comparative analysis revealed a recombination event within the A. thaliana FRI promoter. This relocated A. thaliana FRI to chromosome 4, very close to the nucleolar organizer region, leaving a fragment of FRI in the syntenic location on A. thaliana chromosome 5. Our data show this rearrangement occurred after the divergence from A. lyrata. We explored the allelic variation at BolC.FRI.a within cultivated B. oleracea germplasm and identified two major alleles, which appear equally functional both to each other and A. thaliana FRI, when expressed as fusions in A. thaliana.

CONCLUSIONS

We identify the two Brassica oleracea FRI genes, one of which we show through A. thaliana complementation experiments is functional, and show their genomic location is not syntenic with A. thaliana FRI due to an ancient recombination event. This has complicated previous association analyses of FRI with variation in life history strategy in the Brassica genus.

Evolution and Structural Diversification of PILS Putative Auxin Carriers in Plants

The phytohormone auxin contributes to virtually every aspect of the plant development. The spatiotemporal distribution of auxin depends on a complex interplay between auxin metabolism and intercellular auxin transport. Intracellular auxin compartmentalization provides another link between auxin transport processes and auxin metabolism. The PIN-LIKES (PILS) putative auxin carriers localize to the endoplasmic reticulum (ER) and contribute to cellular auxin homeostasis. PILS proteins regulate intracellular auxin accumulation, the rate of auxin conjugation and, subsequently, affect nuclear auxin signaling. Here, we investigate sequence diversification of the PILS family in Arabidopsis thaliana and provide insights into the evolution of these novel putative auxin carriers in plants. Our data suggest that PILS proteins are conserved throughout the plant lineage and expanded during higher plant evolution. PILS proteins diversified early during plant evolution into three clades. Besides the ancient Clade I encompassing non-land plant species, PILS proteins evolved into two clades. The diversification of Clade II and Clade III occurred already at the level of non-vascular plant evolution and, hence, both clades contain vascular and non-vascular plant species. Nevertheless, Clade III contains fewer non- and increased numbers of vascular plants, indicating higher importance of Clade III for vascular plant evolution. Notably, PILS proteins are distinct and appear evolutionarily older than the prominent PIN-FORMED auxin carriers. Moreover, we revealed particular PILS sequence divergence in Arabidopsis and assume that these alterations could contribute to distinct gene regulations and protein functions.

Genome size in Anthurium evaluated in the context of karyotypes and phenotypes

BACKGROUND AND AIMS

Anthurium is an important horticultural crop from the family Araceae, order Alismatales, a lineage considered to have diverged from other monocots prior to the cereals. Genome size and its distribution in Anthurium were investigated to gain a basic understanding of genome organization in this large genus and to forge a firm foundation for advancement of molecular approaches for the study of Anthurium. Currently, genome size estimates have been reported for only two Anthurium samples.

METHODOLOGY

Bulk nuclear DNA content estimates were obtained by flow cell cytometry using leaf tissue collected from Anthurium species of different subgeneric groups and from commercial cultivars. The most current and well-supported topology of subgeneric, sectional relationships was applied to present genome size estimates in the context of reported chromosome counts, karyotypes, putative phylogenetic relationships, observed phenotypes and pedigree.

PRINCIPAL RESULTS

Genome size estimates based on bulk nuclear DNA content for 77 accessions representing 34 species and 9 cultivars were obtained, including initial estimates for 33 Anthurium species, and both the smallest (Anthurium obtusum; Tetraspermium) and largest (Anthurium roseospadix; Calomystrium) Anthurium genome sizes reported to date. Genome size did not distinguish any subgeneric section, but ranged 5-fold (4.42-20.83 pg/2 C) despite consistent 2N= 30 chromosome counts. Intraspecies genome size variation >20 % is reported for Anthurium ravenii, A. watermaliense and A. gracile.

CONCLUSIONS

Genome size estimates for Anthurium species spanning 13 recognized subgeneric sections indicate that genome size does not generally correlate with chromosome count or phylogenetic relationships. Mechanisms of genome expansion and contraction, including amplification and reduction of repetitive elements, polyploidy, chromosome reorganization/loss, may be involved in genome evolution in Anthurium as in other species. The new information on Anthurium genome sizes provides a platform for molecular studies supporting further research on genome evolution as well as cultivar development.

Proteomic divergence in Arabidopsis autopolyploids and allopolyploids and their progenitors

Autopolyploidy and allopolyploidy are common in many plants and some animals. Rapid changes in genomic composition and gene expression have been observed in both autopolyploids and allopolyploids, but the effects of polyploidy on proteomic divergence are poorly understood. Here, we report quantitative analysis of protein changes in leaves of Arabidopsis autopolyploids and allotetraploids and their progenitors using isobaric tags for relative and absolute quantitation (iTRAQ) coupled with mass spectrometry. In more than 1000 proteins analyzed, the levels of protein divergence were relatively high (~18%) between Arabidopsis thaliana and Arabidopsis arenosa, relatively low (~6.8%) between an A. thaliana diploid and autotetraploid and intermediate (~8.3 and 8.2%) in F(1)- and F(8)-resynthesized allotetraploids relative to mid-parent values, respectively. This pattern of proteomic divergence was consistent with the previously reported gene expression data. In particular, many non-additively accumulated proteins (61-62%) in the F(1) and F(8) allotetraploids were also differentially expressed between the parents. The differentially accumulated proteins in functional categories of abiotic and biotic stresses were overrepresented between an A. thaliana autotetraploid and diploid and between two Arabidopsis species, but not significantly different between allotetraploids and their progenitors. Although the trend of changes is similar, the percentage of differentially accumulated proteins that matched previously reported differentially expressed genes was relatively low. Western blot analysis confirmed several selected proteins with isoforms the cumulative levels of which were differentially expressed. These data suggest high protein divergence between species and rapid changes in post-transcriptional regulation and translational modifications of proteins during polyploidization.

The genome of the mesopolyploid crop species Brassica rapa

We report the annotation and analysis of the draft genome sequence of Brassica rapa accession Chiifu-401-42, a Chinese cabbage. We modeled 41,174 protein coding genes in the B. rapa genome, which has undergone genome triplication. We used Arabidopsis thaliana as an outgroup for investigating the consequences of genome triplication, such as structural and functional evolution. The extent of gene loss (fractionation) among triplicated genome segments varies, with one of the three copies consistently retaining a disproportionately large fraction of the genes expected to have been present in its ancestor. Variation in the number of members of gene families present in the genome may contribute to the remarkable morphological plasticity of Brassica species. The B. rapa genome sequence provides an important resource for studying the evolution of polyploid genomes and underpins the genetic improvement of Brassica oil and vegetable crops.

Extensive chromosome homoeology among Brassiceae species were revealed by comparative genetic mapping with high-density EST-based SNP markers in radish (Raphanus sativus L.)

A linkage map of expressed sequence tag (EST)-based markers in radish (Raphanus sativus L.) was constructed using a low-cost and high-efficiency single-nucleotide polymorphism (SNP) genotyping method named multiplex polymerase chain reaction–mixed probe dot-blot analysis developed in this study. Seven hundred and forty-six SNP markers derived from EST sequences of R. sativus were assigned to nine linkage groups with a total length of 806.7 cM. By BLASTN, 726 markers were found to have homologous genes in Arabidopsis thaliana, and 72 syntenic regions, which have great potential for utilizing genomic information of the model species A. thaliana in basic and applied genetics of R. sativus, were identified. By construction and analysis of the genome structures of R. sativus based on the 24 genomic blocks within the Brassicaceae ancestral karyotype, 23 of the 24 genomic blocks were detected in the genome of R. sativus, and half of them were found to be triplicated. Comparison of the genome structure of R. sativus with those of the A, B, and C genomes of Brassica species and that of Sinapis alba L. revealed extensive chromosome homoeology among Brassiceae species, which would facilitate transfer of the genomic information from one Brassiceae species to another.

The genome sequence of the North-European cucumber (Cucumis sativus L.) unravels evolutionary adaptation mechanisms in plants

Cucumber (Cucumis sativus L.), a widely cultivated crop, has originated from Eastern Himalayas and secondary domestication regions includes highly divergent climate conditions e.g. temperate and subtropical. We wanted to uncover adaptive genome differences between the cucumber cultivars and what sort of evolutionary molecular mechanisms regulate genetic adaptation of plants to different ecosystems and organism biodiversity. Here we present the draft genome sequence of the Cucumis sativus genome of the North-European Borszczagowski cultivar (line B10) and comparative genomics studies with the known genomes of: C. sativus (Chinese cultivar – Chinese Long (line 9930)), Arabidopsis thaliana, Populus trichocarpa and Oryza sativa. Cucumber genomes show extensive chromosomal rearrangements, distinct differences in quantity of the particular genes (e.g. involved in photosynthesis, respiration, sugar metabolism, chlorophyll degradation, regulation of gene expression, photooxidative stress tolerance, higher non-optimal temperatures tolerance and ammonium ion assimilation) as well as in distributions of abscisic acid-, dehydration- and ethylene-responsive cis-regulatory elements (CREs) in promoters of orthologous group of genes, which lead to the specific adaptation features. Abscisic acid treatment of non-acclimated Arabidopsis and C. sativus seedlings induced moderate freezing tolerance in Arabidopsis but not in C. sativus. This experiment together with analysis of abscisic acid-specific CRE distributions give a clue why C. sativus is much more susceptible to moderate freezing stresses than A. thaliana. Comparative analysis of all the five genomes showed that, each species and/or cultivars has a specific profile of CRE content in promoters of orthologous genes. Our results constitute the substantial and original resource for the basic and applied research on environmental adaptations of plants, which could facilitate creation of new crops with improved growth and yield in divergent conditions.

Metal hyperaccumulation and hypertolerance: a model for plant evolutionary genomics

In the course of evolution, plants adapted to widely differing metal availabilities in soils and therefore represent an important source of natural variation of metal homeostasis networks. Research on plant metal homeostasis can thus provide insights into the functioning, regulation and adaptation of biological networks. Here, we describe major recent breakthroughs in the understanding of the genetic and molecular basis of metal hyperaccumulation and associated hypertolerance, a naturally selected complex trait which represents an extreme adaptation of the metal homeostasis network. Investigations in this field reveal further the molecular alterations underlying the evolution of natural phenotypic diversity and provide a highly relevant framework for comparative genomics.

Screening synteny blocks in pairwise genome comparisons through integer programming

Background

It is difficult to accurately interpret chromosomal correspondences such as true orthology and paralogy due to significant divergence of genomes from a common ancestor. Analyses are particularly problematic among lineages that have repeatedly experienced whole genome duplication (WGD) events. To compare multiple "subgenomes" derived from genome duplications, we need to relax the traditional requirements of "one-to-one" syntenic matchings of genomic regions in order to reflect "one-to-many" or more generally "many-to-many" matchings. However this relaxation may result in the identification of synteny blocks that are derived from ancient shared WGDs that are not of interest. For many downstream analyses, we need to eliminate weak, low scoring alignments from pairwise genome comparisons. Our goal is to objectively select subset of synteny blocks whose total scores are maximized while respecting the duplication history of the genomes in comparison. We call this "quota-based" screening of synteny blocks in order to appropriately fill a quota of syntenic relationships within one genome or between two genomes having WGD events.

Results

We have formulated the synteny block screening as an optimization problem known as "Binary Integer Programming" (BIP), which is solved using existing linear programming solvers. The computer program QUOTA-ALIGN performs this task by creating a clear objective function that maximizes the compatible set of synteny blocks under given constraints on overlaps and depths (corresponding to the duplication history in respective genomes). Such a procedure is useful for any pairwise synteny alignments, but is most useful in lineages affected by multiple WGDs, like plants or fish lineages. For example, there should be a 1:2 ploidy relationship between genome A and B if genome B had an independent WGD subsequent to the divergence of the two genomes. We show through simulations and real examples using plant genomes in the rosid superorder that the quota-based screening can eliminate ambiguous synteny blocks and focus on specific genomic evolutionary events, like the divergence of lineages (in cross-species comparisons) and the most recent WGD (in self comparisons).

Conclusions

The QUOTA-ALIGN algorithm screens a set of synteny blocks to retain only those compatible with a user specified ploidy relationship between two genomes. These blocks, in turn, may be used for additional downstream analyses such as identifying true orthologous regions in interspecific comparisons. There are two major contributions of QUOTA-ALIGN: 1) reducing the block screening task to a BIP problem, which is novel; 2) providing an efficient software pipeline starting from all-against-all BLAST to the screened synteny blocks with dot plot visualizations. Python codes and full documentations are publicly available http://github.com/tanghaibao/quota-alignment. QUOTA-ALIGN program is also integrated as a major component in SynMap http://genomevolution.com/CoGe/SynMap.pl, offering easier access to thousands of genomes for non-programmers.

The Arabidopsis lyrata genome sequence and the basis of rapid genome size change

We report the 207-Mb genome sequence of the North American Arabidopsis lyrata strain MN47 based on 8.3× dideoxy sequence coverage. We predict 32,670 genes in this outcrossing species compared to the 27,025 genes in the selfing species Arabidopsis thaliana. The much smaller 125-Mb genome of A. thaliana, which diverged from A. lyrata 10 million years ago, likely constitutes the derived state for the family. We found evidence for DNA loss from large-scale rearrangements, but most of the difference in genome size can be attributed to hundreds of thousands of small deletions, mostly in noncoding DNA and transposons. Analysis of deletions and insertions still segregating in A. thaliana indicates that the process of DNA loss is ongoing, suggesting pervasive selection for a smaller genome. The high-quality reference genome sequence for A. lyrata will be an important resource for functional, evolutionary and ecological studies in the genus Arabidopsis.