Batoids are a diverse clade of flat cartilaginous fishes that occur primarily in benthic marine habitats. The skates and rays typically use their flexible pectoral fins for feeding and propulsion via undulatory swimming. However, two groups of rays have adopted a pelagic or bentho-pelagic lifestyle and utilize oscillatory swimming-the Myliobatidae and Gymnuridae. The myliobatids have evolved cephalic lobes, anteriorly extended appendages that are optimized for feeding, while their pectoral fins exhibit several modifications that likely arose in association with functional optimization of pelagic cruising via oscillatory flight. Here, we examine variation in fin ray distribution and ontogenetic timing of fin ray development in batoid pectoral fins in an evolutionary context using the following methods: radiography, computed tomography, dissections, and cleared and stained specimens. We propose an index for characterizing variation in the distribution of pectoral fin rays. While undulatory swimmers exhibit symmetry or slight anterior bias, we found a posterior shift in the distribution of fin rays that arose in two distinct lineages in association with oscillatory swimming. Undulatory and oscillatory swimmers occupy nonoverlapping morphospace with respect to fin ray distribution illustrating significant remodeling of pectoral fins in oscillatory swimmers. Further, we describe a derived skeletal feature in anterior pectoral fins of the Myliobatidae that is likely associated with optimization of oscillatory swimming. By examining the distribution of fin rays with clearly defined articulation points, we were able to infer evolutionary trends and body plan remodeling associated with invasion of the pelagic environment. Finally, we found that the number and distribution of fin rays is set early in development in the little skate, round stingray, and cownose ray, suggesting that fin ray counts from specimens after birth or hatching are representative of adults and therefore comparable among species.

Adaptive phenotypes are shaped by a combination of genetic and environmental forces, and while the literature is rich with studies focusing on either genetics or environment contributions, those that consider both are rare. Here we utilize the cichlid oral jaw apparatus to fill this knowledge-gap. First, we employed RNA-seq in bony and ligamentous tissues important for jaw opening to identify differentially expressed genes between species and across foraging environments. Our foraging treatments were designed to force animals to employ either suction or biting/scraping, which broadly mimic pelagic or benthic modes of feeding. We found a large number of differentially expressed genes between species, and while we identified relatively few differences between environments, species differences were far more pronounced when reared in the pelagic versus benthic environment. Further, these data carried the signature of genetic assimilation, and implicated cell cycle regulation in shaping the jaw across species and environments. Next, we repeated the foraging experiment and performed ATAC-seq procedures on nuclei harvested from the same tissues. Cross-referencing results from both analyses revealed subsets of genes that were both differentially expressed and differentially accessible in either the pelagic (n=15) or the benthic environment (n=11), as well as loci where differences were robust to foraging environment (n=13). All in all, these data provide novel insights into the epigenetic, genetic, and cellular bases of species- and environment-specific bone shapes, as well as the evolution of phenotypic plasticity in this iconic model system.

A major ongoing challenge in biological research is to understand the origin and maintenance of biodiversity, with broad implications in conservation, ecology and evolutionary biology. Traditionally, these endeavours have involved characterizing the forces and mechanisms operating above the organismal level (e.g. selection, environmental

Manta rays and their relatives of the family Myliobatidae have pectoral fins that have been modified for underwater flight, as well as a pair of fleshy projections at the anterior of the body called cephalic lobes, which are specialized for feeding. As a unique trait with a dedicated function, cephalic lobes offer an excellent opportunity to elucidate the processes by which diverse body plans and features evolve. To shed light on the morphological development and genetic underpinnings of cephalic lobes, we examined paired fin development in cownose rays, which represent the sister taxon to manta rays in the genus Mobula. We find that cephalic lobes develop as anterior pectoral fin domains and lack independent posterior patterning by 5 ′ HoxD genes and Shh, indicating that cephalic lobes are not independent appendages but rather are modified pectoral fin domains. In addition, by leveraging interspecies comparative transcriptomics and domain-specific RNA-sequencing, we identify shared expression of anterior patterning genes, including Alx1, Alx4, Pax9, Hoxa13, Hoxa2, and Hoxd4, in the pectoral fins of cownose ray (Rhinoptera bonasus) and little skate (Leucoraja erinacea), providing evidence supporting homology between the cephalic lobes of myliobatids and the anterior pectoral fins of skates. We also suggest candidate genes that may be involved in development of myliobatid-specific features, including Omd, which is likely associated with development of thick anterior pectoral fin radials of myliobatids, and Dkk1, which may inhibit tissue outgrowth at the posterior boundary of the developing cephalic lobes. Finally, we observe that cephalic lobes share a surprising number of developmental similarities with another paired fin modification: the claspers of male cartilaginous fishes, including enrichment of Hand2, Hoxa13, and androgen receptor. These results suggest that cephalic lobes may have evolved by co-opting developmental pathways that specify novel domains in paired fins. Taken together, these data on morphological development and comparative gene expression patterns illustrate how distinct body plans and seemingly novel features can arise via subtle changes to existing developmental pathways.

Heat Shock Protein 101 (HSP101) in plants, and bacterial and yeast orthologs, is essential for thermotolerance. To investigate thermotolerance mechanisms involving HSP101, we performed a suppressor screen in Arabidopsis thaliana of a missense HSP101 allele (hot1-4). hot1-4 plants are sensitive to acclimation heat treatments that are otherwise permissive for HSP101 null mutants, indicating that the hot1-4 protein is toxic. We report one suppressor (shot2, suppressor of hot1-4 2) has a missense mutation of a conserved residue in CLEAVAGE STIMULATION FACTOR77 (CstF77), a subunit of the polyadenylation complex critical for mRNA 3’ end maturation. We performed ribosomal RNA depletion RNA-Seq and captured transcriptional readthrough with a custom bioinformatics pipeline. Acclimation heat treatment caused transcriptional readthrough in hot1-4 shot2, with more readthrough in heat-induced genes, reducing the levels of toxic hot1-4 protein and suppressing hot1-4 heat sensitivity. Although shot2 mutants develop like the wild type in the absence of stress and survive mild heat stress, reduction of heat-induced genes and decreased HSP accumulation makes shot2 in HSP101 null and wild-type backgrounds sensitive to severe heat stress. Our study reveals the critical function of CstF77 for 3’ end formation of mRNA and the dominant role of HSP101 in dictating the outcome of severe heat stress.