Identifying gene expression programs in single-cell RNA-seq data using linear correlation explanation

OBJECTIVE

Gene expression analysis through single-cell RNA sequencing (scRNA-seq) has revolutionized our understanding of gene regulation in diverse cell types, tissues, and organisms. While existing methods primarily focus on identifying cell type-specific gene expression programs (GEPs), the characterization of GEPs associated with biological processes and stimuli responses remains limited. In this study, we aim to infer biologically meaningful GEPs that are associated with both cellular phenotypes and activity programs directly from scRNA-seq data.

METHODS

We applied linear CorEx, a machine-learning-based approach, to infer GEPs by grouping genes based on total correlation optimization function in simulated and real-world scRNA-seq datasets. Additionally, we utilized a transfer learning approach to project CorEx-inferred GEPs to other scRNA-seq datasets.

RESULTS

By leveraging total correlation optimization, linear CorEx groups genes and demonstrates superior performance in identifying cell types and activity programs compared to similar methods using simulated data. Furthermore, we apply this same approach to real-world scRNA-seq data from the mouse dentate gyrus and embryonic colon development, uncovering biologically relevant GEPs related to cell types, developmental ages, and cell cycle programs. We also demonstrate the potential for transfer learning by evaluating similar datasets, showcasing the cross-species sensitivity of linear CorEx.

CONCLUSION

Our findings validate linear CorEx as a valuable tool for comprehensively analyzing complex signals in scRNA-seq data, leading to deeper insights into gene expression dynamics, cellular heterogeneity, and regulatory mechanisms.

A chromosome-scale fishing cat reference genome for the evaluation of potential germline risk variants

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.

Structurally divergent and recurrently mutated regions of primate genomes

We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or ∼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.

Single-cell insights into development of the bovine placenta†

A central determinant of pregnancy success is proper development of the conceptus (embryo/fetus and associated extraembryonic membranes including the placenta). Although the gross morphology and histology of the bovine placenta have been well studied, the cellular and molecular mechanisms regulating placenta development and trophoblast differentiation and function remain essentially undefined. Here, single-cell transcriptome (scRNA-seq) analysis was performed on the day 17 bovine conceptus and chorion of day 24, 30, and 50 conceptuses (n = 3-4 samples per day) using the 10X Genomics platform. Bioinformatic analyses identified cell types and their ontogeny including trophoblast, mesenchyme, and immune cells. Loss of interferon tau-expressing trophoblast uninucleate cells occurred between days 17 and 30, whereas binucleate cells, identified based on expression of placental lactogen (CSH2) and specific pregnancy-associated glycoprotein genes (PAGs), first appeared on day 24. Several different types of uninucleate cells were present in day 24, 30, and 50 samples, but only one (day 24) or two types of binucleate cells (days 30 and 50). Cell trajectory analyses provided a conceptual framework for uninucleate cell development and binucleate cell differentiation, and bioinformatic analyses identified candidate transcription factors governing differentiation and function of the trophoblasts. The digital atlas of cell types in the developing bovine conceptus reported here serves as a resource to discover key genes and biological pathways regulating its development during the critical periods of implantation and placentation.

Phylogenomics analyses of all species of Swordtails (Genus Xiphophorus ) highlights hybridization precedes speciation

Hybridization has been recognized as an important driving force for evolution, however studies of the genetic consequence and its cause are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the central American genus Xiphophorus were proposed to have evolved with multiple ancient and ongoing hybridization events, and served as a valuable research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genome resource and its annotation of all 26 Xiphophorus species. On this dataset we resolved the so far conflicting phylogeny. Through comparative genomic analyses we investigated the molecular evolution of genes related to melanoma, for a main sexually selected trait and for the genetic control of puberty timing, which are predicted to be involved in pre-and postzygotic isolation and thus to influence the probability of interspecific hybridization in Xiphophorus . We demonstrate dramatic size-variation of some gene families across species, despite the reticulate evolution and short divergence time. Finally, we clarify the hybridization history in the genus Xiphophorus genus, settle the long dispute on the hybridization origin of two Southern swordtails, highlight hybridizations precedes speciation, and reveal the distribution of hybridization ancestry remaining in the fused genome.

A second unveiling: Haplotig masking of the eastern oyster genome improves population-level inference

Genome assembly can be challenging for species that are characterized by high amounts of polymorphism, heterozygosity, and large effective population sizes. High levels of heterozygosity can result in genome mis-assemblies and a larger than expected genome size due to the haplotig versions of a single locus being assembled as separate loci. Here, we describe the first chromosome-level genome for the eastern oyster, Crassostrea virginica. Publicly released and annotated in 2017, the assembly has a scaffold N50 of 54 mb and is over 97.3% complete based on BUSCO analysis. The genome assembly for the eastern oyster is a critical resource for foundational research into molluscan adaptation to a changing environment and for selective breeding for the aquaculture industry. Subsequent resequencing data suggested the presence of haplotigs in the original assembly, and we developed a post hoc method to break up chimeric contigs and mask haplotigs in published heterozygous genomes and evaluated improvements to the accuracy of downstream analysis. Masking haplotigs had a large impact on SNP discovery and estimates of nucleotide diversity and had more subtle and nuanced effects on estimates of heterozygosity, population structure analysis, and outlier detection. We show that haplotig masking can be a powerful tool for improving genomic inference, and we present an open, reproducible resource for the masking of haplotigs in any published genome.

A pangenome graph reference of 30 chicken genomes allows genotyping of large and complex structural variants

BACKGROUND

The red junglefowl, the wild outgroup of domestic chickens, has historically served as a reference for genomic studies of domestic chickens. These studies have provided insight into the etiology of traits of commercial importance. However, the use of a single reference genome does not capture diversity present among modern breeds, many of which have accumulated molecular changes due to drift and selection. While reference-based resequencing is well-suited to cataloging simple variants such as single-nucleotide changes and short insertions and deletions, it is mostly inadequate to discover more complex structural variation in the genome.

METHODS

We present a pangenome for the domestic chicken consisting of thirty assemblies of chickens from different breeds and research lines.

RESULTS

We demonstrate how this pangenome can be used to catalog structural variants present in modern breeds and untangle complex nested variation. We show that alignment of short reads from 100 diverse wild and domestic chickens to this pangenome reduces reference bias by 38%, which affects downstream genotyping results. This approach also allows for the accurate genotyping of a large and complex pair of structural variants at the K feathering locus using short reads, which would not be possible using a linear reference.

CONCLUSIONS

We expect that this new paradigm of genomic reference will allow better pinpointing of exact mutations responsible for specific phenotypes, which will in turn be necessary for breeding chickens that meet new sustainability criteria and are resilient to quickly evolving pathogen threats.

Astyanax mexicanus surface and cavefish chromosome-scale assemblies for trait variation discovery

The ability of organisms to adapt to sudden extreme environmental changes produces some of the most drastic examples of rapid phenotypic evolution. The Mexican Tetra, Astyanax mexicanus, is abundant in the surface waters of northeastern Mexico, but repeated colonizations of cave environments have resulted in the independent evolution of troglomorphic phenotypes in several populations. Here, we present three chromosome-scale assemblies of this species, for one surface and two cave populations, enabling the first whole-genome comparisons between independently evolved cave populations to evaluate the genetic basis for the evolution of adaptation to the cave environment. Our assemblies represent the highest quality of sequence completeness with predicted protein-coding and non-coding gene metrics far surpassing prior resources and, to our knowledge, all long-read assembled teleost genomes, including zebrafish. Whole genome synteny alignments show highly conserved gene order among cave forms in contrast to a higher number of chromosomal rearrangements when compared to other phylogenetically close or distant teleost species. By phylogenetically assessing gene orthology across distant branches of amniotes, we discover gene orthogroups unique to A. mexicanus. When compared to a representative surface fish genome, we find a rich amount of structural sequence diversity, defined here as the number and size of insertions and deletions as well as expanding and contracting repeats across cave forms. These new more complete genomic resources ensure higher trait resolution for comparative, functional, developmental, and genetic studies of drastic trait differences within a species.

Concurrent invasions of European starlings in Australia and North America reveal population-specific differentiation in shared genomic regions

A species' success during the invasion of new areas hinges on an interplay between the demographic processes common to invasions and the specific ecological context of the novel environment. Evolutionary genetic studies of invasive species can investigate how genetic bottlenecks and ecological conditions shape genetic variation in invasions, and our study pairs two invasive populations that are hypothesized to be from the same source population to compare how each population evolved during and after introduction. Invasive European starlings (Sturnus vulgaris) established populations in both Australia and North America in the 19th century. Here, we compare whole-genome sequences among native and independently introduced European starling populations to determine how demographic processes interact with rapid evolution to generate similar genetic patterns in these recent and replicated invasions. Demographic models indicate that both invasive populations experienced genetic bottlenecks as expected based on invasion history, and we find that specific genomic regions have differentiated even on this short evolutionary timescale. Despite genetic bottlenecks, we suggest that genetic drift alone cannot explain differentiation in at least two of these regions. The demographic boom intrinsic to many invasions as well as potential inversions may have led to high population-specific differentiation, although the patterns of genetic variation are also consistent with the hypothesis that this infamous and highly mobile invader adapted to novel selection (e.g., extrinsic factors). We use targeted sampling of replicated invasions to identify and evaluate support for multiple, interacting evolutionary mechanisms that lead to differentiation during the invasion process.

Single-haplotype comparative genomics provides insights into lineage-specific structural variation during cat evolution

The role of structurally dynamic genomic regions in speciation is poorly understood due to challenges inherent in diploid genome assembly. Here we reconstructed the evolutionary dynamics of structural variation in five cat species by phasing the genomes of three interspecies F1 hybrids to generate near-gapless single-haplotype assemblies. We discerned that cat genomes have a paucity of segmental duplications relative to great apes, explaining their remarkable karyotypic stability. X chromosomes were hotspots of structural variation, including enrichment with inversions in a large recombination desert with characteristics of a supergene. The X-linked macrosatellite DXZ4 evolves more rapidly than 99.5% of the genome clarifying its role in felid hybrid incompatibility. Resolved sensory gene repertoires revealed functional copy number changes associated with ecomorphological adaptations, sociality and domestication. This study highlights the value of gapless genomes to reveal structural mechanisms underpinning karyotypic evolution, reproductive isolation and ecological niche adaptation.

Feline precision medicine using whole-exome sequencing identifies a novel frameshift mutation for vitamin D-dependent rickets type 2

OBJECTIVES

A 14-week-old female domestic longhair kitten presented with shifting lameness and disproportionately smaller size compared with a co-housed littermate.

METHODS

Hematology and serum biochemical testing were conducted to investigate causes for delayed growth, and radiographs of the appendicular skeleton were obtained.

RESULTS

The afflicted kitten had marked hypocalcemia, mild hypophosphatemia and substantial elevations in alkaline phosphatase activity, as well as pathognomonic radiographic findings consistent with rickets. Skeletal changes and hypocalcemia prompted testing of concentrations of parathyroid hormone (PTH) and vitamin D metabolites. Endocrine testing demonstrated significant increases in serum concentrations of PTH and 1,25-dihydroxycholecalciferol (calcitriol), supporting a diagnosis of vitamin D-dependent rickets type 2. Provision of analgesia, supraphysiologic doses of calcitriol and calcium carbonate supplementation achieved normalization of the serum calcium concentration and restoration of normal growth, although some skeletal abnormalities persisted. Once skeletally mature, ongoing calcitriol supplementation was not required. Whole-exome sequencing (WES) was conducted to identify the underlying DNA variant. A cytosine deletion at cat chromosome position B4:76777621 in VDR (ENSFCAT00000029466:c.106delC) was identified and predicted to cause a stop codon in exon 2 (p.Arg36Glufs*18), disrupting >90% of the receptor. The variant was unique and homozygous in this patient and absent in the sibling and approximately 400 other cats for which whole-genome and whole-exome data were available.

CONCLUSIONS AND RELEVANCE

A unique, heritable form of rickets was diagnosed in a domestic longhair cat. WES identified a novel frameshift mutation affecting the gene coding for the vitamin D3 receptor, determining the likely causal genetic variant. Precision medicine techniques, including whole-exome and whole-genome sequencing, can be a standard of care in cats to identify disease etiologies, and to target therapeutics and personalize treatment.

Genomic landscape and gene expression profiles of feline oral squamous cell carcinoma

Feline oral squamous cell carcinoma (FOSCC) is a cancer of the squamous cell lining in the oral cavity and represents up to 80% of all oral cancers in cats, with a poor prognosis. We have used whole exome sequencing (WES) and RNA sequencing of the tumor to discover somatic mutations and gene expression changes that may be associated with FOSCC occurrence. FOSCC offers a potential comparative model to study human head and neck squamous cell carcinoma (HNSCC) due to its similar spontaneous formation, and morphological and histological features. In this first study using WES to identify somatic mutations in feline cancer, we have identified tumor-associated gene mutations in six cats with FOSCC and found some overlap with identified recurrently mutated genes observed in HNSCC. Four samples each had mutations in TP53, a common mutation in all cancers, but each was unique. Mutations in other cellular growth control genes were also found such as KAT2B and ARID1A. Enrichment analysis of FOSCC gene expression profiles suggests a molecular similarity to human OSCC as well, including alterations in epithelial to mesenchymal transition and IL6/JAK/STAT pathways. In this preliminary study, we present exome and transcriptome results that further our understanding of FOSCC.

Genomic analysis of two phlebotomine sand fly vectors of leishmania from the new and old World

Phlebotomine sand flies are of global significance as important vectors of human disease, transmitting bacterial, viral, and protozoan pathogens, including the kinetoplastid parasites of the genus Leishmania, the causative agents of devastating diseases collectively termed leishmaniasis. More than 40 pathogenic Leishmania species are transmitted to humans by approximately 35 sand fly species in 98 countries with hundreds of millions of people at risk around the world. No approved efficacious vaccine exists for leishmaniasis and available therapeutic drugs are either toxic and/or expensive, or the parasites are becoming resistant to the more recently developed drugs. Therefore, sand fly and/or reservoir control are currently the most effective strategies to break transmission. To better understand the biology of sand flies, including the mechanisms involved in their vectorial capacity, insecticide resistance, and population structures we sequenced the genomes of two geographically widespread and important sand fly vector species: Phlebotomus papatasi, a vector of Leishmania parasites that cause cutaneous leishmaniasis, (distributed in Europe, the Middle East and North Africa) and Lutzomyia longipalpis, a vector of Leishmania parasites that cause visceral leishmaniasis (distributed across Central and South America). We categorized and curated genes involved in processes important to their roles as disease vectors, including chemosensation, blood feeding, circadian rhythm, immunity, and detoxification, as well as mobile genetic elements. We also defined gene orthology and observed micro-synteny among the genomes. Finally, we present the genetic diversity and population structure of these species in their respective geographical areas. These genomes will be a foundation on which to base future efforts to prevent vector-borne transmission of Leishmania parasites.

The immune cell landscape and response of Marek's disease resistant and susceptible chickens infected with Marek's disease virus

Genetically resistant or susceptible chickens to Marek's disease (MD) have been widely used models to identify the molecular determinants of these phenotypes. However, these prior studies lacked the basic identification and understanding of immune cell types that could be translated toward improved MD control. To gain insights into specific immune cell types and their responses to Marek's disease virus (MDV) infection, we used single-cell RNA sequencing (scRNAseq) on splenic cells from MD resistant and susceptible birds. In total, 14,378 cells formed clusters that identified various immune cell types. Lymphocytes, specifically T cell subtypes, were the most abundant with significant proportional changes in some subtypes upon infection. The largest number of differentially expressed genes (DEG) response was seen in granulocytes, while macrophage DEGs differed in directionality by subtype and line. Among the most DEG in almost all immune cell types were granzyme and granulysin, both associated with cell-perforating processes. Protein interactive network analyses revealed multiple overlapping canonical pathways within both lymphoid and myeloid cell lineages. This initial estimation of the chicken immune cell type landscape and its accompanying response will greatly aid efforts in identifying specific cell types and improving our knowledge of host response to viral infection.

High resolution genomes of multiple Xiphophorus species provide new insights into microevolution, hybrid incompatibility, and epistasis

Because of diverged adaptative phenotypes, fish species of the genus Xiphophorus have contributed to a wide range of research for a century. Existing Xiphophorus genome assemblies are not at the chromosomal level and are prone to sequence gaps, thus hindering advancement of the intra- and inter-species differences for evolutionary, comparative, and translational biomedical studies. Herein, we assembled high-quality chromosome-level genome assemblies for three distantly related Xiphophorus species, namely, X. maculatus, X. couchianus, and X. hellerii Our overall goal is to precisely assess microevolutionary processes in the clade to ascertain molecular events that led to the divergence of the Xiphophorus species and to progress understanding of genetic incompatibility to disease. In particular, we measured intra- and inter-species divergence and assessed gene expression dysregulation in reciprocal interspecies hybrids among the three species. We found expanded gene families and positively selected genes associated with live bearing, a special mode of reproduction. We also found positively selected gene families are significantly enriched in nonpolymorphic transposable elements, suggesting the dispersal of these nonpolymorphic transposable elements has accompanied the evolution of the genes, possibly by incorporating new regulatory elements in support of the Britten-Davidson hypothesis. We characterized inter-specific polymorphisms, structural variants, and polymorphic transposable element insertions and assessed their association to interspecies hybridization-induced gene expression dysregulation related to specific disease states in humans.

Single-nuclei RNA sequencing (snRNA-seq) uncovers trophoblast cell types and lineages in the mature bovine placenta

Ruminants have a semi-invasive placenta, which possess highly vascularized placentomes formed by maternal endometrial caruncles and fetal placental cotyledons and required for fetal development to term. The synepitheliochorial placenta of cattle contains at least two trophoblast cell populations, including uninucleate (UNC) and binucleate (BNC) cells that are most abundant in the cotyledonary chorion of the placentomes. The interplacentomal placenta is more epitheliochorial in nature with the chorion developing specialized areolae over the openings of uterine glands. Of note, the cell types in the placenta and cellular and molecular mechanisms governing trophoblast differentiation and function are little understood in ruminants. To fill this knowledge gap, the cotyledonary and intercotyledonary areas of the mature day 195 bovine placenta were analyzed by single nuclei analysis. Single-nuclei RNA-seq analysis found substantial differences in cell type composition and transcriptional profiles between the two distinct regions of the placenta. Based on clustering and cell marker gene expression, five different trophoblast cell types were identified in the chorion, including proliferating and differentiating UNC and two different types of BNC in the cotyledon. Cell trajectory analyses provided a framework for understanding the differentiation of trophoblast UNC into BNC. The upstream transcription factor binding analysis of differentially expressed genes identified a candidate set of regulator factors and genes regulating trophoblast differentiation. This foundational information is useful to discover essential biological pathways underpinning the development and function of the bovine placenta.

Structurally divergent and recurrently mutated regions of primate genomes

To better understand the pattern of primate genome structural variation, we sequenced and assembled using multiple long-read sequencing technologies the genomes of eight nonhuman primate species, including New World monkeys (owl monkey and marmoset), Old World monkey (macaque), Asian apes (orangutan and gibbon), and African ape lineages (gorilla, bonobo, and chimpanzee). Compared to the human genome, we identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. Across 50 million years of primate evolution, we estimate that 819.47 Mbp or ~27% of the genome has been affected by SVs based on analysis of these primate lineages. We identify 1,607 structurally divergent regions (SDRs) wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (CARDs, ABCD7, OLAH) and new lineage-specific genes are generated (e.g., CKAP2, NEK5) and have become targets of rapid chromosomal diversification and positive selection (e.g., RGPDs). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species for the first time.

Abstract B57: Studying the role of OPN-CD44 interaction in colorectal cancer antitumor immunity suppression using scRNA-seq data

Colorectal cancer (CRC) remains one of the deadliest malignancies worldwide despite recent progress in screening and treatment strategies. Osteopontin (OPN) is a multifunctional protein expressed by multiple cells in the tumor microenvironment (TME) and is associated with immune suppression and poor prognosis in CRC. Tumor-associated macrophages (TAM), are also associated with immune suppression and poor prognosis in CRC, and may express OPN. While both TAM and OPN overlap in the TME, the role of OPN and OPN expressing TAM is poorly understood. In this study, we show the critical interaction and strategy for targeting this axis to overcome therapeutic resistance in CRC. Five publicly available CRC scRNA-seq data (571,818 cells total) from both untreated and treated patients were analyzed to determine the OPN expression and OPN-CD44 interaction in primary CRC and liver metastases using CellChat ligand-receptor interaction analysis. For experimental analyses, patient tissue was collected from the University of Missouri (MU) and the Harry S. Truman VA (HSTVA) on an approved IRB. OPN-positive TAMs were determined by immunofluorescence (IFC) and flow cytometry (FC). Analysis of FC data was completed using analyzed (FlowJo). Wilcoxon rank sum test was performed to find differentially expressed genes between the samples with and without OPN-CD44 interaction in each dataset. Data from scRNA-seq, IFC, and FC all demonstrated that TAMs were the predominant cells expressing OPN in the TME and the majority of TAM are OPN+. Most interestingly, in scRNA-seq data, we found minimal OPN+TAM in primary tumors from treated patients. These patients had a good response to therapy. Conversely, in FC data from our own patients, we found OPN+ TAM were preserved in primary tumors from patients’ post-therapy that had no response. Gene set enrichment analysis (GSEA) in treated patients with minimal OPN+TAM showed decreased exhaustion markers (logFC>0.5, FDR adjusted p-value<0.05). We then found that the strongest OPN-CD44 interactions were predicted between OPN+TAM and other TAM as well as T cells. Separating patients into high and low predicted OPN-CD44 interactors, GSEA demonstrated enrichment of cell death regulation and apoptosis pathways in CD8+ T cells in high interactors. To determine potential targets for modulation of OPN+TAM, we undertook trajectory analysis demonstrating that OPN+TAM were an intermediate phase, suggesting they are malleable. We then found that these OPN+TAM highly express TLR8 in both scRNAseq and FC, suggesting a potential target to modulate OPN+TAM. In this study, we demonstrate that OPN+TAM are the largest proportion of TAM in CRC tumors. Additionally, the presence of OPN+ TAM in treated samples is associated with response to systemic therapy. OPN+TAM appear to represent an intermediate phase of TAMs differentiation and highly express TLR8 suggesting a mechanism for modulation to improve patient response. Further ex vivo studies targeting these cells and their role in therapeutic response are critical.

Citation Format: Yulia I Nussbaum, Yariswamy Manjunath, Elizabeth Shumway, Kaifi T Jussuf, Wesley C Warren, Jonathan B Mitchem. Studying the role of OPN-CD44 interaction in colorectal cancer antitumor immunity suppression using scRNA-seq data [abstract]. In: Proceedings of the AACR Special Conference: Tumor Immunology and Immunotherapy; 2022 Oct 21-24; Boston, MA. Philadelphia (PA): AACR; Cancer Immunol Res 2022;10(12 Suppl):Abstract nr B57.

Somatic mutation variant analysis in rural, resectable non-small cell lung carcinoma patients

Rural non-small cell lung cancer (NSCLC) patients do worse, largely related to lack of access to care. In this study, the mutational characteristics and potential for targeted therapy in rural, resectable NSCLC patients using whole exome sequencing (WES) were analyzed. WES was performed on tumor-adjacent normal pairs from rural patients undergoing resection for NSCLC. Sequencing alignment, variant-calling, annotation, and tumor mutational burden (TMB) calculations were performed using standard methods. cBioportal and OncoKB were used for comparisons of mutational frequencies and actionable targets. Thirty-four NSCLC patients underwent WES after surgical resection. The gene most frequently containing somatic variants was TP53. The median number of somatic variants was 188 (Range 11-1056), and median TMB was 3.30 (0.33-18.56) nonsynonymous mutations per Mb. Tumor stage and survival were not associated with number of variants, TMB or TP53 mutational status. Significant concordance among the most common mutations when cross-referenced to cBioportal (R = 0.78, p < 0.0001) was observed. 24% of patients had variants in actionable genes based on OncoKB annotation. In summary, we demonstrate baseline mutational frequency and establish foundations for targeted adjuvant trials in rural NSCLC patients with specific differences. Future studies must ensure to include rural patients to improve NSCLC patient outcomes.

Transcript- and annotation-guided genome assembly of the European starling

The European starling, Sturnus vulgaris, is an ecologically significant, globally invasive avian species that is also suffering from a major decline in its native range. Here, we present the genome assembly and long-read transcriptome of an Australian-sourced European starling (S. vulgaris vAU), and a second, North American, short-read genome assembly (S. vulgaris vNA), as complementary reference genomes for population genetic and evolutionary characterization. S. vulgaris vAU combined 10× genomics linked-reads, low-coverage Nanopore sequencing, and PacBio Iso-Seq full-length transcript scaffolding to generate a 1050 Mb assembly on 6222 scaffolds (7.6 Mb scaffold N50, 94.6% busco completeness). Further scaffolding against the high-quality zebra finch (Taeniopygia guttata) genome assigned 98.6% of the assembly to 32 putative nuclear chromosome scaffolds. Species-specific transcript mapping and gene annotation revealed good gene-level assembly and high functional completeness. Using S. vulgaris vAU, we demonstrate how the multifunctional use of PacBio Iso-Seq transcript data and complementary homology-based annotation of sequential assembly steps (assessed using a new tool, saaga) can be used to assess, inform, and validate assembly workflow decisions. We also highlight some counterintuitive behaviour in traditional busco metrics, and present buscomp, a complementary tool for assembly comparison designed to be robust to differences in assembly size and base-calling quality. This work expands our knowledge of avian genomes and the available toolkit for assessing and improving genome quality. The new genomic resources presented will facilitate further global genomic and transcriptomic analysis on this ecologically important species.

TBX21 Methylation as a Potential Regulator of Immune Suppression in CMS1 Subtype Colorectal Cancer

Cytotoxic T lymphocyte (CTL) infiltration is associated with survival, recurrence, and therapeutic response in colorectal cancer (CRC). Immune checkpoint inhibitor (ICI) therapy, which requires CTLs for response, does not work for most CRC patients. Therefore, it is critical to improve our understanding of immune resistance in this disease. We utilized 2391 CRC patients and 7 omics datasets, integrating clinical and genomic data to determine how DNA methylation may impact survival and CTL function in CRC. Using comprehensive molecular subtype (CMS) 1 patients as reference, we found TBX21 to be the only gene with altered expression and methylation that was associated with CTL infiltration. We found that CMS1 patients with high TBX21 expression and low methylation had a significant survival advantage. To confirm the role of Tbx21 in CTL function, we utilized scRNAseq data, demonstrating the association of TBX21 with markers of enhanced CTL function. Further analysis using pathway enrichment found that the genes TBX21, MX1, and SP140 had altered expression and methylation, suggesting that the TP53/P53 pathway may modify TBX21 methylation to upregulate TBX21 expression. Together, this suggests that targeting epigenetic modification more specifically for therapy and patient stratification may provide improved outcomes in CRC.

CTSK variant implicated in suspected pyknodysostosis in a domestic cat

Case summary

A 9-month-old entire male domestic longhair cat presented with a history of pathological fractures, chronic musculoskeletal pain and poor growth. Multiple facial and skeletal abnormalities were identified on physical examination and advanced imaging (CT and radiographs). A variant in CTSK was identified in the affected cat following whole-exome sequencing (WES). The cat was managed symptomatically with diet, environmental modifications and analgesia.

Relevance and novel information

This is the first report of a cat with a similar clinical presentation and genetic variant to the hereditary human genetic disorder pyknodysostosis. In this case, WES was performed, which often facilitates the diagnosis of various hereditary disorders (ie, a conceptual framework for practicing feline genomic medicine). Despite the severe skeletal and appendicular abnormalities described, the cat was alive more than 2 years after its initial presentation.

Abstract 267: Xenograft models of non-metastatic non-small cell lung cancer reveals micrometastasis-associated single cell composition

Background: Circulating tumor cells (CTCs) represent micrometastatic disease and may offer unique insights into future recurrences in lethal malignancies, including non-small cell lung cancer (NSCLC). Due to CTC rarity and limited stability, no CTC-derived xenograft (CDX) models have ever been generated from non-metastatic NSCLC patients directly. Alternative strategies are needed to molecularly characterize CTCs in this potentially curable patient group.

Methods: Surgically resected NSCLC primary tumor tissues were implanted in immunodeficient mice to establish ten patient-derived xenografts (PDXs). CTCs from 2/10 PDX models led to generation of two stable metastatic models that were studied by single cell sequencing.

Results: Single cell analysis revealed an additional alveolar epithelial type II (AT2) population in metastatic tumors, besides a common AT2 cluster in PDX/metastatic tumors. This was consistent with an external validation set analysis in primary and metastatic NSCLC patient tumors. Further, AT2 clusters of metastatic tumors expressed higher cancer stemness genes versus primary PDX tumor that was recapitulated in patients primary and metastatic tumors.

Conclusions: Stable metastatic models from early stage NSCLC patients can be generated with CTCs from PDX models. The distinct AT2 population identified in CDX tumors with cancer stemness features might be critical mediator of metastasis that deserves further study to discover personalized strategies against NSCLC micrometastases.

Citation Format: Kanve Nagaraj Suvilesh, Yulia I. Nussbaum, Vijay Radhakrishnan, Yariswamy Manjunath, Diego M. Avella, Kevin F. Staveley-O’Carroll, Eric T. Kimchi, Aadel A. Chaudhuri, Chi-Ren Shyu, Guangfu Li, Klaus Pantel, Wesley C. Warren, Jonathan B. Mitchem, Jussuf T. Kaifi. Xenograft models of non-metastatic non-small cell lung cancer reveals micrometastasis-associated single cell composition [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 267.

Tumorigenic circulating tumor cells from xenograft mouse models of non-metastatic NSCLC patients reveal distinct single cell heterogeneity and drug responses

Background

Circulating tumor cells (CTCs) are liquid biopsies that represent micrometastatic disease and may offer unique insights into future recurrences in non-small cell lung cancer (NSCLC). Due to CTC rarity and limited stability, no stable CTC-derived xenograft (CDX) models have ever been generated from non-metastatic NSCLC patients directly. Alternative strategies are needed to molecularly characterize CTCs and means of potential future metastases in this potentially curable patient group.

Methods

Surgically resected NSCLC primary tumor tissues from non-metastatic patients were implanted subcutaneously in immunodeficient mice to establish primary tumor patient-derived xenograft (ptPDX) models. CTCs were isolated as liquid biopsies from the blood of ptPDX mice and re-implanted subcutaneously into naïve immunodeficient mice to generate liquid biopsy CTC-derived xenograft (CDX) tumor models. Single cell RNA sequencing was performed and validated in an external dataset of non-xenografted human NSCLC primary tumor and metastases tissues. Drug response testing in CDX models was performed with standard of care chemotherapy (carboplatin/paclitaxel). Blockade of MYC, which has a known role in drug resistance, was performed with a MYC/MAX dimerization inhibitor (10058-F4).

Results

Out of ten ptPDX, two (20%) stable liquid biopsy CDX mouse models were generated. Single cell RNA sequencing analysis revealed an additional regenerative alveolar epithelial type II (AT2)-like cell population in CDX tumors that was also identified in non-xenografted NSCLC patients’ metastases tissues. Drug testing using these CDX models revealed different treatment responses to carboplatin/paclitaxel. MYC target genes and c-MYC protein were upregulated in the chemoresistant CDX model, while MYC/MAX dimerization blocking could overcome chemoresistance to carboplatin/paclitaxel.

Conclusions

To overcome the lack of liquid biopsy CDX models from non-metastatic NSCLC patients, CDX models can be generated with CTCs from ptPDX models that were originally established from patients’ primary tumors. Single cell analyses can identify distinct drug responses and cell heterogeneities in CDX tumors that can be validated in NSCLC metastases tissues. CDX models deserve further development and study to discover personalized strategies against micrometastases in non-metastatic NSCLC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-022-01553-5.

Single Circulating-Tumor-Cell-Targeted Sequencing to Identify Somatic Variants in Liquid Biopsies in Non-Small-Cell Lung Cancer Patients

Non-small-cell lung cancer (NSCLC) accounts for most cancer-related deaths worldwide. Liquid biopsy by a blood draw to detect circulating tumor cells (CTCs) is a tool for molecular profiling of cancer using single-cell and next-generation sequencing (NGS) technologies. The aim of the study was to identify somatic variants in single CTCs isolated from NSCLC patients by targeted NGS. Thirty-one subjects (20 NSCLC patients, 11 smokers without cancer) were enrolled for blood draws (7.5 mL). CTCs were identified by immunofluorescence, individually retrieved, and DNA-extracted. Targeted NGS was performed to detect somatic variants (single-nucleotide variants (SNVs) and insertions/deletions (Indels)) across 65 oncogenes and tumor suppressor genes. Cancer-associated variants were classified using OncoKB database. NSCLC patients had significantly higher CTC counts than control smokers (p = 0.0132; Mann–Whitney test). Analyzing 23 CTCs and 13 white blood cells across seven patients revealed a total of 644 somatic variants that occurred in all CTCs within the same subject, ranging from 1 to 137 per patient. The highest number of variants detected in ≥1 CTC within a patient was 441. A total of 18/65 (27.7%) genes were highly mutated. Mutations with oncogenic impact were identified in functional domains of seven oncogenes/tumor suppressor genes (NF1, PTCH1, TP53, SMARCB1, SMAD4, KRAS, and ERBB2). Single CTC-targeted NGS detects heterogeneous and shared mutational signatures within and between NSCLC patients. CTC single-cell genomics have potential for integration in NSCLC precision oncology.

3D genomics across the tree of life reveals condensin II as a determinant of architecture type

We investigated genome folding across the eukaryotic tree of life. We find two types of three-dimensional (3D) genome architectures at the chromosome scale. Each type appears and disappears repeatedly during eukaryotic evolution. The type of genome architecture that an organism exhibits correlates with the absence of condensin II subunits. Moreover, condensin II depletion converts the architecture of the human genome to a state resembling that seen in organisms such as fungi or mosquitoes. In this state, centromeres cluster together at nucleoli, and heterochromatin domains merge. We propose a physical model in which lengthwise compaction of chromosomes by condensin II during mitosis determines chromosome-scale genome architecture, with effects that are retained during the subsequent interphase. This mechanism likely has been conserved since the last common ancestor of all eukaryotes.

Current and Prospective Methods for Assessing Anti-Tumor Immunity in Colorectal Cancer

Colorectal cancer (CRC) remains one of the deadliest malignancies worldwide despite recent progress in treatment strategies. Though immune checkpoint inhibition has proven effective for a number of other tumors, it offers benefits in only a small group of CRC patients with high microsatellite instability. In general, heterogenous cell groups in the tumor microenvironment are considered as the major barrier for unveiling the causes of low immune response. Therefore, deconvolution of cellular components in highly heterogeneous microenvironments is crucial for understanding the immune contexture of cancer. In this review, we assimilate current knowledge and recent studies examining anti-tumor immunity in CRC. We also discuss the utilization of novel immune contexture assessment methods that have not been used in CRC research to date.

Explainable artificial intelligence in high-throughput drug repositioning for subgroup stratifications with interventionable potential

Enabling precision medicine requires developing robust patient stratification methods as well as drugs tailored to homogeneous subgroups of patients from a heterogeneous population. Developing de novo drugs is expensive and time consuming with an ultimately low FDA approval rate. These limitations make developing new drugs for a small portion of a disease population unfeasible. Therefore, drug repositioning is an essential alternative for developing new drugs for a disease subpopulation. This shows the importance of developing data-driven approaches that find druggable homogeneous subgroups within the disease population and reposition the drugs for these subgroups. In this study, we developed an explainable AI approach for patient stratification and drug repositioning. Contrast pattern mining and network analysis were used to discover homogeneous subgroups within a disease population. For each subgroup, a biomedical network analysis was done to find the drugs that are most relevant to a given subgroup of patients. The set of candidate drugs for each subgroup was ranked using an aggregated drug score assigned to each drug. The proposed method represents a human-in-the-loop framework, where medical experts use the data-driven results to generate hypotheses and obtain insights into potential therapeutic candidates for patients who belong to a subgroup. Colorectal cancer (CRC) was used as a case study. Patients' phenotypic and genotypic data was utilized with a heterogeneous knowledge base because it gives a multi-view perspective for finding new indications for drugs outside of their original use. Our analysis of the top candidate drugs for the subgroups identified by medical experts showed that most of these drugs are cancer-related, and most of them have the potential to be a CRC regimen based on studies in the literature.

Towards complete and error-free genome assemblies of all vertebrate species

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

Ultracontinuous Single Haplotype Genome Assemblies for the Domestic Cat (Felis catus) and Asian Leopard Cat (Prionailurus bengalensis)

In addition to including one of the most popular companion animals, species from the cat family Felidae serve as a powerful system for genetic analysis of inherited and infectious disease, as well as for the study of phenotypic evolution and speciation. Previous diploid-based genome assemblies for the domestic cat have served as the primary reference for genomic studies within the cat family. However, these versions suffered from poor resolution of complex and highly repetitive regions, with substantial amounts of unplaced sequence that is polymorphic or copy number variable. We sequenced the genome of a female F1 Bengal hybrid cat, the offspring of a domestic cat (Felis catus) x Asian leopard cat (Prionailurus bengalensis) cross, with PacBio long sequence reads and used Illumina sequence reads from the parents to phase >99.9% of the reads into the 2 species' haplotypes. De novo assembly of the phased reads produced highly continuous haploid genome assemblies for the domestic cat and Asian leopard cat, with contig N50 statistics exceeding 83 Mb for both genomes. Whole-genome alignments reveal the Felis and Prionailurus genomes are colinear, and the cytogenetic differences between the homologous F1 and E4 chromosomes represent a case of centromere repositioning in the absence of a chromosomal inversion. Both assemblies offer significant improvements over the previous domestic cat reference genome, with a 100% increase in contiguity and the capture of the vast majority of chromosome arms in 1 or 2 large contigs. We further demonstrated that comparably accurate F1 haplotype phasing can be achieved with members of the same species when one or both parents of the trio are not available. These novel genome resources will empower studies of feline precision medicine, adaptation, and speciation.

Chromosome-length genome assembly and structural variations of the primal Basenji dog (Canis lupus familiaris) genome

BACKGROUND

Basenjis are considered an ancient dog breed of central African origins that still live and hunt with tribesmen in the African Congo. Nicknamed the barkless dog, Basenjis possess unique phylogeny, geographical origins and traits, making their genome structure of great interest. The increasing number of available canid reference genomes allows us to examine the impact the choice of reference genome makes with regard to reference genome quality and breed relatedness.

RESULTS

Here, we report two high quality de novo Basenji genome assemblies: a female, China (CanFam_Bas), and a male, Wags. We conduct pairwise comparisons and report structural variations between assembled genomes of three dog breeds: Basenji (CanFam_Bas), Boxer (CanFam3.1) and German Shepherd Dog (GSD) (CanFam_GSD). CanFam_Bas is superior to CanFam3.1 in terms of genome contiguity and comparable overall to the high quality CanFam_GSD assembly. By aligning short read data from 58 representative dog breeds to three reference genomes, we demonstrate how the choice of reference genome significantly impacts both read mapping and variant detection.

CONCLUSIONS

The growing number of high-quality canid reference genomes means the choice of reference genome is an increasingly critical decision in subsequent canid variant analyses. The basal position of the Basenji makes it suitable for variant analysis for targeted applications of specific dog breeds. However, we believe more comprehensive analyses across the entire family of canids is more suited to a pangenome approach. Collectively this work highlights the importance the choice of reference genome makes in all variation studies.

The genome of the stable fly, Stomoxys calcitrans, reveals potential mechanisms underlying reproduction, host interactions, and novel targets for pest control

Background

The stable fly, Stomoxys calcitrans, is a major blood-feeding pest of livestock that has near worldwide distribution, causing an annual cost of over $2 billion for control and product loss in the USA alone. Control of these flies has been limited to increased sanitary management practices and insecticide application for suppressing larval stages. Few genetic and molecular resources are available to help in developing novel methods for controlling stable flies.

Results

This study examines stable fly biology by utilizing a combination of high-quality genome sequencing and RNA-Seq analyses targeting multiple developmental stages and tissues. In conjunction, 1600 genes were manually curated to characterize genetic features related to stable fly reproduction, vector host interactions, host-microbe dynamics, and putative targets for control. Most notable was characterization of genes associated with reproduction and identification of expanded gene families with functional associations to vision, chemosensation, immunity, and metabolic detoxification pathways.

Conclusions

The combined sequencing, assembly, and curation of the male stable fly genome followed by RNA-Seq and downstream analyses provide insights necessary to understand the biology of this important pest. These resources and new data will provide the groundwork for expanding the tools available to control stable fly infestations. The close relationship of Stomoxys to other blood-feeding (horn flies and Glossina) and non-blood-feeding flies (house flies, medflies, Drosophila) will facilitate understanding of the evolutionary processes associated with development of blood feeding among the Cyclorrhapha.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00975-9.

The Developmental and Genetic Architecture of the Sexually Selected Male Ornament of Swordtails

Sexual selection results in sex-specific characters like the conspicuously pigmented extension of the ventral tip of the caudal fin-the "sword"-in males of several species of Xiphophorus fishes. To uncover the genetic architecture underlying sword formation and to identify genes that are associated with its development, we characterized the sword transcriptional profile and combined it with genetic mapping approaches. Results showed that the male ornament of swordtails develops from a sexually non-dimorphic prepattern of transcription factors in the caudal fin. Among genes that constitute the exclusive sword transcriptome and are located in the genomic region associated with this trait we identify the potassium channel, Kcnh8, as a sword development gene. In addition to its neural function kcnh8 performs a known role in fin growth. These findings indicate that during evolution of swordtails a brain gene has been co-opted for an additional novel function in establishing a male ornament.

A chromosome-level genome of Astyanax mexicanus surface fish for comparing population-specific genetic differences contributing to trait evolution

Identifying the genetic factors that underlie complex traits is central to understanding the mechanistic underpinnings of evolution. Cave-dwelling Astyanax mexicanus populations are well adapted to subterranean life and many populations appear to have evolved troglomorphic traits independently, while the surface-dwelling populations can be used as a proxy for the ancestral form. Here we present a high-resolution, chromosome-level surface fish genome, enabling the first genome-wide comparison between surface fish and cavefish populations. Using this resource, we performed quantitative trait locus (QTL) mapping analyses and found new candidate genes for eye loss such as dusp26. We used CRISPR gene editing in A. mexicanus to confirm the essential role of a gene within an eye size QTL, rx3, in eye formation. We also generated the first genome-wide evaluation of deletion variability across cavefish populations to gain insight into this potential source of cave adaptation. The surface fish genome reference now provides a more complete resource for comparative, functional and genetic studies of drastic trait differences within a species.

Fixation of allelic gene expression landscapes and expression bias pattern shape the transcriptome of the clonal Amazon molly

The Amazon molly is a unique clonal fish species that originated from an interspecies hybrid between Poecilia species P. mexicana and P. latipinna. It reproduces by gynogenesis, which eliminates paternal genomic contribution to offspring. An earlier study showed that Amazon molly shows biallelic expression for a large portion of the genome, leading to two main questions: (1) Are the allelic expression patterns from the initial hybridization event stabilized or changed during establishment of the asexual species and its further evolution? (2) Is allelic expression biased toward one parental allele a stochastic or adaptive process? To answer these questions, the allelic expression of P. formosa siblings was assessed to investigate intra- and inter-cohort allelic expression variability. For comparison, interspecies hybrids between P. mexicana and P. latipinna were produced in the laboratory to represent the P. formosa ancestor. We have identified inter-cohort and intra-cohort variation in parental allelic expression. The existence of inter-cohort divergence suggests functional P. formosa allelic expression patterns do not simply reflect the atavistic situation of the first interspecies hybrid but potentially result from long-term selection of transcriptional fitness. In addition, clonal fish show a transcriptional trend representing minimal intra-clonal variability in allelic expression patterns compared to the corresponding hybrids. The intra-clonal similarity in gene expression translates to sophisticated genetic functional regulation at the individuum level. These findings suggest the parental alleles inherited by P. formosa form tightly regulated genetic networks that lead to a stable transcriptomic landscape within clonal individuals.

Sequence diversity analyses of an improved rhesus macaque genome enhance its biomedical utility

The rhesus macaque (Macaca mulatta) is the most widely studied nonhuman primate (NHP) in biomedical research. We present an updated reference genome assembly (Mmul_10, contig N50 = 46 Mbp) that increases the sequence contiguity 120-fold and annotate it using 6.5 million full-length transcripts, thus improving our understanding of gene content, isoform diversity, and repeat organization. With the improved assembly of segmental duplications, we discovered new lineage-specific genes and expanded gene families that are potentially informative in studies of evolution and disease susceptibility. Whole-genome sequencing (WGS) data from 853 rhesus macaques identified 85.7 million single-nucleotide variants (SNVs) and 10.5 million indel variants, including potentially damaging variants in genes associated with human autism and developmental delay, providing a framework for developing noninvasive NHP models of human disease.

Sequence analysis in Bos taurus reveals pervasiveness of X-Y arms races in mammalian lineages

Studies of Y Chromosome evolution have focused primarily on gene decay, a consequence of suppression of crossing-over with the X Chromosome. Here, we provide evidence that suppression of X-Y crossing-over unleashed a second dynamic: selfish X-Y arms races that reshaped the sex chromosomes in mammals as different as cattle, mice, and men. Using super-resolution sequencing, we explore the Y Chromosome of Bos taurus (bull) and find it to be dominated by massive, lineage-specific amplification of testis-expressed gene families, making it the most gene-dense Y Chromosome sequenced to date. As in mice, an X-linked homolog of a bull Y-amplified gene has become testis-specific and amplified. This evolutionary convergence implies that lineage-specific X-Y coevolution through gene amplification, and the selfish forces underlying this phenomenon, were dominatingly powerful among diverse mammalian lineages. Together with Y gene decay, X-Y arms races molded mammalian sex chromosomes and influenced the course of mammalian evolution.

Global assessment of organ specific basal gene expression over a diurnal cycle with analyses of gene copies exhibiting cyclic expression patterns

Background

Studying functional divergences between paralogs that originated from genome duplication is a significant topic in investigating molecular evolution. Genes that exhibit basal level cyclic expression patterns including circadian and light responsive genes are important physiological regulators. Temporal shifts in basal gene expression patterns are important factors to be considered when studying genetic functions. However, adequate efforts have not been applied to studying basal gene expression variation on a global scale to establish transcriptional activity baselines for each organ. Furthermore, the investigation of cyclic expression pattern comparisons between genome duplication created paralogs, and potential functional divergence between them has been neglected. To address these questions, we utilized a teleost fish species, Xiphophorus maculatus, and profiled gene expression within 9 organs at 3-h intervals throughout a 24-h diurnal period.

Results

Our results showed 1.3–21.9% of genes in different organs exhibited cyclic expression patterns, with eye showing the highest fraction of cycling genes while gonads yielded the lowest. A majority of the duplicated gene pairs exhibited divergences in their basal level expression patterns wherein only one paralog exhibited an oscillating expression pattern, or both paralogs exhibit oscillating expression patterns, but each gene duplicate showed a different peak expression time, and/or in different organs.

Conclusions

These observations suggest cyclic genes experienced significant sub-, neo-, or non-functionalization following the teleost genome duplication event. In addition, we developed a customized, web-accessible, gene expression browser to facilitate data mining and data visualization for the scientific community.

A new domestic cat genome assembly based on long sequence reads empowers feline genomic medicine and identifies a novel gene for dwarfism

The domestic cat (Felis catus) numbers over 94 million in the USA alone, occupies households as a companion animal, and, like humans, suffers from cancer and common and rare diseases. However, genome-wide sequence variant information is limited for this species. To empower trait analyses, a new cat genome reference assembly was developed from PacBio long sequence reads that significantly improve sequence representation and assembly contiguity. The whole genome sequences of 54 domestic cats were aligned to the reference to identify single nucleotide variants (SNVs) and structural variants (SVs). Across all cats, 16 SNVs predicted to have deleterious impacts and in a singleton state were identified as high priority candidates for causative mutations. One candidate was a stop gain in the tumor suppressor FBXW7. The SNV is found in cats segregating for feline mediastinal lymphoma and is a candidate for inherited cancer susceptibility. SV analysis revealed a complex deletion coupled with a nearby potential duplication event that was shared privately across three unrelated cats with dwarfism and is found within a known dwarfism associated region on cat chromosome B1. This SV interrupted UDP-glucose 6-dehydrogenase (UGDH), a gene involved in the biosynthesis of glycosaminoglycans. Importantly, UGDH has not yet been associated with human dwarfism and should be screened in undiagnosed patients. The new high-quality cat genome reference and the compilation of sequence variation demonstrate the importance of these resources when searching for disease causative alleles in the domestic cat and for identification of feline biomedical models.

The practice of genomic medicine is predicated on the availability of a high quality reference genome and an understanding of the impact of genome variation. Such resources have lead to countless discoveries in humans, however by working exclusively within the framework of human genetics, our potential for understanding diseases biology is limited, as similar analyses in other species have often lead to novel insights. The generation of Felis_catus_9.0, a new high quality reference genome for the domestic cat, helps facilitate the expansion of genomic medicine into the Felis lineage. Using Felis_catus_9.0 we analyze the landscape of genomic variation from a collection of 54 cats within the context of human gene constraint. The distribution of variant impacts in cats is correlated with patterns of gene constraint in humans, indicating the utility of this reference for identifying novel mutations that cause phenotypes relevant to human and cat health. Moreover, structural variant analysis revealed a novel variant for feline dwarfism in UGDH, a gene that has not been associated with dwarfism in any other species, suggesting a role for UGDH in cases of undiagnosed dwarfism in humans.

Neuronal Ceroid Lipofuscinosis in a Domestic Cat Associated with a DNA Sequence Variant That Creates a Premature Stop Codon in CLN6

A neutered male domestic medium-haired cat presented at a veterinary neurology clinic at 20 months of age due to progressive neurological signs that included visual impairment, focal myoclonus, and frequent severe generalized seizures that were refractory to treatment with phenobarbital. Magnetic resonance imaging revealed diffuse global brain atrophy. Due to the severity and frequency of its seizures, the cat was euthanized at 22 months of age. Microscopic examination of the cerebellum, cerebral cortex and brainstem revealed pronounced intracellular accumulations of autofluorescent storage material and inflammation in all 3 brain regions. Ultrastructural examination of the storage material indicated that it consisted almost completely of tightly-packed membrane-like material. The clinical signs and neuropathology strongly suggested that the cat suffered from a form of neuronal ceroid lipofuscinosis (NCL). Whole exome sequence analysis was performed on genomic DNA from the affected cat. Comparison of the sequence data to whole exome sequence data from 39 unaffected cats and whole genome sequence data from an additional 195 unaffected cats revealed a homozygous variant in CLN6 that was unique to the affected cat. This variant was predicted to cause a stop gain in the transcript due to a guanine to adenine transition (ENSFCAT00000025909:c.668G > A; XM_003987007.5:c.668G > A) and was the sole loss of function variant detected. CLN6 variants in other species, including humans, dogs, and sheep, are associated with the CLN6 form of NCL. Based on the affected cat's clinical signs, neuropathology and molecular genetic analysis, we conclude that the cat's disorder resulted from the loss of function of CLN6. This study is only the second to identify the molecular genetic basis of a feline NCL. Other cats exhibiting similar signs can now be screened for the CLN6 variant. This could lead to establishment of a feline model of CLN6 disease that could be used in therapeutic intervention studies.

A novel swine sex-linked marker and its application across different mammalian species

Advances in genome editing tools have reduced barriers to the creation of animal models. Due to their anatomical and physiological similarities to humans, there has been a growing need for pig models to study human diseases, for xenotransplantation and translational research. The ability to determine the sex of genetically modified embryos, cells or fetuses is beneficial for every project involving the production of transgenic animals. This strategy can improve the time-efficiency and lower the production costs. Additionally, sex assessment is very useful for wildlife studies to understand population behavior and structure. Thus, we developed a simple and fast PCR-based protocol for sex determination in pigs by using a unique primer set to amplify either the DDX3X or DDX3Y gene. The sex was 100% correctly assigned when tail genomic DNA, Day-35 fetus and hair samples from pigs were used. For both blastocysts and oocytes (84.6% and 96.5% of efficacy, respectively) the unidentified samples were potentially due to a limitation in sample size. Our assay also worked for domestic sheep (Ovis aries), American bison (Bison bison) and European cattle (Bos taurus) samples and by in silico analysis we confirmed X-Y amplicon length polymorphisms for the DDX3 gene in 12 other mammalian species. This PCR protocol for determining sex in pig tissues and cells showed to be simple, specific, highly reproducible and less time consuming as well as an important tool for other livestock species and wildlife studies.

Mutations in the Kinesin-2 Motor KIF3B Cause an Autosomal-Dominant Ciliopathy

Kinesin-2 enables ciliary assembly and maintenance as an anterograde intraflagellar transport (IFT) motor. Molecular motor activity is driven by a heterotrimeric complex comprised of KIF3A and KIF3B or KIF3C plus one non-motor subunit, KIFAP3. Using exome sequencing, we identified heterozygous KIF3B variants in two unrelated families with hallmark ciliopathy phenotypes. In the first family, the proband presents with hepatic fibrosis, retinitis pigmentosa, and postaxial polydactyly; he harbors a de novo c.748G>C (p.Glu250Gln) variant affecting the kinesin motor domain encoded by KIF3B. The second family is a six-generation pedigree affected predominantly by retinitis pigmentosa. Affected individuals carry a heterozygous c.1568T>C (p.Leu523Pro) KIF3B variant segregating in an autosomal-dominant pattern. We observed a significant increase in primary cilia length in vitro in the context of either of the two mutations while variant KIF3B proteins retained stability indistinguishable from wild type. Furthermore, we tested the effects of KIF3B mutant mRNA expression in the developing zebrafish retina. In the presence of either missense variant, rhodopsin was sequestered to the photoreceptor rod inner segment layer with a concomitant increase in photoreceptor cilia length. Notably, impaired rhodopsin trafficking is also characteristic of recessive KIF3B models as exemplified by an early-onset, autosomal-recessive, progressive retinal degeneration in Bengal cats; we identified a c.1000G>A (p.Ala334Thr) KIF3B variant by genome-wide association study and whole-genome sequencing. Together, our genetic, cell-based, and in vivo modeling data delineate an autosomal-dominant syndromic retinal ciliopathy in humans and suggest that multiple KIF3B pathomechanisms can impair kinesin-driven ciliary transport in the photoreceptor.

The Piranha Genome Provides Molecular Insight Associated to Its Unique Feeding Behavior

The piranha enjoys notoriety due to its infamous predatory behavior but much is still not understood about its evolutionary origins and the underlying molecular mechanisms for its unusual feeding biology. We sequenced and assembled the red-bellied piranha (Pygocentrus nattereri) genome to aid future phenotypic and genetic investigations. The assembled draft genome is similar to other related fishes in repeat composition and gene count. Our evaluation of genes under positive selection suggests candidates for adaptations of piranhas’ feeding behavior in neural functions, behavior, and regulation of energy metabolism. In the fasted brain, we find genes differentially expressed that are involved in lipid metabolism and appetite regulation as well as genes that may control the aggression/boldness behavior of hungry piranhas. Our first analysis of the piranha genome offers new insight and resources for the study of piranha biology and for feeding motivation and starvation in other organisms.

The quail genome: insights into social behaviour, seasonal biology and infectious disease response

BACKGROUND

The Japanese quail (Coturnix japonica) is a popular domestic poultry species and an increasingly significant model species in avian developmental, behavioural and disease research.

RESULTS

We have produced a high-quality quail genome sequence, spanning 0.93 Gb assigned to 33 chromosomes. In terms of contiguity, assembly statistics, gene content and chromosomal organisation, the quail genome shows high similarity to the chicken genome. We demonstrate the utility of this genome through three diverse applications. First, we identify selection signatures and candidate genes associated with social behaviour in the quail genome, an important agricultural and domestication trait. Second, we investigate the effects and interaction of photoperiod and temperature on the transcriptome of the quail medial basal hypothalamus, revealing key mechanisms of photoperiodism. Finally, we investigate the response of quail to H5N1 influenza infection. In quail lung, many critical immune genes and pathways were downregulated after H5N1 infection, and this may be key to the susceptibility of quail to H5N1.

CONCLUSIONS

We have produced a high-quality genome of the quail which will facilitate further studies into diverse research questions using the quail as a model avian species.

Human-specific tandem repeat expansion and differential gene expression during primate evolution

Short tandem repeats (STRs) and variable number tandem repeats (VNTRs) are important sources of natural and disease-causing variation, yet they have been problematic to resolve in reference genomes and genotype with short-read technology. We created a framework to model the evolution and instability of STRs and VNTRs in apes. We phased and assembled 3 ape genomes (chimpanzee, gorilla, and orangutan) using long-read and 10x Genomics linked-read sequence data for 21,442 human tandem repeats discovered in 6 haplotype-resolved assemblies of Yoruban, Chinese, and Puerto Rican origin. We define a set of 1,584 STRs/VNTRs expanded specifically in humans, including large tandem repeats affecting coding and noncoding portions of genes (e.g., MUC3A, CACNA1C). We show that short interspersed nuclear element-VNTR-Alu (SVA) retrotransposition is the main mechanism for distributing GC-rich human-specific tandem repeat expansions throughout the genome but with a bias against genes. In contrast, we observe that VNTRs not originating from retrotransposons have a propensity to cluster near genes, especially in the subtelomere. Using tissue-specific expression from human and chimpanzee brains, we identify genes where transcript isoform usage differs significantly, likely caused by cryptic splicing variation within VNTRs. Using single-cell expression from cerebral organoids, we observe a strong effect for genes associated with transcription profiles analogous to intermediate progenitor cells. Finally, we compare the sequence composition of some of the largest human-specific repeat expansions and identify 52 STRs/VNTRs with at least 40 uninterrupted pure tracts as candidates for genetically unstable regions associated with disease.

Comparative genomic analysis of six Glossina genomes, vectors of African trypanosomes

BACKGROUND

Tsetse flies (Glossina sp.) are the vectors of human and animal trypanosomiasis throughout sub-Saharan Africa. Tsetse flies are distinguished from other Diptera by unique adaptations, including lactation and the birthing of live young (obligate viviparity), a vertebrate blood-specific diet by both sexes, and obligate bacterial symbiosis. This work describes the comparative analysis of six Glossina genomes representing three sub-genera: Morsitans (G. morsitans morsitans, G. pallidipes, G. austeni), Palpalis (G. palpalis, G. fuscipes), and Fusca (G. brevipalpis) which represent different habitats, host preferences, and vectorial capacity.

RESULTS

Genomic analyses validate established evolutionary relationships and sub-genera. Syntenic analysis of Glossina relative to Drosophila melanogaster shows reduced structural conservation across the sex-linked X chromosome. Sex-linked scaffolds show increased rates of female-specific gene expression and lower evolutionary rates relative to autosome associated genes. Tsetse-specific genes are enriched in protease, odorant-binding, and helicase activities. Lactation-associated genes are conserved across all Glossina species while male seminal proteins are rapidly evolving. Olfactory and gustatory genes are reduced across the genus relative to other insects. Vision-associated Rhodopsin genes show conservation of motion detection/tracking functions and variance in the Rhodopsin detecting colors in the blue wavelength ranges.

CONCLUSIONS

Expanded genomic discoveries reveal the genetics underlying Glossina biology and provide a rich body of knowledge for basic science and disease control. They also provide insight into the evolutionary biology underlying novel adaptations and are relevant to applied aspects of vector control such as trap design and discovery of novel pest and disease control strategies.