The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture.

The light brown apple moth, Epiphyas postvittana is an invasive, polyphagous pest of horticultural systems around the world. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, where it has been found on over 500 plants, including many horticultural crops. We have produced a genomic resource, to understand the biological basis of the polyphagous and invasive nature of this and other lepidopteran pests. The assembled genome sequence encompassed 598 Mb and has an N50 of 301.17 kb, with a BUSCO completion rate of 97.9%. Epiphyas postvittana has 34% of its assembled genome represented as repetitive sequences, with the majority of the known elements made up of longer DNA transposable elements (14.07 Mb) and retrotransposons (LINE 17.83 Mb). Of the 31,389 predicted genes, 28,714 (91.5%) were assigned to 11,438 orthogroups across the Lepidoptera, of which 945 were specific to E. postvittana. Twenty gene families showed significant expansions in E. postvittana, including some likely to have a role in its pest status, such as cytochrome p450s, glutathione-S-transferases and UDP-glucuronosyltransferases. Finally, using a RAD-tag approach, we investigated the population genomics of this pest, looking at its likely patterns of invasion.

RAPIDOMICS: rapid genome-wide sequencing in a neonatal intensive care unit-successes and challenges.

Genetic disorders are one of the leading causes of infant mortality and are frequent in neonatal intensive care units (NICUs). Rapid genome-wide sequencing (GWS; whole genome or exome sequencing (ES)), due to its diagnostic capabilities and immediate impacts on medical management, is becoming an appealing testing option in the NICU setting. RAPIDOMICS was a trio-based rapid ES pilot study of 25 babies with suspected genetic disorders in the BC Women's Hospital NICU. ES and bioinformatic analysis were performed after careful patient ascertainment. Trio analysis was performed using an in-house pipeline reporting variants in known disease-causing genes. Variants interpreted by the research team as definitely or possibly causal of the infant's phenotype were Sanger validated in a clinical laboratory. The average time to preliminary diagnosis was 7.2 days. Sanger validation was pursued in 15 patients for 13 autosomal dominant and 2 autosomal recessive disorders, with an overall diagnostic rate (partial or complete) of 60%.Conclusion: In total, 72% of patients enrolled had a genomic diagnosis achieved through ES, multi-gene panel testing or chromosomal microarray analysis. Among these, there was an 83% rate of significant and immediate impact on medical decision-making directly related to new knowledge of the diagnosis. Health service implementation challenges and successes are discussed. What is Known: • Rapid genome-wide sequencing in the neonatal intensive care setting has a greater diagnostic hit rate and impact on medical management than conventional genetic testing. However, the impact of consultation with genetics and patient ascertainment requires further investigation. What is New: • This study demonstrates the importance of genetic consultation and careful patient selection prior to pursuing exome sequencing (ES). • In total, 15/25 (60%) patients achieved a diagnosis through ES and 18/25 (72%) through ES, multi-gene panel testing or chromosomal microarray analysis with 83% of those having immediate effects on medical management.

Estimating the biodiversity of terrestrial invertebrates on a forested island using DNA barcodes and metabarcoding data.

Invertebrates are a major component of terrestrial ecosystems, however, estimating their biodiversity is challenging. We compiled an inventory of invertebrate biodiversity along an elevation gradient on the temperate forested island of Hauturu, New Zealand, by DNA barcoding of specimens obtained from leaf litter samples and pitfall traps. We compared the barcodes and biodiversity estimates from this data set with those from a parallel DNA metabarcoding analysis of soil from the same locations, and with pre-existing sequences in reference databases, before exploring the use of combined data sets as a basis for estimating total invertebrate biodiversity. We obtained 1,282 28S and 1,610 COI barcodes from a total of 1,947 invertebrate specimens, which were clustered into 247 (28S) and 366 (COI) OTUs, of which ≤ 10% were represented in GenBank. Coleoptera were most abundant (730 sequenced specimens), followed by Hymenoptera, Diptera, Lepidoptera, and Amphipoda. The most abundant OTU from both the 28S (153 sequences) and COI (140 sequences) data sets was an undescribed beetle from the family Salpingidae. Based on the occurrences of COI OTUs along the elevation gradient, we estimated there are ~1,000 arthropod species (excluding mites) on Hauturu, including 770 insects, of which 344 are beetles. A DNA metabarcoding analysis of soil DNA from the same sites resulted in the identification of similar numbers of OTUs in most invertebrate groups compared with the DNA barcoding, but less than 10% of the DNA barcoding COI OTUs were also detected by the metabarcoding analysis of soil DNA. A mark-recapture analysis based on the overlap between these data sets estimated the presence of approximately 6,800 arthropod species (excluding mites) on the island, including ~3,900 insects. Estimates of New Zealand-wide biodiversity for selected arthropod groups based on matching of the COI DNA barcodes with pre-existing reference sequences suggested over 13,200 insect species are present, including 4,000 Coleoptera, 2,200 Diptera, and 2,700 Hymenoptera species, and 1,000 arachnid species (excluding mites). These results confirm that metabarcoding analyses of soil DNA tends to recover different components of terrestrial invertebrate biodiversity compared to traditional invertebrate sampling, but the combined methods provide a novel basis for estimating invertebrate biodiversity.

Analysis of the genome of the New Zealand giant collembolan (Holacanthella duospinosa) sheds light on hexapod evolution.

BACKGROUND: The New Zealand collembolan genus Holacanthella contains the largest species of springtails (Collembola) in the world. Using Illumina technology we have sequenced and assembled a draft genome and transcriptome from Holacanthella duospinosa (Salmon). We have used this annotated assembly to investigate the genetic basis of a range of traits critical to the evolution of the Hexapoda, the phylogenetic position of H. duospinosa and potential horizontal gene transfer events. RESULTS: Our genome assembly was ~375 Mbp in size with a scaffold N50 of ~230 Kbp and sequencing coverage of ~180×. DNA elements, LTRs and simple repeats and LINEs formed the largest components and SINEs were very rare. Phylogenomics (370,877 amino acids) placed H. duospinosa within the Neanuridae. We recovered orthologs of the conserved sex determination genes thought to play a role in sex determination. Analysis of CpG content suggested the absence of DNA methylation, and consistent with this we were unable to detect orthologs of the DNA methyltransferase enzymes. The small subunit rRNA gene contained a possible retrotransposon. The Hox gene complex was broken over two scaffolds. For chemosensory ability, at least 15 and 18 ionotropic glutamate and gustatory receptors were identified, respectively. However, we were unable to identify any odorant receptors or their obligate co-receptor Orco. Twenty-three chitinase-like genes were identified from the assembly. Members of this multigene family may play roles in the digestion of fungal cell walls, a common food source for these saproxylic organisms. We also detected 59 and 96 genes that blasted to bacteria and fungi, respectively, but were located on scaffolds that otherwise contained arthropod genes. CONCLUSIONS: The genome of H. duospinosa contains some unusual features including a Hox complex broken over two scaffolds, in a different manner to other arthropod species, a lack of odorant receptor genes and an apparent lack of environmentally responsive DNA methylation, unlike many other arthropods. Our detection of candidate horizontal gene transfer candidates confirms that this phenomenon is occurring across Collembola. These findings allow us to narrow down the regions of the arthropod phylogeny where key innovations have occurred that have facilitated the evolutionary success of Hexapoda.

Not just browsing: an animal that grazes phyllosphere microbes facilitates community heterogeneity.

Although grazers have long been recognized as top-down architects of plant communities, animal roles in determining microbial community composition have seldom been examined, particularly in aboveground systems. To determine the extent to which an animal can shape microbial communities, we conducted a controlled mesocosm study in situ to see if introducing mycophageous tree snails changed phyllosphere fungal community composition relative to matched control mesocosms. Fungal community composition and change was determined by Illumina sequencing of DNA collected from leaf surfaces before snails were introduced, daily for 3 days and weekly for 6 weeks thereafter. Scanning electron microscopy was used to confirm that grazing had occurred, and we recorded 3.5 times more cover of fungal hyphae in control mesocosms compared with those containing snails. Snails do not appear to vector novel microbes and despite grazing, a significant proportion of the initial leaf phyllosphere persisted in the mesocosms. Within-mesocosm diversities of fungi were similar regardless of whether or not snails were added. The greatest differences between the snail-treated and control mesocosms was that grazed mesocosms showed greater infiltration of microbes that were not sampled when the experiment commenced and that the variance in fungal community composition (beta diversity) was greater between leaves in snail-treated mesocosms indicating increased community heterogeneity and ecosystem fragmentation.

Evaluating a multigene environmental DNA approach for biodiversity assessment.

BACKGROUND: There is an increasing demand for rapid biodiversity assessment tools that have a broad taxonomic coverage. Here we evaluate a suite of environmental DNA (eDNA) markers coupled with next generation sequencing (NGS) that span the tree of life, comparing them with traditional biodiversity monitoring tools within ten 20×20 meter plots along a 700 meter elevational gradient. RESULTS: From six eDNA datasets (one from each of 16S, 18S, ITS, trnL and two from COI) we identified sequences from 109 NCBI taxonomy-defined phyla or equivalent, ranging from 31 to 60 for a given eDNA marker. Estimates of alpha and gamma diversity were sensitive to the number of sequence reads, whereas beta diversity estimates were less sensitive. The average within-plot beta diversity was lower than between plots for all markers. The soil beta diversity of COI and 18S markers showed the strongest response to the elevational variation of the eDNA markers (COI: r=0.49, p<0.001; 18S: r=0.48, p<0.001). Furthermore pairwise beta diversities for these two markers were strongly correlated with those calculated from traditional vegetation and invertebrate biodiversity measures. CONCLUSIONS: Using a soil-based eDNA approach, we demonstrate that standard phylogenetic markers are capable of recovering sequences from a broad diversity of eukaryotes, in addition to prokaryotes by 16S. The COI and 18S eDNA markers are the best proxies for aboveground biodiversity based on the high correlation between the pairwise beta diversities of these markers and those obtained using traditional methods.

A novel fatty acyl desaturase from the pheromone glands of Ctenopseustis obliquana and C. herana with specific Z5-desaturase activity on myristic acid.

Sexual communication in the Lepidoptera typically involves a female-produced sex pheromone that attracts males of the same species. The most common type of moth sex pheromone comprises individual or blends of fatty acyl derivatives that are synthesized by a specific enzymatic pathway in the female's pheromone gland, often including a desaturation step. This reaction is catalyzed by fatty acyl desaturases that introduce double bonds at specific locations in the fatty acid precursor backbone. The two tortricid moths, Ctenopseustis obliquana and C. herana (brown-headed leafrollers), which are endemic in New Zealand, both use (Z)-5-tetradecenyl acetate as part of their sex pheromone. In C. herana, (Z)-5-tetradecenyl acetate is the sole component of the pheromone. Labeling experiments have revealed that this compound is produced via an unusual Δ5-desaturation of myristic acid. Previously six desaturases were identified from the pheromone glands of Ctenopseustis and its sibling genus Planotortrix, with one differentially regulated to produce the distinct blends used by individual species. However, none were able to conduct the Δ5-desaturation observed in C. herana, and presumably C. obliquana. We have now identified an additional desaturase gene, desat7, expressed in the pheromone glands of both Ctenopseustis species, which is not closely related to any previously described moth pheromone desaturase. The encoded enzyme displays Δ5-desaturase activity on myristic acid when heterologously expressed in yeast, but is not able to desaturate any other fatty acid (C8-C16). We conclude that desat7 represents a new group of desaturases that has evolved a role in the biosynthesis of sex pheromones in moths.

A Mendelian trait for olfactory sensitivity affects odor experience and food selection.

Humans vary in acuity to many odors [1-4], with variation within olfactory receptor (OR) genes contributing to these differences [5-9]. How such variation also affects odor experience and food selection remains uncertain [10], given that such effects occur for taste [11-15]. Here we investigate ß-ionone, which shows extreme sensitivity differences [4, 16, 17]. ß-ionone is a key aroma in foods and beverages [18-21] and is added to products in order to give a pleasant floral note [22, 23]. Genome-wide and in vitro assays demonstrate rs6591536 as the causal variant for ß-ionone odor sensitivity. rs6591536 encodes a N183D substitution in the second extracellular loop of OR5A1 and explains >96% of the observed phenotypic variation, resembling a monogenic Mendelian trait. Individuals carrying genotypes for ß-ionone sensitivity can more easily differentiate between food and beverage stimuli with and without added ß-ionone. Sensitive individuals typically describe ß-ionone in foods and beverages as "fragrant" and "floral," whereas less-sensitive individuals describe these stimuli differently. rs6591536 genotype also influences emotional associations and explains differences in food and product choices. These studies demonstrate that an OR variant that influences olfactory sensitivity can affect how people experience and respond to foods, beverages, and other products.

Identification of regions associated with variation in sensitivity to food-related odors in the human genome.

Humans vary in their ability to smell numerous odors [1-3], including those associated with food [4-6]. Odor sensitivity is heritable [7-11], with examples linking genetic variation for sensitivity to specific odors typically located near olfactory receptor (OR) genes [12-16]. However, with thousands of aromas and few deorphaned ORs [17, 18], there has been little progress toward linking variation at OR loci to odor sensitivity [19, 20]. We hypothesized that OR genes contain the variation that explains much of the differences in sensitivity for odors, paralleling the genetics of taste [21, 22], which affect the flavor experience of foods [23-25]. We employed a genome-wide association approach for ten food-related odors and identified genetic associations to sensitivity for 2-heptanone (p = 5.1 × 10(-8)), isobutyraldehyde (p = 6.4 × 10(-10)), ß-damascenone (p = 1.6 × 10(-7)), and ß-ionone (p = 1.4 × 10(-31)). Each locus is located in/near distinct clusters of OR genes. These findings increase the number of olfactory sensitivity loci to nine and demonstrate the importance of OR-associated variation in sensory acuity for food-related odors. Analysis of genotype frequencies across human populations implies that variation in sensitivity for these odors is widespread. Furthermore, each participant possessed one of many possible combinations of sensitivities for these odors, supporting the notion that everyone experiences their own unique "flavor world."

Sex pheromone evolution is associated with differential regulation of the same desaturase gene in two genera of leafroller moths.

Chemical signals are prevalent in sexual communication systems. Mate recognition has been extensively studied within the Lepidoptera, where the production and recognition of species-specific sex pheromone signals are typically the defining character. While the specific blend of compounds that makes up the sex pheromones of many species has been characterized, the molecular mechanisms underpinning the evolution of pheromone-based mate recognition systems remain largely unknown. We have focused on two sets of sibling species within the leafroller moth genera Ctenopseustis and Planotortrix that have rapidly evolved the use of distinct sex pheromone blends. The compounds within these blends differ almost exclusively in the relative position of double bonds that are introduced by desaturase enzymes. Of the six desaturase orthologs isolated from all four species, functional analyses in yeast and gene expression in pheromone glands implicate three in pheromone biosynthesis, two Δ9-desaturases, and a Δ10-desaturase, while the remaining three desaturases include a Δ6-desaturase, a terminal desaturase, and a non-functional desaturase. Comparative quantitative real-time PCR reveals that the Δ10-desaturase is differentially expressed in the pheromone glands of the two sets of sibling species, consistent with differences in the pheromone blend in both species pairs. In the pheromone glands of species that utilize (Z)-8-tetradecenyl acetate as sex pheromone component (Ctenopseustis obliquana and Planotortrix octo), the expression levels of the Δ10-desaturase are significantly higher than in the pheromone glands of their respective sibling species (C. herana and P. excessana). Our results demonstrate that interspecific sex pheromone differences are associated with differential regulation of the same desaturase gene in two genera of moths. We suggest that differential gene regulation among members of a multigene family may be an important mechanism of molecular innovation in sex pheromone evolution and speciation.

Patterns of mitochondrial haplotype diversity in the invasive pest Epiphyas postvittana (Lepidoptera: Tortricidae).

The light brown apple moth, Epiphyas postvittana (Walker) (Lepidoptera: Tortricidae), is a horticultural pest of Australia and New Zealand that has more recently invaded Hawaii, Europe, and California. A 2,216-bp region of the mitochondrial genome containing the cytochrome oxidase I and II genes was sequenced from 752 individuals. Haplotype network analyses revealed a major split between a predominantly Western Australian clade and all other samples, suggestive of either a deep genetic divergence or a cryptic species. Nucleotide and haplotype diversity were highest in the country of origin, Australia, and in New Zealand populations, with evidence of haplotype sharing between New Zealand and Tasmania. Nucleotide and haplotype diversity were higher in California than within the British Isles or Hawaii. From the total of 96 haplotypes, seven were found in California, of which four were private. Within California, there have been at least two introductions; based on genetic diversity we were unable to assign a likely source for a single moth found and eradicated in Los Angeles in 2007; however, our data suggest it is unlikely that Hawaii and the British Isles are sources of the major E. postvittana population found throughout the rest of the state since 2006.