Molecular phylogeny and divergence time of Harpalyce (Leguminosae, Papilionoideae), a lineage with amphitropical diversification in seasonally dry forests and savannas

Our knowledge of the systematics of the papilionoid legume tribe Brongniartieae has greatly benefitted from recent advances in molecular phylogenetics. The tribe was initially described to include species marked by a strongly bilabiate calyx and an embryo with a straight radicle, but recent research has placed taxa from the distantly related core Sophoreae and Millettieae within it. Despite these advances, the most species-rich genera within the Brongniartieae are still not well studied, and their morphological and biogeographical evolution remains poorly understood. Comprising 35 species, Harpalyce is one of these poorly studied genera. In this study, we present a comprehensive, multi-locus molecular phylogeny of the Brongniartieae, with an increased sampling of Harpalyce, to investigate morphological and biogeographical evolution within the group. Our results confirm the monophyly of Harpalyce and indicate that peltate glandular trichomes and a strongly bilabiate calyx with a carinal lip and three fused lobes are synapomorphies for the genus, which is internally divided into three distinct ecologically and geographically divergent lineages, corresponding to the previously recognized sections. Our biogeographical reconstructions demonstrate that Brongniartieae originated in South America during the Eocene, with subsequent pulses of diversification in South America, Mesoamerica, and Australia. Harpalyce also originated in South America during the Miocene at around 20 Ma, with almost synchronous later diversification in South America and Mexico/Mesoamerica beginning 10 Ma, but mostly during the Pliocene. Migration of Harpalyce from South to North America was accompanied by a biome and ecological shift from savanna to seasonally dry forest.

The genetic structure and demographic history of Zabelia tyaihyonii, endemic to Korean limestone karst forests, based on genome-wide SNP markers

Similar to the global phenomenon, many plant species endemic to Korean limestone karst forests are at risk of extinction due to human intervention. Zabelia tyaihyonii is a familiar shrub, called "Hardy abelia" and "Fragrant abelia" growing in the karst forests of Korea, where it is one of the most threatened species. We investigated the genetic structure and demographic history of Z. tyaihyonii, which allow us to develop appropriate conservation and management strategies. The genetic structure was evaluated using a total of 187 samples from 14 populations, covering the entire distribution of Z. tyaihyonii in South Korea. We utilized 254 and 1753 SNP loci obtained via MIG-seq (Multiplexed ISSR Genotyping by sequencing) for structure and demographic analyses, respectively. The population demographic modeling was performed with site frequency spectrum. To gain further historical insights, we also employed ENM (Ecological Niche Modeling). We found two distinct clusters (CLI and CLII) of ancient origin (ca. 490 ka). Despite CLII experiencing a more severe bottleneck, both clusters showed similar levels of genetic diversity, indicating mutual historical gene flow. Their historical distribution range seems to have changed very little. We proposed a historical distribution scenario for Z. tyaihyonii, taking into account its intrinsic factors, and emphasized a more complex response to Quaternary climate change beyond simple allopatric speciation models. These findings provide valuable insights for conservation and management strategies for Z. tyaihyonii.

Characterization and comparative analysis of the complete plastid genomes of four Astragalus species

Astragalus is the largest flowering plant genus. We assembled the plastid genomes of four Astragalus species (Astragalus iranicus, A. macropelmatus, A. mesoleios, A. odoratus) using next-generation sequencing and analyzed their plastomes including genome organization, codon usage, nucleotide diversity, prediction of RNA editing and etc. The total length of the newly sequenced Astragalus plastomes ranged from 121,050 bp to 123,622 bp, with 110 genes comprising 76 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosome RNA (rRNA) genes. Comparative analysis of the chloroplast genomes of Astragalus revealed several hypervariable regions comprising three non-coding sites (trnQ(UUG)-accD, rps7 -trnV(GAC) and trnR(ACG)-trnN(GUU)) and four protein-coding genes (ycf1, ycf2, accD and clpP), which have potential as molecular markers. Positive selection signatures were found in five genes in Astragalus species including rps11, rps15, accD, clpP and ycf1. The newly sequenced species, A. macropelmatus, has an approximately 13-kb inversion in IR region. Phylogenetic analysis based on 75 protein-coding gene sequences confirmed that Astragalus form a monophyletic clade within the tribe Galegeae and Oxytropis is sister group to the Coluteoid clade. The results of this study may helpful in elucidating the chloroplast genome structure, understanding the evolutionary dynamics at genus Astragalus and IRLC levels and investigating the phylogenetic relationships. Moreover, the newly plastid genomes sequenced have been increased the plastome data resources on Astragalus that can be useful in further phylogenomic studies.

Plastid phylogenomics uncovers multiple species in Medicago truncatula (Fabaceae) germplasm accessions

Medicago truncatula is a model legume that has been extensively investigated in diverse subdisciplines of plant science. Medicago littoralis can interbreed with M. truncatula and M. italica; these three closely related species form a clade, i.e. TLI clade. Genetic studies have indicated that M. truncatula accessions are heterogeneous but their taxonomic identities have not been verified. To elucidate the phylogenetic position of diverse M. truncatula accessions within the genus, we assembled 54 plastid genomes (plastomes) using publicly available next-generation sequencing data and conducted phylogenetic analyses using maximum likelihood. Five accessions showed high levels of plastid DNA polymorphism. Three of these highly polymorphic accessions contained sequences from both M. truncatula and M. littoralis. Phylogenetic analyses of sequences placed some accessions closer to distantly related species suggesting misidentification of source material. Most accessions were placed within the TLI clade and maximally supported the interrelationships of three subclades. Two Medicago accessions were placed within a M. italica subclade of the TLI clade. Plastomes with a 45-kb (rpl20-ycf1) inversion were placed within the M. littoralis subclade. Our results suggest that the M. truncatula accession genome pool represents more than one species due to possible mistaken identities and gene flow among closely related species.

Remarkably rapid, recent diversification of Cochemiea and Mammillaria in the Baja California, Mexico region

PREMISE

The Cactaceae of northwestern Mexico and the southwestern United States constitute a major component of the angiosperm biodiversity of the region. The Mammilloid clade, (Cactaceae, tribe Cacteae), composed of the genera Cochemiea, Coryphantha, Cumarinia, Mammillaria, and Pelecyphora is especially species rich. We sought to understand the timing, geographical and climate influences correlated with expansion of the Mammilloid clade, through the Sonoran Desert into Baja California.

METHODS

We reconstructed the historical biogeography of the Mammilloid clade, using Bayesian and maximum likelihood methods, based on a strongly supported molecular phylogeny. We also estimated divergence times, the timing of emergence of key characters, and diversification rates and rate shifts of the Mammilloid clade.

RESULTS

We found that the most recent common ancestor of Cochemiea arrived in the Cape region of Baja California from the Sonoran Desert region approximately 5 million years ago, coinciding with the timing of peninsular rifting from the mainland, suggesting dispersal and vicariance as causes of species richness and endemism. The diversification rate for Cochemiea is estimated to be approximately 12 times that of the mean background diversification rate for angiosperms. Divergence time estimation shows that many of the extant taxa in Cochemiea and Baja California Mammillaria emerged from common ancestors 1 million to 200,000 years ago, having a mid-Pleistocene origin.

CONCLUSIONS

Cochemiea and Mammillaria of the Baja California region are examples of recent, rapid diversification. Geological and climatic forces at multiple spatial and temporal scales are correlated with the western distributions of the Mammilloid clade.

Born in the mitochondrion and raised in the nucleus: evolution of a novel tandem repeat family in Medicago polymorpha (Fabaceae)

Plant nuclear genomes harbor sequence elements derived from the organelles (mitochondrion and plastid) through intracellular gene transfer (IGT). Nuclear genomes also show a dramatic range of repeat content, suggesting that any sequence can be readily amplified. These two aspects of plant nuclear genomes are well recognized but have rarely been linked. Through investigation of 31 Medicago taxa we detected exceptionally high post-IGT amplification of mitochondrial (mt) DNA sequences containing rps10 in the nuclear genome of Medicago polymorpha and closely related species. The amplified sequences were characterized as tandem arrays of five distinct repeat motifs (2157, 1064, 987, 971, and 587 bp) that have diverged from the mt genome (mitogenome) in the M. polymorpha nuclear genome. The mt rps10-like arrays were identified in seven loci (six intergenic and one telomeric) of the nuclear chromosome assemblies and were the most abundant tandem repeat family, representing 1.6-3.0% of total genomic DNA, a value approximately three-fold greater than the entire mitogenome in M. polymorpha. Compared to a typical mt gene, the mt rps10-like sequence coverage level was 691.5-7198-fold higher in M. polymorpha and closely related species. In addition to the post-IGT amplification, our analysis identified the canonical telomeric repeat and the species-specific satellite arrays that are likely attributable to an ancestral chromosomal fusion in M. polymorpha. A possible relationship between chromosomal instability and the mt rps10-like tandem repeat family in the M. polymorpha clade is discussed.

Origin and diversification of the saguaro cactus (Carnegiea gigantea): a within-species phylogenomic analysis

Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulf of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro's current range reflects post-glacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects.

Highly Resolved Papilionoid Legume Phylogeny Based on Plastid Phylogenomics

Comprising 501 genera and around 14,000 species, Papilionoideae is not only the largest subfamily of Fabaceae (Leguminosae; legumes), but also one of the most extraordinarily diverse clades among angiosperms. Papilionoids are a major source of food and forage, are ecologically successful in all major biomes, and display dramatic variation in both floral architecture and plastid genome (plastome) structure. Plastid DNA-based phylogenetic analyses have greatly improved our understanding of relationships among the major groups of Papilionoideae, yet the backbone of the subfamily phylogeny remains unresolved. In this study, we sequenced and assembled 39 new plastomes that are covering key genera representing the morphological diversity in the subfamily. From 244 total taxa, we produced eight datasets for maximum likelihood (ML) analyses based on entire plastomes and/or concatenated sequences of 77 protein-coding sequences (CDS) and two datasets for multispecies coalescent (MSC) analyses based on individual gene trees. We additionally produced a combined nucleotide dataset comprising CDS plus matK gene sequences only, in which most papilionoid genera were sampled. A ML tree based on the entire plastome maximally supported all of the deep and most recent divergences of papilionoids (223 out of 236 nodes). The Swartzieae, ADA (Angylocalyceae, Dipterygeae, and Amburaneae), Cladrastis, Andira, and Exostyleae clades formed a grade to the remainder of the Papilionoideae, concordant with nine ML and two MSC trees. Phylogenetic relationships among the remaining five papilionoid lineages (Vataireoid, Dermatophyllum, Genistoid s.l., Dalbergioid s.l., and Baphieae + Non-Protein Amino Acid Accumulating or NPAAA clade) remained uncertain, because of insufficient support and/or conflicting relationships among trees. Our study fully resolved most of the deep nodes of Papilionoideae, however, some relationships require further exploration. More genome-scale data and rigorous analyses are needed to disentangle phylogenetic relationships among the five remaining lineages.

Plastome Structural Evolution and Homoplastic Inversions in Neo-Astragalus (Fabaceae)

The plastid genomes of photosynthetic green plants have largely maintained conserved gene content and order as well as structure over hundreds of millions of years of evolution. Several plant lineages, however, have departed from this conservation and contain many plastome structural rearrangements, which have been associated with an abundance of repeated sequences both overall and near rearrangement endpoints. We sequenced the plastomes of 25 taxa of Astragalus L. (Fabaceae), a large genus in the inverted repeat-lacking clade of legumes, to gain a greater understanding of the connection between repeats and plastome inversions. We found plastome repeat structure has a strong phylogenetic signal among these closely related taxa mostly in the New World clade of Astragalus called Neo-Astragalus. Taxa without inversions also do not differ substantially in their overall repeat structure from four taxa each with one large-scale inversion. For two taxa with inversion endpoints between the same pairs of genes, differences in their exact endpoints indicate the inversions occurred independently. Our proposed mechanism for inversion formation suggests the short inverted repeats now found near the endpoints of the four inversions may be there as a result of these inversions rather than their cause. The longer inverted repeats now near endpoints may have allowed the inversions first mediated by shorter microhomologous sequences to propagate, something that should be considered in explaining how any plastome rearrangement becomes fixed regardless of the mechanism of initial formation.

The chicken or the egg? Plastome evolution and an independent loss of the inverted repeat in papilionoid legumes

The plastid genome (plastome), while surprisingly constant in gene order and content across most photosynthetic angiosperms, exhibits variability in several unrelated lineages. During the diversification history of the legume family Fabaceae, plastomes have undergone many rearrangements, including inversions, expansion, contraction and loss of the typical inverted repeat (IR), gene loss and repeat accumulation in both shared and independent events. While legume plastomes have been the subject of study for some time, most work has focused on agricultural species in the IR-lacking clade (IRLC) and the plant model Medicago truncatula. The subfamily Papilionoideae, which contains virtually all of the agricultural legume species, also comprises most of the plastome variation detected thus far in the family. In this study three non-papilioniods were included among 34 newly sequenced legume plastomes, along with 33 publicly available sequences, to assess plastome structural evolution in the subfamily. In an effort to examine plastome variation across the subfamily, approximately 20% of the sampling represents the IRLC with the remainder selected to represent the early-branching papilionoid clades. A number of IR-related and repeat-mediated changes were identified and examined in a phylogenetic context. Recombination between direct repeats associated with ycf2 resulted in intraindividual plastome heteroplasmy. Although loss of the IR has not been reported in legumes outside of the IRLC, one genistoid taxon was found to completely lack the typical plastome IR. The role of the IR and non-IR repeats in the progression of plastome change is discussed.

In and out: Evolution of viral sequences in the mitochondrial genomes of legumes (Fabaceae)

Plant specific mitoviruses in the 'genus' Mitovirus (Narnaviridae) and their integrated sequences (non-retroviral endogenous RNA viral elements or NERVEs) have been recently identified in various plant lineages. However, the sparse phylogenetic coverage of complete plant mitochondrial genome (mitogenome) sequences and the non-conserved nature of mitochondrial intergenic regions have hindered comparative studies on mitovirus NERVEs in plants. In this study, 10 new mitogenomes were sequenced from legumes (Fabaceae). Based on comparative genomic analysis of 27 total mitogenomes, we identified mitovirus NERVEs and transposable elements across the family. All legume mitogenomes included NERVEs and total NERVE length varied from ca. 2 kb in the papilionoid Trifolium to 35 kb in the mimosoid Acacia. Most of the NERVE integration sites were in highly variable intergenic regions, however, some were positioned in six cis-spliced mitochondrial introns. In the Acacia mitogenome, there were L1-like transposon sequences including an almost full-length copy with target site duplications (TSDs). The integration sites of NERVEs in four introns showed evidence of L1-like retrotransposition events. Phylogenetic analysis revealed that there were multiple instances of precise deletion of NERVEs between TSDs. This study provides clear evidence that a L1-like retrotransposition mechanism has a long history of contributing to the integration of viral RNA into plant mitogenomes while microhomology-mediated deletion can restore the integration site.

Projected climate change threatens significant range contraction of Cochemiea halei (Cactaceae), an island endemic, serpentine-adapted plant species at risk of extinction

AIM

Threats faced by narrowly distributed endemic plant species in the face of the Earth's sixth mass extinction and climate change exposure are especially severe for taxa on islands. We investigated the current and projected distribution and range changes of Cochemiea halei, an endemic island cactus. This taxon is of conservation concern, currently listed as vulnerable on the International Union for the Conservation of Nature Red List and as a species of special concern under Mexican federal law. The goals of this study are to (a) identify the correlations between climate variables and current suitable habitat for C. halei; (b) determine whether the species is a serpentine endemic or has a facultative relationship with ultramafic soils; and (c) predict range changes of the species based on climate change scenarios.

LOCATION

The island archipelago in Bahía Magdalena on the Pacific coast, Baja California Sur, Mexico.

METHODS

We used temperature and precipitation variables at 30-arc second resolution and soil type, employing multiple species distribution modeling methods, to identify important climate and soil conditions driving current habitat suitability. The best model of current suitability is used to predict possible effects of four climate change scenarios based on best-case to worst-case representative concentration pathways, with projected climate data from two general circulation models, over two time periods.

MAIN CONCLUSIONS

The occurrence of the species is found to be strongly correlated with ultramafic soils. The most important climate predictor for habitat suitability is annual temperature range. The species is predicted to undergo range contractions from 21% to 53%, depending on the severity and duration of exposure to climate change. The broader implications for a wide range of narrowly adapted, threatened, and endemic plant species indicate an urgent need for threat assessment based on habitat suitability and climate change modeling.

Diverse genomoviruses representing twenty-nine species identified associated with plants

Genomoviruses (family Genomoviridae) are circular single-stranded DNA viruses that have been mainly identified through metagenomics studies in a wide variety of samples from various environments. Here, we describe 98 genomes of genomoviruses found associated with members of 19 plant families from Australia, Brazil, France, South Africa and the USA. These 98 genomoviruses represent 29 species, 26 of which are new, in the genera Gemykolovirus (n = 37), Gemyduguivirus (n = 9), Gemygorvirus (n = 8), Gemykroznavirus (n = 6), Gemycircularvirus (n = 21) and Gemykibivirus (n = 17).

A Novel Divergent Geminivirus Identified in Asymptomatic New World Cactaceae Plants

Cactaceae comprise a diverse and iconic group of flowering plants which are almost exclusively indigenous to the New World. The wide variety of growth forms found amongst the cacti have led to the trafficking of many species throughout the world as ornamentals. Despite the evolution and physiological properties of these plants having been extensively studied, little research has focused on cactus-associated viral communities. While only single-stranded RNA viruses had ever been reported in cacti, here we report the discovery of cactus-infecting single-stranded DNA viruses. These viruses all apparently belong to a single divergent species of the family Geminiviridae and have been tentatively named Opuntia virus 1 (OpV1). A total of 79 apparently complete OpV1 genomes were recovered from 31 different cactus plants (belonging to 20 different cactus species from both the Cactoideae and Opuntioideae clades) and from nine cactus-feeding cochineal insects (Dactylopius sp.) sampled in the USA and Mexico. These 79 OpV1 genomes all share > 78.4% nucleotide identity with one another and < 64.9% identity with previously characterized geminiviruses. Collectively, the OpV1 genomes display evidence of frequent recombination, with some genomes displaying up to five recombinant regions. In one case, recombinant regions span ~40% of the genome. We demonstrate that an infectious clone of an OpV1 genome can replicate in Nicotiana benthamiana and Opuntia microdasys. In addition to expanding the inventory of viruses that are known to infect cacti, the OpV1 group is so distantly related to other known geminiviruses that it likely represents a new geminivirus genus. It remains to be determined whether, like its cactus hosts, its geographical distribution spans the globe.

Extensive gene tree discordance and hemiplasy shaped the genomes of North American columnar cacti

Few clades of plants have proven as difficult to classify as cacti. One explanation may be an unusually high level of convergent and parallel evolution (homoplasy). To evaluate support for this phylogenetic hypothesis at the molecular level, we sequenced the genomes of four cacti in the especially problematic tribe Pachycereeae, which contains most of the large columnar cacti of Mexico and adjacent areas, including the iconic saguaro cactus (Carnegiea gigantea) of the Sonoran Desert. We assembled a high-coverage draft genome for saguaro and lower coverage genomes for three other genera of tribe Pachycereeae (Pachycereus, Lophocereus, and Stenocereus) and a more distant outgroup cactus, Pereskia We used these to construct 4,436 orthologous gene alignments. Species tree inference consistently returned the same phylogeny, but gene tree discordance was high: 37% of gene trees having at least 90% bootstrap support conflicted with the species tree. Evidently, discordance is a product of long generation times and moderately large effective population sizes, leading to extensive incomplete lineage sorting (ILS). In the best supported gene trees, 58% of apparent homoplasy at amino sites in the species tree is due to gene tree-species tree discordance rather than parallel substitutions in the gene trees themselves, a phenomenon termed "hemiplasy." The high rate of genomic hemiplasy may contribute to apparent parallelisms in phenotypic traits, which could confound understanding of species relationships and character evolution in cacti.

A molecular phylogeny reveals the Cuban enigmatic genus Behaimia as a new piece in the Brongniartieae puzzle of papilionoid legumes

The papilionoid legume tribe Brongniartieae comprises a collection of 15 genera with disparate morphologies that were previously positioned in at least four remotely related tribes. The Brongniartieae displays a wide geographical disjunction between Australia and the New World and previous phylogenetic studies had provided conflicting results about the relationships between the American and Australian genera. We carry out phylogenetic analyses of (1) a plastid matK dataset extensively sampled across legumes to solve the enigmatic relationship of the Cuban-endemic monospecific genus Behaimia; and (2) multilocus datasets with focus on all genera ever referred to Brongniartieae. These analyses resulted in a well-resolved and strongly-supported phylogenetic tree of the Brongniartieae. The monophyly of all American genera of Brongniartieae is strongly supported. The doubtful position of the Australian genus Plagiocarpus is resolved within a clade comprising all Australian genera. Behaimia has been traditionally classified in tribe Millettieae, but our new molecular data and re-assessment of morphological traits have resolved the genus within the early-branching papilionoid tribe Brongniartieae. Characters including the pinnately multifoliolate (vs. unifoliolate) leaves, a sessile (vs. stipitate) ovary, and an indehiscent or late dehiscent one-seeded pod distinguish Behaimia from its closer relatives, the South American genera Cyclolobium and Limadendron.

Exceptional reduction of the plastid genome of saguaro cactus (Carnegiea gigantea): Loss of the ndh gene suite and inverted repeat

•

PREMISE OF THE STUDY

Land-plant plastid genomes have only rarely undergone significant changes in gene content and order. Thus, discovery of additional examples adds power to tests for causes of such genome-scale structural changes.•

METHODS

Using next-generation sequence data, we assembled the plastid genome of saguaro cactus and probed the nuclear genome for transferred plastid genes and functionally related nuclear genes. We combined these results with available data across Cactaceae and seed plants more broadly to infer the history of gene loss and to assess the strength of phylogenetic association between gene loss and loss of the inverted repeat (IR).•

KEY RESULTS

The saguaro plastid genome is the smallest known for an obligately photosynthetic angiosperm (∼113 kb), having lost the IR and plastid ndh genes. This loss supports a statistically strong association across seed plants between the loss of ndh genes and the loss of the IR. Many nonplastid copies of plastid ndh genes were found in the nuclear genome, but none had intact reading frames; nor did three related nuclear-encoded subunits. However, nuclear pgr5, which functions in a partially redundant pathway, was intact.•

CONCLUSIONS

The existence of an alternative pathway redundant with the function of the plastid NADH dehydrogenase-like complex (NDH) complex may permit loss of the plastid ndh gene suite in photoautotrophs like saguaro. Loss of these genes may be a recurring mechanism for overall plastid genome size reduction, especially in combination with loss of the IR.

Oberholzeria (Fabaceae subfam. Faboideae), a new monotypic legume genus from Namibia

Oberholzeria etendekaensis, a succulent biennial or short-lived perennial shrublet is described as a new species, and a new monotypic genus. Discovered in 2012, it is a rare species known only from a single locality in the Kaokoveld Centre of Plant Endemism, north-western Namibia. Phylogenetic analyses of molecular sequence data from the plastid matK gene resolves Oberholzeria as the sister group to the Genisteae clade while data from the nuclear rDNA ITS region showed that it is sister to a clade comprising both the Crotalarieae and Genisteae clades. Morphological characters diagnostic of the new genus include: 1) succulent stems with woody remains; 2) pinnately trifoliolate, fleshy leaves; 3) monadelphous stamens in a sheath that is fused above; 4) dimorphic anthers with five long, basifixed anthers alternating with five short, dorsifixed anthers, and 5) pendent, membranous, one-seeded, laterally flattened, slightly inflated but indehiscent fruits.

Evolution of a horizontally acquired legume gene, albumin 1, in the parasitic plant Phelipanche aegyptiaca and related species

Background

Parasitic plants, represented by several thousand species of angiosperms, use modified structures known as haustoria to tap into photosynthetic host plants and extract nutrients and water. As a result of their direct plant-plant connections with their host plant, parasitic plants have special opportunities for horizontal gene transfer, the nonsexual transmission of genetic material across species boundaries. There is increasing evidence that parasitic plants have served as recipients and donors of horizontal gene transfer (HGT), but the long-term impacts of eukaryotic HGT in parasitic plants are largely unknown.

Results

Here we show that a gene encoding albumin 1 KNOTTIN-like protein, closely related to the albumin 1 genes only known from papilionoid legumes, where they serve dual roles as food storage and insect toxin, was found in Phelipanche aegyptiaca and related parasitic species of family Orobanchaceae, and was likely acquired by a Phelipanche ancestor via HGT from a legume host based on phylogenetic analyses. The KNOTTINs are well known for their unique “disulfide through disulfide knot” structure and have been extensively studied in various contexts, including drug design. Genomic sequences from nine related parasite species were obtained, and 3D protein structure simulation tests and evolutionary constraint analyses were performed. The parasite gene we identified here retains the intron structure, six highly conserved cysteine residues necessary to form a KNOTTIN protein, and displays levels of purifying selection like those seen in legumes. The albumin 1 xenogene has evolved through >150 speciation events over ca. 16 million years, forming a small family of differentially expressed genes that may confer novel functions in the parasites. Moreover, further data show that a distantly related parasitic plant, Cuscuta, obtained two copies of albumin 1 KNOTTIN-like genes from legumes through a separate HGT event, suggesting that legume KNOTTIN structures have been repeatedly co-opted by parasitic plants.

Conclusions

The HGT-derived albumins in Phelipanche represent a novel example of how plants can acquire genes from other plants via HGT that then go on to duplicate, evolve, and retain the specialized features required to perform a unique host-derived function.

Revisiting the phylogeny of papilionoid legumes: New insights from comprehensively sampled early-branching lineages

PREMISE OF STUDY

Phylogenetic relationships of the papilionoid legumes (Papilionoideae) reveal that the early branches are more highly diverse in floral morphology than are other clades of Papilionoideae. This study attempts for the first time to comprehensively sample the early-branching clades of this economically and ecologically important legume subfamily and thus to resolve relationships among them. •

METHODS

Parsimony and Bayesian phylogenetic analyses of the plastid matK and trnL intron sequences included 29 genera not yet sampled in matK phylogenies of the Papilionoideae, 11 of which were sampled for DNA sequence data for the first time. •

KEY RESULTS

The comprehensively sampled matK phylogeny better resolved the deep-branching relationships and increased support for many clades within Papilionoideae. The potentially earliest-branching papilionoid clade does not include any genus traditionally assigned to tribe Swartzieae. Dipterygeae is monophyletic with the inclusion of Monopteryx. The genera Aldina and Amphimas represent two of the nine main but as yet unresolved lineages comprising the large 50-kb inversion clade within papilionoids. The quinolizidine-alkaloid-accumulating genistoid clade is expanded to include a strongly supported subclade containing Ormosia and the previously unplaced Clathrotropis s.s., Panurea, and Spirotropis. Camoensia is the first-branching genus of the core genistoids. •

CONCLUSIONS

The well-resolved phylogeny of the earliest-branching papilionoids generated in this study will greatly facilitate the efforts to redefine and stabilize the classification of this legume subfamily. Many key floral traits did not often predict phylogenetic relationships, so comparative studies on floral evolution and plant-animal interactions, for example, should also benefit from this study.

Keeping it simple: flowering plants tend to retain, and revert to, simple leaves

• A wide range of factors (developmental, physiological, ecological) with unpredictable interactions control variation in leaf form. Here, we examined the distribution of leaf morphologies (simple and complex forms) across angiosperms in a phylogenetic context to detect patterns in the directions of changes in leaf shape. • Seven datasets (diverse angiosperms and six nested clades, Sapindales, Apiales, Papaveraceae, Fabaceae, Lepidium, Solanum) were analysed using maximum likelihood and parsimony methods to estimate asymmetries in rates of change among character states. • Simple leaves are most frequent among angiosperm lineages today, were inferred to be ancestral in angiosperms and tended to be retained in evolution (stasis). Complex leaves slowly originated ('gains') and quickly reverted to simple leaves ('losses') multiple times, with a significantly greater rate of losses than gains. Lobed leaves may be a labile intermediate step between different forms. The nested clades showed mixed trends; Solanum, like the angiosperms in general, had higher rates of losses than gains, but the other clades had higher rates of gains than losses. • The angiosperm-wide pattern could be taken as a null model to test leaf evolution patterns in particular clades, in which patterns of variation suggest clade-specific processes that have yet to be investigated fully.

Phylogeny and character evolution in Medicago (Leguminosae): Evidence from analyses of plastid trnK/matK and nuclear GA3ox1 sequences

PREMISE OF THE STUDY

The genus Medicago, with about 87 species, includes the model legume species M. truncatula, and a number of important forage species such as M. sativa (alfalfa), M. scutellata (snail medic), and M. lupulina (black medic). Relationships within the genus are not yet sufficiently resolved, contributing to difficulty in understanding the evolution of a number of distinguishing characteristics such as aneuploidy and polyploidy, life history, structure of cotyledons, and number of seeds per fruit. •

METHODS

Phylogenetic relationships of 70-73 species of Medicago and its sister genus Trigonella (including Melilotus) were reconstructed from nucleotide sequences of the plastid trnK/matK region and the nuclear-encoded GA3ox1 gene (gibberellin 3-β-hydroxylase) using maximum parsimony and Bayesian inference methods. •

KEY RESULTS

Our results support certain currently recognized taxonomic groups, e.g., sect. Medicago (with M. sativa) and sect. Buceras. However, other strongly supported clades-the "reduced subsection Leptospireae clade" that includes M. lupulina, the "polymorpha clade" that includes M. murex and M. polymorpha and the "subsection Pachyspireae clade" that includes M. truncatula-each of which includes species presently in different subsections of sect. Spirocarpos, contradict the current classification. •

CONCLUSIONS

These results support the hypothesis that some characters considered important in existing taxonomies, for example, single-seeded fruits that have arisen more than once in both Medicago and Trigonella, are indeed homoplastic. Others, such as the 2n = 14 chromosome number, have also arisen independently within the genus. In addition, we demonstrate support for the utility of GA3ox1 sequences for phylogenetic analysis among and within closely related genera of legumes.

The evolution of a capacity to build supra-cellular ropes enabled filamentous cyanobacteria to colonize highly erodible substrates

BACKGROUND

Several motile, filamentous cyanobacteria display the ability to self-assemble into tightly woven or twisted groups of filaments that form macroscopic yarns or ropes, and that are often centimeters long and 50-200 microm in diameter. Traditionally, this trait has been the basis for taxonomic definition of several genera, notably Microcoleus and Hydrocoleum, but the trait has not been associated with any plausible function.

METHOD AND FINDINGS

Through the use of phylogenetic reconstruction, we demonstrate that pedigreed, rope-building cyanobacteria from various habitats do not form a monophyletic group. This is consistent with the hypothesis that rope-building ability was fixed independently in several discrete clades, likely through processes of convergent evolution or lateral transfer. Because rope-building cyanobacteria share the ability to colonize geologically unstable sedimentary substrates, such as subtidal and intertidal marine sediments and non-vegetated soils, it is also likely that this supracellular differentiation capacity imparts a particular fitness advantage in such habitats. The physics of sediment and soil erosion in fact predict that threads in the 50-200 microm size range will attain optimal characteristics to stabilize such substrates on contact.

CONCLUSIONS

Rope building is a supracellular morphological adaptation in filamentous cyanobacteria that allows them to colonize physically unstable sedimentary environments, and to act as successful pioneers in the biostabilization process.

Complete plastid genome sequence of the chickpea (Cicer arietinum) and the phylogenetic distribution of rps12 and clpP intron losses among legumes (Leguminosae)

Chickpea (Cicerarietinum, Leguminosae), an important grain legume, is widely used for food and fodder throughout the world. We sequenced the complete plastid genome of chickpea, which is 125,319bp in size, and contains only one copy of the inverted repeat (IR). The genome encodes 108 genes, including 4 rRNAs, 29 tRNAs, and 75 proteins. The genes rps16, infA, and ycf4 are absent in the chickpea plastid genome, and ndhB has an internal stop codon in the 5'exon, similar to other legumes. Two genes have lost their introns, one in the 3'exon of the transpliced gene rps12, and the one between exons 1 and 2 of clpP; this represents the first documented case of the loss of introns from both of these genes in the same plastid genome. An extensive phylogenetic survey of these intron losses was performed on 302 taxa across legumes and the related family Polygalaceae. The clpP intron has been lost exclusively in taxa from the temperate "IR-lacking clade" (IRLC), whereas the rps12 intron has been lost in most members of the IRLC (with the exception of Wisteria, Callerya, Afgekia, and certain species of Millettia, which represent the earliest diverging lineages of this clade), and in the tribe Desmodieae, which is closely related to the tribes Phaseoleae and Psoraleeae. Data provided here suggest that the loss of the rps12 intron occurred after the loss of the IR. The two new genomic changes identified in the present study provide additional support of the monophyly of the IR-loss clade, and resolution of the pattern of the earliest-branching lineages in this clade. The availability of the complete chickpea plastid genome sequence also provides valuable information on intergenic spacer regions among legumes and endogenous regulatory sequences for plastid genetic engineering.

Compound leaf development and evolution in the legumes

Across vascular plants, Class 1 KNOTTED1-like (KNOX1) genes appear to play a critical role in the development of compound leaves. An exception to this trend is found in the Fabaceae, where pea (Pisum sativum) uses UNIFOLIATA, an ortholog of the floral regulators FLORICAULA (FLO) and LEAFY (LFY), in place of KNOX1 genes to regulate compound leaf development. To assess the phylogenetic distribution of KNOX1-independent compound leaf development, a survey of KNOX1 protein expression across the Fabaceae was undertaken. The majority of compound-leafed Fabaceae have expression of KNOX1 proteins associated with developing compound leaves. However, in a large subclade of the Fabaceae, the inverted repeat-lacking clade (IRLC), of which pea is a member, KNOX1 expression is not associated with compound leaves. These data suggest that the FLO/LFY gene may function in place of KNOX1 genes in generating compound leaves throughout the IRLC. The contribution of FLO/LFY to leaf complexity in a member of the Fabaceae outside of the IRLC was examined by reducing expression of FLO/LFY orthologs in transgenic soybean (Glycine max). Transgenic plants with reduced FLO/LFY expression showed only slight reductions in leaflet number. Overexpression of a KNOX1 gene in alfalfa (Medicago sativa), a member of the IRLC, resulted in an increase in leaflet number. This implies that KNOX1 targets, which promote compound leaf development, are present in alfalfa and are still sensitive to KNOX1 regulation. These data suggest that KNOX1 genes and the FLO/LFY gene may have played partially overlapping roles in compound leaf development in ancestral Fabaceae but that the FLO/LFY gene took over this role in the IRLC.

Evolutionary rates analysis of Leguminosae implicates a rapid diversification of lineages during the tertiary

Tertiary macrofossils of the flowering plant family Leguminosae (legumes) were used as time constraints to estimate ages of the earliest branching clades identified in separate plastid matK and rbcL gene phylogenies. Penalized likelihood rate smoothing was performed on sets of Bayesian likelihood trees generated with the AIC-selected GTR+ Gamma +I substitution model. Unequivocal legume fossils dating from the Recent continuously back to about 56 million years ago were used to fix the family stem clade at 60 million years (Ma), and at 1-Ma intervals back to 70 Ma. Specific fossils that showed distinctive combinations of apomorphic traits were used to constrain the minimum age of 12 specific internal nodes. These constraints were placed on stem rather than respective crown clades in order to bias for younger age estimates. Regardless, the mean age of the legume crown clade differs by only 1.0 to 2.5 Ma from the fixed age of the legume stem clade. Additionally, the oldest caesalpinioid, mimosoid, and papilionoid crown clades show approximately the same age range of 39 to 59 Ma. These findings all point to a rapid family-wide diversification, and predict few if any legume fossils prior to the Cenozoic. The range of the matK substitution rate, 2.1-24.6 x 10(-10) substitutions per site per year, is higher than that of rbcL, 1.6- 8.6 x 10(-10), and is accompanied by more uniform rate variation among codon positions. The matK and rbcL substitution rates are highly correlated across the legume family. For example, both loci have the slowest substitution rates among the mimosoids and the fastest rates among the millettioid legumes. This explains why groups such as the millettioids are amenable to species-level phylogenetic analysis with these loci, whereas other legume groups are not.

Metacommunity process rather than continental tectonic history better explains geographically structured phylogenies in legumes

Penalized likelihood estimated ages of both densely sampled intracontinental and sparsely sampled transcontinental crown clades in the legume family show a mostly Quaternary to Neogene age distribution. The mode ages of the intracontinental crown clades range from 4-6 Myr ago, whereas those of the transcontinental crown clades range from 8-16 Myr ago. Both of these young age estimates are detected despite methodological approaches that bias results toward older ages. Hypotheses that resort to vicariance or continental history to explain continental disjunct distributions are dismissed because they require mostly Palaeogene and older tectonic events. An alternative explanation centring on dispersal that may well explain the geographical as well as the ecological phylogenetic structure of legume phylogenies is Hubbell's unified neutral theory of biodiversity and biogeography. This is the only dispersalist theory that encompasses evolutionary time and makes predictions about phylogenetic structure.

A phylogeny of legumes (Leguminosae) based on analysis of the plastid matK gene resolves many well-supported subclades within the family

Phylogenetic analysis of 330 plastid matK gene sequences, representing 235 genera from 37 of 39 tribes, and four outgroup taxa from eurosids I supports many well-resolved subclades within the Leguminosae. These results are generally consistent with those derived from other plastid sequence data (rbcL and trnL), but show greater resolution and clade support overall. In particular, the monophyly of subfamily Papilionoideae and at least seven major subclades are well-supported by bootstrap and Bayesian credibility values. These subclades are informally recognized as the Cladrastis clade, genistoid sensu lato, dalbergioid sensu lato, mirbelioid, millettioid, and robinioid clades, and the inverted-repeat-lacking clade (IRLC). The genistoid clade is expanded to include genera such as Poecilanthe, Cyclolobium, Bowdichia, and Diplotropis and thus contains the vast majority of papilionoids known to produce quinolizidine alkaloids. The dalbergioid clade is expanded to include the tribe Amorpheae. The mirbelioids include the tribes Bossiaeeae and Mirbelieae, with Hypocalypteae as its sister group. The millettioids comprise two major subclades that roughly correspond to the tribes Millettieae and Phaseoleae and represent the only major papilionoid clade marked by a macromorphological apomorphy, pseudoracemose inflorescences. The robinioids are expanded to include Sesbania and members of the tribe Loteae. The IRLC, the most species-rich subclade, is sister to the robinioids. Analysis of the matK data consistently resolves but modestly supports a clade comprising papilionoid taxa that accumulate canavanine in the seeds. This suggests a single origin for the biosynthesis of this most commonly produced of the nonprotein amino acids in legumes.

Improved bootstrap confidence limits in large-scale phylogenies, with an example from Neo-Astragalus (Leguminosae)

Phylogenetic analyses of large data sets pose special challenges, including the apparent tendency for the bootstrap support for a clade to decline with increased taxon sampling of that clade. We document this decline in data sets with increasing numbers of taxa in Astragalus, the most species-rich angiosperm genus. Support for one subclade, Neo-Astragalus, declined monotonically with increased sampling of taxa inside Neo-Astragalus, irrespective of whether parsimony or neighbor-joining methods were used or of which particular heuristic search algorithm was used (although more stringent algorithms tended to yield higher support). Three possible explanations for this decline were examined, including (1) mistaken assignment of the most recent common ancestor of the taxon sample (and its bootstrap support) with the most recent common ancestor of the clade from which it was sampled; (2) computational limitations of heuristic search strategies; and (3) statistical bias in bootstrap proportions, especially that from random homoplasy distributed among taxa. The best explanation appears to be (3), although computational shortcomings (2) may explain some of the problem. The bootstrap proportion, as currently used in phylogenetic analysis, does not accurately capture the classical notion of confidence assessments on the null hypothesis of nonmonophyly, especially in large data sets. More accurate assessments of confidence as type I error levels (relying on iterated bootstrap methods) remove most of the monotonic decline in confidence with increasing numbers of taxa.

PHYLOGENY OF THE CONJUGATING GREEN ALGAE (ZYGNEMOPHYCEAE) BASED ON rbc L SEQUENCES

Sequences of the gene encoding the large subunit of RUBISCO (rbcL) for 30 genera in the six currently recognized families of conjugating green algae (Desmidiaceae, Gonatozygaceae, Mesotaeniaceae, Peniaceae, and Zygnemataceae) were analyzed using maximum parsimony and maximum likelihood; bootstrap replications were performed as a measure of support for clades. Other Charophyceae sensu Mattox and Stewart and representative land plants were used as outgroups. All analyses supported the monophyly of the conjugating green algae. The Desmidiales, or placoderm desmids, constitute a monophyletic group, with moderate to strong support for the four component families of this assemblage (Closteriaceae, Desmidiaceae, Gonatozygaceae, and Peniaceae). The analyses showed that the two families of Zygnematales (Mesotaeniaceae, Zygnemataceae), which have plesiomorphic, unornamented and unsegmented cell walls, are not monophyletic. However, combined taxa of these two traditional families may constitute a monophyletic group. Partitioning the data by codon position revealed no significant differences across all positions or between partitions of positions one and two versus position three. The trees resulting from parsimony analyses using first plus second positions versus third position differed only in topology of branches with poor bootstrap support. The tree derived from third positions only was more resolved than the tree derived from first and second positions. The rbcL-based phylogeny is largely congruent with published analyses of small subunit rDNA sequences for the Zygnematales. The molecular data do not support hypotheses of monophyly for groups of extant unicellular and filamentous or colonial desmid genera exhibiting a common cell shape. A trend is evident from simple omniradiate cell shapes to taxa with lobed cell and plastid shapes, which supports the hypothesis that chloroplast shape evolved generally from simple to complex. The data imply that multicellular placoderm desmids are monophyletic. Several anomalous placements of genera were found, including the saccoderm desmid Roya in the Gonatozygaceae and the zygnematacean Entransia in the Coleochaetales. The former is strongly supported, although the latter is not, and Entransia's phylogenetic position warrants further study.

Error, bias, and long-branch attraction in data for two chloroplast photosystem genes in seed plants

Sequences of two chloroplast photosystem genes, psaA and psbB, together comprising about 3,500 bp, were obtained for all five major groups of extant seed plants and several outgroups among other vascular plants. Strongly supported, but significantly conflicting, phylogenetic signals were obtained in parsimony analyses from partitions of the data into first and second codon positions versus third positions. In the former, both genes agreed on a monophyletic gymnosperms, with Gnetales closely related to certain conifers. In the latter, Gnetales are inferred to be the sister group of all other seed plants, with gymnosperms paraphyletic. None of the data supported the modern "anthophyte hypothesis," which places Gnetales as the sister group of flowering plants. A series of simulation studies were undertaken to examine the error rate for parsimony inference. Three kinds of errors were examined: random error, systematic bias (both properties of finite data sets), and statistical inconsistency owing to long-branch attraction (an asymptotic property). Parsimony reconstructions were extremely biased for third-position data for psbB. Regardless of the true underlying tree, a tree in which Gnetales are sister to all other seed plants was likely to be reconstructed for these data. None of the combinations of genes or partitions permits the anthophyte tree to be reconstructed with high probability. Simulations of progressively larger data sets indicate the existence of long-branch attraction (statistical inconsistency) for third-position psbB data if either the anthophyte tree or the gymnosperm tree is correct. This is also true for the anthophyte tree using either psaA third positions or psbB first and second positions. A factor contributing to bias and inconsistency is extremely short branches at the base of the seed plant radiation, coupled with extremely high rates in Gnetales and nonseed plant outgroups.

Phylogenetic systematics of the tribe Millettieae (Leguminosae) based on chloroplast trnK/matK sequences and its implications for evolutionary patterns in Papilionoideae

Phylogenetic relationships in the tribe Millettieae and allies in the subfamily Papilionoideae (Leguminosae) were reconstructed from chloroplast trnK/matK sequences. Sixty-two accessions representing 57 traditionally recognized genera of Papilionoideae were sampled, including 27 samples from Millettieae. Phylogenies were constructed using maximum parsimony and are well resolved and supported by high bootstrap values. A well-supported "core Millettieae" clade is recognized, comprising the four large genera Millettia, Lonchocarpus, Derris, and Tephrosia. Several other small genera of Millettieae are not in the core Millettieae clade. Platycyamus is grouped with Phaseoleae (in part). Ostryocarpus, Austrosteenisia, and Dalbergiella are neither in the core Millettieae or Phaseoleae clade. These taxa, along with core Millettieae and Phaseoleae, form a monophyletic sister group to Indigofereae. Cyclolobium and Poecilanthe are close to Brongniartieae. Callerya and Wisteria belong to a large clade that includes all the legumes that lack the inverted repeat in their chloroplast genome, which confirms previous rbcL and phytochrome gene family phylogenies. The evolutionary history of four characters was examined in Millettieae and allies: the presence of canavanine, inflorescence types, the dehiscence of pods, and the presence of winged pods. trnK/matK sequence analysis suggests that the presence of a pseudoraceme or pseudopanicle and the accumulation of nonprotein amino acids are phylogenetically informative for Millettieae and allies with only a few exceptions.

Testing transposable elements as genetic drive mechanisms using Drosophila P element constructs as a model system

The use of transposable elements (TEs) as genetic drive mechanisms was explored using Drosophila melanogaster as a model system. Alternative strategies, employing autonomous and nonautonomous P element constructs were compared for their efficiency in driving the ry+ allele into populations homozygous for a ry- allele at the genomic rosy locus. Transformed flies were introduced at 1%, 5%, and 10% starting frequencies to establish a series of populations that were monitored over the course of 40 generations, using both phenotypic and molecular assays. The transposon-borne ry+ marker allele spread rapidly in almost all populations when introduced at 5% and 10% seed frequencies, but 1% introductions frequently failed to become established. A similar initial rapid increase in frequency of the ry+ transposon occurred in several control populations lacking a source of transposase. Constructs carrying ry+ markers also increased to moderate frequencies in the absence of selection on the marker. The results of Southern and in situ hybridization studies indicated a strong inverse relationship between the degree of conservation of construct integrity and transposition frequency. These finding have relevance to possible future applications of transposons as genetic drive mechanisms.

Persistent nuclear ribosomal DNA sequence polymorphism in the Amelanchier agamic complex (Rosaceae)

Individual plants of several Amelanchier taxa contain many polymorphic nucleotide sites in the internal transcribed spacers (ITS) of nuclear ribosomal DNA (nrDNA). This polymorphism is unusual because it is not recent in origin and thus has resisted homogenization by concerted evolution. Amelanchier ITS sequence polymorphism is hypothesized to be the result of gene flow between two major North American clades resolved by phylogenetic analysis of ITS sequences. Western North American species plus A. humilis and A. sanguinea of eastern North America form one clade (A), and the remaining eastern North American Amelanchier make up clade B. Five eastern North American taxa are polymorphic at many of the nucleotide sites where clades A and B have diverged and are thought to be of hybrid origin, with A. humilis or A. sanguinea as one parent and various members of clade B as the other parent. Morphological evidence suggests that A. humilis is one of the parents of one of the polymorphic taxa, a microspecies that we refer to informally as A. "erecta." Sequences of 21 cloned copies of the ITS1-5.8S gene-ITS2 region from one A. "erecta" individual are identical to A. humilis sequence or to the clade B consensus sequence, or they are apparent recombinants of A. humilis and clade B ITS repeats. Amelanchier "erecta" and another polymorphic taxon are suspected to be relatively old because both grow several hundred kilometers beyond the range of one of their parents. ITS sequence polymorphisms have apparently persisted in these two taxa perhaps because of polyploidy and/or agamospermy (asexual seed production), which are prevalent in the genus.

Paths and patterns: the biology and physics of swimming bacterial populations

The velocity distribution of swimming micro-organisms depends on directional cues supplied by the environment. Directional swimming within a bounded space results in the accumulation of organisms near one or more surfaces. Gravity, gradients of chemical concentration and illumination affect the motile behaviour of individual swimmers. Concentrated populations of organisms scatter and absorb light or consume molecules, such as oxygen. When supply is one-sided, consumption creates gradients; the presence of the population alters the intensity and the symmetry of the environmental cues. Patterns of cues interact dynamically with patterns of the consumer population. In suspensions, spatial variations in the concentration of organisms are equivalent to variations of mean mass density of the fluid. When organisms accumulate in one region whilst moving away from another region, the force of gravity causes convection that translocates both organisms and dissolved substances. The geometry of the resulting concentration-convection patterns has features that are remarkably reproducible. Of interest for biology are (1) the long-range organisation achieved by organisms that do not communicate, and (2) that the entire system, consisting of fluid, cells, directional supply of consumables, boundaries and gravity, generates a dynamic that improves the organisms' habitat by enhancing transport and mixing. Velocity distributions of the bacterium Bacillus subtilis have been measured within the milieu of the spatially and temporally varying oxygen concentration which they themselves create. These distributions of swimming speed and direction are the fundamental ingredients required for a quantitative mathematical treatment of the patterns. The quantitative measurement of swimming behaviour also contributes to our understanding of aerotaxis of individual cells.

Regulation of the SOS response in Bacillus subtilis: evidence for a LexA repressor homolog

The inducible SOS response for DNA repair and mutagenesis in the bacterium Bacillus subtilis resembles the extensively characterized SOS system of Escherichia coli. In this report, we demonstrate that the cellular repressor of the E. coli SOS system, the LexA protein, is specifically cleaved in B. subtilis following exposure of the cells to DNA-damaging treatments that induce the SOS response. The in vivo cleavage of LexA is dependent upon the functions of the E. coli RecA protein homolog in B. subtilis (B. subtilis RecA) and results in the same two cleavage fragments as produced in E. coli cells following the induction of the SOS response. We also show that a mutant form of the E. coli RecA protein (RecA430) can partially substitute for the nonfunctional cellular RecA protein in the B. subtilis recA4 mutant, in a manner consistent with its known activities and deficiencies in E. coli. RecA430 protein, which has impaired repressor cleaving (LexA, UmuD, and bacteriophage lambda cI) functions in E.coli, partially restores genetic exchange to B. subtilis recA4 strains but, unlike wild-type E. coli RecA protein, is not capable of inducing SOS functions (expression of DNA damage-inducible [din::Tn917-lacZ] operons or RecA synthesis) in B. subtilis in response to DNA-damaging agents or those functions that normally accompany the development of physiological competence. Our results provide support for the existence of a cellular repressor in B. subtilis that is functionally homologous to the E. coli LexA repressor and suggest that the mechanism by which B. subtilis RecA protein (like RecA of E. coli) becomes activated to promote the induction of the SOS response is also conserved.

DNA repair and the evolution of transformation in Bacillus subtilis. II. Role of inducible repair

In Bacillus subtilis, DNA repair and recombination are intimately associated with competence, the physiological state in which the bacterium can bind, take up and recombine exogenous DNA. Previously, we have shown that the homologous DNA transformation rate (ratio of transformants to total cells) increases with increasing UV dosage if cells are transformed after exposure to UV radiation (UV-DNA), whereas the transformation rate decreases if cells are transformed before exposure to UV (DNA-UV). In this report, by using different DNA repair-deficient mutants, we show that the greater increase in transformation rate in UV-DNA experiments than in DNA-UV experiments does not depend upon excision repair or inducible SOS-like repair, although certain quantitative aspects of the response do depend upon these repair systems. We also show that there is no increase in the transformation rate in a UV-DNA experiment when repair and recombination proficient cells are transformed with nonhomologous plasmid DNA, although the results in a DNA-UV experiment are essentially unchanged by using plasmid DNA. We have used din operon fusions as a sensitive means of assaying for the expression of genes under the control of the SOS-like regulon in both competent and noncompetent cell subpopulations as a consequence of competence development and our subsequent experimental treatments. Results indicate that the SOS-like system is induced in both competent and noncompetent subpopulations in our treatments and so should not be a major factor in the differential response in transformation rate observed in UV-DNA and DNA-UV treatments. These results provide further support to the hypothesis that the evolutionary function of competence is to bring DNA into the cell for use as template in the repair of DNA damage.

DNA repair and the evolution of transformation in the bacterium Bacillus subtilis

The purpose of the work reported here is to test the hypothesis that natural genetic transformation in the bacterium Bacillus subtilis has evolved as a DNA repair system. Specifically, tests were made to determine whether transformation functions to provide DNA template for the bacteria] cell to use in recombinational repair. The survivorship and the homologous transformation rate as a function of dose of ultraviolet irradiation (UV) was studied in two experimental treatments, in which cells were either transformed before (DNA-UV), or after (UV-DNA), treatment with UV. The results show that there is a qualitative difference in the relationship between the survival of transformed cells (sexual cells) and total cells (primarily asexual cells) in the two treatments. As predicted by the repair hypothesis, in the UV-DNA treatment, transformed cells had greater average survivorship than total cells, while in the DNA-UV treatment this relationship was reversed. There was also a consistent and qualitative difference between the UV-DNA and DNA-UV treatments in the relationship between the homologous transformation rate (transformed cells/total cells) and UV dosage. As predicted by the repair hypothesis, the homologous transformation rate increases with UV dose in the UV-DNA experiments but decreases with UV dose in the DNA-UV treatments. However, the transformation rate for plasmid DNA does not increase in a UV-DNA treatment. These results support the DNA repair hypothesis for the evolution of transformation in particular, and sex generally.

Inhibition of DNA methyltransferases in vitro by benzo[a]pyrene diol epoxide-modified substrates

Covalent adducts formed from the ultimate carcinogen 7 beta,8 alpha-dihydroxy-9 alpha, 10 alpha-epoxy-7,8,9,10-tetrahydrobenzo[ a]pyrene inhibit the enzyme-catalyzed transfer of methyl groups from S-adenosylmethionine to cytosine residues in DNA. Two DNA methyltransferase enzymes, isolated from the bacterium Haemophilus and mouse spleen nuclei, were tested for their ability to methylate carcinogen-modified substrates in vitro. These model enzymes possess the known methylation activities found in mammalian cells, de novo, and maintenance methylation of CpG-containing nucleotide sequences. The in vitro alkylation of DNA substrates by the carcinogen effectively decreases the methyltransferase reaction of both enzymes in a manner that is directly dependent upon the level of covalent modification of the DNA. Inhibition of de novo methylation activity can be detected at very low levels of carcinogen modification, 1 hydrocarbon residue per 20,000-40,000 nucleotides. Adduct levels in this range are capable of initiating transformation. Both enzymes are inactivated by direct reaction with the carcinogen in the absence of DNA. We also find that carcinogen adducts are capable of inhibiting DNA methylation at CpG sites removed from the primary lesion. These results support the proposal that carcinogen-induced DNA damage can cause alterations in methylation patterns that may eventually lead to heritable changes in gene expression.