Search

UNLABELLED: Biodegradable plastics are urgently needed to replace petroleum-derived polymeric materials and prevent their accumulation in the environment. To this end, we isolated and characterized a halophilic and alkaliphilic bacterium from the Great Salt Lake in Utah. The isolate was identified as a Halomonas species and designated CUBES01. Full-genome sequencing and genomic reconstruction revealed the unique genetic traits and metabolic capabilities of the strain, including the common polyhydroxyalkanoate (PHA) biosynthesis pathway. Fluorescence staining identified intracellular polyester granules that accumulated predominantly during the strains exponential growth, a feature rarely found among natural PHA producers. CUBES01 was found to metabolize a range of renewable carbon feedstocks, including glucosamine and acetyl-glucosamine, as well as sucrose, glucose, fructose, and further glycerol, propionate, and acetate. Depending on the substrate, the strain accumulated up to ~60% of its biomass (dry wt/wt) in poly(3-hydroxybutyrate), while reaching a doubling time of 1.7 h at 30°C and an optimum osmolarity of 1 M sodium chloride and a pH of 8.8. The physiological preferences of the strain may not only enable long-term aseptic cultivation but also facilitate the release of intracellular products through osmolysis. The development of a minimal medium also allowed the estimation of maximum polyhydroxybutyrate production rates, which were projected to exceed 5 g/h. Finally, also, the genetic tractability of the strain was assessed in conjugation experiments: two orthogonal plasmid vectors were stable in the heterologous host, thereby opening the possibility of genetic engineering through the introduction of foreign genes. IMPORTANCE: The urgent need for renewable replacements for synthetic materials may be addressed through microbial biotechnology. To simplify the large-scale implementation of such bio-processes, robust cell factories that can utilize sustainable and widely available feedstocks are pivotal. To this end, non-axenic growth-associated production could reduce operational costs and enhance biomass productivity, thereby improving commercial competitiveness. Another major cost factor is downstream processing, especially in the case of intracellular products, such as bio-polyesters. Simplified cell-lysis strategies could also further improve economic viability.

The American black bear, Ursus americanus, is a widespread and ecologically important species in North America. In California, the black bear plays an important role in a variety of ecosystems and serves as an important species for recreational hunting. While research suggests that the populations in California are currently healthy, continued monitoring is critical, with genomic analyses providing an important surveillance tool. Here we report a high-quality, near chromosome-level genome assembly from a U. americanus sample from California. The primary assembly has a total length of 2.5 Gb contained in 316 scaffolds, a contig N50 of 58.9 Mb, a scaffold N50 of 67.6 Mb, and a BUSCO completeness score of 96%. This U. americanus genome assembly will provide an important resource for the targeted management of black bear populations in California, with the goal of achieving an appropriate balance between the recreational value of black bears and the maintenance of viable populations. The high quality of this genome assembly will also make it a valuable resource for comparative genomic analyses among black bear populations and among bear species.

As demonstrated by the SARS-CoV-2 pandemic, the emergence of novel viral strains with increased transmission rates poses a serious threat to global health. Statistical models of genome sequence evolution may provide a critical tool for early detection of these strains. Using a novel stochastic model that links transmission rates to the entire viral genome sequence, we study the utility of phylogenetic methods that use a phylogenetic tree relating viral samples versus count-based methods that use case counts of variants over time exclusively to detect increased transmission rates and identify candidate causative mutations. We find that phylogenies in particular can detect novel transmission-enhancing variants very soon after their origin and may facilitate the development of early detection systems for outbreak surveillance.

The higher classification of termites requires substantial revision as the Neoisoptera, the most diverse termite lineage, comprise many paraphyletic and polyphyletic higher taxa. Here, we produce an updated termite classification using genomic-scale analyses. We reconstruct phylogenies under diverse substitution models with ultraconserved elements analyzed as concatenated matrices or within the multi-species coalescence framework. Our classification is further supported by analyses controlling for rogue loci and taxa, and topological tests. We show that the Neoisoptera are composed of seven family-level monophyletic lineages, including the Heterotermitidae Froggatt, Psammotermitidae Holmgren, and Termitogetonidae Holmgren, raised from subfamilial rank. The species-rich Termitidae are composed of 18 subfamily-level monophyletic lineages, including the new subfamilies Crepititermitinae, Cylindrotermitinae, Forficulitermitinae, Neocapritermitinae, Protohamitermitinae, and Promirotermitinae; and the revived Amitermitinae Kemner, Microcerotermitinae Holmgren, and Mirocapritermitinae Kemner. Building an updated taxonomic classification on the foundation of unambiguously supported monophyletic lineages makes it highly resilient to potential destabilization caused by the future availability of novel phylogenetic markers and methods. The taxonomic stability is further guaranteed by the modularity of the new termite classification, designed to accommodate as-yet undescribed species with uncertain affinities to the herein delimited monophyletic lineages in the form of new families or subfamilies.

Papua New Guinea (PNG) hosts distinct environments mainly represented by the ecoregions of the Highlands and Lowlands that display increased altitude and a predominance of pathogens, respectively. Since its initial peopling approximately 50,000 years ago, inhabitants of these ecoregions might have differentially adapted to the environmental pressures exerted by each of them. However, the genetic basis of adaptation in populations from these areas remains understudied. Here, we investigated signals of positive selection in 62 highlanders and 43 lowlanders across 14 locations in the main island of PNG using whole-genome genotype data from the Oceanian Genome Variation Project (OGVP) and searched for signals of positive selection through population differentiation and haplotype-based selection scans. Additionally, we performed archaic ancestry estimation to detect selection signals in highlanders within introgressed regions of the genome. Among highland populations we identified candidate genes representing known biomarkers for mountain sickness (SAA4, SAA1, PRDX1, LDHA) as well as candidate genes of the Notch signaling pathway (PSEN1, NUMB, RBPJ, MAML3), a novel proposed pathway for high altitude adaptation in multiple organisms. We also identified candidate genes involved in oxidative stress, inflammation, and angiogenesis, processes inducible by hypoxia, as well as in components of the eye lens and the immune response. In contrast, candidate genes in the lowlands are mainly related to the immune response (HLA-DQB1, HLA-DQA2, TAAR6, TAAR9, TAAR8, RNASE4, RNASE6, ANG). Moreover, we find two candidate regions to be also enriched with archaic introgressed segments, suggesting that archaic admixture has played a role in the local adaptation of PNG populations.

MOTIVATION: Somatic mosaicism has been implicated in several developmental disorders, cancers, and other diseases. Short tandem repeats (STRs) consist of repeated sequences of 1-6 bp and comprise >1 million loci in the human genome. Somatic mosaicism at STRs is known to play a key role in the pathogenicity of loci implicated in repeat expansion disorders and is highly prevalent in cancers exhibiting microsatellite instability. While a variety of tools have been developed to genotype germline variation at STRs, a method for systematically identifying mosaic STRs is lacking. RESULTS: We introduce prancSTR, a novel method for detecting mosaic STRs from individual high-throughput sequencing datasets. prancSTR is designed to detect loci characterized by a single high-frequency mosaic allele, but can also detect loci with multiple mosaic alleles. Unlike many existing mosaicism detection methods for other variant types, prancSTR does not require a matched control sample as input. We show that prancSTR accurately identifies mosaic STRs in simulated data, demonstrate its feasibility by identifying candidate mosaic STRs in Illumina whole genome sequencing data derived from lymphoblastoid cell lines for individuals sequenced by the 1000 Genomes Project, and evaluate the use of prancSTR on Element and PacBio data. In addition to prancSTR, we present simTR, a novel simulation framework which simulates raw sequencing reads with realistic error profiles at STRs. AVAILABILITY AND IMPLEMENTATION: prancSTR and simTR are freely available at https://github.com/gymrek-lab/trtools. Detailed documentation is available at https://trtools.readthedocs.io/.

Schizophyllum commune is a mushroom-forming fungus notable for its distinctive fruiting bodies with split gills. It is used as a model organism to study mushroom development, lignocellulose degradation and mating type loci. It is a hypervariable species with considerable genetic and phenotypic diversity between the strains. In this study, we systematically phenotyped 16 dikaryotic strains for aspects of mushroom development and 18 monokaryotic strains for lignocellulose degradation. There was considerable heterogeneity among the strains regarding these phenotypes. The majority of the strains developed mushrooms with varying morphologies, although some strains only grew vegetatively under the tested conditions. Growth on various carbon sources showed strain-specific profiles. The genomes of seven monokaryotic strains were sequenced and analyzed together with six previously published genome sequences. Moreover, the related species Schizophyllum fasciatum was sequenced. Although there was considerable genetic variation between the genome assemblies, the genes related to mushroom formation and lignocellulose degradation were well conserved. These sequenced genomes, in combination with the high phenotypic diversity, will provide a solid basis for functional genomics analyses of the strains of S. commune.

The impact of variations in the three-dimensional structure of the genome has been recognized, but solid cancer tissue studies are limited. Here, we performed integrated deep Hi-C sequencing with matched whole-genome sequencing, whole-genome bisulfite sequencing, 5-hydroxymethylcytosine (5hmC) sequencing and RNA sequencing across a cohort of 80 biopsy samples from patients with metastatic castration-resistant prostate cancer. Dramatic differences were present in gene expression, 5-methylcytosine/5hmC methylation and in structural variation versus mutation rate between A and B (open and closed) chromatin compartments. A subset of tumors exhibited depleted regional chromatin contacts at the AR locus, linked to extrachromosomal circular DNA (ecDNA) and worse response to AR signaling inhibitors. We also identified topological subtypes associated with stark differences in methylation structure, gene expression and prognosis. Our data suggested that DNA interactions may predispose to structural variant formation, exemplified by the recurrent TMPRSS2-ERG fusion. This comprehensive integrated sequencing effort represents a unique clinical tumor resource.

With over 16 million submitted genomic sequences, the SARS-CoV-2 (SC2) virus, the cause of the most recent worldwide COVID-19 pandemic, has become the most sequenced genome of all known viruses, revealing, for example, a vast number of expanding viral lineages. Since the pandemic phase appears to be over, we performed a retrospective re-examination of the demographic grouping pattern and their genomic characteristics during the entire pandemic period up to the peak of the last pandemic wave. For our study, we extracted from the NCBI only unique viral sequences and converted each sequence data to a relational vector, indicating the presence/absence of each variational event compared to a reference sequence. Our study revealed several genomic features that are unexpected or different from those of previous studies. For example, approximately 44,000 variants with unique sequences emerged during the pandemic period; they group into only four major viral-genomic groups and each has a set of mostly unique highly-conserved variant-genotypes (HCVGs); and a small set from the first (ancestral) group was inherited by the three (descendant) groups, suggesting that HCVGs in the next group may be predictable from the current group(s). Such a concept may be potentially important in designing panvalent vaccines against the current and future waves of viral infections.

Identifying and characterizing mobile genetic elements in sequencing data is essential for understanding their diversity, ecology, biotechnological applications and impact on public health. Here we introduce geNomad, a classification and annotation framework that combines information from gene content and a deep neural network to identify sequences of plasmids and viruses. geNomad uses a dataset of more than 200,000 marker protein profiles to provide functional gene annotation and taxonomic assignment of viral genomes. Using a conditional random field model, geNomad also detects proviruses integrated into host genomes with high precision. In benchmarks, geNomad achieved high classification performance for diverse plasmids and viruses (Matthews correlation coefficient of 77.8% and 95.3%, respectively), substantially outperforming other tools. Leveraging geNomad's speed and scalability, we processed over 2.7 trillion base pairs of sequencing data, leading to the discovery of millions of viruses and plasmids that are available through the IMG/VR and IMG/PR databases. geNomad is available at https://portal.nersc.gov/genomad .

The roles of Asgard archaea in eukaryogenesis and marine biogeochemical cycles are well studied, yet their contributions in soil ecosystems remain unknown. Of particular interest are Asgard archaeal contributions to methane cycling in wetland soils. To investigate this, we reconstructed two complete genomes for soil-associated Atabeyarchaeia, a new Asgard lineage, and a complete genome of Freyarchaeia, and predicted their metabolism in situ. Metatranscriptomics reveals expression of genes for [NiFe]-hydrogenases, pyruvate oxidation and carbon fixation via the Wood-Ljungdahl pathway. Also expressed are genes encoding enzymes for amino acid metabolism, anaerobic aldehyde oxidation, hydrogen peroxide detoxification and carbohydrate breakdown to acetate and formate. Overall, soil-associated Asgard archaea are predicted to include non-methanogenic acetogens, highlighting their potential role in carbon cycling in terrestrial environments.

BACKGROUND: The obligate intracellular bacterial family Chlamydiaceae comprises a number of different species that cause disease in various vertebrate hosts including humans. Chlamydia suis, primarily found in the gastrointestinal tract of pigs, is the only species of the Chlamydiaceae family to have naturally gained tetracycline resistance (TetR), through a genomic island (Tet-island), integrated into the middle of chromosomal invasin-like gene inv. Previous studies have hypothesised that the uptake of the Tet-island from a host outside the Chlamydiaceae family was a unique event, followed by spread among C. suis through homologous recombination. In vitro recombination studies have shown that Tet-island exchange between C. suis strains is possible. Our aim in this study was to gain a deeper understanding of the interclade recombination of the Tet-island, among currently circulating C. suis field strains compared to in vitro-generated recombinants, using published whole genome sequences of C. suis field strains (n = 35) and in vitro-generated recombinants (n = 63). RESULTS: We found that the phylogeny of inv better reflected the phylogeny of the Tet-island than that of the whole genome, supporting recombination rather than site-specific insertion as the means of transfer. There were considerable differences between the distribution of recombinations within in vitro-generated strains compared to that within the field strains. These differences are likely because in vitro-generated recombinants were selected for a tetracycline and rifamycin/rifampicin resistant background, leading to the largest peak of recombination across the Tet-island. Finally, we found that interclade recombinations across the Tet-island were more variable in length downstream of the Tet-island than upstream. CONCLUSIONS: Our study supports the hypothesis that the occurrence of TetR strains in both clades of C. suis came about through interclade recombination after a single ancestral horizontal gene transfer event.

Patterns of transcriptional activity are encoded in our genome through regulatory elements such as promoters or enhancers that, paradoxically, contain similar assortments of sequence-specific transcription factor (TF) binding sites1-3. Knowledge of how these sequence motifs encode multiple, often overlapping, gene expression programs is central to understanding gene regulation and how mutations in non-coding DNA manifest in disease4,5. Here, by studying gene regulation from the perspective of individual transcription start sites (TSSs), using natural genetic variation, perturbation of endogenous TF protein levels and massively parallel analysis of natural and synthetic regulatory elements, we show that the effect of TF binding on transcription initiation is position dependent. Analysing TF-binding-site occurrences relative to the TSS, we identified several motifs with highly preferential positioning. We show that these patterns are a combination of a TF's distinct functional profiles-many TFs, including canonical activators such as NRF1, NFY and Sp1, activate or repress transcription initiation depending on their precise position relative to the TSS. As such, TFs and their spacing collectively guide the site and frequency of transcription initiation. More broadly, these findings reveal how similar assortments of TF binding sites can generate distinct gene regulatory outcomes depending on their spatial configuration and how DNA sequence polymorphisms may contribute to transcription variation and disease and underscore a critical role for TSS data in decoding the regulatory information of our genome.

Strains across the Lactobacillaceae family form the basis for a trillion-dollar industry. Our understanding of the genomic basis for their key traits is fragmented, however, including the metabolism that is foundational to their industrial uses. Pangenome analysis of publicly available Lactobacillaceae genomes allowed us to generate genome-scale metabolic network reconstructions for 26 species of industrial importance. Their manual curation led to more than 75,000 gene-protein-reaction associations that were deployed to generate 2,446 genome-scale metabolic models. Cross-referencing genomes and known metabolic traits allowed for manual metabolic network curation and validation of the metabolic models. As a result, we provide the first pangenomic basis for metabolism in the Lactobacillaceae family and a collection of predictive computational metabolic models that enable a variety of practical uses.IMPORTANCELactobacillaceae, a bacterial family foundational to a trillion-dollar industry, is increasingly relevant to biosustainability initiatives. Our study, leveraging approximately 2,400 genome sequences, provides a pangenomic analysis of Lactobacillaceae metabolism, creating over 2,400 curated and validated genome-scale models (GEMs). These GEMs successfully predict (i) unique, species-specific metabolic reactions; (ii) niche-enriched reactions that increase organism fitness; (iii) essential media components, offering insights into the global amino acid essentiality of Lactobacillaceae; and (iv) fermentation capabilities across the family, shedding light on the metabolic basis of Lactobacillaceae-based commercial products. This quantitative understanding of Lactobacillaceae metabolic properties and their genomic basis will have profound implications for the food industry and biosustainability, offering new insights and tools for strain selection and manipulation.

Bacteroides fragilis is a Gram-negative commensal bacterium commonly found in the human colon, which differentiates into two genomospecies termed divisions I and II. Through a comprehensive collection of 694 B. fragilis whole genome sequences, we identify novel features distinguishing these divisions. Our study reveals a distinct geographic distribution with division I strains predominantly found in North America and division II strains in Asia. Additionally, division II strains are more frequently associated with bloodstream infections, suggesting a distinct pathogenic potential. We report differences between the two divisions in gene abundance related to metabolism, virulence, stress response, and colonization strategies. Notably, division II strains harbor more antimicrobial resistance (AMR) genes than division I strains. These findings offer new insights into the functional roles of division I and II strains, indicating specialized niches within the intestine and potential pathogenic roles in extraintestinal sites.ImportanceUnderstanding the distinct functions of microbial species in the gut microbiome is crucial for deciphering their impact on human health. Classifying division II strains as Bacteroides fragilis can lead to erroneous associations, as researchers may mistakenly attribute characteristics observed in division II strains to the more extensively studied division I B. fragilis. Our findings underscore the necessity of recognizing these divisions as separate species with distinct functions. We unveil new findings of differential gene prevalence between division I and II strains in genes associated with intestinal colonization and survival strategies, potentially influencing their role as gut commensals and their pathogenicity in extraintestinal sites. Despite the significant niche overlap and colonization patterns between these groups, our study highlights the complex dynamics that govern strain distribution and behavior, emphasizing the need for a nuanced understanding of these microorganisms.

Drosophila prolongata is a member of the melanogaster species group and rhopaloa subgroup native to the subtropical highlands of southeast Asia. This species exhibits an array of recently evolved male-specific morphological, physiological, and behavioral traits that distinguish it from its closest relatives, making it an attractive model for studying the evolution of sexual dimorphism and testing theories of sexual selection. The lack of genomic resources has impeded the dissection of the molecular basis of sex-specific development and behavior in this species. To address this, we assembled the genome of D. prolongata using long-read sequencing and Hi-C scaffolding, resulting in a highly complete and contiguous (scaffold N50 2.2Mb) genome assembly of 220Mb. The repetitive content of the genome is 24.6%, the plurality of which are LTR retrotransposons (33.2%). Annotations based on RNA-seq data and homology to related species revealed a total of 19,330 genes, of which 16,170 are protein-coding. The assembly includes 98.5% of Diptera BUSCO genes, including 93.8% present as a single copy. Despite some likely regional duplications, the completeness of this genome suggests that it can be readily used for gene expression, GWAS, and other genomic analyses.

Mycena s.s. is a ubiquitous mushroom genus whose members degrade multiple dead plant substrates and opportunistically invade living plant roots. Having sequenced the nuclear genomes of 24 Mycena species, we find them to defy the expected patterns for fungi based on both their traditionally perceived saprotrophic ecology and substrate specializations. Mycena displayed massive genome expansions overall affecting all gene families, driven by novel gene family emergence, gene duplications, enlarged secretomes encoding polysaccharide degradation enzymes, transposable element (TE) proliferation, and horizontal gene transfers. Mainly due to TE proliferation, Arctic Mycena species display genomes of up to 502 Mbp (2-8× the temperate Mycena), the largest among mushroom-forming Agaricomycetes, indicating a possible evolutionary convergence to genomic expansions sometimes seen in Arctic plants. Overall, Mycena show highly unusual, varied mosaic-like genomic structures adaptable to multiple lifestyles, providing genomic illustration for the growing realization that fungal niche adaptations can be far more fluid than traditionally believed.

The jaguar (Panthera onca) is the largest living cat species native to the Americas and one of few large American carnivorans to have survived into the Holocene. However, the extent to which jaguar diversity declined during the end-Pleistocene extinction event remains unclear. For example, Pleistocene jaguar fossils from North America are notably larger than the average extant jaguar, leading to hypotheses that jaguars from this continent represent a now-extinct subspecies (Panthera onca augusta) or species (Panthera augusta). Here, we used a hybridization capture approach to recover an ancient mitochondrial genome from a large, late Pleistocene jaguar from Kingston Saltpeter Cave, Georgia, United States, which we sequenced to 26-fold coverage. We then estimated the evolutionary relationship between the ancient jaguar mitogenome and those from other extinct and living large felids, including multiple jaguars sampled across the species' current range. The ancient mitogenome falls within the diversity of living jaguars. All sampled jaguar mitogenomes share a common mitochondrial ancestor ~400 thousand years ago, indicating that the lineage represented by the ancient specimen dispersed into North America from the south at least once during the late Pleistocene. While genomic data from additional and older specimens will continue to improve understanding of Pleistocene jaguar diversity in the Americas, our results suggest that this specimen falls within the variation of extant jaguars despite the relatively larger size and geographic location and does not represent a distinct taxon.

During the last decade, the spotted wing drosophila, Drosophila suzukii, has spread from eastern Asia to the Americas, Europe, and Africa. This fly attacks many species of cultivated and wild fruits with soft, thin skins, where its serrated ovipositor allows it to lay eggs in undamaged fruit. Parasitoids from the native range of D. suzukii may provide sustainable management of this polyphagous pest. Among these parasitoids, host-specificity testing has revealed a lineage of Ganaspis near brasiliensis, referred to in this paper as G1, that appears to be a cryptic species more host-specific to D. suzukii than other parasitoids. Differentiation among cryptic species is critical for introduction and subsequent evaluation of their impact on D. suzukii. Here, we present results on divergence in genomic sequences and architecture and reproductive isolation between lineages of Ganaspis near brasiliensis that appear to be cryptic species. We studied five populations, two from China, two from Japan, and one from Canada, identified as the G1 vs G3 lineages based on differences in cytochrome oxidase l sequences. We assembled and annotated the genomes of these populations and analyzed divergences in sequence and genome architecture between them. We also report results from crosses to test reproductive compatibility between the G3 lineage from China and the G1 lineage from Japan. The combined results on sequence divergence, differences in genome architectures, ortholog divergence, reproductive incompatibility, differences in host ranges and microhabitat preferences, and differences in morphology show that these lineages are different species. Thus, the decision to evaluate the lineages separately and only import and introduce the more host-specific lineage to North America and Europe was appropriate.

Tandem repeats are frequent across the human genome, and variation in repeat length has been linked to a variety of traits. Recent improvements in long read sequencing technologies have the potential to greatly improve tandem repeat analysis, especially for long or complex repeats. Here, we introduce LongTR, which accurately genotypes tandem repeats from high-fidelity long reads available from both PacBio and Oxford Nanopore Technologies. LongTR is freely available at https://github.com/gymrek-lab/longtr and https://zenodo.org/doi/10.5281/zenodo.11403979 .

Historically neglected by microbial ecologists, soil viruses are now thought to be critical to global biogeochemical cycles. However, our understanding of their global distribution, activities and interactions with the soil microbiome remains limited. Here we present the Global Soil Virus Atlas, a comprehensive dataset compiled from 2,953 previously sequenced soil metagenomes and composed of 616,935 uncultivated viral genomes and 38,508 unique viral operational taxonomic units. Rarefaction curves from the Global Soil Virus Atlas indicate that most soil viral diversity remains unexplored, further underscored by high spatial turnover and low rates of shared viral operational taxonomic units across samples. By examining genes associated with biogeochemical functions, we also demonstrate the viral potential to impact soil carbon and nutrient cycling. This study represents an extensive characterization of soil viral diversity and provides a foundation for developing testable hypotheses regarding the role of the virosphere in the soil microbiome and global biogeochemistry.

Ninu (greater bilby, Macrotis lagotis) are desert-dwelling, culturally and ecologically important marsupials. In collaboration with Indigenous rangers and conservation managers, we generated the Ninu chromosome-level genome assembly (3.66 Gbp) and genome sequences for the extinct Yallara (lesser bilby, Macrotis leucura). We developed and tested a scat single-nucleotide polymorphism panel to inform current and future conservation actions, undertake ecological assessments and improve our understanding of Ninu genetic diversity in managed and wild populations. We also assessed the beneficial impact of translocations in the metapopulation (N = 363 Ninu). Resequenced genomes (temperate Ninu, 6; semi-arid Ninu, 6; and Yallara, 4) revealed two major population crashes during global cooling events for both species and differences in Ninu genes involved in anatomical and metabolic pathways. Despite their 45-year captive history, Ninu have fewer long runs of homozygosity than other larger mammals, which may be attributable to their boom-bust life history. Here we investigated the unique Ninu biology using 12 tissue transcriptomes revealing expression of all 115 conserved eutherian chorioallantoic placentation genes in the uterus, an XY1Y2 sex chromosome system and olfactory receptor gene expansions. Together, we demonstrate the holistic value of genomics in improving key conservation actions, understanding unique biological traits and developing tools for Indigenous rangers to monitor remote wild populations.

Consortia of multicellular magnetotactic bacteria (MMB) are currently the only known example of bacteria without a unicellular stage in their life cycle. Because of their recalcitrance to cultivation, most previous studies of MMB have been limited to microscopic observations. To study the biology of these unique organisms in more detail, we use multiple culture-independent approaches to analyze the genomics and physiology of MMB consortia at single-cell resolution. We separately sequenced the metagenomes of 22 individual MMB consortia, representing 8 new species, and quantified the genetic diversity within each MMB consortium. This revealed that, counter to conventional views, cells within MMB consortia are not clonal. Single consortia metagenomes were then used to reconstruct the species-specific metabolic potential and infer the physiological capabilities of MMB. To validate genomic predictions, we performed stable isotope probing (SIP) experiments and interrogated MMB consortia using fluorescence in situ hybridization (FISH) combined with nanoscale secondary ion mass spectrometry (NanoSIMS). By coupling FISH with bioorthogonal noncanonical amino acid tagging (BONCAT), we explored their in situ activity as well as variation of protein synthesis within cells. We demonstrate that MMB consortia are mixotrophic sulfate reducers and that they exhibit metabolic differentiation between individual cells, suggesting that MMB consortia are more complex than previously thought. These findings expand our understanding of MMB diversity, ecology, genomics, and physiology, as well as offer insights into the mechanisms underpinning the multicellular nature of their unique lifestyle.

The decrease in sequencing expenses has facilitated the creation of reference genomes and proteomes for an expanding array of organisms. Nevertheless, no established repository that details organism-specific genomic and proteomic sequences of specific lengths, referred to as kmers, exists to our knowledge. In this article, we present kmerDB, a database accessible through an interactive web interface that provides kmer-based information from genomic and proteomic sequences in a systematic way. kmerDB currently contains 202,340,859,107 base pairs and 19,304,903,356 amino acids, spanning 54,039 and 21,865 reference genomes and proteomes, respectively, as well as 6,905,362 and 149,305,183 genomic and proteomic species-specific sequences, termed quasi-primes. Additionally, we provide access to 5,186,757 nucleic and 214,904,089 peptide sequences absent from every genome and proteome, termed primes. kmerDB features a user-friendly interface offering various search options and filters for easy parsing and searching. The service is available at: www.kmerdb.com.

Escherichia coli has been attributed to playing a major role in a cascade of events that affect the prevalence and severity of uterine disease in cattle. The objectives of this project were to (i) define the association between the prevalence of specific antimicrobial resistance and virulence factor genes in E. coli with the clinical status related to uterine infection, (ii) identify the genetic relationship between E. coli isolates from cows with diarrhea, with mastitis, and with and without metritis, and (iii) determine the association between the phenotypic and genotypic antimicrobial resistance identified on the E. coli isolated from postpartum cattle. Bacterial isolates (n = 148) were obtained from a larger cross-sectional study. Cows were categorized into one of three clinical groups before enrollment: metritis, cows with purulent discharge, and control cows. For genomic comparison, public genomes (n = 130) from cows with diarrhea, mastitis, and metritis were included in a genome-wide association study, to evaluate differences between the drug classes or the virulence factor category among clinical groups. A distinct E. coli genotype associated with metritis could not be identified. Instead, a high genetic diversity among the isolates from uterine sources was present. A virulence factor previously associated with metritis (fimH) using PCR was not associated with metritis. There was moderate accuracy for whole-genome sequencing to predict phenotypic resistance, which varied depending on the antimicrobial tested. Findings from this study contradict the traditional pathotype classification and the unique intrauterine E. coli genotype associated with metritis in dairy cows.IMPORTANCEMetritis is a common infectious disease in dairy cattle and the second most common reason for treating a cow with antimicrobials. The pathophysiology of the disease is complex and is not completely understood. Specific endometrial pathogenic Escherichia coli have been reported to be adapted to the endometrium and sometimes lead to uterine disease. Unfortunately, the specific genomic details of the endometrial-adapted isolates have not been investigated using enough genomes to represent the genomic diversity of this organism to identify specific virulence genes that are consistently associated with disease development and severity. Results from this study provide key microbial ecological advances by elucidating and challenging accepted concepts for the role of Intrauterine E. coli in metritis in dairy cattle, especially contradicting the existence of a unique intrauterine E. coli genotype associated with metritis in dairy cows, which was not found in our study.

The genome of Candidatus Xenohaliotis californiensis was assembled from shotgun metagenomic sequencing of experimentally infected white abalone. Ninety-one percent genome completeness was achieved with low contamination. Sequencing this genome provides the opportunity to track pathogen evolution over time, conduct gene expression experiments, and study dynamics between this pathogen and its phage.

Cowpea is a highly drought-adapted leguminous crop with great promise for improving agricultural sustainability and food security. Here, we report analyses derived from array-based genotyping of 2,021 accessions constituting a core subset of the worlds largest cowpea collection, held at the International Institute of Tropical Agriculture (IITA) in Ibadan, Nigeria. We used this dataset to examine genetic variation and population structure in worldwide cowpea. We confirm that the primary pattern of population structure is two geographically defined subpopulations originating in West and East Africa, respectively, and that population structure is associated with shifts in phenotypic distribution. Furthermore, we establish the cowpea core collection as a resource for genome-wide association studies by mapping the genetic basis of several phenotypes, with a focus on seed coat pigmentation patterning and color. We anticipate that the genotyped IITA Cowpea Core Collection will serve as a powerful tool for mapping complex traits, facilitating the acceleration of breeding programs to enhance the resilience of this crop in the face of rapid global climate change.

Cotton (Gossypium hirsutum L.) is the key renewable fibre crop worldwide, yet its yield and fibre quality show high variability due to genotype-specific traits and complex interactions among cultivars, management practices and environmental factors. Modern breeding practices may limit future yield gains due to a narrow founding gene pool. Precision breeding and biotechnological approaches offer potential solutions, contingent on accurate cultivar-specific data. Here we address this need by generating high-quality reference genomes for three modern cotton cultivars (UGA230, UA48 and CSX8308) and updating the TM-1 cotton genetic standard reference. Despite hypothesized genetic uniformity, considerable sequence and structural variation was observed among the four genomes, which overlap with ancient and ongoing genomic introgressions from Pima cotton, gene regulatory mechanisms and phenotypic trait divergence. Differentially expressed genes across fibre development correlate with fibre production, potentially contributing to the distinctive fibre quality traits observed in modern cotton cultivars. These genomes and comparative analyses provide a valuable foundation for future genetic endeavours to enhance global cotton yield and sustainability.

We describe the genome of an Eremiobacterota population from tundra soil that contains the minimal set of nif genes needed to fix atmospheric N2. This putative diazotroph population, which we name Candidatus Lamibacter sapmiensis, links for the first time Eremiobacterota and N2 fixation. The integrity of the genome and its nif genes are well supported by both environmental and taxonomic signals. Ca. Lamibacter sapmiensis contains three nifH homologues and the complementary set of nifDKENB genes that are needed to assemble a functional nitrogenase. The putative diazotrophic role of Ca. Lamibacter sapmiensis is supported by the presence of genes that regulate N2 fixation and other genes involved in downstream processes such as ammonia assimilation. Similar to other Eremiobacterota, Ca. Lamibacter sapmiensis encodes the potential for atmospheric chemosynthesis via CO2 fixation coupled with H2 and CO oxidation. Interestingly, the presence of a N2O reductase indicates that this population could play a role as a N2O sink in tundra soils. Due to the lack of activity data, it remains uncertain if Ca. Lamibacter sapmiensis is able to assemble a functional nitrogenase and participate in N2 fixation. Confirmation of this ability would be a testament to the great metabolic versatility of Eremiobacterota, which appears to underlie their ecological success in cold and oligotrophic environments.

BACKGROUND: The colonization of land and the diversification of terrestrial plants is intimately linked to the evolutionary history of their symbiotic fungal partners. Extant representatives of these fungal lineages include mutualistic plant symbionts, the arbuscular mycorrhizal (AM) fungi in Glomeromycota and fine root endophytes in Endogonales (Mucoromycota), as well as fungi with saprotrophic, pathogenic and endophytic lifestyles. These fungal groups separate into three monophyletic lineages but their evolutionary relationships remain enigmatic confounding ancestral reconstructions. Their taxonomic ranks are currently fluid. RESULTS: In this study, we recognize these three monophyletic linages as phyla, and use a balanced taxon sampling and broad taxonomic representation for phylogenomic analysis that rejects a hard polytomy and resolves Glomeromycota as sister to a clade composed of Mucoromycota and Mortierellomycota. Low copy numbers of genes associated with plant cell wall degradation could not be assigned to the transition to a plant symbiotic lifestyle but appears to be an ancestral phylogenetic signal. Both plant symbiotic lineages, Glomeromycota and Endogonales, lack numerous thiamine metabolism genes but the lack of fatty acid synthesis genes is specific to AM fungi. Many genes previously thought to be missing specifically in Glomeromycota are either missing in all analyzed phyla, or in some cases, are actually present in some of the analyzed AM fungal lineages, e.g. the high affinity phosphorus transporter Pho89. CONCLUSION: Based on a broad taxon sampling of fungal genomes we present a well-supported phylogeny for AM fungi and their sister lineages. We show that among these lineages, two independent evolutionary transitions to mutualistic plant symbiosis happened in a genomic background profoundly different from that known from the emergence of ectomycorrhizal fungi in Dikarya. These results call for further reevaluation of genomic signatures associated with plant symbiosis.

Variants in cis-regulatory elements link the noncoding genome to human pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS), enhances noncoding variant analysis by integrating both whole-genome sequencing (WGS) and user-provided functional data. With simplified parameter settings and an efficient multiple testing correction method, CWAS-Plus conducts the CWAS workflow 50 times faster than CWAS, making it more accessible and user-friendly for researchers. Here, we used a single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type-specific enhancers and promoters. Examining autism spectrum disorder WGS data (n = 7280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimers disease WGS data (n = 1087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Pluss utility in genomic disorders and scalability for processing large-scale WGS data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/, respectively.

BackgroundInfection prevention (IP) measures are designed to mitigate the transmission of pathogens in healthcare. Using large-scale viral genomic and social network analyses, we determined if IP measures used during the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic were adequate in protecting healthcare workers (HCWs) and patients from acquiring SARS-CoV-2.MethodsWe performed retrospective cross-sectional analyses of viral genomics from all available SARS-CoV-2 viral samples collected at UC San Diego Health and social network analysis using the electronic medical record to derive temporospatial overlap of infections among related viromes and supplemented with contact tracing data. The outcome measure was any instance of healthcare transmission, defined as cases with closely related viral genomes and epidemiological connection within the healthcare setting during the infection window. Between November 2020 through January 2022, 12 933 viral genomes were obtained from 35 666 patients and HCWs.ResultsAmong 5112 SARS-CoV-2 viral samples sequenced from the second and third waves of SARS-CoV-2 (pre-Omicron), 291 pairs were derived from persons with a plausible healthcare overlap. Of these, 34 pairs (12%) were phylogenetically linked: 19 attributable to household and 14 to healthcare transmission. During the Omicron wave, 2106 contact pairs among 7821 sequences resulted in 120 (6%) related pairs among 32 clusters, of which 10 were consistent with healthcare transmission. Transmission was more likely to occur in shared spaces in the older hospital compared with the newer hospital (2.54 vs 0.63 transmission events per 1000 admissions, P < .001).ConclusionsIP strategies were effective at identifying and preventing healthcare SARS-CoV-2 transmission.

The Yellow Warbler (Setophaga petechia) is a small songbird in the wood-warbler family (Parulidae) that exhibits phenotypic and ecological differences across a widespread distribution and is important to California's riparian habitat conservation. Here, we present a high-quality de novo genome assembly of a vouchered female Yellow Warbler from southern California. Using HiFi long-read and Omni-C proximity sequencing technologies, we generated a 1.22 Gb assembly including 687 scaffolds with a contig N50 of 6.80 Mb, scaffold N50 of 21.18 Mb, and a BUSCO completeness score of 96.0%. This highly contiguous genome assembly provides an essential resource for understanding the history of gene flow, divergence, and local adaptation in Yellow Warblers and can inform conservation management of this charismatic bird species.

Large-genome bacteriophages (jumbo phages) of the proposed family Chimalliviridae assemble a nucleus-like compartment bounded by a protein shell that protects the replicating phage genome from host-encoded restriction enzymes and DNA-targeting CRISPR-Cas nucleases. While the nuclear shell provides broad protection against host nucleases, it necessitates transport of mRNA out of the nucleus-like compartment for translation by host ribosomes, and transport of specific proteins into the nucleus-like compartment to support DNA replication and mRNA transcription. Here, we identify a conserved phage nuclear shell-associated protein that we term Chimallin C (ChmC), which adopts a nucleic acid-binding fold, binds RNA with high affinity in vitro, and binds phage mRNAs in infected cells. ChmC also forms phase-separated condensates with RNA in vitro. Targeted knockdown of ChmC using mRNA-targeting dCas13d results in accumulation of phage-encoded mRNAs in the phage nucleus, reduces phage protein production, and compromises virion assembly. Taken together, our data show that the conserved ChmC protein plays crucial roles in the viral life cycle, potentially by facilitating phage mRNA translocation through the nuclear shell to promote protein production and virion development.

Tanoak (Notholithocarpus densiflorus) is an evergreen tree in the Fagaceae family found in California and southern Oregon. Historically, tanoak acorns were an important food source for Native American tribes, and the bark was used extensively in the leather tanning process. Long considered a disjunct relictual element of the Asian stone oaks (Lithocarpus spp.), phylogenetic analysis has determined that the tanoak is an example of convergent evolution. Tanoaks are deeply divergent from oaks (Quercus) of the Pacific Northwest and comprise a new genus with a single species. These trees are highly susceptible to sudden oak death (SOD), a plant pathogen (Phytophthora ramorum) that has caused widespread deaths of tanoaks. In this study, we set out to assemble the genome and perform comparative studies among a number of individuals that demonstrated varying levels of susceptibility to SOD. First, we sequenced and de novo assembled a draft reference genome of N. densiflorus using cobarcoded library processing methods and an MGI DNBSEQ-G400 sequencer. To increase the contiguity of the final assembly, we also sequenced Oxford Nanopore long reads to 30× coverage. To our knowledge, the draft genome reported here is one of the more contiguous and complete genomes of a tree species published to date, with a contig N50 of ∼1.2 Mb, a scaffold N50 of ∼2.1 Mb, and a complete gene score of 95.5% through BUSCO analysis. In addition, we sequenced 11 genetically distinct individuals and mapped these onto the draft reference genome, enabling the discovery of almost 25 million single nucleotide polymorphisms and ∼4.4 million small insertions and deletions. Finally, using cobarcoded data, we were able to generate a complete haplotype coverage of all 11 genomes.

Whitebark pine (WBP, Pinus albicaulis) is a white pine of subalpine regions in the Western contiguous United States and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola) and additional threats from mountain pine beetle (Dendroctonus ponderosae), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short reads of haploid megagametophyte tissue and Oxford Nanopore long reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gb of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gb). Approximately 87.2% (24.0 Gb) of total sequence was placed on the 12 WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the 3 subclasses of NLRs. Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo-assembled transcriptomes.

The Alliance of Genome Resources (Alliance) is an extensible coalition of knowledgebases focused on the genetics and genomics of intensively studied model organisms. The Alliance is organized as individual knowledge centers with strong connections to their research communities and a centralized software infrastructure, discussed here. Model organisms currently represented in the Alliance are budding yeast, Caenorhabditis elegans, Drosophila, zebrafish, frog, laboratory mouse, laboratory rat, and the Gene Ontology Consortium. The project is in a rapid development phase to harmonize knowledge, store it, analyze it, and present it to the community through a web portal, direct downloads, and application programming interfaces (APIs). Here, we focus on developments over the last 2 years. Specifically, we added and enhanced tools for browsing the genome (JBrowse), downloading sequences, mining complex data (AllianceMine), visualizing pathways, full-text searching of the literature (Textpresso), and sequence similarity searching (SequenceServer). We enhanced existing interactive data tables and added an interactive table of paralogs to complement our representation of orthology. To support individual model organism communities, we implemented species-specific landing pages and will add disease-specific portals soon; in addition, we support a common community forum implemented in Discourse software. We describe our progress toward a central persistent database to support curation, the data modeling that underpins harmonization, and progress toward a state-of-the-art literature curation system with integrated artificial intelligence and machine learning (AI/ML).

Vitamin A deficiency remains prevalent on a global scale, including in regions where maize constitutes a high percentage of human diets. One solution for alleviating this deficiency has been to increase grain concentrations of provitamin A carotenoids in maize (Zea mays ssp. mays L.)-an example of biofortification. The International Maize and Wheat Improvement Center (CIMMYT) developed a Carotenoid Association Mapping panel of 380 inbred lines adapted to tropical and subtropical environments that have varying grain concentrations of provitamin A and other health-beneficial carotenoids. Several major genes have been identified for these traits, 2 of which have particularly been leveraged in marker-assisted selection. This project assesses the predictive ability of several genomic prediction strategies for maize grain carotenoid traits within and between 4 environments in Mexico. Ridge Regression-Best Linear Unbiased Prediction, Elastic Net, and Reproducing Kernel Hilbert Spaces had high predictive abilities for all tested traits (β-carotene, β-cryptoxanthin, provitamin A, lutein, and zeaxanthin) and outperformed Least Absolute Shrinkage and Selection Operator. Furthermore, predictive abilities were higher when using genome-wide markers rather than only the markers proximal to 2 or 13 genes. These findings suggest that genomic prediction models using genome-wide markers (and assuming equal variance of marker effects) are worthwhile for these traits even though key genes have already been identified, especially if breeding for additional grain carotenoid traits alongside β-carotene. Predictive ability was maintained for all traits except lutein in between-environment prediction. The TASSEL (Trait Analysis by aSSociation, Evolution, and Linkage) Genomic Selection plugin performed as well as other more computationally intensive methods for within-environment prediction. The findings observed herein indicate the utility of genomic prediction methods for these traits and could inform their resource-efficient implementation in biofortification breeding programs.

The coppery titi monkey (Plecturocebus cupreus) is an emerging nonhuman primate model system for behavioral and neurobiological research. At the same time, the almost entire absence of genomic resources for the species has hampered insights into the genetic underpinnings of the phenotypic traits of interest. To facilitate future genotype-to-phenotype studies, we here present a high-quality, fully annotated de novo genome assembly for the species with chromosome-length scaffolds spanning the autosomes and chromosome X (scaffold N50 = 130.8 Mb), constructed using data obtained from several orthologous short- and long-read sequencing and scaffolding techniques. With a base-level accuracy of ∼99.99% in chromosome-length scaffolds as well as benchmarking universal single-copy ortholog and k-mer completeness scores of >99.0% and 95.1% at the genome level, this assembly represents one of the most complete Pitheciidae genomes to date, making it an invaluable resource for comparative evolutionary genomics research to improve our understanding of lineage-specific changes underlying adaptive traits as well as deleterious mutations associated with disease.

Health equity is the state in which everyone has fair and just opportunities to attain their highest level of health. The field of human genomics has fallen short in increasing health equity, largely because the diversity of the human population has been inadequately reflected among participants of genomics research. This lack of diversity leads to disparities that can have scientific and clinical consequences. Achieving health equity related to genomics will require greater effort in addressing inequities within the field. As part of the commitment of the National Human Genome Research Institute (NHGRI) to advancing health equity, it convened experts in genomics and health equity research to make recommendations and performed a review of current literature to identify the landscape of gaps and opportunities at the interface between human genomics and health equity research. This Perspective describes these findings and examines health equity within the context of human genomics and genomic medicine.

Zygnematophyceae are the algal sisters of land plants. Here we sequenced four genomes of filamentous Zygnematophyceae, including chromosome-scale assemblies for three strains of Zygnema circumcarinatum. We inferred traits in the ancestor of Zygnematophyceae and land plants that might have ushered in the conquest of land by plants: expanded genes for signaling cascades, environmental response, and multicellular growth. Zygnematophyceae and land plants share all the major enzymes for cell wall synthesis and remodifications, and gene gains shaped this toolkit. Co-expression network analyses uncover gene cohorts that unite environmental signaling with multicellular developmental programs. Our data shed light on a molecular chassis that balances environmental response and growth modulation across more than 600 million years of streptophyte evolution.

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mitochondrial genome in different tissues throughout ageing. We used ultrasensitive duplex sequencing to profile ~2.5 million mitochondrial genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloguing ~1.2 million mitochondrial somatic and ultralow-frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the light strand origin of replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared with primates with a surfeit of reactive oxygen species-associated G > T/C > A mutations, and that somatic mutations in protein-coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that re-align mito-nuclear ancestry within an organisms lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

Plants rely on Nucleotide-binding, Leucine-rich repeat Receptors (NLRs) for pathogen recognition. Highly variable NLRs (hvNLRs) show remarkable intraspecies diversity, while their low-variability paralogs (non-hvNLRs) are conserved between ecotypes. At a population level, hvNLRs provide new pathogen-recognition specificities, but the association between allelic diversity and genomic and epigenomic features has not been established. Our investigation of NLRs in Arabidopsis Col-0 has revealed that hvNLRs show higher expression, less gene body cytosine methylation, and closer proximity to transposable elements than non-hvNLRs. hvNLRs show elevated synonymous and nonsynonymous nucleotide diversity and are in chromatin states associated with an increased probability of mutation. Diversifying selection maintains variability at a subset of codons of hvNLRs, while purifying selection maintains conservation at non-hvNLRs. How these features are established and maintained, and whether they contribute to the observed diversity of hvNLRs is key to understanding the evolution of plant innate immune receptors.

Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium-an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

The genus Clostridium is a large and diverse group within the Bacillota (formerly Firmicutes), whose members can encode useful complex traits such as solvent production, gas-fermentation, and lignocellulose breakdown. We describe 270 genome sequences of solventogenic clostridia from a comprehensive industrial strain collection assembled by Professor David Jones that includes 194 C. beijerinckii, 57 C. saccharobutylicum, 4 C. saccharoperbutylacetonicum, 5 C. butyricum, 7 C. acetobutylicum, and 3 C. tetanomorphum genomes. We report methods, analyses and characterization for phylogeny, key attributes, core biosynthetic genes, secondary metabolites, plasmids, prophage/CRISPR diversity, cellulosomes and quorum sensing for the 6 species. The expanded genomic data described here will facilitate engineering of solvent-producing clostridia as well as non-model microorganisms with innately desirable traits. Sequences could be applied in conventional platform biocatalysts such as yeast or Escherichia coli for enhanced chemical production. Recently, gene sequences from this collection were used to engineer Clostridium autoethanogenum, a gas-fermenting autotrophic acetogen, for continuous acetone or isopropanol production, as well as butanol, butanoic acid, hexanol and hexanoic acid production.

Sugarcane, the world's most harvested crop by tonnage, has shaped global history, trade and geopolitics, and is currently responsible for 80% of sugar production worldwide1. While traditional sugarcane breeding methods have effectively generated cultivars adapted to new environments and pathogens, sugar yield improvements have recently plateaued2. The cessation of yield gains may be due to limited genetic diversity within breeding populations, long breeding cycles and the complexity of its genome, the latter preventing breeders from taking advantage of the recent explosion of whole-genome sequencing that has benefited many other crops. Thus, modern sugarcane hybrids are the last remaining major crop without a reference-quality genome. Here we take a major step towards advancing sugarcane biotechnology by generating a polyploid reference genome for R570, a typical modern cultivar derived from interspecific hybridization between the domesticated species (Saccharum officinarum) and the wild species (Saccharum spontaneum). In contrast to the existing single haplotype ('monoploid') representation of R570, our 8.7 billion base assembly contains a complete representation of unique DNA sequences across the approximately 12 chromosome copies in this polyploid genome. Using this highly contiguous genome assembly, we filled a previously unsized gap within an R570 physical genetic map to describe the likely causal genes underlying the single-copy Bru1 brown rust resistance locus. This polyploid genome assembly with fine-grain descriptions of genome architecture and molecular targets for biotechnology will help accelerate molecular and transgenic breeding and adaptation of sugarcane to future environmental conditions.

The fishing cat, Prionailurus viverrinus, faces a population decline, increasing the importance of maintaining healthy zoo populations. Unfortunately, zoo-managed individuals currently face a high prevalence of transitional cell carcinoma (TCC), a form of bladder cancer. To investigate the genetics of inherited diseases among captive fishing cats, we present a chromosome-scale assembly, generate the pedigree of the zoo-managed population, reaffirm the close genetic relationship with the Asian leopard cat (Prionailurus bengalensis), and identify 7.4 million single nucleotide variants (SNVs) and 23,432 structural variants (SVs) from whole genome sequencing (WGS) data of healthy and TCC cats. Only BRCA2 was found to have a high recurrent number of missense mutations in fishing cats diagnosed with TCC when compared to inherited human cancer risk variants. These new fishing cat genomic resources will aid conservation efforts to improve their genetic fitness and enhance the comparative study of feline genomes.

Reference genome assemblies have been created from multiple lineages within the Canidae family; however, despite its phylogenetic relevance as a basal genus within the clade, there is currently no reference genome for the gray fox (Urocyon cinereoargenteus). Here, we present a chromosome-level assembly for the gray fox (U. cinereoargenteus), which represents the most contiguous, non-domestic canid reference genome available to date, with 90% of the genome contained in just 34 scaffolds and a contig N50 and scaffold N50 of 59.4 and 72.9 Megabases, respectively. Repeat analyses identified an increased number of simple repeats relative to other canids. Based on mitochondrial DNA, our Vermont sample clusters with other gray fox samples from the northeastern United States and contains slightly lower levels of heterozygosity than gray foxes on the west coast of California. This new assembly lays the groundwork for future studies to describe past and present population dynamics, including the delineation of evolutionarily significant units of management relevance. Importantly, the phylogenetic position of Urocyon allows us to verify the loss of PRDM9 functionality in the basal canid lineage, confirming that pseudogenization occurred at least 10 million years ago.

Duckweeds are among the fastest reproducing plants, able to clonally divide at exponential rates. However, the genetic and epigenetic impact of clonality on plant genomes is poorly understood. 5-methylcytosine (5mC) is a modified base often described as necessary for the proper regulation of certain genes and transposons and for the maintenance of genome integrity in plants. However, the extent of this dogma is limited by the current phylogenetic sampling of land plant species diversity. Here we analyzed DNA methylomes, small RNAs, mRNA-seq, and H3K9me2 histone modification for Spirodela polyrhiza. S. polyrhiza has lost highly conserved genes involved in de novo methylation of DNA at sites often associated with repetitive DNA, and within genes, however, symmetrical DNA methylation and heterochromatin are maintained during cell division at certain transposons and repeats. Consequently, small RNAs that normally guide methylation to silence repetitive DNA like retrotransposons are diminished. Despite the loss of a highly conserved methylation pathway, and the reduction of small RNAs that normally target repetitive DNA, transposons have not proliferated in the genome, perhaps due in part to the rapid, clonal growth lifestyle of duckweeds.

The genus Xanthomonas has been primarily studied for pathogenic interactions with plants. However, besides host and tissue-specific pathogenic strains, this genus also comprises nonpathogenic strains isolated from a broad range of hosts, sometimes in association with pathogenic strains, and other environments, including rainwater. Based on their incapacity or limited capacity to cause symptoms on the host of isolation, nonpathogenic xanthomonads can be further characterized as commensal and weakly pathogenic. This study aimed to understand the diversity and evolution of nonpathogenic xanthomonads compared to their pathogenic counterparts based on their cooccurrence and phylogenetic relationship and to identify genomic traits that form the basis of a life history framework that groups xanthomonads by ecological strategies. We sequenced genomes of 83 strains spanning the genus phylogeny and identified eight novel species, indicating unexplored diversity. While some nonpathogenic species have experienced a recent loss of a type III secretion system, specifically the hrp2 cluster, we observed an apparent lack of association of the hrp2 cluster with lifestyles of diverse species. We performed association analysis on a large data set of 337 Xanthomonas strains to explain how xanthomonads may have established association with the plants across the continuum of lifestyles from commensals to weak pathogens to pathogens. Presence of distinct transcriptional regulators, distinct nutrient utilization and assimilation genes, transcriptional regulators, and chemotaxis genes may explain lifestyle-specific adaptations of xanthomonads.