Haplotyping germline and cancer genomes with high-throughput linked-read sequencing

Single-cell RNA sequencing reveals heterogeneity and differential expression of the maternal-fetal interface during ICP and normal pregnancy

OBJECTIVE

Intrahepatic cholestasis of pregnancy (ICP) can cause adverse perinatal outcomes. Previous studies have demonstrated that the placenta of an ICP pregnancy differs in morphology and gene expression from the placenta of a normal pregnancy. To date, however, the genetic mechanism by which ICP affects the placenta is poorly understood. Therefore, the aim of this study was to investigate the differences in main cell types, gene signatures, cell ratio, and functional changes in the placenta between ICP and normal pregnancy.

METHODS

Single-cell RNA sequencing (scRNA-seq) technology was used to detect the gene expression of all cells at the placental maternal-fetal interface. Two individuals were analyzed - one with ICP and one without ICP. The classification of cell types was determined by a graph-based clustering algorithm. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the R software phyper () function and DAVID website. The differentially expressed genes (DEGs) encoding transcription factors (TFs) were identified using getorf and DIAMOND software.

RESULTS

We identified 14 cell types and 22 distinct cell subtypes that showed unique functional properties. Additionally, we found differences in the proportions of fibroblasts 1, helper T (Th) cells, extravillous trophoblasts, and villous cytotrophoblasts, and we observed heterogeneity of gene expression between ICP and control placentas. Furthermore, we identified 263 DEGs that belonged to TF families, including zf-C2H2, HMGI/HMGY, and Homeobox. In addition, 28 imprinted genes were preferentially expressed in specific cell types, such as PEG3 and PEG10 in trophoblasts as well as DLK1 and DIO3 in fibroblasts.

CONCLUSIONS

Our results revealed the differences in cell-type ratios, gene expression, and functional changes between ICP and normal placentas, and heterogeneity was found among cell subgroups. Hence, the imbalance of various cell types affects placental activity to varying degrees, indicating the complexity of the cell networks that form the placental tissue system, and this alteration of placental function is associated with adverse events in the perinatal period.

ABCA7-dependent induction of neuropeptide Y is required for synaptic resilience in Alzheimer's disease through BDNF/NGFR signaling

Genetic variants in ABCA7, an Alzheimer's disease (AD)-associated gene, elevate AD risk, yet its functional relevance to the etiology is unclear. We generated a CRISPR-Cas9-mediated abca7 knockout zebrafish to explore ABCA7's role in AD. Single-cell transcriptomics in heterozygous abca7+/- knockout combined with Aβ42 toxicity revealed that ABCA7 is crucial for neuropeptide Y (NPY), brain-derived neurotrophic factor (BDNF), and nerve growth factor receptor (NGFR) expressions, which are crucial for synaptic integrity, astroglial proliferation, and microglial prevalence. Impaired NPY induction decreased BDNF and synaptic density, which are rescuable with ectopic NPY. In induced pluripotent stem cell-derived human neurons exposed to Aβ42, ABCA7-/- suppresses NPY. Clinical data showed reduced NPY in AD correlated with elevated Braak stages, genetic variants in NPY associated with AD, and epigenetic changes in NPY, NGFR, and BDNF promoters linked to ABCA7 variants. Therefore, ABCA7-dependent NPY signaling via BDNF-NGFR maintains synaptic integrity, implicating its impairment in increased AD risk through reduced brain resilience.

Whole-genome de novo sequencing reveals genomic variants associated with differences of sex development in SRY negative pigs

BACKGROUND

Differences of sex development (DSD) are congenital conditions in which chromosomal, gonadal, or phenotypic sex is atypical. In more than 50% of human DSD cases, a molecular diagnosis is not available. In intensively farmed pig populations, the incidence of XX DSD pigs is relatively high, leading to economic losses for pig breeders. Interestingly, in the majority of 38, XX DSD pigs, gonads still develop into testis-like structures or ovotestes despite the absence of the testis-determining gene (SRY). However, the current understanding of the molecular background of XX DSD pigs remains limited.

METHODS

Anatomical and histological characteristics of XX DSD pigs were analysed using necropsy and HE staining. We employed whole-genome sequencing (WGS) with 10× Genomics technology and used de novo assembly methodology to study normal female and XX DSD pigs. Finally, the identified variants were validated in 32 XX DSD pigs, and the expression levels of the candidate variants in the gonads of XX DSD pigs were further examined.

RESULTS

XX DSD pigs are characterised by the intersex reproductive organs and the absence of germ cells in the seminiferous tubules of the gonads. We identified 4,950 single-nucleotide polymorphisms (SNPs) from non-synonymous mutations in XX DSD pigs. Cohort validation results highlighted two specific SNPs, "c.218T > C" in the "Interferon-induced transmembrane protein 1 gene (IFITM1)" and "c.1043C > G" in the "Newborn ovary homeobox gene (NOBOX)", which were found exclusively in XX DSD pigs. Moreover, we verified 14 candidate structural variants (SVs) from 1,474 SVs, identifying a 70 bp deletion fragment in intron 5 of the WW domain-containing oxidoreductase gene (WWOX) in 62.5% of XX DSD pigs. The expression levels of these three candidate genes in the gonads of XX DSD pigs were significantly different from those of normal female pigs.

CONCLUSION

The nucleotide changes of IFITM1 (c.218T > C), NOBOX (c.1043 C > G), and a 70 bp deletion fragment of the WWOX were the most dominant variants among XX DSD pigs. This study provides a theoretical basis for better understanding the molecular background of XX DSD pigs. DSD are conditions affecting development of the gonads or genitalia. These disorders can happen in many different types of animals, including pigs, goats, dogs, and people. In people, DSD happens in about 0.02-0.13% of births, and in pigs, the rate is between 0.08% and 0.75%. Pigs have a common type of DSD where the animal has female chromosomes (38, XX) but no SRY gene, which is usually found on the Y chromosome in males. XX DSD pigs may look like both males and females on the outside and have testis-like or ovotestis (a mix of ovary and testis) gonads inside. XX DSD pigs often lead to not being able to have piglets, slower growth, lower chance of survival, and poorer meat quality. Here, we used a method called whole-genome de novo sequencing to look for variants in the DNA of XX DSD pigs. We then checked these differences in a larger group of pigs. Our results reveal the nucleotide changes in IFITM1 (c.218T > C), NOBOX (c.1043 C > G), and a 70 bp deletion fragment in intron 5 of the WWOX, all linked to XX DSD pigs. The expression levels of these three genes were also different in the gonads of XX DSD pigs compared to normal female pigs. These variants are expected to serve as valuable molecular markers for XX DSD pigs. Because pigs are a lot like humans in their genes, physiology, and body structure, this research could help us learn more about what causes DSD in people.

An overview of recent technological developments in bovine genomics

Cattle are regarded as highly valuable animals because of their milk, beef, dung, fur, and ability to draft. The scientific community has tried a number of strategies to improve the genetic makeup of bovine germplasm. To ensure higher returns for the dairy and beef industries, researchers face their greatest challenge in improving commercially important traits. One of the biggest developments in the last few decades in the creation of instruments for cattle genetic improvement is the discovery of the genome. Breeding livestock is being revolutionized by genomic selection made possible by the availability of medium- and high-density single nucleotide polymorphism (SNP) arrays coupled with sophisticated statistical techniques. It is becoming easier to access high-dimensional genomic data in cattle. Continuously declining genotyping costs and an increase in services that use genomic data to increase return on investment have both made a significant contribution to this. The field of genomics has come a long way thanks to groundbreaking discoveries such as radiation-hybrid mapping, in situ hybridization, synteny analysis, somatic cell genetics, cytogenetic maps, molecular markers, association studies for quantitative trait loci, high-throughput SNP genotyping, whole-genome shotgun sequencing to whole-genome mapping, and genome editing. These advancements have had a significant positive impact on the field of cattle genomics. This manuscript aimed to review recent advances in genomic technologies for cattle breeding and future prospects in this field.

Implementing next-generation sequencing for diagnosis and management of hereditary hearing impairment: a comprehensive review

INTRODUCTION

Sensorineural hearing impairment (SNHI), a common childhood disorder with heterogeneous genetic causes, can lead to delayed language development and psychosocial problems. Next-generation sequencing (NGS) offers high-throughput screening and high-sensitivity detection of genetic etiologies of SNHI, enabling clinicians to make informed medical decisions, provide tailored treatments, and improve prognostic outcomes.

AREAS COVERED

This review covers the diverse etiologies of HHI and the utility of different NGS modalities (targeted sequencing and whole exome/genome sequencing), and includes HHI-related studies on newborn screening, genetic counseling, prognostic prediction, and personalized treatment. Challenges such as the trade-off between cost and diagnostic yield, detection of structural variants, and exploration of the non-coding genome are also highlighted.

EXPERT OPINION

In the current landscape of NGS-based diagnostics for HHI, there are both challenges (e.g. detection of structural variants and non-coding genome variants) and opportunities (e.g. the emergence of medical artificial intelligence tools). The authors advocate the use of technological advances such as long-read sequencing for structural variant detection, multi-omics analysis for non-coding variant exploration, and medical artificial intelligence for pathogenicity assessment and outcome prediction. By integrating these innovations into clinical practice, precision medicine in the diagnosis and management of HHI can be further improved.

Somatic mutation phasing and haplotype extension using linked-reads in multiple myeloma

Somatic mutation phasing informs our understanding of cancer-related events, like driver mutations. We generated linked-read whole genome sequencing data for 23 samples across disease stages from 14 multiple myeloma (MM) patients and systematically assigned somatic mutations to haplotypes using linked-reads. Here, we report the reconstructed cancer haplotypes and phase blocks from several MM samples and show how phase block length can be extended by integrating samples from the same individual. We also uncover phasing information in genes frequently mutated in MM, including DIS3, HIST1H1E, KRAS, NRAS, and TP53, phasing 79.4% of 20,705 high-confidence somatic mutations. In some cases, this enabled us to interpret clonal evolution models at higher resolution using pairs of phased somatic mutations. For example, our analysis of one patient suggested that two NRAS hotspot mutations occurred on the same haplotype but were independent events in different subclones. Given sufficient tumor purity and data quality, our framework illustrates how haplotype-aware analysis of somatic mutations in cancer can be beneficial for some cancer cases.

Resolving the 22q11.2 deletion using CTLR-Seq reveals chromosomal rearrangement mechanisms and individual variance in breakpoints

We developed a generally applicable method, CRISPR/Cas9-targeted long-read sequencing (CTLR-Seq), to resolve, haplotype-specifically, the large and complex regions in the human genome that had been previously impenetrable to sequencing analysis, such as large segmental duplications (SegDups) and their associated genome rearrangements. CTLR-Seq combines in vitro Cas9-mediated cutting of the genome and pulse-field gel electrophoresis to isolate intact large (i.e., up to 2,000 kb) genomic regions that encompass previously unresolvable genomic sequences. These targets are then sequenced (amplification-free) at high on-target coverage using long-read sequencing, allowing for their complete sequence assembly. We applied CTLR-Seq to the SegDup-mediated rearrangements that constitute the boundaries of, and give rise to, the 22q11.2 Deletion Syndrome (22q11DS), the most common human microdeletion disorder. We then performed de novo assembly to resolve, at base-pair resolution, the full sequence rearrangements and exact chromosomal breakpoints of 22q11.2DS (including all common subtypes). Across multiple patients, we found a high degree of variability for both the rearranged SegDup sequences and the exact chromosomal breakpoint locations, which coincide with various transposons within the 22q11.2 SegDups, suggesting that 22q11DS can be driven by transposon-mediated genome recombination. Guided by CTLR-Seq results from two 22q11DS patients, we performed three-dimensional chromosomal folding analysis for the 22q11.2 SegDups from patient-derived neurons and astrocytes and found chromosome interactions anchored within the SegDups to be both cell type-specific and patient-specific. Lastly, we demonstrated that CTLR-Seq enables cell-type specific analysis of DNA methylation patterns within the deletion haplotype of 22q11DS.

AsmMix: an efficient haplotype-resolved hybrid de novo genome assembling pipeline

Accurate haplotyping facilitates distinguishing allele-specific expression, identifying cis-regulatory elements, and characterizing genomic variations, which enables more precise investigations into the relationship between genotype and phenotype. Recent advances in third-generation single-molecule long read and synthetic co-barcoded read sequencing techniques have harnessed long-range information to simplify the assembly graph and improve assembly genomic sequence. However, it remains methodologically challenging to reconstruct the complete haplotypes due to high sequencing error rates of long reads and limited capturing efficiency of co-barcoded reads. We here present a pipeline, AsmMix, for generating both contiguous and accurate diploid genomes. It first assembles co-barcoded reads to generate accurate haplotype-resolved assemblies that may contain many gaps, while the long-read assembly is contiguous but susceptible to errors. Then two assembly sets are integrated into haplotype-resolved assemblies with reduced misassembles. Through extensive evaluation on multiple synthetic datasets, AsmMix consistently demonstrates high precision and recall rates for haplotyping across diverse sequencing platforms, coverage depths, read lengths, and read accuracies, significantly outperforming other existing tools in the field. Furthermore, we validate the effectiveness of our pipeline using a human whole genome dataset (HG002), and produce highly contiguous, accurate, and haplotype-resolved assemblies. These assemblies are evaluated using the GIAB benchmarks, confirming the accuracy of variant calling. Our results demonstrate that AsmMix offers a straightforward yet highly efficient approach that effectively leverages both long reads and co-barcoded reads for haplotype-resolved assembly.

Neotelomeres and telomere-spanning chromosomal arm fusions in cancer genomes revealed by long-read sequencing

Alterations in the structure and location of telomeres are pivotal in cancer genome evolution. Here, we applied both long-read and short-read genome sequencing to assess telomere repeat-containing structures in cancers and cancer cell lines. Using long-read genome sequences that span telomeric repeats, we defined four types of telomere repeat variations in cancer cells: neotelomeres where telomere addition heals chromosome breaks, chromosomal arm fusions spanning telomere repeats, fusions of neotelomeres, and peri-centromeric fusions with adjoined telomere and centromere repeats. These results provide a framework for the systematic study of telomeric repeats in cancer genomes, which could serve as a model for understanding the somatic evolution of other repetitive genomic elements.

Simultaneous de novo calling and phasing of genetic variants at chromosome-scale using NanoStrand-seq

The successful accomplishment of the first telomere-to-telomere human genome assembly, T2T-CHM13, marked a milestone in achieving completeness of the human reference genome. The upcoming era of genome study will focus on fully phased diploid genome assembly, with an emphasis on genetic differences between individual haplotypes. Most existing sequencing approaches only achieved localized haplotype phasing and relied on additional pedigree information for further whole-chromosome scale phasing. The short-read-based Strand-seq method is able to directly phase single nucleotide polymorphisms (SNPs) at whole-chromosome scale but falls short when it comes to phasing structural variations (SVs). To shed light on this issue, we developed a Nanopore sequencing platform-based Strand-seq approach, which we named NanoStrand-seq. This method allowed for de novo SNP calling with high precision (99.52%) and acheived a superior phasing accuracy (0.02% Hamming error rate) at whole-chromosome scale, a level of performance comparable to Strand-seq for haplotype phasing of the GM12878 genome. Importantly, we demonstrated that NanoStrand-seq can efficiently resolve the MHC locus, a highly polymorphic genomic region. Moreover, NanoStrand-seq enabled independent direct calling and phasing of deletions and insertions at whole-chromosome level; when applied to long genomic regions of SNP homozygosity, it outperformed the strategy that combined Strand-seq with bulk long-read sequencing. Finally, we showed that, like Strand-seq, NanoStrand-seq was also applicable to primary cultured cells. Together, here we provided a novel methodology that enabled interrogation of a full spectrum of haplotype-resolved SNPs and SVs at whole-chromosome scale, with broad applications for species with diploid or even potentially polypoid genomes.

Targeted Linked-Read Sequencing for Direct Haplotype Phasing of Parental GJB2/SLC26A4 Alleles: A Universal and Dependable Noninvasive Prenatal Diagnosis Method Applied to Autosomal Recessive Nonsyndromic Hearing Loss in At-Risk Families

Noninvasive prenatal diagnosis (NIPD) for autosomal recessive nonsyndromic hearing loss (ARNSHL) has been rarely reported until recent years. Additionally, the existing method can not be used for challenging genome loci (eg, copy number variations, deletions, inversions, or gene recombinants) or on families without proband genotype. This study assessed the performance of relative haplotype dosage analysis (RHDO)-based NIPD for identifying fetal genotyping in pregnancies at risk of ARNSHL. Fifty couples carrying pathogenic variants associated with ARNSHL in either GJB2 or SLC26A4 were recruited. The RHDO-based targeted linked-read sequencing combined with whole gene coverage probes was used to genotype the fetal cell-free DNA of 49 families who met the quality control standard. Fetal amniocyte samples were genotyped using invasive prenatal diagnosis (IPD) to assess the performance of NIPD. The NIPD results showed 100% (49/49) concordance with those obtained through IPD. Two families with copy number variation and recombination were also successfully identified. Sufficient specific informative single-nucleotide polymorphisms for haplotyping, as well as the fetal cell-free DNA concentration and sequencing depth, are prerequisites for RHDO-based NIPD. This method has the merits of covering the entire genes of GJB2 and SLC26A4, qualifying for copy number variation and recombination analysis with remarkable sensitivity and specificity. Therefore, it has clinical potential as an alternative to traditional IPD for ARNSHL.

Extant and extinct bilby genomes combined with Indigenous knowledge improve conservation of a unique Australian marsupial

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

Fly Me to the Micron: Microtechnologies for Drosophila Research

Multicellular model organisms, such as Drosophila melanogaster (fruit fly), are frequently used in a myriad of biological research studies due to their biological significance and global standardization. However, traditional tools used in these studies generally require manual handling, subjective phenotyping, and bulk treatment of the organisms, resulting in laborious experimental protocols with limited accuracy. Advancements in microtechnology over the course of the last two decades have allowed researchers to develop automated, high-throughput, and multifunctional experimental tools that enable novel experimental paradigms that would not be possible otherwise. We discuss recent advances in microtechnological systems developed for small model organisms using D. melanogaster as an example. We critically analyze the state of the field by comparing the systems produced for different applications. Additionally, we suggest design guidelines, operational tips, and new research directions based on the technical and knowledge gaps in the literature. This review aims to foster interdisciplinary work by helping engineers to familiarize themselves with model organisms while presenting the most recent advances in microengineering strategies to biologists.

Vole genomics links determinate and indeterminate growth of teeth

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.

Genome-wide association analysis uncovers rice blast resistance alleles of Ptr and Pia

A critical step to maximize the usefulness of genome-wide association studies (GWAS) in plant breeding is the identification and validation of candidate genes underlying genetic associations. This is of particular importance in disease resistance breeding where allelic variants of resistance genes often confer resistance to distinct populations, or races, of a pathogen. Here, we perform a genome-wide association analysis of rice blast resistance in 500 genetically diverse rice accessions. To facilitate candidate gene identification, we produce de-novo genome assemblies of ten rice accessions with various rice blast resistance associations. These genome assemblies facilitate the identification and functional validation of novel alleles of the rice blast resistance genes Ptr and Pia. We uncover an allelic series for the unusual Ptr rice blast resistance gene, and additional alleles of the Pia resistance genes RGA4 and RGA5. By linking these associations to three thousand rice genomes we provide a useful tool to inform future rice blast breeding efforts. Our work shows that GWAS in combination with whole-genome sequencing is a powerful tool for gene cloning and to facilitate selection of specific resistance alleles for plant breeding.

De novo assembly and analysis of Sonneratia ovata genome and population analysis

Mangroves are an important part of coastal and estuarine ecosystems where they serve as nurseries for marine species and prevent coastal erosion. Here we report the genome of Sonneratia ovata, which is a true mangrove that grows in estuarine environments and can tolerate moderate salt exposure. We sequenced the S. ovata genome and assembled it into chromosome-level scaffolds through the use of Hi-C. The genome is 212.3 Mb and contains 12 chromosomes that range in size from 12.2 to 23.2 Mb. Annotation identified 29,829 genes with a BUSCO completeness of 95.9%. We identified salt genes and found copy number expansion of salt genes such as ADP-ribosylation factor 1, and elongation factor 1-alpha. Population analysis identified a low level of genetic variation and a lack of population structure within S. ovata.

Targeted phasing of 2-200 kilobase DNA fragments with a short-read sequencer and a single-tube linked-read library method

In the human genome, heterozygous sites refer to genomic positions with a different allele or nucleotide variant on the maternal and paternal chromosomes. Resolving these allelic differences by chromosomal copy, also known as phasing, is achievable on a short-read sequencer when using a library preparation method that captures long-range genomic information. TELL-Seq is a library preparation that captures long-range genomic information with the aid of molecular identifiers (barcodes). The same barcode is used to tag the reads derived from the same long DNA fragment within a range of up to 200 kilobases (kb), generating linked-reads. This strategy can be used to phase an entire genome. Here, we introduce a TELL-Seq protocol developed for targeted applications, enabling the phasing of enriched loci of varying sizes, purity levels, and heterozygosity. To validate this protocol, we phased 2-200 kb loci enriched with different methods: CRISPR/Cas9-mediated excision coupled with pulse-field electrophoresis for the longest fragments, CRISPR/Cas9-mediated protection from exonuclease digestion for mid-size fragments, and long PCR for the shortest fragments. All selected loci have known clinical relevance: BRCA1, BRCA2, MLH1, MSH2, MSH6, APC, PMS2, SCN5A-SCN10A, and PKI3CA. Collectively, the analyses show that TELL-Seq can accurately phase 2-200 kb targets using a short-read sequencer.

Multiomics Analyses Provide New Insight into Genetic Variation of Reproductive Adaptability in Tibetan Sheep

Domestication and artificial selection during production-oriented breeding have greatly shaped the level of genomic variability in sheep. However, the genetic variation associated with increased reproduction remains elusive. Here, two groups of samples from consecutively monotocous and polytocous sheep were collected for genome-wide association, transcriptomic, proteomic, and metabolomic analyses to explore the genetic variation in fecundity in Tibetan sheep. Genome-wide association study revealed strong associations between BMPR1B (p.Q249R) and litter size, as well as between PAPPA and lambing interval; these findings were validated in 1,130 individuals. Furthermore, we constructed the first single-cell atlas of Tibetan sheep ovary tissues and identified a specific mural granulosa cell subtype with PAPPA-specific expression and differential expression of BMPR1B between the two groups. Bulk RNA-seq indicated that BMPR1B and PAPPA expressions were similar between the two groups of sheep. 3D protein structure prediction and coimmunoprecipitation analysis indicated that mutation and mutually exclusive exons of BMPR1B are the main mechanisms for prolific Tibetan sheep. We propose that PAPPA is a key gene for stimulating ovarian follicular growth and development, and steroidogenesis. Our work reveals the genetic variation in reproductive performance in Tibetan sheep, providing insights and valuable genetic resources for the discovery of genes and regulatory mechanisms that improve reproductive success.

TREM2-IGF1 Mediated Glucometabolic Enhancement Underlies Microglial Neuroprotective Properties During Ischemic Stroke

Microglia, the major resident immune cells in the central nervous system, serve as the frontline soldiers against cerebral ischemic injuries, possibly along with metabolic alterations. However, signaling pathways involved in the regulation of microglial immunometabolism in ischemic stroke remain to be further elucidated. In this study, using single-nuclei RNA sequencing, a microglial subcluster up-regulated in ischemic brain tissues is identified, with high expression of Igf1 and Trem2, neuroprotective transcriptional signature and enhanced oxidative phosphorylation. Microglial depletion by PLX3397 exacerbates ischemic brain damage, which is reversed by repopulating the microglia with high Igf1 and Trem2 phenotype. Mechanistically, Igf1 serves as one of the major down-stream molecules of Trem2, and Trem2-Igf1 signaling axis regulates microglial functional and metabolic profiles, exerting neuroprotective effects on ischemic stroke. Overexpression of Igf1 and supplementation of cyclocreatine restore microglial glucometabolic levels and cellular functions even in the absence of Trem2. These findings suggest that Trem2-Igf1 signaling axis reprograms microglial immunometabolic profiles and shifts microglia toward a neuroprotective phenotype, which has promising therapeutic potential in treating ischemic stroke.

Chromosome-scale genome assembly of Lepus oiostolus (Lepus, Leporidae)

Lepus oiostolus (L. oiostolus) is a species endemic to the Qinghai-Tibet Plateau. However, the absence of a reference genome limits genetic studies. Here, we reported a high-quality L. oiostolus genome assembly, with scaffolds anchored to 24 chromosomes and a total assembled length of 2.80 Gb (contig N50 = 64.25 Mb). Genomic annotation uncovered 22,295 protein-coding genes and identified 49.84% of the sequences as transposable elements. Long interspersed nuclear elements (LINEs) constitute a high proportion of the genome. Our study is at the first time to report the chromosome-scale genome for the species of the L. oiostolus. It provides a valuable genomic resource for future research on the evolution of the Leporidae.

Pasa: leveraging population pangenome graph to scaffold prokaryote genome assemblies

Whole genome sequencing has increasingly become the essential method for studying the genetic mechanisms of antimicrobial resistance and for surveillance of drug-resistant bacterial pathogens. The majority of bacterial genomes sequenced to date have been sequenced with Illumina sequencing technology, owing to its high-throughput, excellent sequence accuracy, and low cost. However, because of the short-read nature of the technology, these assemblies are fragmented into large numbers of contigs, hindering the obtaining of full information of the genome. We develop Pasa, a graph-based algorithm that utilizes the pangenome graph and the assembly graph information to improve scaffolding quality. By leveraging the population information of the bacteria species, Pasa is able to utilize the linkage information of the gene families of the species to resolve the contig graph of the assembly. We show that our method outperforms the current state of the arts in terms of accuracy, and at the same time, is computationally efficient to be applied to a large number of existing draft assemblies.

Heteromultivalency enables enhanced detection of nucleic acid mutations

Detecting genetic mutations such as single nucleotide polymorphisms (SNPs) is necessary to prescribe effective cancer therapies, perform genetic analyses and distinguish similar viral strains. Traditionally, SNP sensing uses short oligonucleotide probes that differentially bind the SNP and wild-type targets. However, DNA hybridization-based techniques require precise tuning of the probe's binding affinity to manage the inherent trade-off between specificity and sensitivity. As conventional hybridization offers limited control over binding affinity, here we generate heteromultivalent DNA-functionalized particles and demonstrate optimized hybridization specificity for targets containing one or two mutations. By investigating the role of oligo lengths, spacer lengths and binding orientation, we reveal that heteromultivalent hybridization enables fine-tuned specificity for a single SNP and dramatic enhancements in specificity for two non-proximal SNPs empowered by highly cooperative binding. Capitalizing on these abilities, we demonstrate straightforward discrimination between heterozygous cis and trans mutations and between different strains of the SARS-CoV-2 virus. Our findings indicate that heteromultivalent hybridization offers substantial improvements over conventional monovalent hybridization-based methods.

Recent advances in the culture-independent discovery of natural products using metagenomic approaches

Natural products derived from bacterial sources have long been pivotal in the discovery of drug leads. However, the cultivation of only about 1% of bacteria in laboratory settings has left a significant portion of biosynthetic diversity hidden within the genomes of uncultured bacteria. Advances in sequencing technologies now enable the exploration of genetic material from these metagenomes through culture-independent methods. This approach involves extracting genetic sequences from environmental DNA and applying a hybrid methodology that combines functional screening, sequence tag-based homology screening, and bioinformatic-assisted chemical synthesis. Through this process, numerous valuable natural products have been identified and synthesized from previously uncharted metagenomic territories. This paper provides an overview of the recent advancements in the utilization of culture-independent techniques for the discovery of novel biosynthetic gene clusters and bioactive small molecules within metagenomic libraries.

Pharmacogenomics Beyond Single Common Genetic Variants: The Way Forward

Interindividual variability in genes encoding drug-metabolizing enzymes, transporters, receptors, and human leukocyte antigens has a major impact on a patient's response to drugs with regard to efficacy and safety. Enabled by both technological and conceptual advances, the field of pharmacogenomics is developing rapidly. Major progress in omics profiling methods has enabled novel genotypic and phenotypic characterization of patients and biobanks. These developments are paralleled by advances in machine learning, which have allowed us to parse the immense wealth of data and establish novel genetic markers and polygenic models for drug selection and dosing. Pharmacogenomics has recently become more widespread in clinical practice to personalize treatment and to develop new drugs tailored to specific patient populations. In this review, we provide an overview of the latest developments in the field and discuss the way forward, including how to address the missing heritability, develop novel polygenic models, and further improve the clinical implementation of pharmacogenomics.

LRTK: a platform agnostic toolkit for linked-read analysis of both human genome and metagenome

BACKGROUND

Linked-read sequencing technologies generate high-base quality short reads that contain extrapolative information on long-range DNA connectedness. These advantages of linked-read technologies are well known and have been demonstrated in many human genomic and metagenomic studies. However, existing linked-read analysis pipelines (e.g., Long Ranger) were primarily developed to process sequencing data from the human genome and are not suited for analyzing metagenomic sequencing data. Moreover, linked-read analysis pipelines are typically limited to 1 specific sequencing platform.

FINDINGS

To address these limitations, we present the Linked-Read ToolKit (LRTK), a unified and versatile toolkit for platform agnostic processing of linked-read sequencing data from both human genome and metagenome. LRTK provides functions to perform linked-read simulation, barcode sequencing error correction, barcode-aware read alignment and metagenome assembly, reconstruction of long DNA fragments, taxonomic classification and quantification, and barcode-assisted genomic variant calling and phasing. LRTK has the ability to process multiple samples automatically and provides users with the option to generate reproducible reports during processing of raw sequencing data and at multiple checkpoints throughout downstream analysis. We applied LRTK on linked reads from simulation, mock community, and real datasets for both human genome and metagenome. We showcased LRTK's ability to generate comparative performance results from preceding benchmark studies and to report these results in publication-ready HTML document plots.

CONCLUSIONS

LRTK provides comprehensive and flexible modules along with an easy-to-use Python-based workflow for processing linked-read sequencing datasets, thereby filling the current gap in the field caused by platform-centric genome-specific linked-read data analysis tools.

Exploring the cobia (Rachycentron canadum) genome: unveiling putative male heterogametic regions and identification of sex-specific markers

BACKGROUND

Cobia (Rachycentron canadum) is the only member of the Rachycentridae family and exhibits considerable sexual dimorphism in growth rate. Sex determination in teleosts has been a long-standing basic biological question, and the molecular mechanisms of sex determination/differentiation in cobia are completely unknown.

RESULTS

Here, we reported 2 high-quality, chromosome-level annotated male and female cobia genomes with assembly sizes of 586.51 Mb (contig/scaffold N50: 86.0 kb/24.3 Mb) and 583.88 Mb (79.9 kb/22.5 Mb), respectively. Synteny inference among perciform genomes revealed that cobia and the remora Echeneis naucrates were sister groups. Further, whole-genome resequencing of 31 males and 60 females, genome-wide association study, and sequencing depth analysis identified 3 short male-specific regions within a 10.7-kb continuous genomic region on male chromosome 18, which hinted at an undifferentiated sex chromosome system with a putative XX/XY mode of sex determination in cobia. Importantly, the only 2 genes within/between the male-specific regions, epoxide hydrolase 1 (ephx1, renamed cephx1y) and transcription factor 24 (tcf24, renamed ctcf24y), showed testis-specific/biased gene expression, whereas their counterparts cephx1x and ctf24x, located in female chromosome 18, were similarly expressed in both sexes. In addition, male-specific PCR targeting the cephx1y gene revealed that this genomic feature is conserved in cobia populations from Panama, Brazil, Australia, and Japan.

CONCLUSION

The first comprehensive genomic survey presented here is a valuable resource for future studies on cobia population structure and dynamics, conservation, and evolutionary history. Furthermore, it establishes evidence of putative male heterogametic regions with 2 genes playing a potential role in the sex determination of the species, and it provides further support for the rapid evolution of sex-determining mechanisms in teleost fish.

ABCA7-dependent Neuropeptide-Y signalling is a resilience mechanism required for synaptic integrity in Alzheimer's disease

Alzheimer's disease (AD) remains a complex challenge characterized by cognitive decline and memory loss. Genetic variations have emerged as crucial players in the etiology of AD, enabling hope for a better understanding of the disease mechanisms; yet the specific mechanism of action for those genetic variants remain uncertain. Animal models with reminiscent disease pathology could uncover previously uncharacterized roles of these genes. Using CRISPR/Cas9 gene editing, we generated a knockout model for abca7, orthologous to human ABCA7 - an established AD-risk gene. The abca7 +/- zebrafish showed reduced astroglial proliferation, synaptic density, and microglial abundance in response to amyloid beta 42 (Aβ42). Single-cell transcriptomics revealed abca7 -dependent neuronal and glial cellular crosstalk through neuropeptide Y (NPY) signaling. The abca7 knockout reduced the expression of npy, bdnf and ngfra , which are required for synaptic integrity and astroglial proliferation. With clinical data in humans, we showed reduced NPY in AD correlates with elevated Braak stage, predicted regulatory interaction between NPY and BDNF , identified genetic variants in NPY associated with AD, found segregation of variants in ABCA7, BDNF and NGFR in AD families, and discovered epigenetic changes in the promoter regions of NPY, NGFR and BDNF in humans with specific single nucleotide polymorphisms in ABCA7 . These results suggest that ABCA7-dependent NPY signaling is required for synaptic integrity, the impairment of which generates a risk factor for AD through compromised brain resilience.

Abstract Figure

Comparative Evolutionary Genomics in Insects

Genome sequencing quality, in terms of both read length and accuracy, is constantly improving. By combining long-read sequencing technologies with various scaffolding techniques, chromosome-level genome assemblies are now achievable at an affordable price for non-model organisms. Insects represent an exciting taxon for studying the genomic underpinnings of evolutionary innovations, due to ancient origins, immense species-richness, and broad phenotypic diversity. Here we summarize some of the most important methods for carrying out a comparative genomics study on insects. We describe available tools and offer concrete tips on all stages of such an endeavor from DNA extraction through genome sequencing, annotation, and several evolutionary analyses. Along the way we describe important insect-specific aspects, such as DNA extraction difficulties or gene families that are particularly difficult to annotate, and offer solutions. We describe results from several examples of comparative genomics analyses on insects to illustrate the fascinating questions that can now be addressed in this new age of genomics research.

BLR: a flexible pipeline for haplotype analysis of multiple linked-read technologies

Linked-read sequencing promises a one-method approach for genome-wide insights including single nucleotide variants (SNVs), structural variants, and haplotyping. We introduce Barcode Linked Reads (BLR), an open-source haplotyping pipeline capable of handling millions of barcodes and data from multiple linked-read technologies including DBS, 10× Genomics, TELL-seq and stLFR. Running BLR on DBS linked-reads yielded megabase-scale phasing with low (<0.2%) switch error rates. Of 13616 protein-coding genes phased in the GIAB benchmark set (v4.2.1), 98.6% matched the BLR phasing. In addition, large structural variants showed concordance with HPRC-HG002 reference assembly calls. Compared to diploid assembly with PacBio HiFi reads, BLR phasing was more continuous when considering switch errors. We further show that integrating long reads at low coverage (∼10×) can improve phasing contiguity and reduce switch errors in tandem repeats. When compared to Long Ranger on 10× Genomics data, BLR showed an increase in phase block N50 with low switch-error rates. For TELL-Seq and stLFR linked reads, BLR generated longer or similar phase block lengths and low switch error rates compared to results presented in the original publications. In conclusion, BLR provides a flexible workflow for comprehensive haplotype analysis of linked reads from multiple platforms.

Genetic variants associated with hantavirus infection in a reservoir host are related to regulation of inflammation and immune surveillance

The immunogenetics of wildlife populations influence the epidemiology and evolutionary dynamic of the host-pathogen system. Profiling immune gene diversity present in wildlife may be especially important for those species that, while not at risk of disease or extinction themselves, are host to diseases that are a threat to humans, other wildlife, or livestock. Hantaviruses (genus: Orthohantavirus) are globally distributed zoonotic RNA viruses with pathogenic strains carried by a diverse group of rodent hosts. The marsh rice rat (Oryzomys palustris) is the reservoir host of Orthohantavirus bayoui, a hantavirus that causes fatal cases of hantavirus cardiopulmonary syndrome in humans. We performed a genome wide association study (GWAS) using the rice rat "immunome" (i.e., all exons related to the immune response) to identify genetic variants associated with infection status in wild-caught rice rats naturally infected with their endemic strain of hantavirus. First, we created an annotated reference genome using 10× Chromium Linked Reads sequencing technology. This reference genome was used to create custom baits which were then used to target enrich prepared rice rat libraries (n = 128) and isolate their immunomes prior to sequencing. Top SNPs in the association test were present in four genes (Socs5, Eprs, Mrc1, and Il1f8) which have not been previously implicated in hantavirus infections. However, these genes correspond with other loci or pathways with established importance in hantavirus susceptibility or infection tolerance in reservoir hosts: the JAK/STAT, MHC, and NFκB. These results serve as informative markers for future exploration and highlight the importance of immune pathways that repeatedly emerge across hantavirus systems. Our work aids in creating cross-species comparisons for better understanding mechanisms of genetic susceptibility and host-pathogen coevolution in hantavirus systems.

Neotelomeres and Telomere-Spanning Chromosomal Arm Fusions in Cancer Genomes Revealed by Long-Read Sequencing

Alterations in the structure and location of telomeres are key events in cancer genome evolution. However, previous genomic approaches, unable to span long telomeric repeat arrays, could not characterize the nature of these alterations. Here, we applied both long-read and short-read genome sequencing to assess telomere repeat-containing structures in cancers and cancer cell lines. Using long-read genome sequences that span telomeric repeat arrays, we defined four types of telomere repeat variations in cancer cells: neotelomeres where telomere addition heals chromosome breaks, chromosomal arm fusions spanning telomere repeats, fusions of neotelomeres, and peri-centromeric fusions with adjoined telomere and centromere repeats. Analysis of lung adenocarcinoma genome sequences identified somatic neotelomere and telomere-spanning fusion alterations. These results provide a framework for systematic study of telomeric repeat arrays in cancer genomes, that could serve as a model for understanding the somatic evolution of other repetitive genomic elements.

Simultaneous deep transcriptome and proteome profiling in a single mouse oocyte

Although single-cell multi-omics technologies are undergoing rapid development, simultaneous transcriptome and proteome analysis of a single-cell individual still faces great challenges. Here, we developed a single-cell simultaneous transcriptome and proteome (scSTAP) analysis platform based on microfluidics, high-throughput sequencing, and mass spectrometry technology to achieve deep and joint quantitative analysis of transcriptome and proteome at the single-cell level, providing an important resource for understanding the relationship between transcription and translation in cells. This platform was applied to analyze single mouse oocytes at different meiotic maturation stages, reaching an average quantification depth of 19,948 genes and 2,663 protein groups in single mouse oocytes. In particular, we analyzed the correlation of individual RNA and protein pairs, as well as the meiosis regulatory network with unprecedented depth, and identified 30 transcript-protein pairs as specific oocyte maturational signatures, which could be productive for exploring transcriptional and translational regulatory features during oocyte meiosis.

IL-6 translation is a therapeutic target of human cytokine release syndrome

Chimeric antigen receptor (CAR) T therapies have achieved remarkable success for treating hematologic malignancies, yet are often accompanied by severe cytokine release syndrome (CRS). Here, an accidental clinical observation raised the possibility that metoprolol, an FDA-approved β1 adrenergic receptor blocker widely used for cardiovascular conditions, may alleviate CAR T-induced CRS. Metoprolol effectively blocked IL-6 production in human monocytes through unexpected mechanisms of action of targeting IL-6 protein translation but not IL6 mRNA expression. Mechanistically, metoprolol diminished IL-6 protein synthesis via attenuating eEF2K-eEF2 axis-regulated translation elongation. Furthermore, an investigator-initiated phase I/II clinical trial demonstrated a favorable safety profile of metoprolol in CRS management and showed that metoprolol significantly alleviated CAR T-induced CRS without compromising CAR T efficacy. These results repurposed metoprolol, a WHO essential drug, as a potential therapeutic for CRS and implicated IL-6 translation as a mechanistic target of metoprolol, opening venues for protein translation-oriented drug developments for human inflammatory diseases.

The Impact of Chromosomal Rearrangements in Speciation: From Micro- to Macroevolution

Chromosomal rearrangements (CRs) have been known since almost the beginning of genetics. While an important role for CRs in speciation has been suggested, evidence primarily stems from theoretical and empirical studies focusing on the microevolutionary level (i.e., on taxon pairs where speciation is often incomplete). Although the role of CRs in eukaryotic speciation at a macroevolutionary level has been supported by associations between species diversity and rates of evolution of CRs across phylogenies, these findings are limited to a restricted range of CRs and taxa. Now that more broadly applicable and precise CR detection approaches have become available, we address the challenges in filling some of the conceptual and empirical gaps between micro- and macroevolutionary studies on the role of CRs in speciation. We synthesize what is known about the macroevolutionary impact of CRs and suggest new research avenues to overcome the pitfalls of previous studies to gain a more comprehensive understanding of the evolutionary significance of CRs in speciation across the tree of life.

Disfunction of communication among immune cells in minimal-deviation adenocarcinoma of the cervix as an immunotherapeutic opportunity

Minimal deviation adenocarcinoma (MDA) of the uterine cervix, also referred to as malignant adenoma, is a rare subtype of cervical adenocarcinoma that exhibits histological characteristics resembling those of benign tumors, resulting in a low diagnostic rate and a lack of effective treatment options. The transcriptomic features of MDA at the single-cell resolution and within the tumor microenvironment (TME) remain unclear. In this study, we conducted single-cell transcriptomic analyses of MDA samples (Ca) and adjacent normal tissues (PCa). The present study reveals the prevalence of dendritic cells (DCs) and T cells in the carcinoma (Ca) of mammary ductal adenocarcinoma (MDA), with DCs undergoing significant metabolic reprogramming and immune stress. Additionally, our findings demonstrate the crucial involvement of DCs and T cells in the pathogenesis and metastatic progression of MDA, as evidenced by single-cell transcriptomic profiling of MDA and HPV samples. This resource provides a more profound understanding of the indolent nature of MDA and may prove useful in the development of MDA immunotherapy.

A genomic panel for studying C3-C4 intermediate photosynthesis in the Brassiceae tribe

Research on C4 and C3-C4 photosynthesis has attracted significant attention because the understanding of the genetic underpinnings of these traits will support the introduction of its characteristics into commercially relevant crop species. We used a panel of 19 taxa of 18 Brassiceae species with different photosynthesis characteristics (C3 and C3-C4) with the following objectives: (i) create draft genome assemblies and annotations, (ii) quantify orthology levels using synteny maps between all pairs of taxa, (iii) ⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠⁠describe the phylogenetic relatedness across all the species, and (iv) track the evolution of C3-C4 intermediate photosynthesis in the Brassiceae tribe. Our results indicate that the draft de novo genome assemblies are of high quality and cover at least 90% of the gene space. Therewith we more than doubled the sampling depth of genomes of the Brassiceae tribe that comprises commercially important as well as biologically interesting species. The gene annotation generated high-quality gene models, and for most genes extensive upstream sequences are available for all taxa, yielding potential to explore variants in regulatory sequences. The genome-based phylogenetic tree of the Brassiceae contained two main clades and indicated that the C3-C4 intermediate photosynthesis has evolved five times independently. Furthermore, our study provides the first genomic support of the hypothesis that Diplotaxis muralis is a natural hybrid of D. tenuifolia and D. viminea. Altogether, the de novo genome assemblies and the annotations reported in this study are a valuable resource for research on the evolution of C3-C4 intermediate photosynthesis.

MGA-seq: robust identification of extrachromosomal DNA and genetic variants using multiple genetic abnormality sequencing

Genomic abnormalities are strongly associated with cancer and infertility. In this study, we develop a simple and efficient method - multiple genetic abnormality sequencing (MGA-Seq) - to simultaneously detect structural variation, copy number variation, single-nucleotide polymorphism, homogeneously staining regions, and extrachromosomal DNA (ecDNA) from a single tube. MGA-Seq directly sequences proximity-ligated genomic fragments, yielding a dataset with concurrent genome three-dimensional and whole-genome sequencing information, enabling approximate localization of genomic structural variations and facilitating breakpoint identification. Additionally, by utilizing MGA-Seq, we map focal amplification and oncogene coamplification, thus facilitating the exploration of ecDNA's transcriptional regulatory function.

Genomic Variation, Population History, and Long-Term Genetic Adaptation to High Altitudes in Tibetan Partridge (Perdix hodgsoniae)

Species residing across elevational gradients display adaptations in response to environmental changes such as oxygen availability, ultraviolet radiation, and temperature. Here, we study genomic variation, gene expression, and long-term adaptation in Tibetan Partridge (Perdix hodgsoniae) populations residing across the elevational gradient of the Tibetan Plateau. We generated a high-quality draft genome and used it to carry out downstream population genomic and transcriptomic analysis. The P. hodgsoniae populations residing across various elevations were genetically distinct, and their phylogenetic clustering was consistent with their geographic distribution. We identified possible evidence of gene flow between populations residing in <3,000 and >4,200 m elevation that is consistent with known habitat expansion of high-altitude populations of P. hodgsoniae to a lower elevation. We identified a 60 kb haplotype encompassing the Estrogen Receptor 1 (ESR1) gene, showing strong genetic divergence between populations of P. hodgsoniae. We identified six single nucleotide polymorphisms within the ESR1 gene fixed for derived alleles in high-altitude populations that are strongly conserved across vertebrates. We also compared blood transcriptome profiles and identified differentially expressed genes (such as GAPDH, LDHA, and ALDOC) that correlated with differences in altitude among populations of P. hodgsoniae. These candidate genes from population genomics and transcriptomics analysis were enriched for neutrophil degranulation and glycolysis pathways, which are known to respond to hypoxia and hence may contribute to long-term adaptation to high altitudes in P. hodgsoniae. Our results highlight Tibetan Partridges as a useful model to study molecular mechanisms underlying long-term adaptation to high altitudes.

SpecHLA enables full-resolution HLA typing from sequencing data

Reconstructing diploid sequences of human leukocyte antigen (HLA) genes, i.e., full-resolution HLA typing, from sequencing data is challenging. The high homogeneity across HLA genes and the high heterogeneity within HLA alleles complicate the identification of genomic source loci for sequencing reads. Here, we present SpecHLA, which utilizes fine-tuned reads binning and local assembly to achieve accurate full-resolution HLA typing. SpecHLA accepts sequencing data from paired-end, 10×-linked-reads, high-throughput chromosome conformation capture (Hi-C), Pacific Biosciences (PacBio), and Oxford Nanopore Technology (ONT). It can also incorporate pedigree data and genotype frequency to refine typing. In 32 Human Genome Structural Variation Consortium, Phase 2 (HGSVC2) samples, SpecHLA achieved 98.6% accuracy for G-group-resolution HLA typing, inferring entire HLA alleles with an average of three mismatches fewer, ten gaps fewer, and 590 bp less edit distance than HISAT-genotype per allele. Additionally, SpecHLA exhibited a 2-field typing accuracy of 98.6% in 875 real samples. Finally, SpecHLA detected HLA loss of heterozygosity with 99.7% specificity and 96.8% sensitivity in simulated samples of cancer cell lines.

The effect of hyperthyroidism on cognitive function, neuroinflammation, and necroptosis in APP/PS1 mice

BACKGROUND

Increasing evidence has linked the thyroid dysfunction to the pathogenesis of dementia. Evidence from clinical studies has demonstrated that hypothyroidism is related to an increased risk of dementia. But the association of hyperthyroidism with dementia is largely unknown.

METHODS

We used the adenovirus containing thyrotropin receptor (TSHR) amino acid residues 1-289 (Ad-TSHR289)-induced Graves' disease (GD) phenotype in Alzheimer's disease (AD) model mice (APP/PS1 mice) to evaluate the effect of hyperthyroidism on the cognitive function and β-amyloid (Aβ) accumulation.

RESULTS

GD mice exhibited a stable long-term hyperthyroidism and cognitive deficits. Single Cell RNA-sequencing analysis indicated that microglia function played a critical role in the pathophysiological processes in GD mice. Neuroinflammation and polarization of microglia (M1/M2 phenotype) and activated receptor-interacting serine/threonine protein kinase 3 (RIPK3)/mixed lineage kinase domain-like pseudo-kinase (MLKL)-mediated necroptosis contributed to the pathological process, including Aβ deposition and neuronal loss. RIPK3 inhibitor could inhibit GD-mediated Aβ accumulation and neuronal loss.

CONCLUSIONS

Our findings reveal that GD hyperthyroidism aggravates cognitive deficits in AD mice and induces Aβ deposition and neuronal loss by inducing neuroinflammation and RIPK3/MLKL-mediated necroptosis.

Evolutionary dynamics of genome size and content during the adaptive radiation of Heliconiini butterflies

Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve.

Deciphering complex breakage-fusion-bridge genome rearrangements with Ambigram

Breakage-fusion-bridge (BFB) is a complex rearrangement that leads to tumor malignancy. Existing models for detecting BFBs rely on the ideal BFB hypothesis, ruling out the possibility of BFBs entangled with other structural variations, that is, complex BFBs. We propose an algorithm Ambigram to identify complex BFB and reconstruct the rearranged structure of the local genome during the cancer subclone evolution process. Ambigram handles data from short, linked, long, and single-cell sequences, and optical mapping technologies. Ambigram successfully deciphers the gold- or silver-standard complex BFBs against the state-of-the-art in multiple cancers. Ambigram dissects the intratumor heterogeneity of complex BFB events with single-cell reads from melanoma and gastric cancer. Furthermore, applying Ambigram to liver and cervical cancer data suggests that the BFB mechanism may mediate oncovirus integrations. BFB also exists in noncancer genomics. Investigating the complete human genome reference with Ambigram suggests that the BFB mechanism may be involved in two genome reorganizations of Homo Sapiens during evolution. Moreover, Ambigram discovers the signals of recurrent foldback inversions and complex BFBs in whole genome data from the 1000 genome project, and congenital heart diseases, respectively.

Microfluidic delivery of cutting enzymes for fragmentation of surface-adsorbed DNA molecules

We describe a method for fragmenting, in-situ, surface-adsorbed and immobilized DNAs on polymethylmethacrylate(PMMA)-coated silicon substrates using microfluidic delivery of the cutting enzyme DNase I. Soft lithography is used to produce silicone elastomer (Sylgard 184) gratings which form microfluidic channels for delivery of the enzyme. Bovine serum albumin (BSA) is used to reduce DNase I adsorption to the walls of the microchannels and enable diffusion of the cutting enzyme to a distance of 10mm. Due to the DNAs being immobilized, the fragment order is maintained on the surface. Possible methods of preserving the order for application to sequencing are discussed.

Genomics in the long-read sequencing era

Long-read sequencing (LRS) technologies have provided extremely powerful tools to explore genomes. While in the early years these methods suffered technical limitations, they have recently made significant progress in terms of read length, throughput, and accuracy and bioinformatics tools have strongly improved. Here, we aim to review the current status of LRS technologies, the development of novel methods, and the impact on genomics research. We will explore the most impactful recent findings made possible by these technologies focusing on high-resolution sequencing of genomes and transcriptomes and the direct detection of DNA and RNA modifications. We will also discuss how LRS methods promise a more comprehensive understanding of human genetic variation, transcriptomics, and epigenetics for the coming years.

Co-evolution of AR gene copy number and structural complexity in endocrine therapy resistant prostate cancer

Androgen receptor (AR) inhibition is standard of care for advanced prostate cancer (PC). However, efficacy is limited by progression to castration-resistant PC (CRPC), usually due to AR re-activation via mechanisms that include AR amplification and structural rearrangement. These two classes of AR alterations often co-occur in CRPC tumors, but it is unclear whether this reflects intercellular or intracellular heterogeneity of AR. Resolving this is important for developing new therapies and predictive biomarkers. Here, we analyzed 41 CRPC tumors and 6 patient-derived xenografts (PDXs) using linked-read DNA-sequencing, and identified 7 tumors that developed complex, multiply-rearranged AR gene structures in conjunction with very high AR copy number. Analysis of PDX models by optical genome mapping and fluorescence in situ hybridization showed that AR residing on extrachromosomal DNA (ecDNA) was an underlying mechanism, and was associated with elevated levels and diversity of AR expression. This study identifies co-evolution of AR gene copy number and structural complexity via ecDNA as a mechanism associated with endocrine therapy resistance.

Ariadne: synthetic long read deconvolution using assembly graphs

Synthetic long read sequencing techniques such as UST's TELL-Seq and Loop Genomics' LoopSeq combine 3[Formula: see text] barcoding with standard short-read sequencing to expand the range of linkage resolution from hundreds to tens of thousands of base-pairs. However, the lack of a 1:1 correspondence between a long fragment and a 3[Formula: see text] unique molecular identifier confounds the assignment of linkage between short reads. We introduce Ariadne, a novel assembly graph-based synthetic long read deconvolution algorithm, that can be used to extract single-species read-clouds from synthetic long read datasets to improve the taxonomic classification and de novo assembly of complex populations, such as metagenomes.

Rapid gene content turnover on the germline-restricted chromosome in songbirds

The germline-restricted chromosome (GRC) of songbirds represents a taxonomically widespread example of programmed DNA elimination. Despite its apparent indispensability, we still know very little about the GRC's genetic composition, function, and evolutionary significance. Here we assemble the GRC in two closely related species, the common and thrush nightingale. In total we identify 192 genes across the two GRCs, with many of them present in multiple copies. Interestingly, the GRC appears to be under little selective pressure, with the genetic content differing dramatically between the two species and many GRC genes appearing to be pseudogenized fragments. Only one gene, cpeb1, has a complete coding region in all examined individuals of the two species and shows no copy number variation. The acquisition of this gene by the GRC corresponds with the earliest estimates of the GRC origin, making it a good candidate for the functional indispensability of the GRC in songbirds.

MTG-Link: leveraging barcode information from linked-reads to assemble specific loci

BACKGROUND

Local assembly with short and long reads has proven to be very useful in many applications: reconstruction of the sequence of a locus of interest, gap-filling in draft assemblies, as well as alternative allele reconstruction of large Structural Variants. Whereas linked-read technologies have a great potential to assemble specific loci as they provide long-range information while maintaining the power and accuracy of short-read sequencing, there is a lack of local assembly tools for linked-read data.

RESULTS

We present MTG-Link, a novel local assembly tool dedicated to linked-reads. The originality of the method lies in its read subsampling step which takes advantage of the barcode information contained in linked-reads mapped in flanking regions. We validated our approach on several datasets from different linked-read technologies. We show that MTG-Link is able to assemble successfully large sequences, up to dozens of Kb. We also demonstrate that the read subsampling step of MTG-Link considerably improves the local assembly of specific loci compared to other existing short-read local assembly tools. Furthermore, MTG-Link was able to fully characterize large insertion variants and deletion breakpoints in a human genome and to reconstruct dark regions in clinically-relevant human genes. It also improved the contiguity of a 1.3 Mb locus of biological interest in several individual genomes of the mimetic butterfly Heliconius numata.

CONCLUSIONS

MTG-Link is an efficient local assembly tool designed for different linked-read sequencing technologies. MTG-Link source code is available at https://github.com/anne-gcd/MTG-Link and as a Bioconda package.

Hybridisation has shaped a recent radiation of grass-feeding aphids

BACKGROUND

Aphids are common crop pests. These insects reproduce by facultative parthenogenesis involving several rounds of clonal reproduction interspersed with an occasional sexual cycle. Furthermore, clonal aphids give birth to live young that are already pregnant. These qualities enable rapid population growth and have facilitated the colonisation of crops globally. In several cases, so-called "super clones" have come to dominate agricultural systems. However, the extent to which the sexual stage of the aphid life cycle has shaped global pest populations has remained unclear, as have the origins of successful lineages. Here, we used chromosome-scale genome assemblies to disentangle the evolution of two global pests of cereals-the English (Sitobion avenae) and Indian (Sitobion miscanthi) grain aphids.

RESULTS

Genome-wide divergence between S. avenae and S. miscanthi is low. Moreover, comparison of haplotype-resolved assemblies revealed that the S. miscanthi isolate used for genome sequencing is likely a hybrid, with one of its diploid genome copies closely related to S. avenae (~ 0.5% divergence) and the other substantially more divergent (> 1%). Population genomics analyses of UK and China grain aphids showed that S. avenae and S. miscanthi are part of a cryptic species complex with many highly differentiated lineages that predate the origins of agriculture. The complex consists of hybrid lineages that display a tangled history of hybridisation and genetic introgression.

CONCLUSIONS

Our analyses reveal that hybridisation has substantially contributed to grain aphid diversity, and hence, to the evolutionary potential of this important pest species. Furthermore, we propose that aphids are particularly well placed to exploit hybridisation events via the rapid propagation of live-born "frozen hybrids" via asexual reproduction, increasing the likelihood of hybrid lineage formation.

Pairwise comparative analysis of six haplotype assembly methods based on users' experience

BACKGROUND

A haplotype is a set of DNA variants inherited together from one parent or chromosome. Haplotype information is useful for studying genetic variation and disease association. Haplotype assembly (HA) is a process of obtaining haplotypes using DNA sequencing data. Currently, there are many HA methods with their own strengths and weaknesses. This study focused on comparing six HA methods or algorithms: HapCUT2, MixSIH, PEATH, WhatsHap, SDhaP, and MAtCHap using two NA12878 datasets named hg19 and hg38. The 6 HA algorithms were run on chromosome 10 of these two datasets, each with 3 filtering levels based on sequencing depth (DP1, DP15, and DP30). Their outputs were then compared.

RESULT

Run time (CPU time) was compared to assess the efficiency of 6 HA methods. HapCUT2 was the fastest HA for 6 datasets, with run time consistently under 2 min. In addition, WhatsHap was relatively fast, and its run time was 21 min or less for all 6 datasets. The other 4 HA algorithms' run time varied across different datasets and coverage levels. To assess their accuracy, pairwise comparisons were conducted for each pair of the six packages by generating their disagreement rates for both haplotype blocks and Single Nucleotide Variants (SNVs). The authors also compared them using switch distance (error), i.e., the number of positions where two chromosomes of a certain phase must be switched to match with the known haplotype. HapCUT2, PEATH, MixSIH, and MAtCHap generated output files with similar numbers of blocks and SNVs, and they had relatively similar performance. WhatsHap generated a much larger number of SNVs in the hg19 DP1 output, which caused it to have high disagreement percentages with other methods. However, for the hg38 data, WhatsHap had similar performance as the other 4 algorithms, except SDhaP. The comparison analysis showed that SDhaP had a much larger disagreement rate when it was compared with the other algorithms in all 6 datasets.

CONCLUSION

The comparative analysis is important because each algorithm is different. The findings of this study provide a deeper understanding of the performance of currently available HA algorithms and useful input for other users.