Evolutionary developmental biology and the problem of variation

Abstract. One of the oldest problems in evolutionary biology remains largely unsolved. Which mutations generate evolutionarily relevant phenotypic variation? What kinds of molecular changes do they entail? What are the phenotypic magnitudes, frequencies of origin, and pleiotropic effects of such mutations? How is the genome constructed to allow the observed abundance of phenotypic diversity? Historically, the neo‐Darwinian synthesizers stressed the predominance of micromutations in evolution, whereas others noted the similarities between some dramatic mutations and evolutionary transitions to argue for macromutationism. Arguments on both sides have been biased by misconceptions of the developmental effects of mutations. For example, the traditional view that mutations of important developmental genes always have large pleiotropic effects can now be seen to be a conclusion drawn from observations of a small class of mutations with dramatic effects. It is possible that some mutations, for example, those in cis‐regulatory DNA, have few or no pleiotropic effects and may be the predominant source of morphological evolution. In contrast, mutations causing dramatic phenotypic effects, although superficially similar to hypothesized evolutionary transitions, are unlikely to fairly represent the true path of evolution. Recent developmental studies of gene function provide a new way of conceptualizing and studying variation that contrasts with the traditional genetic view that was incorporated into neo‐Darwinian theory and population genetics. This new approach in developmental biology is as important for micro‐evolutionary studies as the actual results from recent evolutionary developmental studies. In particular, this approach will assist in the task of identifying the specific mutations generating phenotypic variation and elucidating how they alter gene function. These data will provide the current missing link between molecular and phenotypic variation in natural populations.     

Sensitivity to the two‐peptide bacteriocin lactococcin G is dependent on UppP,an enzyme involved in cell‐wall synthesis

Most bacterially produced antimicrobial peptides (bacteriocins) are thought to kill target cells by a receptor‐mediated mechanism. However, for most bacteriocins the receptor is unknown. For instance, no target receptor has been identified for the two‐peptide bacteriocins (class IIb), whose activity requires the combined action of two individual peptides. To identify the receptor for the class IIb bacteriocin lactococcin G, which targets strains of Lactococcus lactis, we generated 12 lactococcin G‐resistant mutants and performed whole‐genome sequencing to identify mutations causing the resistant phenotype. Remarkably, all had a mutation in or near the gene uppP (bacA), encoding an undecaprenyl pyrophosphate phosphatase; a membrane protein involved in peptidoglycan synthesis. Nine mutants had stop codons or frameshifts in the uppP gene, two had point mutations in putative regulatory regions and one caused an amino acid substitution in UppP. To verify the receptor function of UppP, it was shown that growth of non‐sensitive Streptococcus pneumoniae could be inhibited by lactococcin G when L. lactis uppP was expressed in this bacterium. Furthermore, we show that the related class IIb bacteriocin enterocin 1071 also uses UppP as receptor. The approach used here should be broadly applicable to identify receptors for other bacteriocins as well.     

Variegation mutants and mechanisms of chloroplast biogenesis

Variegated plants typically have green‐ and white‐sectored leaves. Cells in the green sectors contain normal‐appearing chloroplasts, whereas cells in the white sectors lack pigments and appear to be blocked at various stages of chloroplast biogenesis. Variegations can be caused by mutations in nuclear, chloroplast or mitochondrial genes. In some plants, the green and white sectors have different genotypes, but in others they have the same (mutant) genotype. One advantage of variegations is that they provide a means of studying genes for proteins that are important for chloroplast development, but for which mutant analysis is difficult, either because mutations in a gene of interest are lethal or because they do not show a readily distinguishable phenotype. This paper focuses on Arabidopsis variegations, for which the most information is available at the molecular level. Perhaps the most interesting of these are variegations caused by defective nuclear gene products in which the cells of the mutant have a uniform genotype. Two questions are of paramount interest: (1) What is the gene product and how does it function in chloroplast biogenesis? (2) What is the mechanism of variegation and why do green sectors arise in plants with a uniform (mutant) genotype? Two paradigms of variegation mechanism are described: immutans (im) and variegated2 (var2). Both mechanisms emphasize compensating activities and the notion of plastid autonomy, but redundant gene products are proposed to play a role in var2, but not in im. It is hypothesized that threshold levels of certain activities are necessary for normal chloroplast development.     

Novel mutants of the aubergine gene

Aubergine is an RNA-binding protein of the Piwi clade, functioning in germline in the piRNA pathway that silences transposons and repetitive sequences. Several mutations of this gene exist, but they mostly result in truncated proteins or correspond to mutations that also affect neighboring genes. We have generated complete aubergine knock-out mutants that do not disrupt the neighboring genes. These novel mutants are characterized by PCR and sequencing. Their nature is confirmed by female sterility and by the presence of crystals in testes, common to the aubergine loss of function mutations. These mutants provide novel and more appropriate tools for the study of the piRNA pathway that controls genome stability.     

Forward and Reverse Genetic Analysis of Microtubule Motors in Chlamydomonas

The ability to integrate biochemical, cell biological, and genetic approaches makes Chlamydomonas reinhardtii the premier model organism for studies of the eukaryotic flagellum and its associated molecular motors. Hundreds of motility mutations have been identified in Chlamydomonas, including many that affect dyneins and kinesins. These mutations have yielded much information on the structure and function of the motors as well as the roles of individual subunits within the motors. The development of insertional mutagenesis has opened the door to powerful new approaches for genetic analysis in Chlamydomonas. Insertional mutants are created by transforming cells with DNA-containing selectable markers. The DNA is randomly integrated throughout the genome and usually deletes part of the chromosome at the site of insertion, thereby creating mutations that are marked by the integrated DNA. These mutations can be used for forward genetic approaches where one characterizes a mutant phenotype and then clones the relevant gene using the integrated DNA as a tag. The insertional mutants also may be used in a reverse genetic approach in which mutants lacking a gene of interest are identified by DNA hybridization. We describe methods to generate and characterize insertional mutants, using mutations that affect the outer dynein arm as examples.     

Rapid identification of lettuce seed germination mutants by bulked segregant analysis and whole genome sequencing

Lettuce (Lactuca sativa) seeds exhibit thermoinhibition, or failure to complete germination when imbibed at warm temperatures. Chemical mutagenesis was employed to develop lettuce lines that exhibit germination thermotolerance. Two independent thermotolerant lettuce seed mutant lines, TG01 and TG10, were generated through ethyl methanesulfonate mutagenesis. Genetic and physiological analyses indicated that these two mutations were allelic and recessive. To identify the causal gene(s), we applied bulked segregant analysis by whole genome sequencing. For each mutant, bulked DNA samples of segregating thermotolerant (mutant) seeds were sequenced and analyzed for homozygous single‐nucleotide polymorphisms. Two independent candidate mutations were identified at different physical positions in the zeaxanthin epoxidase gene (ABSCISIC ACID DEFICIENT 1/ZEAXANTHIN EPOXIDASE, or ABA1/ZEP) in TG01 and TG10. The mutation in TG01 caused an amino acid replacement, whereas the mutation in TG10 resulted in alternative mRNA splicing. Endogenous abscisic acid contents were reduced in both mutants, and expression of the ABA1 gene from wild‐type lettuce under its own promoter fully complemented the TG01 mutant. Conventional genetic mapping confirmed that the causal mutations were located near the ZEP/ABA1 gene, but the bulked segregant whole genome sequencing approach more efficiently identified the specific gene responsible for the phenotype.     

Dbest1, a drosophila homolog of human bestrophin,is not required for viability or photoreceptor integrity

Summary: Best macular dystrophy (BMD) is an autosomal dominant human disease characterized by macular degeneration with juvenile onset (OMIM 153700). The disease is most often associated with mutations in Bestrophin, which encodes a novel protein with four putative transmembrane domains. However, complete loss‐of‐function mutations in Bestrophin have not been reported in humans or mice. We have identified three homologs of human Bestrophin in the Drosophila genome (dbest1‐3). The protein products of these three genes share significant homology to a 364 amino acid N‐terminal domain of human Bestrophin. We used P‐element mutagenesis to delete dbest1, which encodes a protein with the highest amino acid similarity to Bestrophin. Three independent dbest1 mutants were recovered from the mutagenesis screen. Homozygous null mutations in dbest1 do not significantly alter the viability or fertility of mutant flies. Moreover, dbest1 mutants have normal photoreceptor morphology and function. genesis 31:130–136, 2001. © 2001 Wiley‐Liss, Inc.     

The ANGULATA7 gene encodes a DnaJ‐like zinc finger‐domain protein involved in chloroplast function and leaf development in Arabidopsis

The characterization of mutants with altered leaf shape and pigmentation has previously allowed the identification of nuclear genes that encode plastid‐localized proteins that perform essential functions in leaf growth and development. A large‐scale screen previously allowed us to isolate ethyl methanesulfonate‐induced mutants with small rosettes and pale green leaves with prominent marginal teeth, which were assigned to a phenotypic class that we dubbed Angulata. The molecular characterization of the 12 genes assigned to this phenotypic class should help us to advance our understanding of the still poorly understood relationship between chloroplast biogenesis and leaf morphogenesis. In this article, we report the phenotypic and molecular characterization of the angulata7‐1 (anu7‐1) mutant of Arabidopsis thaliana, which we found to be a hypomorphic allele of the EMB2737 gene, which was previously known only for its embryonic‐lethal mutations. ANU7 encodes a plant‐specific protein that contains a domain similar to the central cysteine‐rich domain of DnaJ proteins. The observed genetic interaction of anu7‐1 with a loss‐of‐function allele of GENOMES UNCOUPLED1 suggests that the anu7‐1 mutation triggers a retrograde signal that leads to changes in the expression of many genes that normally function in the chloroplasts. Many such genes are expressed at higher levels in anu7‐1 rosettes, with a significant overrepresentation of those required for the expression of plastid genome genes. Like in other mutants with altered expression of plastid‐encoded genes, we found that anu7‐1 exhibits defects in the arrangement of thylakoidal membranes, which appear locally unappressed.     

Structural aspects of HDAC8 mechanism and dysfunction in Cornelia de Lange syndrome spectrum disorders

Cornelia de Lange Syndrome (CdLS) encompasses a broad spectrum of phenotypes characterized by distinctive craniofacial abnormalities, limb malformations, growth retardation, and intellectual disability. CdLS spectrum disorders are referred to as cohesinopathies, with ～70% of patients having a mutation in a gene encoding a core cohesin protein (SMC1A, SMC3, or RAD21) or a cohesin regulatory protein (NIPBL or HDAC8). Notably, the regulatory function of HDAC8 in cohesin biology has only recently been discovered. This Zn²⁺‐dependent hydrolase catalyzes the deacetylation of SMC3, a necessary step for cohesin recycling during the cell cycle. To date, 23 different missense mutants in the gene encoding HDAC8 have been identified in children with developmental features that overlap those of CdLS. Enzymological, biophysical, and structural studies of CdLS HDAC8 protein mutants have yielded critical insight on compromised catalysis in vitro. Most CdLS HDAC8 mutations trigger structural changes that directly or indirectly impact substrate binding and catalysis. Additionally, several mutations significantly compromise protein thermostability. Intriguingly, catalytic activity in many HDAC8 mutants can be partially or fully restored by an N‐acylthiourea activator, suggesting a plausible strategy for the chemical rescue of compromised HDAC8 catalysis in vivo.     

Saturating the genetic map of Arabidopsis thaliana with embryonic mutations

One goal of the Arabidopsis genome project is to identify every gene with an essential function in growth and development. Towards that end, the results are reported here of a mapping project designed to enhance the linkage map of Arabidopsis and establish a valuable resource of mutations in essential genes with known map locations. Embryo-defective (emb) mutations were chosen because they represent the most common heritable defect identified following mutagenesis in Arabidopsis. Multiple marker lines with easily scored phenotypes were constructed to facilitate mapping efforts. Recombination data were obtained for 169 mutants defective in embryo-genesis. The chromosomal locations of 110 emb genes are presented in this report. Twenty-six of these genes are tagged with T-DNA. Nine other mutants isolated following seed transformation appear to contain chromosomal translocations. Another 31 mutant genes in the collectiohave been assigned to a linkage group but not yet placed on the map. Nineteen examples of duplicate alleles have also been found. This is consistent with the estimate that approximately 500 genes readily mutate to give an embryo-defective phenotype in Arabidopsis. With continued progress, it may therefore be possible to approach saturation for this important class of mutations. Molecular cloning of these genes should be facilitated by identifying cDNAs and genomic sequences that map to similar locations.     

FUScinating insights into motor neuron degeneration

Point mutations in FUS cause amyotrophic lateral sclerosis (ALS), a devastating neurodegenerative disease—but do they do that by a loss of the protein's normal function, or by endowing it with novel toxic functions, or both? In this issue of The EMBO Journal, Scekic‐Zahirovic et al ( 2016 ) report that mutant FUS, but not the complete loss of FUS, triggers motor neuron degeneration in mice, arguing for a toxic gain‐of‐function mechanism.     

Molecular Population Genetics of Rice Domestication

Domestication is a selection process that genetically modifies species to meet human needs. A most intriguing feature of domestication is the extreme phenotypic diversification among breeds. What could be the ultimate source of such genetic variations? Another notable outcome of artificial selection is the reduction in the fitness of domesticated species when they live in the wild without human assistance. The complete sequences of the two subspecies of rice cultivars provide an opportunity to address these questions. Between the two subspecies, we found much higher rates of non‐synonymous (N) than synonymous (S) substitutions and the N/S ratios are higher between cultivars than between wild species. Most interestingly, substitutions of highly dissimilar amino acids that are deleterious and uncommon between natural species are disproportionately common between the two subspecies of rice. We suggest strong selection in the absence of effective recombination may be the driving force, which we called the domestication‐associated Hill‐Robertson effect. These hitchhiking mutations may contribute to some fitness reduction in cultivars. Comparisons of the two genomes also reveal the existence of highly divergent regions in the genomes. Haplotypes in these regions often form highly polymorphic linkage blocks that are much older than speciation between wild species. Genes from such regions could contribute to the differences between indica and japonica and are likely to be involved in the diversifying selection under domestication. Their existence suggests that the amount of genetic variation within the single progenitor species Oryza rufipogon may be insufficient to account for the variation among rice cultivars, which may come from a more inclusive gene pool comprising most of the A‐genome wild species. Genes from the highly polymorphic regions also provide strong support for the independent domestication of the two subspecies. The genomic variation in rice has revealing implications for studying the genetic basis of indica‐japonica differentiation under rice domestication and subsequent improvement.     

Heterogeneous genetic invasions of three insecticide resistance mutations in Indo‐Pacific populations of Aedes aegypti (L.)

Nations throughout the Indo‐Pacific region use pyrethroid insecticides to control Aedes aegypti, the mosquito vector of dengue, often without knowledge of pyrethroid resistance status of the pest or origin of resistance. Two mutations (V1016G + F1534C) in the sodium channel gene (Vssc) of Ae. aegypti modify ion channel function and cause target‐site resistance to pyrethroid insecticides, with a third mutation (S989P) having a potential additive effect. Of 27 possible genotypes involving these mutations, some allelic combinations are never seen whereas others predominate. Here, five allelic combinations common in Ae. aegypti from the Indo‐Pacific region are described and their geographical distributions investigated using genome‐wide SNP markers. We tested the hypothesis that resistance allele combinations evolved de novo in populations versus the alternative that dispersal of Ae. aegypti between populations facilitated genetic invasions of allele combinations. We used latent factor mixed‐models to detect SNPs throughout the genome that showed structuring in line with resistance allele combinations and compared variation at SNPs within the Vssc gene with genome‐wide variation. Mixed‐models detected an array of SNPs linked to resistance allele combinations, all located within or in close proximity to the Vssc gene. Variation at SNPs within the Vssc gene was structured by resistance profile, whereas genome‐wide SNPs were structured by population. These results demonstrate that alleles near to resistance mutations have been transferred between populations via linked selection. This indicates that genetic invasions have contributed to the widespread occurrence of Vssc allele combinations in Ae. aegypti in the Indo‐Pacific region, pointing to undocumented mosquito invasions between countries.     

Mutations affecting the mitochondrial genes encoding the cytochrome oxidase subunit I and apocytochrome b of Chlamydomonas reinhardtii

Mitochondrial mutants of the green alga Chlamydomonas reinhardtii that are inactivated in the cytochrome pathway of respiration have previously been isolated. Despite the fact that the alternative oxidase pathway is still active the mutants have lost the capacity to grow heterotrophically (dark + acetate) and display reduced growth under mixotrophic conditions (light + acetate). In crosses between wild-type and mutant cells, the meiotic progeny only inherit the character transmitted by the mt ^– parent, which indicates that the mutations are located in the 15.8 kb linear mitochondrial genome. Two new mutants (dum-18 and dum-19) have now been isolated and characterized genetically, biochemically and at the molecular level. In addition, two previously isolated mutants (dum-11 and dum-15) were characterized in more detail. dum-11 contains two types of deleted mitochondrial DNA molecules: 15.1 kb monomers lacking the subterminal part of the genome, downstream of codon 147 of the apocytochrome b (COB) gene, and dimers resulting from head-to-head fusion of asymmetrically deleted monomers (15.1 and 9.5 kb DNA molecules, respectively). As in the wild type, the three other mutants contain only 15.8 kb mitochondrial DNA molecules. dum-15 is mutated at codon 140 of the COB gene, a serine (TCT) being changed into a tyrosine (TAC). dum-18 and dum-19 both inactivate cytochrome c oxidase, as a result of frameshift mutations (addition or deletion of 1 bp) at codons 145 and 152, respectively, of the COX1 gene encoding subunit I of cytochrome c oxidase. In a total of ten respiratory deficient mitochondrial mutants characterized thus far, only mutations located in COB or COXI have been isolated. The possibility that the inactivation of the other mitochondrial genes is lethal for the cells is discussed.     

Strong-polar mutations in the transferase gene of the galactose operon in E.coli

Summary A class of mutations in the transferase gene of the galactose operon in E. coli is described, which is strongly polar for the synthesis of kinase. The latter enzyme is made only to the extent of about 0.1% of the amount made in the induced wildtype. This amount is not dependent on the map position of the mutations and the residual synthesis is non-inducible. The mutants thus resemble 0° mutants in the same operon.Epimerase, which is coded for by the gene proximal to the transferase gene with respect to the operator, is made in normal amounts and its synthesis is normally inducible.The mutants do not seem to belong either to the nonsense or to the frameshift class on the basis of reversion pattern, suppressibility, and degree of polarity. The possible nature of the mutations is discussed.     

Genetic analysis of artificial pigmentation mutants in Porphyra yezoensis Ueda (Bangiales, Rhodophyta)

Porphyra yezoensis Ueda artificial pigmentation mutants, yel (green), fre (red‐orange) and bop (pink), obtained by treatment with /V‐methyl‐/V′‐nitro‐N‐nitrosoguanidine, were genetically analysed. The mutations associated with color phenotypes are recessive because all of the heterozygous conchocelis resembled the wild type color when they were crossed with the wild type (wt). In the reciprocal crosses of yel × wt, both parental colors and eight types of blades appeared in the F₁ gametophytic blades from the heterozygous conchocelis. Both colors segregated in the sectored F₁ blades in a 1:1 ratio, indicating that the color pheno‐type of yel resulted from a single mutation in the nuclear gene. In the reciprocal crosses of fre × wt, however, four colors and more than 40 types of blades appeared in the F₁ blades from the heterozygous conchocelis, indicating that the color phenotype of fre resulted from two mutations in different genes. In the reciprocal crosses of bop×wt, three colors and 12 types of blades were observed in the F₁ blades from the heterozygous conchocelis. Both parental colors appeared far more frequently than the third new color. These results indicated that the color phenotype of bop resulted from two closely linked mutations in different genes, and the epistasis occurred in the F₁ blades. The mutants, yel, fre and bop, differ from the spontaneous green (C‐O), the red (H‐25) and the violet (V‐O) mutants of P. yezoensis, respectively.     

Development of an Agrobacterium‐delivered CRISPR/Cas9 system for wheat genome editing

CRISPR/Cas9 has been widely used for genome editing in many organisms, including important crops like wheat. Despite the tractability in designing CRISPR/Cas9, efficacy in the application of this powerful genome editing tool also depends on DNA delivery methods. In wheat, the biolistics based transformation is the most used method for delivery of the CRISPR/Cas9 complex. Due to the high frequency of gene silencing associated with co‐transferred plasmid backbone and low edit rate in wheat, a large T₀ transgenic plant population are required for recovery of desired mutations, which poses a bottleneck for many genome editing projects. Here, we report an Agrobacterium‐delivered CRISPR/Cas9 system in wheat, which includes a wheat codon optimized Cas9 driven by a maize ubiquitin gene promoter and a guide RNA cassette driven by wheat U6 promoters in a single binary vector. Using this CRISPR/Cas9 system, we have developed 68 edit mutants for four grain‐regulatory genes, TaCKX2‐1, TaGLW7, TaGW2, and TaGW8, in T₀, T₁, and T₂ generation plants at an average edit rate of 10% without detecting off‐target mutations in the most Cas9‐active plants. Homozygous mutations can be recovered from a large population in a single generation. Different from most plant species, deletions over 10 bp are the dominant mutation types in wheat. Plants homozygous of 1160‐bp deletion in TaCKX2‐D1 significantly increased grain number per spikelet. In conclusion, our Agrobacterium‐delivered CRISPR/Cas9 system provides an alternative option for wheat genome editing, which requires a small number of transformation events because CRISPR/Cas9 remains active for novel mutations through generations.