Organismal duplication, inclusive fitness theory, and altruism: understanding the evolution of endosperm and the angiosperm reproductive syndrome.

For almost a century, events relating to the evolutionary origin of endosperm, a unique embryo-nourishing tissue that is essential to the reproductive process in flowering plants, have remained a mystery. Integration of recent advances in phylogenetic reconstruction, comparative reproductive biology, and genetic theory can be used to elucidate the evolutionary events and forces associated with the establishment of endosperm. Endosperm is shown to be derived from one of two embryos formed during a rudimentary process of "double fertilization" that evolved in the ancestors of angiosperms. Acquisition of embryo-nourishing behavior (with accompanying loss of individual fitness) by this supernumerary fertilization product was dependent upon compensatory gains in the inclusive fitness of related embryos. The result of the loss of individual fitness by one of the two original products of double fertilization was the establishment of endosperm, a highly modified embryo/organism that reproduces cryptically through behavior that enhances the fitness of its associated embryo within a seed. Finally, although triploid endosperm remains a synapomorphy of angiosperms, inclusive fitness analysis demonstrates that the embryo-nourishing properties of endosperm initially evolved in a diploid condition.

Insight into the evolution of the Solanaceae from the parental genomes of Petunia hybrida

Petunia hybrida is a popular bedding plant that has a long history as a genetic model system. We report the whole-genome sequencing and assembly of inbred derivatives of its two wild parents, P. axillaris N and P. inflata S6. The assemblies include 91.3% and 90.2% coverage of their diploid genomes (1.4 Gb; 2n = 14) containing 32,928 and 36,697 protein-coding genes, respectively. The genomes reveal that the Petunia lineage has experienced at least two rounds of hexaploidization: the older gamma event, which is shared with most Eudicots, and a more recent Solanaceae event that is shared with tomato and other solanaceous species. Transcription factors involved in the shift from bee to moth pollination reside in particularly dynamic regions of the genome, which may have been key to the remarkable diversity of floral colour patterns and pollination systems. The high-quality genome sequences will enhance the value of Petunia as a model system for research on unique biological phenomena such as small RNAs, symbiosis, self-incompatibility and circadian rhythms.

Genetic variation for carbon isotope discrimination in sunflower: Association with transpiration efficiency and evidence for cytoplasmic inheritance

Plants accumulate isotopes of carbon at different rates because of discrimination against C-13 relative to C-12. In plants that fix carbon by the C-3 pathway, the amount of discrimination correlates negatively with transpiration efficiency (TE) where TE is the amount of dry matter accumulated per unit water transpired. Therefore, carbon isotope discrimination (Delta) has become a useful tool for selecting genotypes with improved TE and performance in dry environments. Surveys of 161 sunflower (Helianthus spp.) genotypes of diverse origin revealed a large and unprecedented range of genetic variation for Delta (19.5-23.8parts per thousand). A strong negative genetic correlation (r(g)) between TE and Delta (r(g) = -0.87, P < 0.001) was observed in glasshouse studies. Gas exchange measurements of field grown plants indicated that Delta was strongly correlated with stomatal conductance to water vapor (g), (r(g) 0.64, P < 0.01), and the ratio of net assimilation rate (A) to g, (r(g) = 0.86, P < 0.001), an instantaneous measure of TE. Genotype CMSHA89MAX1 had the lowest TE (and highest Delta) of all genotypes tested in these studies and low yields in hybrid combination. Backcrossing studies showed that the TE of this genotype was due to an adverse effect of the MAX1 cytoplasm, which was inherited from the diploid perennial H. maximiliani Schrader. Overall, these studies suggested that there is an excellent opportunity for breeders to develop sunflower germplasm with improved TE. This can be achieved, in part, by avoiding cytoplasms such as the MAX1 cytoplasm.

A genetic linkage map in autotetraploid lucerne adapted to northern Australia, and use of the map to identify DNA markers linked to resistance to Phytophthora medicaginis

Phytophthora root rot, caused by Phytophthora medicaginis, is a major limitation to lucerne ( Medicago sativa L.) production in Australia and North America. Quantitative trait loci (QTLs) involved in resistance to P. medicaginis were identified in a lucerne backcross population of 120 individuals. A genetic linkage map was constructed for tetraploid lucerne using 50 RAPD ( randomly amplified polymorphic DNA), 104 AFLP (amplified fragment length polymorphism) markers, and one SSR ( simple sequence repeat or microsatellite) marker, which originated from the resistant parent (W116); 13 markers remain unlinked. The linkage map contains 18 linkage groups covering 2136.5 cM, with an average distance of 15.0 cM between markers. Four of the linkage groups contained only either 2 or 3 markers. Using duplex markers and repulsion phase linkages the map condensed to 7 homology groups and 2 unassigned linkage groups. Three regions located on linkage groups 2, 14, and 18, were identified as associated with root reaction and the QTLs explained 6 - 15% of the phenotypic variation. The research also indicates that different resistance QTLs are involved in conferring resistance in different organs. Two QTLs were identified as associated with disease resistance expressed after inoculation of detached leaves. The marker, W11-2 on group 18, identified as associated with root reaction, contributed 7% of the phenotypic variation in leaf response in our population. This marker appears to be linked to a QTL encoding a resistance factor contributing to both root and leaf reaction. One other QTL, not identified as associated with root reaction, was positioned on group 1 and contributed to 6% of the variation. This genetic linkage map provides an entry point for future molecular-based improvement of lucerne in Australia, and markers linked to the QTLs we have reported should be useful for marker-assisted selection for partial resistance to P. medicaginis in lucerne.

Meiotic and epigenetic defects in Dnmt3L-knockout mouse spermatogenesis

The production of mature germ cells capable of generating totipotent zygotes is a highly specialized and sexually dimorphic process. The transition from diploid primordial germ cell to haploid spermatozoa requires genome-wide reprogramming of DNA methylation, stage- and testis-specific gene expression, mitotic and meiotic division, and the histone-protamine transition, all requiring unique epigenetic control. Dnmt3L, a DNA methyltransferase regulator, is expressed during gametogenesis, and its deletion results in sterility. We found that during spermatogenesis, Dnmt3L contributes to the acquisition of DNA methylation at paternally imprinted regions, unique nonpericentric heterochromatic sequences, and interspersed repeats, including autonomous transposable elements. We observed retrotransposition of an LTR-ERV1 element in the DNA from Dnmt3L(-/-) germ cells, presumably as a result of hypomethylation. Later in development, in Dnmt3L(-/-) meiotic spermatocytes, we detected abnormalities in the status of biochemical markers of heterochromatin, implying aberrant chromatin packaging. Coincidentally, homologous chromosomes fail to align and form synaptonemal complexes, spermatogenesis arrests, and spermatocytes are lost by apoptosis and sloughing. Because Dnmt3L expression is restricted to gonocytes, the presence of defects in later stages reveals a mechanism whereby early genome reprogramming is linked inextricably to changes in chromatin structure required for completion of spermatogenesis.

DNA markers linked to yield, yield components, and morphological traits in autotetraploid lucerne (Medicago sativa L.)

We have mapped and identifed DNA markers linked to morphology, yield, and yield components of lucerne, using a backcross population derived from winter-active parents. The high-yielding and recurrent parent, D, produced individual markers that accounted for up to 18% of total yield over 6 harvests, at Gatton, south-eastern Queensland. The same marker, AC/TT8, was consistently identified at each individual harvest, and in individual harvests accounted for up to 26% of the phenotypic variation for yield. This marker was located in linkage group 2 of the D map, and several other markers positively associated with yield were consistently identified in this linkage group. Similarly, markers negatively associated with yield were consistently identified in the W116 map, W116 being the low-yielding parent. Highly significant positive correlations were observed between total yield and yield for harvests 1-6, and between total yield and stem length, tiller number, leaf yield/plant, leaf yield/5 stems, stem yield/plant, and stem yield/5 stems. Highly significant QTL were located for all these characters as well as for leaf shape and pubescence.

Production of triploid Kuruma shrimp, Marsupenaeus (Penaeus) japonicus (Bate) nauplii through inhibition of polar body I, or polar body I and II extrusion using 6-dimethylaminopurine

This study investigated the chromosome ploidy level of Marsupenaeus (Penaeus) japonicus (Bate) non-viable (unhatched) embryos and nauplii after exposure to 6-dimethylaminopurine (6-DMAP), timed to stop either polar body (PB) I, or PBI and II extrusion. Embryos from eight separate families or spawnings were exposed to 150 or 200 mu M 6-DMAP from 1- to 3-min post-spawning detection (psd) for a 4- to 5-min duration (timed to stop PBI extrusion). Separate aliquots of embryos from five of the same spawnings were also exposed to 200 mu M of 6-DMAP from 1- to 3-min psd for a 16-min duration (timed to stop both PBI and II extrusion). For one spawning, a third aliquot of embryos was exposed to 400 p M of 6-DMAP from 1- to 3-min psd for a 16-min duration (timed to stop both PBI and II extrusion). At 18-h psd, non-viable embryo and nauplii samples were taken separately for fluorescent activated cell sorting (FACS). FACS revealed that there were diploids and triploids among all treated non-viable embryos and nauplii. All control non-viable embryos and nauplii were diploid. Percentages of triploid induction for the 4- to 5-min and 16-min durations were not significantly different (P > 0.05). Additionally, no difference was found in the triploidy level of nonviable embryos compared to nauplii in these treatments. The percentage of triploid embryos and nauplii when exposed to 6-DMAP for a 4- to 5-min duration ranged from 29.57% to 99.23% (average 55.28 +/- 5.45%) and from 5.60% to 98.85% (average 46.70 +/- 7.20%), respectively. The percentage of triploid embryos and nauplii when exposed to 6-DMAP for a 16-min duration ranged from 11.71% to 98.96% (average 52.49 +/- 11.00%) and from 47.5% to 99.24% (average 79.38 +/- 5.24%), respectively. To our knowledge, this is the first documentation of successful PBI or PBI and II inhibition in shrimp. This study conclusively shows that treatment of M. japonicus embryos with 6-DMAP at 1- to 3-min pscl for either a 4- to 5-min duration (timed to stop PBl extrusion) or 16-min duration (timed to stop both PBI and II extrusion) results in viable triploid nauplii. (c) 2006 Elsevier B.V. All rights reserved.

Preliminary evaluation of the common bean core collection at CIAT

Common bean (Phaseolus vulgaris L.), an annual diploid (2n=2x=22) species (Maréchal 1970; Delgado Salinas 1985), is adapted to mild temperatures (18°C to 35°C) and grown worldwide in a broad range of environments and in diverse production systems. Common bean is grown for its green leaves, green pods, and green and dry seeds. Dry leaves, threshed pods and stalks are fed to animals and used as fuel for cooking, especially in the developing countries of Africa and Asia (Singh 1991).

Genetic structure and male-mediated gene flow in the ghost bat (Macroderma gigas)

The Australian ghost bat is a large, opportunistic carnivorous species that has undergone a marked range contraction toward more mesic, tropical sites over the past century. Comparison of mitochondrial DNA (mtDNA) control region sequences and six nuclear microsatellite loci in 217 ghost bats from nine populations across subtropical and tropical Australia revealed strong population subdivision (mtDNA phi(ST) = 0.80; microsatellites URST = 0.337). Low-latitude (tropical) populations had higher heterozygosity and less marked phylogeographic structure and lower subdivision among sites within regions (within Northern Territory [NT] and within North Queensland [NQ]) than did populations at higher latitudes (subtropical sites; central Queensland [CQ]), although sampling of geographically proximal breeding sites is unavoidably restricted for the latter. Gene flow among populations within each of the northern regions appears to be male biased in that the difference in population subdivision for mtDNA and microsatellites (NT phi(ST) = 0.39, URST = 0.02; NQ phi(ST) = 0.60, URST = -0.03) is greater than expected from differences in the effective population size of haploid versus diploid loci. The high level of population subdivision across the range of the ghost bat contrasts with evidence for high gene flow in other chiropteran species and may be due to narrow physiological tolerances and consequent limited availability of roosts for ghost bats, particularly across the subtropical and relatively arid regions. This observation is consistent with the hypothesis that the contraction of the species' range is associated with late Holocene climate change. The extreme isolation among higher-latitude populations may predispose them to additional local extinctions if the processes responsible for the range contraction continue to operate.

Modelling the dynamics of genetic algorithms using statistical mechanics

A formalism for modelling the dynamics of Genetic Algorithms (GAs) using methods from statistical mechanics, originally due to Prugel-Bennett and Shapiro, is reviewed, generalized and improved upon. This formalism can be used to predict the averaged trajectory of macroscopic statistics describing the GA's population. These macroscopics are chosen to average well between runs, so that fluctuations from mean behaviour can often be neglected. Where necessary, non-trivial terms are determined by assuming maximum entropy with constraints on known macroscopics. Problems of realistic size are described in compact form and finite population effects are included, often proving to be of fundamental importance. The macroscopics used here are cumulants of an appropriate quantity within the population and the mean correlation (Hamming distance) within the population. Including the correlation as an explicit macroscopic provides a significant improvement over the original formulation. The formalism is applied to a number of simple optimization problems in order to determine its predictive power and to gain insight into GA dynamics. Problems which are most amenable to analysis come from the class where alleles within the genotype contribute additively to the phenotype. This class can be treated with some generality, including problems with inhomogeneous contributions from each site, non-linear or noisy fitness measures, simple diploid representations and temporally varying fitness. The results can also be applied to a simple learning problem, generalization in a binary perceptron, and a limit is identified for which the optimal training batch size can be determined for this problem. The theory is compared to averaged results from a real GA in each case, showing excellent agreement if the maximum entropy principle holds. Some situations where this approximation brakes down are identified. In order to fully test the formalism, an attempt is made on the strong sc np-hard problem of storing random patterns in a binary perceptron. Here, the relationship between the genotype and phenotype (training error) is strongly non-linear. Mutation is modelled under the assumption that perceptron configurations are typical of perceptrons with a given training error. Unfortunately, this assumption does not provide a good approximation in general. It is conjectured that perceptron configurations would have to be constrained by other statistics in order to accurately model mutation for this problem. Issues arising from this study are discussed in conclusion and some possible areas of further research are outlined.

Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.

Genetic Analysis of Mitotic Recombination in Saccharomyces cerevisiae

Mitotic genome instability can occur during the repair of double-strand breaks (DSBs) in DNA, which arise from endogenous and exogenous sources. Studying the mechanisms of DNA repair in the budding yeast, Saccharomyces cerevisiae has shown that Homologous Recombination (HR) is a vital repair mechanism for DSBs. HR can result in a crossover event, in which the broken molecule reciprocally exchanges information with a homologous repair template. The current model of double-strand break repair (DSBR) also allows for a tract of information to non-reciprocally transfer from the template molecule to the broken molecule. These “gene conversion” events can vary in size and can occur in conjunction with a crossover event or in isolation. The frequency and size of gene conversions in isolation and gene conversions associated with crossing over has been a source of debate due to the variation in systems used to detect gene conversions and the context in which the gene conversions are measured.

In Chapter 2, I use an unbiased system that measures the frequency and size of gene conversion events, as well as the association of gene conversion events with crossing over between homologs in diploid yeast. We show mitotic gene conversions occur at a rate of 1.3x10-6 per cell division, are either large (median 54.0kb) or small (median 6.4kb), and are associated with crossing over 43% of the time.

DSBs can arise from endogenous cellular processes such as replication and transcription. Two important RNA/DNA hybrids are involved in replication and transcription: R-loops, which form when an RNA transcript base pairs with the DNA template and displaces the non-template DNA strand, and ribonucleotides embedded into DNA (rNMPs), which arise when replicative polymerase errors insert ribonucleotide instead of deoxyribonucleotide triphosphates. RNaseH1 (encoded by RNH1) and RNaseH2 (whose catalytic subunit is encoded by RNH201) both recognize and degrade the RNA in within R-loops while RNaseH2 alone recognizes, nicks, and initiates removal of rNMPs embedded into DNA. Due to their redundant abilities to act on RNA:DNA hybrids, aberrant removal of rNMPs from DNA has been thought to lead to genome instability in an rnh201Δ background.

In Chapter 3, I characterize (1) non-selective genome-wide homologous recombination events and (2) crossing over on chromosome IV in mutants defective in RNaseH1, RNaseH2, or RNaseH1 and RNaseH2. Using a mutant DNA polymerase that incorporates 4-fold fewer rNMPs than wild type, I demonstrate that the primary recombinogenic lesion in the RNaseH2-defective genome is not rNMPs, but rather R-loops. This work suggests different in-vivo roles for RNaseH1 and RNaseH2 in resolving R-loops in yeast and is consistent with R-loops, not rNMPs, being the the likely source of pathology in Aicardi-Goutières Syndrome patients defective in RNaseH2.

Aneuploidy Tolerance in a Polyploid Organ

Endopolyploid cells (hereafter - polyploid cells), which contain whole genome duplications in an otherwise diploid organism, play vital roles in development and physiology of diverse organs such as our heart and liver. Polyploidy is also observed with high frequency in many tumors, and division of such cells frequently creates aneuploidy (chromosomal imbalances), a hallmark of cancer. Despite its frequent occurrence and association with aneuploidy, little is known about the specific role that polyploidy plays in diverse contexts. Using a new model tissue, the Drosophila rectal papilla, we sought to uncover connections between polyploidy and aneuploidy during organ development. Our lab previously discovered that the papillar cells of the Drosophila hindgut undergo developmentally programmed polyploid cell divisions, and that these polyploid cell divisions are highly error-prone. Time-lapse studies of polyploid mitosis revealed that the papillar cells undergo a high percentage of tripolar anaphase, which causes extreme aneuploidy. Despite this massive chromosome imbalance, we found the tripolar daughter cells are viable and support normal organ development and function, suggesting acquiring extra genome sets enables a cell to tolerate the genomic alterations incurred by aneuploidy. We further extended these findings by seeking mechanisms by which the papillar cells tolerated this resultant aneuploidy.

Physical and Genetic Analysis of the CUP1 Tandem Array in the Yeast Saccharomyces cerevisiae

The genomes of many strains of baker’s yeast, Saccharomyces cerevisiae, contain multiple repeats of the copper-binding protein Cup1. Cup1 is a member of the metallothionein family, and is found in a tandem array on chromosome VIII. In this thesis, I describe studies that characterized these tandem arrays and their mechanism of formation across diverse strains of yeast. I show that CUP1 arrays are an illuminating model system for observing recombination in eukaryotes, and describe insights derived from these observations.

In our first study, we analyzed 101 natural isolates of S. cerevisiae in order to examine the diversity of CUP1-containing repeats across different strains. We identified five distinct classes of repeats that contain CUP1. We also showed that some strains have only a single copy of CUP1. By comparing the sequences of all the strains, we were able to elucidate the mechanism of formation of the CUP1 tandem arrays, which involved unequal non-homologous recombination events starting from a strain that had only a single CUP1 gene. Our observation of CUP1 repeat formation allows more general insights about the formation of tandem repeats from single-copy genes in eukaryotes, which is one of the most important mechanisms by which organisms evolve.

In our second study, we delved deeper into our mechanistic investigations by measuring the relative rates of inter-homolog and intra-/inter-sister chromatid recombination in CUP1 tandem arrays. We used a diploid strain that is heterozygous both for insertion of a selectable marker (URA3) inside the tandem array, and also for markers at either end of the array. The intra-/inter-sister chromatid recombination rate turned out to be more than ten-fold greater than the inter-homolog rate. Moreover, we found that loss of the proteins Rad51 and Rad52, which are required for most inter-homolog recombination, did not greatly reduce recombination in the CUP1 tandem repeats. Additionally, we investigated the effects of elevated copper levels on the rate of each type of recombination at the CUP1 locus. Both types of recombination are increased at high concentrations of copper (as is known to be the case for CUP1 transcription). Furthermore, the inter-homolog recombination rate at the CUP1 locus is higher than the average over the genome during mitosis, but is lower than the average during meiosis.

The research described in Chapter 2 is published in 2014.

Seawater carbonate chemistry and effects of CO2 and their modulation by light in the life-cycle stages of the coccolithophore Emiliania huxleyi strains RCC 1216 and 1217 during experiments, 2012

The effects of ocean acidification on the life-cycle stages of the coccolithophore Emiliania huxleyi and their by light were examined. Calcifying diploid and noncalcifying haploid cells (Roscoff culture collection 1216 and 1217) were acclimated to present-day and elevated CO2 partial pressures (PCO2; 38.5 vs. 101.3 Pa, ., 380 vs. 1000 matm) under low and high light (50 vs. 300 mmol photons m-2 s-1). Growth rates as well as quotas and production rates of C and N were measured. Sources of inorganic C for biomass buildup were using a 14C disequilibrium assay. Photosynthetic O2 evolution was measured as a function of dissolved inorganic C and light by means of membrane-inlet mass spectrometry. The diploid stage responded to elevated PCO2 by shunting resources from the production of particulate inorganic C toward organic C yet keeping the production of total particulate C constant. As the effect of ocean acidification was stronger under low light, the diploid stage might be less affected by increased acidity when energy availability is high. The haploid stage maintained elemental composition and production rates under elevated PCO2. Although both life-cycle stages involve different ways of dealing with elevated PCO2, the responses were generally modulated by energy availability, being typically most pronounced under low light. Additionally, PCO2 responses resembled those induced by high irradiances, indicating that ocean acidification affects the interplay between energy-generating processes (photosynthetic light reactions) and processes competing for energy (biomass buildup and calcification). A conceptual model is put forward explaining why the magnitude of single responses is determined by energy availability.