Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly.

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.

Alternative Splicing of Opioid Receptor Genes Shows a Conserved Pattern for 6TM Receptor Variants.

The opioid receptor (OPR) family comprises the mu-, delta-, and kappa-opioid, and nociceptin receptors that belong to the superfamily of 7-transmembrane spanning G protein-coupled receptors (GPCRs). The mu-opioid receptor is the main target for clinically used opioid analgesics, and its biology has been extensively studied. The N-terminally truncated 6TM receptors isoform produced through alternative splicing of the OPRM1 gene displays unique signaling and analgesic properties, but it is unclear if other OPRs have the same ability. In this study, we have built a comprehensive map of alternative splicing events that produce 6TM receptor variants in all the OPRs and demonstrated their evolutionary conservation. We then obtained evidence for their translation through ribosomal footprint analysis. We discovered that N-terminally truncated 6TM GPCRs are rare in the human genome and OPRs are overrepresented in this group. Finally, we also observed a significant enrichment of 6TM GPCR genes among genes associated with pain, psychiatric disorders, and addiction. Understanding the biology of 6TM receptors and leveraging this knowledge for drug development should pave the way for novel therapies.

De novo assembly of the olive fruit fly (Bactrocera oleae) genome with linked-reads and long-read technologies minimizes gaps and provides exceptional Y chromosome assembly.

BACKGROUND: The olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly's biology and proposing alternative control methods to pesticide use. RESULTS: Genomic DNA was sequenced from male and female Demokritos strain flies, maintained in the laboratory for over 45 years. We used short-, mate-pair-, and long-read sequencing technologies to generate a combined male-female genome assembly (GenBank accession GCA_001188975.2). Genomic DNA sequencing from male insects using 10x Genomics linked-reads technology followed by mate-pair and long-read scaffolding and gap-closing generated a highly contiguous 489 Mb genome with a scaffold N50 of 4.69 Mb and L50 of 30 scaffolds (GenBank accession GCA_001188975.4). RNA-seq data generated from 12 tissues and/or developmental stages allowed for genome annotation. Short reads from both males and females and the chromosome quotient method enabled identification of Y-chromosome scaffolds which were extensively validated by PCR. CONCLUSIONS: The high-quality genome generated represents a critical tool in olive fruit fly research. We provide an extensive RNA-seq data set, and genome annotation, critical towards gaining an insight into the biology of the olive fruit fly. In addition, elucidation of Y-chromosome sequences will advance our understanding of the Y-chromosome's organization, function and evolution and is poised to provide avenues for sterile insect technique approaches.

SARS-CoV-2 rapidly evolves lineage-specific phenotypic differences when passaged repeatedly in immune-naïve mice.

The persistence of SARS-CoV-2 despite the development of vaccines and a degree of herd immunity is partly due to viral evolution reducing vaccine and treatment efficacy. Serial infections of wild-type (WT) SARS-CoV-2 in Balb/c mice yield mouse-adapted strains with greater infectivity and mortality. We investigate if passaging unmodified B.1.351 (Beta) and B.1.617.2 (Delta) 20 times in K18-ACE2 mice, expressing the human ACE2 receptor, in a BSL-3 laboratory without selective pressures, drives human health-relevant evolution and if evolution is lineage-dependent. Late-passage virus causes more severe disease, at organism and lung tissue scales, with late-passage Delta demonstrating antibody resistance and interferon suppression. This resistance co-occurs with a de novo spike S371F mutation, linked with both traits. S371F, an Omicron-characteristic mutation, is co-inherited at times with spike E1182G per Nanopore sequencing, existing in different within-sample viral variants at others. Both S371F and E1182G are linked to mammalian GOLGA7 and ZDHHC5 interactions, which mediate viral-cell entry and antiviral response. This study demonstrates SARS-CoV-2's tendency to evolve with phenotypic consequences, its evolution varying by lineage, and suggests non-dominant quasi-species contribution.

Genomic and cytogenetic analysis of the Ceratitis capitata temperature-sensitive lethal region.

Genetic sexing strains (GSS) are an important tool in support of sterile insect technique (SIT) applications against insect pests and disease vectors. The yet unknown temperature-sensitive lethal (tsl) gene and the recently identified white pupae (wp) gene have been used as selectable markers in the most successful GSS developed so far, the Ceratitis capitata (medfly) VIENNA 8 GSS. The molecular identification of the tsl gene may open the way for its use as a marker for the development of GSS in other insect pests and disease vectors of SIT importance. Prior studies have already shown that the tsl gene is located on the right arm of chromosome 5, between the wp and Zw loci (tsl genomic region). In the present study, we used genomic, transcriptomic, bioinformatic, and cytogenetic approaches to characterize and analyze this genomic region in wild-type and tsl mutant medfly strains. Our results suggested the presence of 561 genes, with 322 of them carrying SNPs and/or insertion-deletion (indel) mutations in the tsl genomic region. Furthermore, comparative transcriptomic analysis indicated the presence of 32 differentially expressed genes, and bioinformatic analysis revealed the presence of 33 orthologs with a described heat-sensitive phenotype of Drosophila melanogaster in this region. These data can be used in functional genetic studies to identify the tsl gene(s) and the causal mutation(s) responsible for the temperature-sensitive lethal phenotype in medfly, and potentially additional genes causing a similar phenotype.

Panmixia in the American eel extends to its tropical range of distribution: Biological implications and policymaking challenges.

The American eel (Anguilla rostrata) has long been regarded as a panmictic fish and has been confirmed as such in the northern part of its range. In this paper, we tested for the first time whether panmixia extends to the tropical range of the species. To do so, we first assembled a reference genome (975 Mbp, 19 chromosomes) combining long (PacBio and Nanopore and short (Illumina paired-end) reads technologies to support both this study and future research. To test for population structure, we estimated genotype likelihoods from low-coverage whole-genome sequencing of 460 American eels, collected at 21 sampling sites (in seven geographic regions) ranging from Canada to Trinidad and Tobago. We estimated genetic distance between regions, performed ADMIXTURE-like clustering analysis and multivariate analysis, and found no evidence of population structure, thus confirming that panmixia extends to the tropical range of the species. In addition, two genomic regions with putative inversions were observed, both geographically widespread and present at similar frequencies in all regions. We discuss the implications of lack of genetic population structure for the species. Our results are key for the future genomic research in the American eel and the implementation of conservation measures throughout its geographic range. Additionally, our results can be applied to fisheries management and aquaculture of the species.

Concerted epithelial and stromal changes during progression of Barrett's Esophagus to invasive adenocarcinoma exposed by multi-scale, multi-omics analysis.

Esophageal adenocarcinoma arises from Barrett's esophagus, a precancerous metaplastic replacement of squamous by columnar epithelium in response to chronic inflammation. Multi-omics profiling, integrating single-cell transcriptomics, extracellular matrix proteomics, tissue-mechanics and spatial proteomics of 64 samples from 12 patients' paths of progression from squamous epithelium through metaplasia, dysplasia to adenocarcinoma, revealed shared and patient-specific progression characteristics. The classic metaplastic replacement of epithelial cells was paralleled by metaplastic changes in stromal cells, ECM and tissue stiffness. Strikingly, this change in tissue state at metaplasia was already accompanied by appearance of fibroblasts with characteristics of carcinoma-associated fibroblasts and of an NK cell-associated immunosuppressive microenvironment. Thus, Barrett's esophagus progresses as a coordinated multi-component system, supporting treatment paradigms that go beyond targeting cancerous cells to incorporating stromal reprogramming.

Analysis of early C2C12 myogenesis identifies stably and differentially expressed transcriptional regulators whose knock-down inhibits myoblast differentiation.

Myogenesis is a tightly controlled process involving the transcriptional activation and repression of thousands of genes. Although many components of the transcriptional network regulating the later phases of myogenesis have been identified, relatively few studies have described the transcriptional landscape during the first 24 h, when myoblasts commit to differentiate. Through dense temporal profiling of differentiating C2C12 myoblasts, we identify 193 transcriptional regulators (TRs) whose expression is significantly altered within the first 24 h of myogenesis. A high-content shRNA screen of 77 TRs involving 427 stable lines identified 42 genes whose knockdown significantly inhibits differentiation of C2C12 myoblasts. Of the TRs that were differentially expressed within the first 24 h, over half inhibited differentiation when knocked down, including known regulators of myogenesis (Myod1, Myog, and Myf5), as well as 19 TRs not previously associated with this process. Surprisingly, a similar proportion (55%) of shRNAs targeting TRs whose expression did not change also inhibited C2C12 myogenesis. We further show that a subset of these TRs inhibits myogenesis by downregulating expression of known regulatory and structural proteins. Our findings clearly illustrate that several TRs critical for C2C12 myogenesis are not differentially regulated, suggesting that approaches that focus functional studies on differentially-expressed transcripts will fail to provide a comprehensive view of this complex process.

Chromosome-level assembly reveals a putative Y-autosomal fusion in the sex determination system of the Greenland Halibut (Reinhardtius hippoglossoides).

Despite the commercial importance of Greenland Halibut (Reinhardtius hippoglossoides), important gaps still persist in our knowledge of this species, including its reproductive biology and sex determination mechanism. Here, we combined single-molecule sequencing of long reads (Pacific Sciences) with chromatin conformation capture sequencing (Hi-C) data to assemble the first chromosome-level reference genome for this species. The high-quality assembly encompassed more than 598 Megabases (Mb) assigned to 1594 scaffolds (scaffold N50 = 25 Mb) with 96% of its total length distributed among 24 chromosomes. Investigation of the syntenic relationship with other economically important flatfish species revealed a high conservation of synteny blocks among members of this phylogenetic clade. Sex determination analysis revealed that similar to other teleost fishes, flatfishes also exhibit a high level of plasticity and turnover in sex determination mechanisms. A low-coverage whole-genome sequence analysis of 198 individuals revealed that Greenland Halibut possesses a male heterogametic XY system and several putative candidate genes implied in the sex determination of this species. Our study also suggests for the first time in flatfishes that a putative Y-autosomal fusion could be associated with a reduction of recombination typical of the early steps of sex chromosome evolution.

A chromosome-anchored genome assembly for Lake Trout (Salvelinus namaycush).

Here, we present an annotated, chromosome-anchored, genome assembly for Lake Trout (Salvelinus namaycush) - a highly diverse salmonid species of notable conservation concern and an excellent model for research on adaptation and speciation. We leveraged Pacific Biosciences long-read sequencing, paired-end Illumina sequencing, proximity ligation (Hi-C) sequencing, and a previously published linkage map to produce a highly contiguous assembly composed of 7378 contigs (contig N50 = 1.8 Mb) assigned to 4120 scaffolds (scaffold N50 = 44.975 Mb). Long read sequencing data were generated using DNA from a female double haploid individual. 84.7% of the genome was assigned to 42 chromosome-sized scaffolds and 93.2% of Benchmarking Universal Single Copy Orthologues were recovered, putting this assembly on par with the best currently available salmonid genomes. Estimates of genome size based on k-mer frequency analysis were highly similar to the total size of the finished genome, suggesting that the entirety of the genome was recovered. A mitochondrial genome assembly was also produced. Self-versus-self synteny analysis allowed us to identify homeologs resulting from the salmonid specific autotetraploid event (Ss4R) as well as regions exhibiting delayed rediploidization. Alignment with three other salmonid genomes and the Northern Pike (Esox lucius) genome also allowed us to identify homologous chromosomes in related taxa. We also generated multiple resources useful for future genomic research on Lake Trout, including a repeat library and a sex-averaged recombination map. A novel RNA sequencing data set for liver tissue was also generated in order to produce a publicly available set of annotations for 49,668 genes and pseudogenes. Potential applications of these resources to population genetics and the conservation of native populations are discussed.

The living microarray: a high-throughput platform for measuring transcription dynamics in single cells.

BACKGROUND: Current methods of measuring transcription in high-throughput have led to significant improvements in our knowledge of transcriptional regulation and Systems Biology. However, endpoint measurements obtained from methods that pool populations of cells are not amenable to studying time-dependent processes that show cell heterogeneity. RESULTS: Here we describe a high-throughput platform for measuring transcriptional changes in real time in single mammalian cells. By using reverse transfection microarrays we are able to transfect fluorescent reporter plasmids into 600 independent clusters of cells plated on a single microscope slide and image these clusters every 20 minutes. We use a fast-maturing, destabilized and nuclear-localized reporter that is suitable for automated segmentation to accurately measure promoter activity in single cells. We tested this platform with synthetic drug-inducible promoters that showed robust induction over 24 hours. Automated segmentation and tracking of over 11 million cell images during this period revealed that cells display substantial heterogeneity in their responses to the applied treatment, including a large proportion of transfected cells that do not respond at all. CONCLUSIONS: The results from our single-cell analysis suggest that methods that measure average cellular responses, such as DNA microarrays, RT-PCR and chromatin immunoprecipitation, characterize a response skewed by a subset of cells in the population. Our method is scalable and readily adaptable to studying complex systems, including cell proliferation, differentiation and apoptosis.

Haplotype-resolved de novo assembly of the Vero cell line genome.

The Vero cell line is the most used continuous cell line for viral vaccine manufacturing with more than 40 years of accumulated experience in the vaccine industry. Additionally, the Vero cell line has shown a high affinity for infection by MERS-CoV, SARS-CoV, and recently SARS-CoV-2, emerging as an important discovery and screening tool to support the global research and development efforts in this COVID-19 pandemic. However, the lack of a reference genome for the Vero cell line has limited our understanding of host-virus interactions underlying such affinity of the Vero cell towards key emerging pathogens, and more importantly our ability to redesign high-yield vaccine production processes using Vero genome editing. In this paper, we present an annotated highly contiguous 2.9 Gb assembly of the Vero cell genome. In addition, several viral genome insertions, including Adeno-associated virus serotypes 3, 4, 7, and 8, have been identified, giving valuable insights into quality control considerations for cell-based vaccine production systems. Variant calling revealed that, in addition to interferon, chemokines, and caspases-related genes lost their functions. Surprisingly, the ACE2 gene, which was previously identified as the host cell entry receptor for SARS-CoV and SARS-CoV-2, also lost function in the Vero genome due to structural variations.

White pupae phenotype of tephritids is caused by parallel mutations of a MFS transporter.

Mass releases of sterilized male insects, in the frame of sterile insect technique programs, have helped suppress insect pest populations since the 1950s. In the major horticultural pests Bactrocera dorsalis, Ceratitis capitata, and Zeugodacus cucurbitae, a key phenotype white pupae (wp) has been used for decades to selectively remove females before releases, yet the gene responsible remained unknown. Here, we use classical and modern genetic approaches to identify and functionally characterize causal wp- mutations in these distantly related fruit fly species. We find that the wp phenotype is produced by parallel mutations in a single, conserved gene. CRISPR/Cas9-mediated knockout of the wp gene leads to the rapid generation of white pupae strains in C. capitata and B. tryoni. The conserved phenotype and independent nature of wp- mutations suggest this technique can provide a generic approach to produce sexing strains in other major medical and agricultural insect pests.

Invasive growth associated with cold-inducible RNA-binding protein expression drives recurrence of surgically resected brain metastases.

BACKGROUND: Sixty percent of surgically resected brain metastases (BrM) recur within 1 year. These recurrences have long been thought to result from the dispersion of cancer cells during surgery. We tested the alternative hypothesis that invasion of cancer cells into the adjacent brain plays a significant role in local recurrence and shortened overall survival. METHODS: We determined the invasion pattern of 164 surgically resected BrM and correlated with local recurrence and overall survival. We performed single-cell RNA sequencing (scRNAseq) of >15,000 cells from BrM and adjacent brain tissue. Validation of targets was performed with a novel cohort of BrM patient-derived xenografts (PDX) and patient tissues. RESULTS: We demonstrate that invasion of metastatic cancer cells into the adjacent brain is associated with local recurrence and shortened overall survival. scRNAseq of paired tumor and adjacent brain samples confirmed the existence of invasive cancer cells in the tumor-adjacent brain. Analysis of these cells identified cold-inducible RNA-binding protein (CIRBP) overexpression in invasive cancer cells compared to cancer cells located within the metastases. Applying PDX models that recapitulate the invasion pattern observed in patients, we show that CIRBP is overexpressed in highly invasive BrM and is required for efficient invasive growth in the brain. CONCLUSIONS: These data demonstrate peritumoral invasion as a driver of treatment failure in BrM that is functionally mediated by CIRBP. These findings improve our understanding of the biology underlying postoperative treatment failure and lay the groundwork for rational clinical trial development based upon invasion pattern in surgically resected BrM.

The transcriptional landscape of Shh medulloblastoma.

Sonic hedgehog medulloblastoma encompasses a clinically and molecularly diverse group of cancers of the developing central nervous system. Here, we use unbiased sequencing of the transcriptome across a large cohort of 250 tumors to reveal differences among molecular subtypes of the disease, and demonstrate the previously unappreciated importance of non-coding RNA transcripts. We identify alterations within the cAMP dependent pathway (GNAS, PRKAR1A) which converge on GLI2 activity and show that 18% of tumors have a genetic event that directly targets the abundance and/or stability of MYCN. Furthermore, we discover an extensive network of fusions in focally amplified regions encompassing GLI2, and several loss-of-function fusions in tumor suppressor genes PTCH1, SUFU and NCOR1. Molecular convergence on a subset of genes by nucleotide variants, copy number aberrations, and gene fusions highlight the key roles of specific pathways in the pathogenesis of Sonic hedgehog medulloblastoma and open up opportunities for therapeutic intervention.

Addressable nanowell arrays formed using reversibly sealable hybrid elastomer-metal stencils.

There are two major array formats used in life science research and biomedical analysis. The first is the microwell plate format with millimeter-sized wells each with microliter capacity addressed individually and repeatedly during experiments. The second is the microarray format with micrometer-sized spots that are patterned initially but not addressable individually thereafter. Here, we present an addressable nanoliter-well plate with micrometer sized wells that combines the advantages of the two array formats. The nanowells are formed by reversibly sealing a steel stencil featuring an array of micrometer-scale openings to an optically transparent substrate. The nanowells have a capacity of approximately 1 nL, are approximately 140 microm in diameter, and are arrayed at a density of 1600 wells cm(-2). A soft polymer is patterned photolithographically around each opening so as to form a microgasket for pressure sensitive, liquid tight, and reversible sealing to any type of smooth substrate, either hydrophilic or hydrophobic. The rigidity of the steel prevents the distortion that occurs in soft, all-polymeric stencils and permits accurate registration across the entire array, which in turn allows for repeated, individual addressing of wells using an inkjet spotter. The stencils are used to pattern cells, make protein microarrays, and create nanowells on surfaces to study reverse transfection by first spotting plasmids encoding fluorescent proteins into the wells, seeding cells, and monitoring the transfection of the cells in real time using time-lapse imaging. The hybrid elastomer-metal stencils (HEMSs) are versatile and useful for multiplexed analysis of drugs, biomolecules, and cells with microarray density.

Methodologies for Transcript Profiling Using Long-Read Technologies.

RNA sequencing using next-generation sequencing technologies (NGS) is currently the standard approach for gene expression profiling, particularly for large-scale high-throughput studies. NGS technologies comprise high throughput, cost efficient short-read RNA-Seq, while emerging single molecule, long-read RNA-Seq technologies have enabled new approaches to study the transcriptome and its function. The emerging single molecule, long-read technologies are currently commercially available by Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT), while new methodologies based on short-read sequencing approaches are also being developed in order to provide long range single molecule level information-for example, the ones represented by the 10x Genomics linked read methodology. The shift toward long-read sequencing technologies for transcriptome characterization is based on current increases in throughput and decreases in cost, making these attractive for de novo transcriptome assembly, isoform expression quantification, and in-depth RNA species analysis. These types of analyses were challenging with standard short sequencing approaches, due to the complex nature of the transcriptome, which consists of variable lengths of transcripts and multiple alternatively spliced isoforms for most genes, as well as the high sequence similarity of highly abundant species of RNA, such as rRNAs. Here we aim to focus on single molecule level sequencing technologies and single-cell technologies that, combined with perturbation tools, allow the analysis of complete RNA species, whether short or long, at high resolution. In parallel, these tools have opened new ways in understanding gene functions at the tissue, network, and pathway levels, as well as their detailed functional characterization. Analysis of the epi-transcriptome, including RNA methylation and modification and the effects of such modifications on biological systems is now enabled through direct RNA sequencing instead of classical indirect approaches. However, many difficulties and challenges remain, such as methodologies to generate full-length RNA or cDNA libraries from all different species of RNAs, not only poly-A containing transcripts, and the identification of allele-specific transcripts due to current error rates of single molecule technologies, while the bioinformatics analysis on long-read data for accurate identification of 5' and 3' UTRs is still in development.

Model-based analysis of sample index hopping reveals its widespread artifacts in multiplexed single-cell RNA-sequencing.

Index hopping is the main cause of incorrect sample assignment of sequencing reads in multiplexed pooled libraries. We introduce a statistical model for estimating the sample index-hopping rate in multiplexed droplet-based single-cell RNA-seq data and for probabilistic inference of the true sample of origin of hopped reads. We analyze several datasets and estimate the sample index hopping probability to range between 0.003-0.009, a small number that counter-intuitively gives rise to a large fraction of phantom molecules - the fraction of phantom molecules exceeds 8% in more than 25% of samples and reaches as high as 85% in low-complexity samples. Phantom molecules lead to widespread complications in downstream analyses, including transcriptome mixing across cells, emergence of phantom copies of cells from other samples, and misclassification of empty droplets as cells. We demonstrate that our approach can correct for these artifacts by accurately purging the majority of phantom molecules from the data.

Single-cell analysis of human adipose tissue identifies depot and disease specific cell types.

The complex relationship between metabolic disease risk and body fat distribution in humans involves cellular characteristics which are specific to body fat compartments. Here we show depot-specific differences in the stromal vascual fraction of visceral and subcutaneous adipose tissue by performing single-cell RNA sequencing of tissue specimen from obese individuals. We characterize multiple immune cells, endothelial cells, fibroblasts, adipose and hematopoietic stem cell progenitors. Subpopulations of adipose-resident immune cells are metabolically active and associated with metabolic disease status and those include a population of potential dysfunctional CD8+ T cells expressing metallothioneins. We identify multiple types of adipocyte progenitors that are common across depots, including a subtype enriched in individuals with type 2 diabetes. Depot-specific analysis reveals a class of adipocyte progenitors unique to visceral adipose tissue, which shares common features with beige preadipocytes. Our human single-cell transcriptome atlas across fat depots provides a resource to dissect functional genomics of metabolic disease.