Unprecedented reorganization of holocentric chromosomes provides insights into the enigma of lepidopteran chromosome evolution

Chromosome evolution presents an enigma in the mega-diverse Lepidoptera. Most species exhibit constrained chromosome evolution with nearly identical haploid chromosome counts and chromosome-level gene collinearity among species more than 140 million years divergent. However, a few species possess radically inflated chromosomal counts due to extensive fission and fusion events. To address this enigma of constraint in the face of an exceptional ability to change, we investigated an unprecedented reorganization of the standard lepidopteran chromosome structure in the green-veined white butterfly (Pieris napi). We find that gene content in P. napi has been extensively rearranged in large collinear blocks, which until now have been masked by a haploid chromosome number close to the lepidopteran average. We observe that ancient chromosome ends have been maintained and collinear blocks are enriched for functionally related genes suggesting both a mechanism and a possible role for selection in determining the boundaries of these genome-wide rearrangements.Peer reviewe

Advances in studying the role of genetic divergence and recombination in adaptation in non-model species

Understanding the role of genetic divergence and recombination in adaptation is crucial to understanding the evolutionary potential of species since they can directly affect the levels of genetic variation present within populations or species. Genetic variation in the functional parts of the genome such as exons or regulatory regions is the raw material for evolution, because natural selection can only operate on phenotypic variation already present in the population. When natural selection acts on a phenotype, it usually results in reduction in the levels of genetic variation at the causal loci, and the surrounding linked loci, due to recombination dynamics (i.e. linkage); the degree to which natural selection influences the genetic differentiation in the linked regions depends on the local recombination rates. Studies investigating the role of genetic divergence and recombination are common in model species such as Drosophila melanogaster. Only recently have genomic tools allowed us to start investigating their role in shaping genetic variation in non-model species. This thesis adds to the growing research in that domain. In this thesis, I have asked a diverse set of questions to understand the role of genetic divergence and recombination in adaptation in non-model species, with a focus on Lepidoptera. First, how do we identify causal genetic variation causing adaptive phenotypes? This question is fundamental to evolutionary biology and addressing it requires a well-assembled genome, the generation of which is a cost, labor, and time intensive task. In paper I, I present a tool, MESPA, that stitches together exonic sequences in fragmented assemblies to produce high-quality gene models. These high-quality gene models can be used by researchers in the downstream analyses, providing genomic insights for a fraction of cost of a high quality genome.  Second, what does the pattern of recombination rate look like in chromosomes that lack centromeres (i.e.holocentric chromosomes)? In paper II, I compare the recombination landscape and the patterns of nucleotide diversity in three Lepidotera with holocentric chromosomes, Pieris napi, Bombyx mandarina, and Bombyx mori, with a monocentric species. Our results show that on average these three Lepidoptera have high rates of recombination across the vast majority of their genome. Our results also suggest that given similar effective population sizes, these species are likely to harbor more genetic diversity compared to monocentric species, which has important evolutionary consequences for these species. Third, what is the potential for parallelism at the genetic level in convergent melanic phenotypes? In paper III, I investigated the genetic basis of the female-limited melanic phenotype in the green-veined white (Pieris napi) butterfly, and found a 20kb region, approximately 50kb from the gene cortex, associated with this trait. This gene has been implicated in melanic phenotypes in other Lepidoptera that diverged from Pieris approximately 100my, indicating very high predictability for this trait. Finally, what is the role of cis-regulatory variation in local adaptation? In paper IV, I analyzed the relationship between allele specific expression (ASE) and genetic divergence (FST) in the F1 hybrids of Pieris napi napi and Pieris napi adalwinda. I show that intersecting results from ASE with FST is a powerful approach to identify genes involved in local adaptation.At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 2: Manuscript. Paper 3: Manuscript. Paper 4: Manuscript.</p

Advances in studying the role of genetic divergence and recombination in adaptation in non-model species

Understanding the role of genetic divergence and recombination in adaptation is crucial to understanding the evolutionary potential of species since they can directly affect the levels of genetic variation present within populations or species. Genetic variation in the functional parts of the genome such as exons or regulatory regions is the raw material for evolution, because natural selection can only operate on phenotypic variation already present in the population. When natural selection acts on a phenotype, it usually results in reduction in the levels of genetic variation at the causal loci, and the surrounding linked loci, due to recombination dynamics (i.e. linkage); the degree to which natural selection influences the genetic differentiation in the linked regions depends on the local recombination rates. Studies investigating the role of genetic divergence and recombination are common in model species such as Drosophila melanogaster. Only recently have genomic tools allowed us to start investigating their role in shaping genetic variation in non-model species. This thesis adds to the growing research in that domain. In this thesis, I have asked a diverse set of questions to understand the role of genetic divergence and recombination in adaptation in non-model species, with a focus on Lepidoptera. First, how do we identify causal genetic variation causing adaptive phenotypes? This question is fundamental to evolutionary biology and addressing it requires a well-assembled genome, the generation of which is a cost, labor, and time intensive task. In paper I, I present a tool, MESPA, that stitches together exonic sequences in fragmented assemblies to produce high-quality gene models. These high-quality gene models can be used by researchers in the downstream analyses, providing genomic insights for a fraction of cost of a high quality genome.  Second, what does the pattern of recombination rate look like in chromosomes that lack centromeres (i.e.holocentric chromosomes)? In paper II, I compare the recombination landscape and the patterns of nucleotide diversity in three Lepidotera with holocentric chromosomes, Pieris napi, Bombyx mandarina, and Bombyx mori, with a monocentric species. Our results show that on average these three Lepidoptera have high rates of recombination across the vast majority of their genome. Our results also suggest that given similar effective population sizes, these species are likely to harbor more genetic diversity compared to monocentric species, which has important evolutionary consequences for these species. Third, what is the potential for parallelism at the genetic level in convergent melanic phenotypes? In paper III, I investigated the genetic basis of the female-limited melanic phenotype in the green-veined white (Pieris napi) butterfly, and found a 20kb region, approximately 50kb from the gene cortex, associated with this trait. This gene has been implicated in melanic phenotypes in other Lepidoptera that diverged from Pieris approximately 100my, indicating very high predictability for this trait. Finally, what is the role of cis-regulatory variation in local adaptation? In paper IV, I analyzed the relationship between allele specific expression (ASE) and genetic divergence (FST) in the F1 hybrids of Pieris napi napi and Pieris napi adalwinda. I show that intersecting results from ASE with FST is a powerful approach to identify genes involved in local adaptation.At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 2: Manuscript. Paper 3: Manuscript. Paper 4: Manuscript.</p

Scripts for analyzing Papilio and Corvus data

Contains the commands used for the <i>Papilio</i> and <i>Corvus</i> analyses including the R commands that was used for generating the plots

Single copy orthologs in the Drosophila validation work

<p>1) dmel_70_100_80cov.fa – Contains the single copy orthologs from <i>Drosophila melanogaster</i> that had atleast 80% coverage and 70% identity in a reciprocal best blast search with single copy orthologs from <i>Drosophila virilis</i>.</p> <p>2) dviri_70_100_80cov.fa – Contains the single copy orthologs from <i>Drosophila</i> <i>virilis </i>that had atleast 80% coverage and 70% identity in a reciprocal best blast search with single copy orthologs from <i>Drosophila melanogaster</i>.</p

A transposable element insertion is associated with an alternative life history strategy

Tradeoffs affect resource allocation during development and result in fitness consequences that drive the evolution of life history strategies. Yet despite their importance, we know little about the mechanisms underlying life history tradeoffs. Many species of Colias butterflies exhibit an alternative life history strategy (ALHS) where females divert resources from wing pigment synthesis to reproductive and somatic development. Due to this reallocation, a wing color polymorphism is associated with the ALHS: either yellow/orange or white. Here we map the locus associated with this ALHS in Colias crocea to a transposable element insertion located downstream of the Colias homolog of BarH-1, a homeobox transcription factor. Using CRISPR/Cas9 gene editing, antibody staining, and electron microscopy we find white-specific expression of BarH-1 suppresses the formation of pigment granules in wing scales and gives rise to white wing color. Lipid and transcriptome analyses reveal physiological differences associated with the ALHS. Together, these findings characterize a mechanism for a female-limited ALHS.Peer reviewe

Evolutionary history of host use, rather than plant phylogeny, determines gene expression in a generalist butterfly

Background: Although most insect species are specialized on one or few groups of plants, there are phytophagous insects that seem to use virtually any kind of plant as food. Understanding the nature of this ability to feed on a wide repertoire of plants is crucial for the control of pest species and for the elucidation of the macroevolutionary mechanisms of speciation and diversification of insect herbivores. Here we studied Vanessa cardui, the species with the widest diet breadth among butterflies and a potential insect pest, by comparing tissue-specific transcriptomes from caterpillars that were reared on different host plants. We tested whether the similarities of gene-expression response reflect the evolutionary history of adaptation to these plants in the Vanessa and related genera, against the null hypothesis of transcriptional profiles reflecting plant phylogenetic relatedness. Result: Using both unsupervised and supervised methods of data analysis, we found that the tissue-specific patterns of caterpillar gene expression are better explained by the evolutionary history of adaptation of the insects to the plants than by plant phylogeny. Conclusion: Our findings suggest that V. cardui may use two sets of expressed genes to achieve polyphagy, one associated with the ancestral capability to consume Rosids and Asterids, and another allowing the caterpillar to incorporate a wide range of novel host-plants

COSIM - Coating Simulation

In the given coating chamber is a system of rotating holders, which holds rotating tools. The chamber also contains electrodes (of various materials), from which the material evaporates and is coating the tools. The goal of this thesis is to simulate the process of coating and visualize the resulting layer in given places of coated tools depending on their placement, number, rotation speed and other optional settings

A transposable element insertion is associated with an&nbsp;alternative life history strategy.

Tradeoffs affect resource allocation during development and result in fitness consequences that drive the evolution of life history strategies. Yet despite their importance, we know little about the mechanisms underlying life history tradeoffs. Many species of Colias butterflies exhibit an alternative life history strategy (ALHS) where females divert resources from wing pigment synthesis to reproductive and somatic development. Due to this reallocation, a wing color polymorphism is associated with the ALHS: either yellow/orange or white. Here we map the locus associated with this ALHS in Colias crocea to a transposable element insertion located downstream of the Colias homolog of BarH-1, a homeobox transcription factor. Using CRISPR/Cas9 gene editing, antibody staining, and electron microscopy we find white-specific expression of BarH-1 suppresses the formation of pigment granules in wing scales and gives rise to white wing color. Lipid and transcriptome analyses reveal physiological differences associated with the ALHS. Together, these findings characterize a mechanism for a female-limited ALHS

<i>De Novo</i> Transcriptome Assembly (NGS) of <i>Curcuma longa</i> L. Rhizome Reveals Novel Transcripts Related to Anticancer and Antimalarial Terpenoids

<div><p>Herbal remedies are increasingly being recognised in recent years as alternative medicine for a number of diseases including cancer. <i>Curcuma longa</i> L., commonly known as turmeric is used as a culinary spice in India and in many Asian countries has been attributed to lower incidences of gastrointestinal cancers. Curcumin, a secondary metabolite isolated from the rhizomes of this plant has been shown to have significant anticancer properties, in addition to antimalarial and antioxidant effects. We sequenced the transcriptome of the rhizome of the 3 varieties of <i>Curcuma longa</i> L. using Illumina reversible dye terminator sequencing followed by <i>de novo</i> transcriptome assembly. Multiple databases were used to obtain a comprehensive annotation and the transcripts were functionally classified using GO, KOG and PlantCyc. Special emphasis was given for annotating the secondary metabolite pathways and terpenoid biosynthesis pathways. We report for the first time, the presence of transcripts related to biosynthetic pathways of several anti-cancer compounds like taxol, curcumin, and vinblastine in addition to anti-malarial compounds like artemisinin and acridone alkaloids, emphasizing turmeric's importance as a highly potent phytochemical. Our data not only provides molecular signatures for several terpenoids but also a comprehensive molecular resource for facilitating deeper insights into the transcriptome of <i>C. longa</i>.</p></div