Improved draft reference genome for the Glassy-winged Sharpshooter (Homalodisca vitripennis), a vector for Pierce's disease.

Homalodisca vitripennis (Hemiptera: Cicadellidae), known as the glassy-winged sharpshooter, is a xylem feeding leafhopper and an important agricultural pest as a vector of Xylella fastidiosa, which causes Pierce's disease in grapes and a variety of other scorch diseases. The current H. vitripennis reference genome from the Baylor College of Medicine's i5k pilot project is a 1.4-Gb assembly with 110,000 scaffolds, which still has significant gaps making identification of genes difficult. To improve on this effort, we used a combination of Oxford Nanopore long-read sequencing technology combined with Illumina sequencing reads to generate a better assembly and first-pass annotation of the whole genome sequence of a wild-caught Californian (Tulare County) individual of H. vitripennis. The improved reference genome assembly for H. vitripennis is 1.93-Gb in length (21,254 scaffolds, N50 = 650 Mb, BUSCO completeness = 94.3%), with 33.06% of the genome masked as repetitive. In total, 108,762 gene models were predicted including 98,296 protein-coding genes and 10,466 tRNA genes. As an additional community resource, we identified 27 orthologous candidate genes of interest for future experimental work including phenotypic marker genes like white. Furthermore, as part of the assembly process, we generated four endosymbiont metagenome-assembled genomes, including a high-quality near complete 1.7-Mb Wolbachia sp. genome (1 scaffold, CheckM completeness = 99.4%). The improved genome assembly and annotation for H. vitripennis, curated set of candidate genes, and endosymbiont MAGs will be invaluable resources for future research of H. vitripennis.

Genomic and anatomical comparisons of skin support independent adaptation to life in water by cetaceans and hippos.

The macroevolutionary transition from terra firma to obligatory inhabitance of the marine hydrosphere has occurred twice in the history of Mammalia: Cetacea and Sirenia. In the case of Cetacea (whales, dolphins, and porpoises), molecular phylogenies provide unambiguous evidence that fully aquatic cetaceans and semiaquatic hippopotamids (hippos) are each other's closest living relatives. Ancestral reconstructions suggest that some adaptations to the aquatic realm evolved in the common ancestor of Cetancodonta (Cetacea + Hippopotamidae). An alternative hypothesis is that these adaptations evolved independently in cetaceans and hippos. Here, we focus on the integumentary system and evaluate these hypotheses by integrating new histological data for cetaceans and hippos, the first genome-scale data for pygmy hippopotamus, and comprehensive genomic screens and molecular evolutionary analyses for protein-coding genes that have been inactivated in hippos and cetaceans. We identified eight skin-related genes that are inactivated in both cetaceans and hippos, including genes that are related to sebaceous glands, hair follicles, and epidermal differentiation. However, none of these genes exhibit inactivating mutations that are shared by cetaceans and hippos. Mean dates for the inactivation of skin genes in these two clades serve as proxies for phenotypic changes and suggest that hair reduction/loss, the loss of sebaceous glands, and changes to the keratinization program occurred ∼16 Ma earlier in cetaceans (∼46.5 Ma) than in hippos (∼30.5 Ma). These results, together with histological differences in the integument and prior analyses of oxygen isotopes from stem hippopotamids ("anthracotheres"), support the hypothesis that aquatic skin adaptations evolved independently in hippos and cetaceans.

Odontogenic ameloblast-associated (ODAM) is inactivated in toothless/enamelless placental mammals and toothed whales.

BACKGROUND: The gene for odontogenic ameloblast-associated (ODAM) is a member of the secretory calcium-binding phosphoprotein gene family. ODAM is primarily expressed in dental tissues including the enamel organ and the junctional epithelium, and may also have pleiotropic functions that are unrelated to teeth. Here, we leverage the power of natural selection to test competing hypotheses that ODAM is tooth-specific versus pleiotropic. Specifically, we compiled and screened complete protein-coding sequences, plus sequences for flanking intronic regions, for ODAM in 165 placental mammals to determine if this gene contains inactivating mutations in lineages that either lack teeth (baleen whales, pangolins, anteaters) or lack enamel on their teeth (aardvarks, sloths, armadillos), as would be expected if the only essential functions of ODAM are related to tooth development and the adhesion of the gingival junctional epithelium to the enamel tooth surface. RESULTS: We discovered inactivating mutations in all species of placental mammals that either lack teeth or lack enamel on their teeth. A surprising result is that ODAM is also inactivated in a few additional lineages including all toothed whales that were examined. We hypothesize that ODAM inactivation is related to the simplified outer enamel surface of toothed whales. An alternate hypothesis is that ODAM inactivation in toothed whales may be related to altered antimicrobial functions of the junctional epithelium in aquatic habitats. Selection analyses on ODAM sequences revealed that the composite dN/dS value for pseudogenic branches is close to 1.0 as expected for a neutrally evolving pseudogene. DN/dS values on transitional branches were used to estimate ODAM inactivation times. In the case of pangolins, ODAM was inactivated ~ 65 million years ago, which is older than the oldest pangolin fossil (Eomanis, 47 Ma) and suggests an even more ancient loss or simplification of teeth in this lineage. CONCLUSION: Our results validate the hypothesis that the only essential functions of ODAM that are maintained by natural selection are related to tooth development and/or the maintenance of a healthy junctional epithelium that attaches to the enamel surface of teeth.

Genomic perspectives of spider silk genes through target capture sequencing: Conservation of stabilization mechanisms and homology-based structural models of spidroin terminal regions.

A powerful system for studying protein aggregation, particularly rapid self-assembly, is spider silk. Spider silks are proteinaceous and silk proteins are synthesized and stored within silk glands as liquid dope. As needed, liquid dope is near-instantaneously transformed into solid fibers or viscous adhesives. The dominant constituents of silks are spidroins (spider fibroins) and their terminal domains are vital for the tight control of silk self-assembly. To better understand spidroin termini, we used target capture and deep sequencing to identify spidroin gene sequences from six species representing the araneoid families of Araneidae, Nephilidae, and Theridiidae. We obtained 145 terminal regions, of which 103 are newly annotated here, as well as novel variants within nine diverse spidroin types. Our comparative analyses demonstrated the conservation of acidic, basic, and cysteine amino acid residues across spidroin types that had been proposed to be important for monomer stability, dimer formation, and self-assembly from a limited sampling of spidroins. Computational, protein homology modeling revealed areas of spidroin terminal regions that are highly conserved in three-dimensions despite sequence divergence across spidroin types. Analyses of our dense sampling of terminal regions suggest that most spidroins share stabilization mechanisms, dimer formation, and tertiary structure, despite producing functionally distinct materials.

Duplication and concerted evolution of MiSp-encoding genes underlie the material properties of minor ampullate silks of cobweb weaving spiders.

BACKGROUND: Orb-web weaving spiders and their relatives use multiple types of task-specific silks. The majority of spider silk studies have focused on the ultra-tough dragline silk synthesized in major ampullate glands, but other silk types have impressive material properties. For instance, minor ampullate silks of orb-web weaving spiders are as tough as draglines, due to their higher extensibility despite lower strength. Differences in material properties between silk types result from differences in their component proteins, particularly members of the spidroin (spider fibroin) gene family. However, the extent to which variation in material properties within a single silk type can be explained by variation in spidroin sequences is unknown. Here, we compare the minor ampullate spidroins (MiSp) of orb-weavers and cobweb weavers. Orb-web weavers use minor ampullate silk to form the auxiliary spiral of the orb-web while cobweb weavers use it to wrap prey, suggesting that selection pressures on minor ampullate spidroins (MiSp) may differ between the two groups. RESULTS: We report complete or nearly complete MiSp sequences from five cobweb weaving spider species and measure material properties of minor ampullate silks in a subset of these species. We also compare MiSp sequences and silk properties of our cobweb weavers to published data for orb-web weavers. We demonstrate that all our cobweb weavers possess multiple MiSp loci and that one locus is more highly expressed in at least two species. We also find that the proportion of ß-spiral-forming amino acid motifs in MiSp positively correlates with minor ampullate silk extensibility across orb-web and cobweb weavers. CONCLUSIONS: MiSp sequences vary dramatically within and among spider species, and have likely been subject to multiple rounds of gene duplication and concerted evolution, which have contributed to the diverse material properties of minor ampullate silks. Our sequences also provide templates for recombinant silk proteins with tailored properties.

Evidence from Multiple Species that Spider Silk Glue Component ASG2 is a Spidroin.

Spiders in the superfamily Araneoidea produce viscous glue from aggregate silk glands. Aggregate glue coats prey-capture threads and hampers the escape of prey from webs, thereby increasing the foraging success of spiders. cDNAs for Aggregate Spider Glue 1 (ASG1) and 2 (ASG2) have been previously described from the golden orb-weaver, Nephila clavipes, and Western black widow, Latrodectus hesperus. To further investigate aggregate glues, we assembled ASG1 and ASG2 from genomic target capture libraries constructed from three species of cob-web weavers and three species of orb-web weavers, all araneoids. We show that ASG1 is unlikely to be a glue, but rather is part of a widespread arthropod gene family, the peritrophic matrix proteins. For ASG2, we demonstrate its remarkable architectural and sequence similarities to spider silk fibroins, indicating that ASG2 is a member of the spidroin gene family. Thus, spidroins have diversified into glues in addition to task-specific, high performance fibers.

Comparison of fibroin cDNAs from webspinning insects: insight into silk formation and function.

Embiopterans (webspinning insects) are renowned for their prolific use of silk. These organisms spin silk to construct elaborate networks of tubes in which they live, forage, and reproduce. The silken galleries are essential for protecting these soft-bodied insects from predators and other environmental hazards. Despite the ecological importance of embiopteran silk, very little is known about its constituent proteins. Here, we characterize the silk protein cDNAs from four embiopteran species to better understand the function and evolution of these adaptive molecules. We show that webspinner fibroins (silk proteins) are highly repetitive in sequence and possess several conserved characteristics, despite differences in habitat preferences across species. The most striking similarities are in the codon usage biases of the fibroin genes, particularly in the repetitive regions, as well as sequence conservation of the carboxyl-terminal regions of the fibroins. Based on analyses of the silk genes, we propose hypotheses regarding codon bias and its effect on the translation and replication of these unusual genes. Furthermore, we discuss the significance of specific fibroin motifs to the mechanical and structural characteristics of silk fibers. Lastly, we report that the conservation of webspinner fibroin carboxyl-terminal regions suggests that fiber formation may occur through a mechanism analogous to that found in Lepidoptera. From these results, insight is gained into the tempo and mode of evolution that has shaped embiopteran fibroins.

Molecular evolution of lepidopteran silk proteins: insights from the ghost moth, Hepialus californicus.

Silk production has independently evolved in numerous arthropod lineages, such as Lepidoptera, the moths and butterflies. Lepidopteran larvae (caterpillars) synthesize silk proteins in modified salivary glands and spin silk fibers into protective tunnels, escape lines, and pupation cocoons. Molecular sequence data for these proteins are necessary to determine critical features of their function and evolution. To this end, we constructed an expression library from the silk glands of the ghost moth, Hepialus californicus, and characterized light chain fibroin and heavy chain fibroin gene transcripts. The predicted H. californicus silk fibroins share many elements with other lepidopteran and trichopteran fibroins, such as conserved placements of cysteine, aromatic, and polar amino acid residues. Further comparative analyses were performed to determine site-specific signatures of selection and to assess whether fibroin genes are informative as phylogenetic markers. We found that purifying selection has constrained mutation within the fibroins and that light chain fibroin is a promising molecular marker. Thus, by characterizing the H. californicus fibroins, we identified key functional amino acids and gained insight into the evolutionary processes that have shaped these adaptive molecules.

Comparison of embiopteran silks reveals tensile and structural similarities across Taxa.

Embioptera is a little studied order of widely distributed, but rarely seen, insects. Members of this group, also called embiids or webspinners, all heavily rely on silken tunnels in which they live and reproduce. However, embiids vary in their substrate preferences and these differences may result in divergent silk mechanical properties. Here, we present diameter measurements, tensile tests, and protein secondary structural analyses of silks spun by several embiid species. Despite their diverse habitats and phylogenetic relationships, these species have remarkably similar silk diameters and ultimate stress values. Yet, ultimate strain, Young's modulus, and toughness vary considerably. To better understand these tensile properties, Fourier transformed infrared spectroscopy was used to quantify secondary structural components. Compared to other arthropod silks, embiid silks are shown to have consistent secondary structures, suggesting that commonality of amino acid sequence motifs and small differences in structural composition can lead to significant changes in tensile properties.

Characterization of silk spun by the embiopteran, Antipaluria urichi.

Silks are renowned for being lightweight materials with impressive mechanical properties. Though moth and spider silks have received the most study, silk production has evolved in many other arthropods. One insect group that has been little investigated is Embioptera (webspinners). Embiopterans produce silk from unique tarsal spinning structures during all life stages. We characterize the molecular and mechanical properties of Antipaluria urichi (Embioptera) silk through multiple approaches. First, we quantify the number of silk secretory structures on their forelimbs and the tensile properties of Antipaluria silk. Second, we present silk protein (fibroin) transcripts from an embiopteran forelimb protarsomere cDNA library. We describe a fibroin that shares several features with other arthropod silks, including a subrepetitive core region, a non-repetitive carboxyl-terminal sequence, and a composition rich in glycine, alanine, and serine. Despite these shared attributes, embiopteran silk has several different tensile properties compared to previously measured silks. For example, the tensile strength of Antipaluria silk is much lower than that of Bombyx mori silk. We discuss the observed mechanical properties in relation to the fibroin sequence, spinning system, and embiopteran silk use.

Blueprint for a high-performance biomaterial: full-length spider dragline silk genes.

Spider dragline (major ampullate) silk outperforms virtually all other natural and manmade materials in terms of tensile strength and toughness. For this reason, the mass-production of artificial spider silks through transgenic technologies has been a major goal of biomimetics research. Although all known arthropod silk proteins are extremely large (>200 kiloDaltons), recombinant spider silks have been designed from short and incomplete cDNAs, the only available sequences. Here we describe the first full-length spider silk gene sequences and their flanking regions. These genes encode the MaSp1 and MaSp2 proteins that compose the black widow's high-performance dragline silk. Each gene includes a single enormous exon (>9000 base pairs) that translates into a highly repetitive polypeptide. Patterns of variation among sequence repeats at the amino acid and nucleotide levels indicate that the interaction of selection, intergenic recombination, and intragenic recombination governs the evolution of these highly unusual, modular proteins. Phylogenetic footprinting revealed putative regulatory elements in non-coding flanking sequences. Conservation of both upstream and downstream flanking sequences was especially striking between the two paralogous black widow major ampullate silk genes. Because these genes are co-expressed within the same silk gland, there may have been selection for similarity in regulatory regions. Our new data provide complete templates for synthesis of recombinant silk proteins that significantly improve the degree to which artificial silks mimic natural spider dragline fibers.