Genome-wide association studies (GWAS) have been highly informative in discovering disease-associated loci but are not designed to capture all structural variations in the human genome. Using long-read sequencing data, we discovered widespread structural variation within SINE-VNTR-Alu (SVA) elements, a class of great ape-specific transposable elements with gene-regulatory roles, which represents a major source of structural variability in the human population. We highlight the presence of structurally variable SVAs (SV-SVAs) in neurological disease-associated loci, and we further associate SV-SVAs to disease-associated SNPs and differential gene expression using luciferase assays and expression quantitative trait loci data. Finally, we genetically deleted SV-SVAs in the BIN1 and CD2AP Alzheimer’s disease-associated risk loci and in the BCKDK Parkinson’s disease-associated risk locus and assessed multiple aspects of their gene-regulatory influence in a human neuronal context. Together, this study reveals a novel layer of genetic variation in transposable elements that may contribute to identification of the structural variants that are the actual drivers of disease associations of GWAS loci.

The genus Nomada Scopoli (Hymenoptera: Apidae) is the largest genus of brood parasitic bees with nearly 800 species found across the globe and in nearly all biogeographic realms except Antarctica. There is no previous molecular phylogeny focused on Nomada despite their high species abundance nor is there an existing comprehensive biogeography for the genus. Using ultraconserved element (UCE) phylogenomic data, we constructed the first molecular phylogeny for the genus Nomada and tested the monophyly of 16 morphologically established species groups. We also estimated divergence dates using fossil calibration points and inferred the origin and of this genus around the globe. Our phylogeny recovered 14 of the 16 previously established species groups as monophyletic. The superba and ruficornis groups, however, were recovered as non-monophyletic and need to be re-evaluated using morphology. Divergence dating and historic biogeographic analyses performed on the phylogenetic reconstruction indicates that Nomada most likely originated in the Holarctic ∼65 Mya. Geodispersal into the southern hemisphere occurred three times; once during the Eocene into the Afrotropics, once during the Oligocene into the Neotropics, and once during the Miocene into Australasia. Geodispersal across the Holarctic was most frequent and occurred repeatedly throughout the Cenozoic era, using the De Geer, Thulean, and the Bering Land Bridges. This is the first instance of a bee using both the Thulean and De Geer land bridges and has implications of how early bee species dispersed throughout the Palearctic in the late Cretaceous and early Paleogene.

Archerfishes (Toxotidae) are variously found in the fresh- and brackish-water environments of Asia Pacific and are well known for their ability to shoot water at terrestrial prey. These shots of water are intended to strike their prey and cause it to fall into the water for capture and consumption. While this behavior is well known, there are competing hypotheses (blowpipe vs. pressure tank hypothesis) of how archerfishes shoot and which oral structures are involved. Current understanding of archerfish shooting structures is largely based on two species, Toxotes chatareus and T. jaculatrix. We do not know if all archerfishes possess the same oral structures to shoot water, if anatomical variation is present within these oral structures, or how these features have evolved. Additionally, there is little information on the evolution of the Toxotidae as a whole, with all previous systematic works focusing on the interrelationships of the family. We first investigate the limits of archerfish species using new and previously published genetic data. Our analyses highlight that the current taxonomy of archerfishes does not conform to the relationships we recover. Toxotes mekongensis and T. siamensis are placed in synonymy of T. chatareus, Toxotes carpentariensis is recognized as a species and removed from synonymy of T. chatareus, and the genus Protoxotes is recognized for T. lorentzi based on the results of our analyses. We then take an integrative approach, using a combined analysis of discrete hard- and soft-tissue morphological characters with genetic data, to construct a phylogeny of the Toxotidae. Using the resulting phylogenetic hypothesis, we then characterize the evolutionary history and anatomical variation within the archerfishes. We discuss the variation in the oral structures and the evolution of the mechanism with respect to the interrelationships of archerfishes, and find that the oral structures of archerfishes support the blowpipe hypothesis but soft-tissue oral structures may also play a role in shooting. Finally, by comparing the morphology of archerfishes to their sister group, we find that the Leptobramidae has relevant shooting features in the oral cavity, suggesting that some components of the archerfish shooting mechanism are examples of co-opted or exapted traits.

Abstract Besides the consequences of retrotransposition, long interspersed element 1 (L1) retrotransposons can affect the host genome through their antisense promoter. In addition to the sense promoter, the evolutionarily recent L1 retrotransposons, which are present in several thousand copies, also possess an anti-sense promoter that can produce L1 chimeric transcripts (LCT) composed of the L1 5’ UTR followed by the adjacent genomic sequence. The full extent to which LCT expression occurs in a given tissue and whether disruption of the defense mechanisms that normally control L1 retrotransposons affects their expression and function in cancer cells, remain to be established. By using CLIFinder, a dedicated bioinformatics tool, we found that LCT expression was widespread in normal brain and aggressive glioma samples, and that approximately 17% of recent L1 retrotransposons, from the L1PA1 to L1PA7 subfamilies, were involved in their production. Importantly, the transcriptional activities of the L1 antisense promoters and of their host loci were coupled. Accordingly, we detected LCT-producing L1 retrotransposons mainly in transcriptionally active genes and genomic loci. Moreover, changes in the host genomic locus expression level in glioma were associated with a similar change in LCT expression level, regardless of the L1 promoter methylation status. Our findings support a model in which the host genomic locus transcriptional activity is the main driving force of LCT expression. We hypothesize that this model is more applicable when host gene and LCT are transcribed from the same strand.

The parasitoid lifestyle is largely regarded as a key innovation that contributed to the evolutionary success and extreme species richness of the order Hymenoptera. Understanding the phylogenetic history of hyperdiverse parasitoid groups is a fundamental step in elucidating the evolution of biological traits linked to parasitoidism. We used a genomic-scale dataset based on ultra-conserved elements and the most comprehensive taxon sampling to date to estimate the evolutionary relationships of Braconidae, the second largest family of Hymenoptera. Based on our results, we propose the existence of 41 braconid subfamilies, confirmed a number of subfamilial placements and proposed subfamily-level taxonomic changes, notably the restoration of Trachypetinae stat. rev. and Masoninae stat. rev. as subfamilies of Braconidae, confirmation that Apozyx penyai Mason belongs in Braconidae placed in the subfamily Apozyginae and the recognition of Ichneutinae sensu stricto and Proteropinae as non-cyclostome subfamilies robustly supported in a phylogenetic context. The correlation between koinobiosis with endoparasitoidism and idiobiosis with ectoparasitoidism, long thought to be an important aspect in parasitoid life history, was formally tested and confirmed in a phylogenetic framework. Using ancestral reconstruction methods based on both parsimony and maximum likelihood, we suggest that the ancestor of the braconoid complex was a koinobiont endoparasitoid, as was that of the cyclostome sensu lato clade. Our results also provide strong evidence for one transition from endo- to ectoparasitoidism and three reversals back to endoparasitoidism within the cyclostome sensu stricto lineage. Transitions of koino- and idiobiosis were identical to that inferred for endo- versus ectoparasitoidism, except with one additional reversal back to koinobiosis in the small subfamily Rhysipolinae.

We investigated the genetic diversity and historic relationships among southern African sheep as well as the relationships between them and sheep outside the continent by sourcing both archaeological and modern sheep samples. Archaeological sheep samples derived from the site Die Kelders 1, near Cape Town, date to approximately 1500 years ago. The modern samples were taken as ear snips from Damara, Namaqua Afrikaner, and Ronderib Afrikaner sheep on a farm in Prieska in the Northern Cape. Illumina sequencing libraries were constructed for both ancient and modern specimens. Ancient specimens were enriched for the mitochondrial genome using an in-solution hybridization protocol and modern specimens were subjected to shotgun sequencing. Sequences were mapped to the Ovis aries reference genome, assigned to haplogroups and subhaplogroups, and used to calculate a phylogenetic tree using previously published, geographically dispersed mitochondrial genome sheep sequences. Genetic diversity statistics show that southern African sheep have lower diversity than sheep in other regions. Phylogenetic analysis reveals that many modern southern African sheep are likely descended from prehistoric indigenous sheep populations and not from sheep imported from Europe during the historic period.

Recognition of the untapped potential of photosynthesis to improve crop yields has spurred research to identify targets for breeding. The CO2-fixing enzyme Rubisco is characterised by a number of inefficiencies and frequently limits carbon assimilation at the top of the canopy, representing a clear target for wheat improvement. Two bread wheat lines with similar genetic backgrounds and contrasting in vivo maximum carboxylation activity of Rubisco per unit leaf nitrogen (Vc,max,25/Narea) determined using high throughput phenotyping methods were selected for detailed study from a panel of 80 spring wheat lines. Detailed phenotyping of photosynthetic traits in the two lines using glasshouse-grown plants showed no difference in Vc,max,25/Narea determined directly via in vivo and in vitro methods. Detailed phenotyping of glasshouse-grown plants of the 80 wheat lines also showed no correlation between photosynthetic traits measured via high throughput phenotyping of field-grown plants. Our findings suggest that the complex interplay between traits determining crop productivity and the dynamic environments experienced by field-grown plants needs to be considered when designing strategies for effective wheat crop yield improvement when breeding for particular environments.

Is there a universal genetically programmed defense providing tolerance to antibiotics when bacteria grow as biofilms? A comparison between biofilms of three different bacterial species by transcriptomic and metabolomic approaches uncovered no evidence of one. Single-species biofilms of three bacterial species ( Pseudomonas aeruginosa , Staphylococcus aureus , and Acinetobacter baumannii ) were grown in vitro for 3 days and then challenged with respective antibiotics (ciprofloxacin, daptomycin, and tigecycline) for an additional 24 h. , ABSTRACT Is there a universal genetically programmed defense providing tolerance to antibiotics when bacteria grow as biofilms? A comparison between biofilms of three different bacterial species by transcriptomic and metabolomic approaches uncovered no evidence of one. Single-species biofilms of three bacterial species ( Pseudomonas aeruginosa , Staphylococcus aureus , and Acinetobacter baumannii ) were grown in vitro for 3 days and then challenged with respective antibiotics (ciprofloxacin, daptomycin, and tigecycline) for an additional 24 h. All three microorganisms displayed reduced susceptibility in biofilms compared to planktonic cultures. Global transcriptomic profiling of gene expression comparing biofilm to planktonic and antibiotic-treated biofilm to untreated biofilm was performed. Extracellular metabolites were measured to characterize the utilization of carbon sources between biofilms, treated biofilms, and planktonic cells. While all three bacteria exhibited a species-specific signature of stationary phase, no conserved gene, gene set, or common functional pathway could be identified that changed consistently across the three microorganisms. Across the three species, glucose consumption was increased in biofilms compared to planktonic cells, and alanine and aspartic acid utilization were decreased in biofilms compared to planktonic cells. The reasons for these changes were not readily apparent in the transcriptomes. No common shift in the utilization pattern of carbon sources was discerned when comparing untreated to antibiotic-exposed biofilms. Overall, our measurements do not support the existence of a common genetic or biochemical basis for biofilm tolerance against antibiotics. Rather, there are likely myriad genes, proteins, and metabolic pathways that influence the physiological state of individual microorganisms in biofilms and contribute to antibiotic tolerance.

Present-day Tibetans have adapted both genetically and culturally to the high altitude environment of the Tibetan Plateau, but fundamental questions about their origins remain unanswered. Recent archaeological and genetic research suggests the presence of an early population on the Plateau within the past 40 thousand years, followed by the arrival of subsequent groups within the past 10 thousand years. Here, we obtain new genome-wide data for 33 ancient individuals from high elevation sites on the southern fringe of the Tibetan Plateau in Nepal, who we show are most closely related to present-day Tibetans. They derive most of their ancestry from groups related to Late Neolithic populations at the northeastern edge of the Tibetan Plateau but also harbor a minor genetic component from a distinct and deep Paleolithic Eurasian ancestry. In contrast to their Tibetan neighbors, present-day non-Tibetan Tibeto-Burman speakers living at mid-elevations along the southern and eastern margins of the Plateau form a genetic cline that reflects a distinct genetic history. Finally, a comparison between ancient and present-day highlanders confirms ongoing positive selection of high altitude adaptive alleles.