Comparative linkage maps suggest that fission, not polyploidy, underlies near-doubling of chromosome number within monkeyflowers (Mimulus; Phrymaceae).

Changes in chromosome number and structure are important contributors to adaptation, speciation and macroevolution. In flowering plants, polyploidy and subsequent reductions in chromosome number by fusion are major sources of chromosomal evolution, but chromosome number increase by fission has been relatively unexplored. Here, we use comparative linkage mapping with gene-based markers to reconstruct chromosomal synteny within the model flowering plant genus Mimulus (monkeyflowers). Two sections of the genus with haploid numbers ≥ 14 have been inferred to be relatively recent polyploids because they are phylogenetically nested within numerous taxa with low base numbers (n=8-10). We combined multiple data sets to build integrated genetic maps of the M. guttatus species complex (section Simiolus, n=14) and the M. lewisii group (section Erythranthe; n=8), and then aligned the two integrated maps using >100 shared markers. We observed strong segmental synteny between M. lewisii and M. guttatus maps, with essentially 1-to-1 correspondence across each of 16 chromosomal blocks. Assuming that the M. lewisii (and widespread) base number of 8 is ancestral, reconstruction of 14 M. guttatus chromosomes requires at least eight fission events (likely shared by Simiolus and sister section Paradanthus (n=16)), plus two fusion events. This apparent burst of fission in the yellow monkeyflower lineages raises new questions about mechanisms and consequences of chromosomal fission in plants. Our comparative maps also provide insight into the origins of a chromosome exhibiting centromere-associated female meiotic drive and create a framework for transferring M. guttatus genome resources across the entire genus.

Is local adaptation in Mimulus guttatus caused by trade-offs at individual loci?

Local adaptation is considered to be the result of fitness trade-offs for particular phenotypes across different habitats. However, it is unclear whether such phenotypic trade-offs exist at the level of individual genetic loci. Local adaptation could arise from trade-offs of alternative alleles at individual loci or by complementary sets of loci with different fitness effects of alleles in one habitat but selective neutrality in the alternative habitat. To evaluate the genome-wide basis of local adaptation, we performed a field-based quantitative trait locus (QTL) mapping experiment on recombinant inbred lines (RILs) created from coastal perennial and inland annual races of the yellow monkeyflower (Mimulus guttatus) grown reciprocally in native parental habitats. Overall, we detected 19 QTLs affecting one or more of 16 traits measured in two environments, most of small effect. We identified 15 additional QTL effects at two previously identified candidate QTLs [DIVERGENCE (DIV)]. Significant QTL by environment interactions were detected at the DIV loci, which was largely attributable to genotypic differences at a single field site. We found no detectable evidence for trade-offs for any one component of fitness, although DIV2 showed a trade-off involving different fitness traits between sites, suggesting that local adaptation is largely controlled by non-overlapping loci. This is surprising for an outcrosser, implying that reduced gene flow prevents the evolution of individuals adapted to multiple environments. We also determined that native genotypes were not uniformly adaptive, possibly reflecting fixed mutational load in one of the populations.

Annotation and analysis of low-complexity protein families of Anopheles gambiae that are associated with cuticle.

We have characterized four new families of homologous genes of the mosquito, Anopheles gambiae, all of which include members shown by previous work to be cuticular in nature. The CPLCG, CPLCW, CPLCP, and CPLCA families (where CPLC is 'cuticular protein of low complexity') encode proteins with a high proportion of low-complexity sequence. We have also annotated the An. gambiae Tweedle genes, a family of cuticular protein genes first described in Drosophila, and additional ungrouped An. gambiae cuticular proteins identified by proteomics. Our annotations reveal multiple gene-family expansions that are specific to Diptera or Culicidae. The CPLCG and CPLCW families occur within a large and dynamic tandem array on chromosome 3R that includes sets of concertedly evolving genes. Most gene families exhibit two or more different expression profiles during development.

Isolated deficiencies of OXPHOS complexes I and IV are identified accurately and quickly by simple enzyme activity immunocapture assays.

OXPHOS deficits are associated with most reported cases of inherited, degenerative and acquired mitochondrial disease. Traditional methods of measuring OXPHOS activities in patients provide valuable clinical information but require fifty to hundreds of milligrams of biopsy tissue samples in order to isolate mitochondria for analysis. We have worked to develop assays that require less sample and here report novel immunocapture assays (lateral flow dipstick immunoassays) to determine the activities of complexes I and IV, which are far and away the most commonly affected complexes in the class of OXPHOS diseases. These assays are extremely simple to perform, rapid (1-1.5 h) and reproducible with low intra-assay and inter-assay coefficients of variability (CVs) s (<10%). Importantly, there is no need to purify mitochondria as crude extracts of whole cells or tissues are suitable samples. Therefore, the assays allow use of samples obtained non-invasively such as cheek swabs and whole blood, which are not amenable to traditional mitochondrial purification and OXPHOS enzyme analysis. As a first step to assess clinical utility of these novel assays, they were used to screen a panel of cultured fibroblasts derived from patients with isolated deficiencies in complex I or IV caused by identified genetic defects. All patients (5/5) with isolated complex IV deficiencies were identified in this population. Similarly, almost all (22/24) patients with isolated complex I deficiencies were identified. We believe that this assay approach should find widespread utility in initial screening of patients suspected of having mitochondrial disease.

Patterns of nucleotide variation and reproductive isolation between a Mimulus allotetraploid and its progenitor species.

Here we report our characterization of a widespread, highly selfing Mimulus allotetraploid formed by interspecific hybridization between M. nasutus and M. guttatus. Nucleotide variation at two nuclear loci (mCYCA and mAP3) within and among tetraploid populations resolves two haplotype clusters for each locus: one shares near identity with sequences from M. nasutus and the other group shares substantial variation with M. guttatus. With respect to the two loci studied, each allotetraploid individual is a 'fixed heterozygote' carrying sequences from both clusters. Moreover, mCYCA variation is consistent with at least two evolutionary origins for the Mimulus allotetraploid. We show that the allotetraploid is strongly reproductively isolated from M. nasutus and M. guttatus; interploidy crosses produce almost no viable seeds. By extension, we infer strong triploid block and argue that Mimulus allotetraploid formation might proceed in one step via the union of unreduced gametes in an M. nasutus-M. guttatus F(1) hybrid. We also discuss the potential roles of mating system and flowering asynchrony in allotetraploid establishment.

Is floral diversification associated with pollinator divergence? Flower shape, flower colour and pollinator preference in Chilean Mimulus.

BACKGROUND AND AIMS: Adaptation to different pollinators is thought to drive divergence in flower colour and morphology, and may lead to interspecific reproductive isolation. Floral diversity was tested for association with divergent pollinator preferences in a group of four closely related wildflower species: the yellow-flowered Mimulus luteus var. luteus and the red-pigmented M. l. variegatus, M. naiandinus and M. cupreus. METHODS: Patterns of pollinator visitation were evaluated in natural plant populations in central Chile, including both single-species and mixed-species sites. Floral anthocyanin pigments were identified, and floral morphology and nectar variation were quantified in a common garden experiment using seeds collected from the study sites. KEY RESULTS: Mimulus l. luteus, M. l. variegatus and M. naiandinus are morphologically similar and share a single generalist bumblebee pollinator, Bombus dahlbomii. Mimulus cupreus differs significantly from the first three taxa in corolla shape as well as nectar characteristics, and had far fewer pollinator visits. CONCLUSIONS: This system shows limited potential for pollinator-mediated restriction of gene flow as a function of flower colour, and no evidence of transition to a novel pollinator. Mimulus cupreus may experience reduced interspecific gene flow due to a lack of bumblebee visitation, but not because of its red pigmentation: rare yellow morphs are equally undervisited by pollinators. Overall, the results suggest that factors other than pollinator shifts may contribute to the maintenance of floral diversity in these Chilean Mimulus species.

Mimulus is an emerging model system for the integration of ecological and genomic studies.

The plant genus Mimulus is rapidly emerging as a model system for studies of evolutionary and ecological functional genomics. Mimulus contains a wide array of phenotypic, ecological and genomic diversity. Numerous studies have proven the experimental tractability of Mimulus in laboratory and field studies. Genomic resources currently under development are making Mimulus an excellent system for determining the genetic and genomic basis of adaptation and speciation. Here, we introduce some of the phenotypic and genetic diversity in the genus Mimulus and highlight how direct genetic studies with Mimulus can address a wide spectrum of ecological and evolutionary questions. In addition, we present the genomic resources currently available for Mimulus and discuss future directions for research. The integration of ecology and genetics with bioinformatics and genome technology offers great promise for exploring the mechanistic basis of adaptive evolution and the genetics of speciation.

A conserved domain in arthropod cuticular proteins binds chitin.

Many insect cuticular proteins include a 35-36 amino acid motif known as the R&R consensus. The extensive conservation of this region led to the suggestion that it functions to bind chitin. Provocatively, it has no sequence similarity to the well-known cysteine-containing chitin-binding domain found in chitinases and some peritrophic membrane proteins. Using fusion proteins expressed in E. coli, we show that an extended form of the R&R consensus from proteins of hard cuticles is necessary and sufficient for chitin binding. Recombinant AGCP2b, a putative cuticular protein from the mosquito Anopheles gambiae, was expressed in E. coli and the purified protein shown to bind to chitin beads. A stretch of 65 amino acids from AGCP2b, including the R&R consensus, conferred chitin binding to glutathione-S-transferase (GST). Directed mutagenesis of some conserved amino acids within this extended R&R consensus from hard cuticle eliminated chitin binding. Thus arthropods have two distinct classes of chitin binding proteins, those with the chitin-binding domain found in lectins, chitinases and peritrophic membranes (cysCBD) and those with the cuticular protein chitin-binding domain (non-cysCBD).

"Soft"-cuticle protein secondary structure as revealed by FT-Raman, ATR FT-IR and CD spectroscopy.

The nature of the interaction of insect cuticular proteins and chitin is unknown even though about half of the cuticular proteins sequenced thus far share a consensus region that has been predicted to be the site of chitin binding. We previously predicted the preponderance of a beta-pleated sheet in the consensus region and proposed its responsibility for the formation of helicoidal cuticle (Iconomidou et al., Insect Biochem. Mol. Biol. 29 (1999) 285). In this study, we examined experimentally the secondary structure of intact and guanidine hydrochloride extracted cuticle and the cuticular protein extract. The studied cuticle came from the larval dorsal abdomen of the lepidopteran Hyalophora cecropia, a classical example of "soft" cuticle. Analysis with FT-Raman, ATR FT-IR and CD spectroscopy indicates that antiparallel beta-pleated sheet is the predominant molecular conformation of "soft-cuticle" proteins both in situ in the cuticle and following extraction. It seems that this conformation dictates the modes of chitin-protein interaction in cuticle, in agreement with earlier proposals (Atkins, J. Biosci. 8 (1985) 375).

Deleterious mutations and genetic variation for flower size in Mimulus guttatus.

Mimulus guttatus is a wildflower that exhibits substantial genetic variation in flower size. Here, we test the hypothesis that this variation is caused by deleterious mutations maintained through mutation-selection balance. The deleterious-mutation model predicts that rare, partially recessive alleles will be the primary source of variation. We test this prediction by measuring the change in the mean flower size (deltaM) and the directional dominance of flower size (deltaB) within a selection experiment. If variation is due to rare (partially) recessive alleles, deltaB/deltaM is expected to be positive and exceed one. However, we obtain negative values for deltaB/deltaM from three independent selection lines. This result is statistically inconsistent with the deleterious-mutation model.

A genetic map in the Mimulus guttatus species complex reveals transmission ratio distortion due to heterospecific interactions.

As part of a study of the genetics of floral adaptation and speciation in the Mimulus guttatus species complex, we constructed a genetic linkage map of an interspecific cross between M. guttatus and M. nasutus. We genotyped an F(2) mapping population (N = 526) at 255 AFLP, microsatellite, and gene-based markers and derived a framework map through repeated rounds of ordering and marker elimination. The final framework map consists of 174 marker loci on 14 linkage groups with a total map length of 1780 cM Kosambi. Genome length estimates (2011-2096 cM) indicate that this map provides thorough coverage of the hybrid genome, an important consideration for QTL mapping. Nearly half of the markers in the full data set (49%) and on the framework map (48%) exhibited significant transmission ratio distortion (alpha = 0.05). We localized a minimum of 11 transmission ratio distorting loci (TRDLs) throughout the genome, 9 of which generate an excess of M. guttatus alleles and a deficit of M. nasutus alleles. This pattern indicates that the transmission ratio distortion results from particular interactions between the heterospecific genomes and suggests that substantial genetic divergence has occurred between these Mimulus species. We discuss possible causes of the unequal representation of parental genomes in the F(2) generation.

Evidence for Dobzhansky-Muller incompatibilites contributing to the sterility of hybrids between Mimulus guttatus and M. nasutus.

Both chromosomal rearrangements and negative interactions among loci (Dobzhansky-Muller incompatibilities) have been advanced as the genetic mechanism underlying the sterility of interspecific hybrids. These alternatives invoke very different evolutionary histories during speciation and also predict different patterns of sterility in artificial hybrids. Chromosomal rearrangements require drift, inbreeding, or other special conditions for initial fixation and, because heterozygosity per se generates any problems with gamete formation, F1 hybrids will be most infertile. In contrast, Dobzhansky-Muller incompatibilities may arise as byproducts of adaptive evolution and often affect the segregating F2 generation most severely. To distinguish the effects of these two mechanisms early in divergence, we investigated the quantitative genetics of hybrid sterility in a line cross between two members of the Mimulus guttatus species complex (M. guttatus and M. nasutus). Hybrids showed partial male and female sterility, and the patterns of infertility were not consistent with the action of chromosomal rearrangements alone. F2 and F1 hybrids exhibited equal decreases in pollen viability (> 40%) relative to the highly fertile parental lines. A large excess of completely pollen-sterile F2 genotypes also pointed to the segregation of Dobzhansky-Muller incompatibility factors affecting male fertility. Female fertility showed a pattern similarly consistent with epistatic interactions: F2 hybrids produced far fewer seeds per flower than F1 hybrids (88.0 +/- 2.8 vs. 162.9 +/- 8.5 SE, respectively) and either parental line, and many F2 genotypes were completely female sterile. Dobzhansky-Muller interactions also resulted in the breakdown of several nonreproductive characters and appear to contribute to correlations between male and female fertility in the F2 generation. These results parallel and contrast with the genetics of postzygotic isolation in model animal systems and are a first step toward understanding the process of speciation in this well-studied group of flowering plants.

Inbreeding load, average dominance and the mutation rate for mildly deleterious alleles in Mimulus guttatus.

The goal of this study is to provide information on the genetics of inbreeding depression in a primarily outcrossing population of Mimulus guttatus. Previous studies of this population indicate that there is tremendous inbreeding depression for nearly every fitness component and that almost all of this inbreeding depression is due to mildly deleterious alleles rather than recessive lethals or steriles. In this article I assayed the homozygous and heterozygous fitnesses of 184 highly inbred lines extracted from a natural population. Natural selection during the five generations of selfing involved in line formation essentially eliminated major deleterious alleles but was ineffective in purging alleles with minor fitness effects and did not appreciably diminish overall levels of inbreeding depression. Estimates of the average degree of dominance of these mildly deleterious alleles, obtained from the regression of heterozygous fitness on the sum of parental homozygous fitness, indicate that the detrimental alleles are partially recessive for most fitness traits, with h approximately 0.15 for cumulative measures of fitness. The inbreeding load, B, for total fitness is approximately 1.0 in this experiment. These results are consistent with the hypothesis that spontaneous mildly deleterious mutations occur at a rate >0.1 mutation per genome per generation.

The contribution of male-sterility mutations to inbreeding depression in Mimulus guttatus.

The magnitude of inbreeding depression can influence many aspects of a population's ecology and evolution, including the nature of selection acting on the mating system and the chances that the population will go extinct during periods of small population size. If inbreeding depression is caused primarily by recessive mutations of large effect on fitness, such as lethals or steriles, then it is expected to be purged rapidly by selection with moderate amounts of inbreeding. In contrast, inbreeding depression primarily caused by many genes with mild effects on fitness will not be rapidly purged with inbreeding, so it should be a more resilient barrier to the evolution of self-fertilization and a more significant threat to the survival of endangered species. Here I show that recessive male-sterility alleles at individual loci are common in a primarily outcrossing population of the plant Mimulus guttatus. Despite the high frequency of these major mutations, most of the inbreeding depression for male fertility and cumulative measures of lifetime fitness results from more mildly deleterious alleles. Male-sterility alleles contribute to 31% of the inbreeding depression for the fraction of viable pollen grains, and to 26% of the inbreeding depression for total fitness. These results suggest that most of the inbreeding depression for male fertility in this population would not be purged, in the short term, with moderate inbreeding.

Spontaneous deleterious mutation in Arabidopsis thaliana.

The frequency and selective impact of deleterious mutations are fundamental parameters in evolutionary theory, yet they have not been directly measured in a plant species. To estimate these quantities, we allowed spontaneous mutations to accumulate for 10 generations in 1,000 inbred lines of the annual, self-fertilizing plant Arabidopsis thaliana and assayed fitness differences between generations 0 and 10 in a common garden. Germination rate, fruit set, and number of seeds per fruit each declined by less than 1% per generation in the mutation lines, and total fitness declined by 0.9% per generation. Among-line variances increased in the mutation lines for all traits. Application of an equal-effects model suggests a downwardly biased genomic deleterious mutation rate of 0.1 and a upwardly biased effect of individual mutations on total fitness of 20%. This genomic deleterious mutation rate is consistent with estimates of nucleotide substitution rates in flowering plants, the genome size of Arabidopsis, and the equilibrium inbreeding depression observed in this highly selfing plant species.

Is beta-pleated sheet the molecular conformation which dictates formation of helicoidal cuticle?

Over 100 sequences for cuticular proteins are now available, but there have been no formal analyses of how these sequences might contribute to the helicoidal architecture of cuticle or to the interaction of these proteins with chitin. A secondary structure prediction scheme (Hamodrakas, S.J., 1988. A protein secondary structure prediction scheme for the IBM PC and compatibles. CABIOS 4, 473-477) that combines six different algorithms predicting alpha-helix, beta-strands and beta-turn/loops/coil has been used to predict the secondary structure of chorion proteins and experimental confirmation has established its utility (Hamodrakas, S.J., 1992. Molecular architecture of helicoidal proteinaceous eggshells. In: Case, S.T. (Ed.), Results and Problems in Cell Differentiation, Vol. 19, Berlin-Heidelberg, Springer Verlag, pp. 116-186 and references therein). We have used this same scheme with eight cuticular protein sequences associated with hard cuticles and nineteen from soft cuticles. Secondary structure predictions were restricted to a conserved 68 amino acid region that begins with a preponderance of hydrophilic residues and ends with a 33 amino acid consensus region, first identified by Rebers and Riddiford (Rebers, J.F., Riddiford, L.M., 1988. Structure and expression of a Manduca sexta larval cuticle gene homologous to Drosophila cuticle genes. J. Mol. Biol. 203, 411-423). Both classes of sequences showed a preponderance of beta-pleated sheet, with four distinct strands in the proteins from 'hard' cuticles and three from 'soft'. In both cases, tyrosine and phenylalanine were found on one face within a sheet, an optimal location for interaction with chitin. We propose that this beta-sheet dictates formation of helicoidal cuticle.

Recombination, RNA evolution, and bifunctional RNA molecules isolated through chimeric SELEX.

Exchange of RNA structural domains through recombination can be used to engineer RNAs with novel functions and may have played an important role in the early evolution of life. The degree of function an RNA element retains upon recombination into a new sequence context is a measure of how deleterious or beneficial recombination will be. When we fused pairs of aptamers previously selected to bind coenzyme A, chloramphenicol, or adenosine, the chimerae retained some ability to bind both targets, but with reduced binding activity both in solution and on affinity resins, probably due to misfolding. Complex populations of recombined RNAs gave similar results. Applying dual selection pressure to recombined populations yielded the combinations that were best suited to binding both targets. Most reselected RNAs folded into the active conformation more readily than chimerae built from arbitrarily chosen aptamers, as indicated both by solution Kd measurements and affinity resin binding activity. Deletion/selection experiments confirmed that the sequences required for binding are fully contained within the respective domains and not derived from interaction between the domains, consistent with the modular architecture of their original design. The combinatorial nature of the recombination methods presented here takes advantage of the full sequence diversity of the starting populations and yields large numbers of bifunctional molecules (10(6) to more than 1012). The method can be easily generalized and should be applicable to engineering dual-function RNAs for a wide variety of applications, including catalysis, novel therapeutics, and studies of long-range RNA structure.

A family of pupal-specific cuticular protein genes in the mosquito Anopheles gambiae.

We have cloned and sequenced members of a cuticular protein multi-gene family from the mosquito Anopheles gambiae. Three genes (agcp2a-c), each approximately 1 kb in length, were found in a 17.4 kb genomic phage clone. Analysis of ten cDNAs revealed that at least four related genes are present. The open reading frame of the genes and cDNAs showed 95% sequence identity. Divergence was observed in the sequence of the 3' ends and the number of copies of two repeated coding sequences. In situ hybridizations with a probe prepared from one of these circular protein genes physically mapped to two loci, 26B on chromosome 2L and 37A on 3R. Transcription of these An. gambiae cuticular protein genes appears to be limited to pharate pupae and the expressed protein(s) is found in early pupae. The deduced amino acid sequence of these proteins contains a hydrophilic region with significant similarity to other cuticular proteins including the pupal-specific cuticular protein, EDG84, of Drosophila melanogaster (Apple and Fristrom).

Proteins of crustacean exoskeleton: IV. Partial amino acid sequences of exoskeletal proteins from the Bermuda land crab, Gecarcinus lateralis, and comparisons to certain insect proteins.

As in all decapod Crustacea, the exoskeleton of the land crab Gecarcinus lateralis consists of four layers. Prior electrophoretic analysis of proteins extracted from these layers revealed an abundance of small M(r) proteins with acidic pIs are found in insect cuticle (O'Brien et al. [1991 Biol. Bull., 181:427-441). Further, immunological cross-reactivity between crab exoskeletal proteins and insect cuticular proteins has been demonstrated (Kumari and Skinner [1993] J. Exp. Zool., 265:195-210). Partial amino acid sequences of a number of proteins from the four exoskeletal layers are described here. Proteins were electrophoresed on two-dimensional (2D) gels, transferred to polyvinylidene difluoride (PVDF) membranes, and stained; individual spots were recovered and their N-termini were sequenced. In addition, a 14-kDa protein (pI = 5.4) from membranous layer (ML14) was eluted from 2D gels and digested with endoproteinase Lys-C; N-termini of its constituent peptides were sequenced. The two epicuticular proteins differed from each other. Three proteins with identical electrophoretic mobility isolated from exocuticle, endocuticle, and membranous layer appeared to have identical N termini, while another electrophoretically identical set from the three layers appeared identical with each other but differed in three positions from the first set. Two proteins from the membranous layer both had a mass of 25 kDa but different isoelectric points. Their sequences were indistinguishable from each other but clearly distinct from another membranous layer protein. Another distinct sequence was found in a 14-kDa protein from endocuticle, while a less acidic pair of 14-kDa proteins from endocuticle and membranous layer were quite similar to one another. The three internal peptide fragments from ML14 were distinct, but one had regions similar to the ML14 N terminus. One crab exoskeletal protein sequence was similar to some structural proteins of vertebrates, whereas others had motifs found in insect cuticular proteins. The sequence similarities identified did not account for the antibody cross-reactivity.

Individual variation in inbreeding depression: the roles of inbreeding history and mutation.

We use mutation-selection recursion models to evaluate the relative contributions of mutation and inbreeding history to variation among individuals in inbreeding depression and the ability of experiments to detect associations between individual inbreeding depression and mating system genotypes within populations. Poisson mutation to deleterious additive or recessive alleles generally produces far more variation among individuals in inbreeding depression than variation in history of inbreeding, regardless of selfing rate. Moreover, variation in inbreeding depression can be higher in a completely outcrossing or selfing population than in a mixed-mating population. In an initially random mating population, the spread of a dominant selfing modifier with no pleiotropic effects on male outcross success causes a measurable increase in inbreeding depression variation if its selfing rate is large and inbreeding depression is caused by recessive lethals. This increase is observable during a short period as the modifier spreads rapidly to fixation. If the modifier alters selfing rate only slightly, it fails to spread or causes no measurable increase in inbreeding depression variance. These results suggest that genetic associations between mating loci and inbreeding depression loci could be difficult to demonstrate within populations and observable only transiently during rapid evolution to a substantially new selfing rate.