The genetic consequences of species-wide declines are rarely quantified because the timing and extent of the decline varies across the species’ range. The sea otter (Enhydra lutris) is a unique model in this regard. Their dramatic decline from thousands to fewer than 100 individuals per population occurred range-wide and nearly simultaneously due to the 18th-19th century fur trade. Consequently, each sea otter population represents an independent natural experiment of recovery after extreme population decline. We designed sequence capture probes for 50 megabases of sea otter exonic and neutral genomic regions. We sequenced 107 sea otters from five populations that span the species range to high coverage (18-76x) and three historic Californian samples from 1500 and 200 years ago to low coverage (1.5-3.5X). We observe distinct population structure and find that sea otters in California are the last survivors of a divergent lineage isolated for thousands of years and therefore warrant special conservation concern. We detect signals of extreme population decline in every surviving sea otter population and use this demographic history to design forward-in-time simulations of coding sequence. Our simulations indicate that this decline could lower the fitness of recovering populations for generations. However, the simulations also demonstrate how historically low effective population sizes prior to the fur trade may have mitigated the effects of population decline on genetic health. Our comprehensive approach shows how demographic inference from genomic data, coupled with simulations, allows assessment of extinction risk and different models of recovery.

Until recently many historical museum specimens were largely inaccessible to genomic inquiry, but high-throughput sequencing (HTS) approaches have allowed researchers to successfully sequence genomic DNA from dried and fluid-preserved museum specimens. In addition to preserved specimens, many museums contain large series of allozyme supernatant samples but the amenability of these samples to HTS has not yet been assessed. Here, we compared the performance of a target-capture approach using alternative sources of genomic DNA from ten specimens of spring salamanders (Plethodontidae: Gyrinophilus porphyriticus) collected between 1985 and 1990: allozyme supernatants, allozyme homogenate pellets, and formalin-fixed tissues. We designed capture probes based on double-digest restriction-site associated sequencing (RADseq) derived loci from frozen blood samples available for seven of the specimens and assessed the success and consistency of capture and RADseq approaches. This study design enabled direct comparisons of data quality and potential biases among the different datasets for phylogenomic and population genomic analyses. We found that in phylogenetic analyses, all enrichment types for a given specimen clustered together. In principal component space all capture-based samples clustered together, but RADseq samples did not cluster with corresponding capture-based samples. SNP calls were on average 18.3% different between enrichment types for a given individual, but these discrepancies were primarily due to differences in heterozygous/homozygous SNP calls. We demonstrate that both allozyme supernatant and formalin-fixed samples can be successfully used for population genomic analyses and we discuss ways to identify and reduce biases associated with combining capture and RADseq data.

The octocoral genus Chrysogorgia Duchassaing & Michelotti, 1864 contains about 80 nominal species that are ecologically important components of benthic communities. A review of morphological variation within Chrysogorgia identified nine groups that may represent distinct lineages; which would help clarify their taxonomy and guide genus revision. Here, we applied a recently developed universal target enrichment bait method for octocoral exons and ultraconserved elements (UCEs) on 96 specimens varying in morphology, collection ages and DNA quality and quantity to determine whether there was genetic support for morphologically defined groups within Chrysogorgia. Following Illumina sequencing and SPAdes assembly we recovered 1682 of 1700 targeted exon loci and 1333 of 1340 targeted UCE loci. Loci recovery per sample was highly variable and significantly correlated with time since specimen collection (2-60 years) and DNA quantity and quality. Phylogenetically informative sites in UCE and exon loci ranged from 34.75 – 36.10% for 50% and 75% taxon-occupancy matrices, respectively. Maximum likelihood analyses recovered highly resolved trees with topologies supporting the recognition of 11 candidate genera, nine of which are novel. Our results also demonstrate that this target-enrichment approach can be applied to degraded museum specimens of up to 60 years old. This study shows that an integrative approach consisting of molecular and morphological methods is essential in a proper revision of Chrysogorgiath taxonomy and regional diversity of these ecologically important corals.

Crop brassicas include three diploid [ Brassica rapa (AA; 2 n = 2 x = 16), B. nigra (BB; 2 n = 2 x = 18), and B. oleracea (CC; 2 n = 2 x = 20)] and three derived allotetraploid species. It is difficult to distinguish Brassica chromosomes as they are small and morphologically similar. We aimed to develop a genome-sequence based cytogenetic toolkit for reproducible identification of Brassica chromosomes and their structural variations. A bioinformatic pipeline was used to extract repeat-free sequences from the whole genome assembly of B. rapa . Identified sequences were subsequently used to develop four c. 47-mer oligonucleotide libraries comprising 27,100, 11,084, 9,291, and 16,312 oligonucleotides. We selected these oligonucleotides after removing repeats from 18 identified sites (500–1,000 kb) with 1,997–5,420 oligonucleotides localized at each site in B. rapa . For one set of probes, a new method for amplification or immortalization of the library is described. oligonucleotide probes produced specific and reproducible in situ hybridization patterns for all chromosomes belonging to A, B, C, and R ( Raphanus sativu s) genomes. The probes were able to identify structural changes between the genomes, including translocations, fusions, and deletions. Furthermore, the probes were able to identify a structural translocation between a pak choi and turnip cultivar of B. rapa. Overall, the comparative chromosomal mapping helps understand the role of chromosome structural changes during genome evolution and speciation in the family Brassicaceae. The probes can also be used to identify chromosomes in aneuploids such as addition lines used for gene mapping, and to track transfer of chromosomes in hybridization and breeding programs.

The systematics and taxonomy of the tropical Asian jumping spiders of the tribe Baviini is reviewed, with a molecular phylogenetic study (UCE sequence capture, traditional Sanger sequencing) guiding a reclassification of the group’s genera. The well-studied members of the group are placed into six genera: Bavia Simon, 1877, Indopadilla Caleb & Sankaran, 2019, Padillothorax Simon, 1901, Piranthus Thorell, 1895, Stagetillus Simon, 1885, and one new genus, Maripanthus Maddison, gen. nov. The identity of Padillothorax is clarified, and Bavirecta Kanesharatnam & Benjamin, 2018 synonymized with it. Hyctiota Strand, 1911 is synonymized with Stagetillus. The molecular phylogeny divides the baviines into three clades, the Piranthus clade with a long embolus (Piranthus, Maripanthus), the genus Padillothorax with a flat body and short embolus, and the Bavia clade with a higher body and (usually) short embolus (remaining genera). In general, morphological synapomorphies support or extend the molecularly delimited groups. Eighteen new species are described: Bavia nessagyna, Indopadilla bamilin, I. kodagura, I. nesinor, I. redunca, I. redynis, I. sabivia, I. vimedaba, Maripanthus draconis (type species of Maripanthus), M. jubatus, M. reinholdae, Padillothorax badut, P. mulu, Piranthus api, P. bakau, P. kohi, P. mandai, and Stagetillus irri, all sp. nov., with taxonomic authority W. Maddison. The distinctions between baviines and the astioid Nungia Żabka, 1985 are reviewed, leading to four species being moved into Nungia from Bavia and other genera. Fifteen new combinations are established: Bavia maurerae (Freudenschuss & Seiter, 2016), Indopadilla annamita (Simon, 1903), I. kahariana (Prószyński & Deeleman-Reinhold, 2013), I. sonsorol (Berry, Beatty & Prószyński, 1997), I. suhartoi (Prószyński & Deeleman-Reinhold, 2013), Maripanthus menghaiensis (Cao & Li, 2016), M. smedleyi (Reimoser, 1929), Nungia hatamensis (Thorell, 1881), N. modesta (Keyserling, 1883), N. papakula (Strand, 1911), N. xiaolonghaensis (Cao & Li, 2016), Padillothorax casteti (Simon, 1900), P. exilis (Cao & Li, 2016), P. flavopunctus (Kanesharatnam & Benjamin, 2018), Stagetillus banda (Strand, 1911), all comb. nov. One combination is restored, Bavia capistrata (C. L. Koch, 1846). Five of these new or restored combinations correct previous errors of placing species in genera that have superficially similar palps but extremely different body forms, in fact belonging in distantly related tribes, emphasizing that the general shape of male palps should be used with caution in determining relationships. A little-studied genus, Padillothorus Prószyński, 2018, is tentatively assigned to the Baviini. Ligdus Thorell, 1895 is assigned to the Ballini.

Repeatable, convergent outcomes are prima facie evidence for determinism in evolutionary processes. Among fishes, well-known examples include microevolutionary habitat transitions into the water column, where freshwater populations (e.g., sticklebacks, cichlids, and whitefishes) recurrently diverge toward slender-bodied pelagic forms and deep-bodied benthic forms. However, the consequences of such processes at deeper macroevolutionary scales in the marine environment are less clear. We applied a phylogenomics-based integrative, comparative approach to test hypotheses about the scope and strength of convergence in a marine fish clade with a worldwide distribution (snappers and fusiliers, family Lutjanidae) featuring multiple water-column transitions over the past 45 million years. We collected genome-wide exon data for 110 (∼80%) species in the group and aggregated data layers for body shape, habitat occupancy, geographic distribution, and paleontological and geological information. We also implemented approaches using genomic subsets to account for phylogenetic uncertainty in comparative analyses. Our results show independent incursions into the water column by ancestral benthic lineages in all major oceanic basins. These evolutionary transitions are persistently associated with convergent phenotypes, where deep-bodied benthic forms with truncate caudal fins repeatedly evolve into slender midwater species with furcate caudal fins. Lineage diversification and transition dynamics vary asymmetrically between habitats, with benthic lineages diversifying faster and colonizing midwater habitats more often than the reverse. Convergent ecological and functional phenotypes along the benthic–pelagic axis are pervasive among different lineages and across vastly different evolutionary scales, achieving predictable high-fitness solutions for similar environmental challenges, ultimately demonstrating strong determinism in fish body-shape evolution.

Parmeliaceae is the largest family of lichen-forming fungi with a worldwide distribution. We used a target enrichment data set and a qualitative selection method for 250 out of 350 genes to infer the phylogeny of the major clades in this family including 81 taxa, with both subfamilies and all seven major clades previously recognized in the subfamily Parmelioideae. The reduced genome-scale data set was analyzed using concatenated-based Bayesian inference and two different Maximum Likelihood analyses, and a coalescent-based species tree method. The resulting topology was strongly supported with the majority of nodes being fully supported in all three concatenated-based analyses. The two subfamilies and each of the seven major clades in Parmelioideae were strongly supported as monophyletic. In addition, most backbone relationships in the topology were recovered with high nodal support. The genus Parmotrema was found to be polyphyletic and consequently, it is suggested to accept the genus Crespoa to accommodate the species previously placed in Parmotrema subgen. Crespoa. This study demonstrates the power of reduced genome-scale data sets to resolve phylogenetic relationships with high support. Due to lower costs, target enrichment methods provide a promising avenue for phylogenetic studies including larger taxonomic/specimen sampling than whole genome data would allow.

Modified-live herpesvirus vaccines are widely used in humans and animals, but field strains can emerge that have a higher virulence and break vaccinal protection. Since the introduction of the first vaccine in the 1970s, Marek’s disease virus overcame the vaccine barrier by the acquisition of numerous genomic mutations. However, the evolutionary adaptations in the herpesvirus genome responsible for the vaccine breaks have remained elusive. Here, we demonstrate that point mutations in the multifunctional meq gene acquired during evolution can significantly alter virulence. Defined mutations found in highly virulent strains also allowed the virus to overcome innate cellular responses and vaccinal protection. Concomitantly, the adaptations in meq enhanced virus shedding into the environment, likely providing a selective advantage for the virus. Our study provides the first experimental evidence that few point mutations in a single herpesviral gene result in drastically increased virulence, enhanced shedding, and escape from vaccinal protection.