Population genomics of the critically endangered kākāpō

Summary The kākāpō is a flightless parrot endemic to New Zealand. Once common in the archipelago, only 201 individuals remain today, most of them descending from an isolated island population. We report the first genome-wide analyses of the species, including a high-quality genome assembly for kākāpō, one of the first chromosome-level reference genomes sequenced by the Vertebrate Genomes Project (VGP). We also sequenced and analyzed 35 modern genomes from the sole surviving island population and 14 genomes from the extinct mainland population. While theory suggests that such a small population is likely to have accumulated deleterious mutations through genetic drift, our analyses on the impact of the long-term small population size in kākāpō indicate that present-day island kākāpō have a reduced number of harmful mutations compared to mainland individuals. We hypothesize that this reduced mutational load is due to the island population having been subjected to a combination of genetic drift and purging of deleterious mutations, through increased inbreeding and purifying selection, since its isolation from the mainland ∼10,000 years ago. Our results provide evidence that small populations can survive even when isolated for hundreds of generations. This work provides key insights into kākāpō breeding and recovery and more generally into the application of genetic tools in conservation efforts for endangered species

Influence of management practice on the microbiota of a critically endangered species: A longitudinal study of kākāpō chick faeces and associated nest litter

Background: The critically endangered kākāpō is a flightless, nocturnal parrot endemic to Aotearoa New Zealand. Recent efforts to describe the gastrointestinal microbial community of this threatened herbivore revealed a low-diversity microbiota that is often dominated by Escherichia-Shigella bacteria. Given the importance of associated microbial communities to animal health, and increasing appreciation of their potential relevance to threatened species conservation, we sought to better understand the development of this unusual gut microbiota profile. To this end, we conducted a longitudinal analysis of faecal material collected from kākāpō chicks during the 2019 breeding season, in addition to associated nest litter material. Results: Using an experimental approach rarely seen in studies of threatened species microbiota, we evaluated the impact of a regular conservation practice on the developing kākāpō microbiota, namely the removal of faecal material from nests. Artificially removing chick faeces from nests had negligible impact on bacterial community diversity for either chicks or nests (p > 0.05). However, the gut microbiota did change significantly over time as chick age increased (p < 0.01), with an increasing relative abundance of Escherichia-Shigella coli over the study period and similar observations for the associated nest litter microbiota (p < 0.01). Supplementary feeding substantially altered gut bacterial diversity of kākāpō chicks (p < 0.01), characterised by a significant increase in Lactobacillus bacteria. Conclusions: Overall, chick age and hand rearing conditions had the most marked impact on faecal bacterial communities. Similarly, the surrounding nest litter microbiota changed significantly over time since a kākāpō chick was first placed in the nest, though we found no evidence that removal of faecal material influenced the bacterial communities of either litter or faecal samples. Taken together, these observations will inform ongoing conservation and management of this most enigmatic of bird species

Draft Genome Sequences of Four Novel Thermal- and Alkaline-Tolerant EgyptianRhizobiumStrains Nodulating Berseem Clover

exterior, view of Roman minaret (1164), 198

The potential for a CRISPR gene drive to eradicate or suppress globally invasive social wasps

CRISPR gene drives have potential for widespread and cost-efficient pest control, but are highly controversial. We examined a potential gene drive targeting spermatogenesis to control the invasive common wasp (Vespula vulgaris) in New Zealand. Vespula wasps are haplodiploid. Their life cycle makes gene drive production challenging, as nests are initiated by single fertilized queens in spring followed by several cohorts of sterile female workers and the production of reproductives in autumn. We show that different spermatogenesis genes have different levels of variation between introduced and native ranges, enabling a potential ‘precision drive’ that could target the reduced genetic diversity and genotypes within the invaded range. In vitro testing showed guide-RNA target specificity and efficacy that was dependent on the gene target within Vespula, but no cross-reactivity in other Hymenoptera. Mathematical modelling incorporating the genetic and life history traits of Vespula wasps identified characteristics for a male sterility drive to achieve population control. There was a trade-off between drive infiltration and impact: a drive causing complete male sterility would not spread, while partial sterility could be effective in limiting population size if the homing rate is high. Our results indicate that gene drives may offer viable suppression for wasps and other haplodiploid pests

Data from: Transcriptomic basis of genome by genome variation in a legume-rhizobia mutualism

In the legume-rhizobia mutualism, the benefit each partner derives from the other depends on the genetic identity of both host and rhizobial symbiont. To gain insight into the extent of genome x genome interactions on hosts at the molecular level and to identify potential mechanisms responsible for the variation, we examined host gene expression within nodules (the plant organ where the symbiosis occurs) of four genotypes of Medicago truncatula grown with either Ensifer meliloti or E. medicae symbionts. These host x symbiont combinations show significant variation in nodule and biomass phenotypes. Likewise, combinations differ in their transcriptomes:  host, symbiont, and host x symbiont affected the expression of 70%, 27% and 21%, respectively, of the approximately 27,000 host genes expressed in nodules. Genes with the highest levels of expression often varied between hosts and/or symbiont strain and include leghemoglobins that modulate oxygen availability and hundreds of Nodule Cysteine-Rich (NCR) peptides involved in symbiont differentiation and viability in nodules. Genes with host x symbiont dependent expression were enriched for functions related to resource exchange between partners (sugar/sulfate/iron/amino acid transport and dicarboxylate/amino acid synthesis). These enrichments suggest mechanisms for host control of the currencies of the mutualism. The transcriptome of M. truncatula accession HM101 (A17), the reference genome used for most molecular research, was less affected by symbiont identity than the other hosts. These findings underscore the importance of assessing the molecular basis of variation in ecologically important traits, particularly those involved in biotic interactions, in multiple genetic contexts

Genome-wide PVE analyses

Code that used GEMMA to get BSLMM and LMM chip-heritability estimates and summaries of results

Master README

Lists location of most important data and results and describes the analysis workflow

Single-variant Association Analyses

Code for and results from the single-variant association analyses with GEMMA

Galba: genome annotation with miniprot and AUGUSTUS

Abstract Background The Earth Biogenome Project has rapidly increased the number of available eukaryotic genomes, but most released genomes continue to lack annotation of protein-coding genes. In addition, no transcriptome data is available for some genomes. Results Various gene annotation tools have been developed but each has its limitations. Here, we introduce GALBA, a fully automated pipeline that utilizes miniprot, a rapid protein-to-genome aligner, in combination with AUGUSTUS to predict genes with high accuracy. Accuracy results indicate that GALBA is particularly strong in the annotation of large vertebrate genomes. We also present use cases in insects, vertebrates, and a land plant. GALBA is fully open source and available as a docker image for easy execution with Singularity in high-performance computing environments. Conclusions Our pipeline addresses the critical need for accurate gene annotation in newly sequenced genomes, and we believe that GALBA will greatly facilitate genome annotation for diverse organisms

Data from: Select and resequence reveals relative fitness of bacteria in symbiotic and free-living environments

Assays to accurately estimate relative fitness of bacteria growing in multistrain communities can advance our understanding of how selection shapes diversity within a lineage. Here, we present a variant of the “evolve and resequence” approach both to estimate relative fitness and to identify genetic variants responsible for fitness variation of symbiotic bacteria in free-living and host environments. We demonstrate the utility of this approach by characterizing selection by two plant hosts and in two free-living environments (sterilized soil and liquid media) acting on synthetic communities of the facultatively symbiotic bacterium Ensifer meliloti. We find (i) selection that hosts exert on rhizobial communities depends on competition among strains, (ii) selection is stronger inside hosts than in either free-living environment, and (iii) a positive host-dependent relationship between relative strain fitness in multistrain communities and host benefits provided by strains in single-strain experiments. The greatest changes in allele frequencies in response to plant hosts are in genes associated with motility, regulation of nitrogen fixation, and host/rhizobia signaling. The approach we present provides a powerful complement to experimental evolution and forward genetic screens for characterizing selection in bacterial populations, identifying gene function, and surveying the functional importance of naturally occurring genomic variation