Modernizing reference genome assemblies

Using a novel panel of drug-resistant triple-negative breast cancer cell lines to identify candidate therapeutic targets and biomarkers

Here, we introduce a novel set of triple-negative breast cancer (TNBC) cell lines consisting of MDA-MB-468, HCC38, and HCC1806 and their sublines adapted to cisplatin, doxorubicin, eribulin, paclitaxel, gemcitabine, or 5-fluorouracil. Whole exome sequencing combined with TCGA-derived patient data resulted in the identification of 682 biomarker candidates in a pan-cancer analysis. Thirty-five genes were considered the most promising candidates because they harbored resistance-associated variants in at least two resistant sublines, and their expression correlated with TNBC patient survival. Exome sequencing and response profiles to cytotoxic drugs and DNA damage response inhibitors identified revealed remarkably little overlap between the resistant sublines, suggesting that each resistance formation process follows a unique route. This reflects recent findings on cancer cell evolution in patients, supporting the relevance of drug-adapted cancer cell lines as preclinical models of acquired resistance. Moreover, all of the drug-resistant TNBC sublines remained sensitive or even displayed collateral sensitivity to a range of tested compounds. Cross-resistance levels were lowest for the CHK2 inhibitor CCT241533, the PLK1 inhibitor SBE13, and the RAD51 recombinase inhibitor B02, suggesting that CHK2, PLK1, and RAD51 are potential drug targets for therapy-refractory TNBC. In conclusion, we present novel preclinical models of acquired drug resistance in TNBC and the identification of novel candidate therapeutic targets and biomarkers for this disease.

SNP associations in the L-citrulline metabolic pathway and vascular aging in the Japanese population

Decreased nitric oxide (NO) production from the vascular endothelium is a major factor for vascular aging. Because vascular aging has few specific subjective symptoms, assessing the susceptibility to vascular aging is beneficial for its early detection and improvement. Therefore, this study evaluated the associations between single nucleotide polymorphisms (SNPs) on L-citrulline (Cit) pathways essential for NO production and the health characteristics involved in vascular aging using a candidate gene approach. Associations with a significance level included those between the KCNMB4 rs17108108 C allele and tendency to gain weight, the ADCY8 rs6470860 G allele and numbness of limbs, the NOS1 rs2271987 T allele and lower back pain, and the PDE9A rs2284972 G allele and body pain with negative mood states. A genome-wide association study was also conducted to analyze SNPs more extensively across genes related to Cit and NO metabolism, which revealed genome-wide significant associations between the PRMT6 rs12028323 C allele and mood disturbance. These significant associations could be explained by the change in downstream NO signaling, supporting the relationship between the investigated traits and vascular function. These traits can be manifested as subjective symptoms of vascular aging. Therefore, the identified SNPs could predict susceptibility to the subjective symptoms of vascular aging, which could lead to genotype-based personalized interventions of Cit for an efficient improvement of vascular health.

Navigating complexities of polymorphic microbiomes in endometrial cancer

The microbiome is key to understanding endometrial cancer (EC) etiology and prevention strategies, implicated in the regulation of estrogen in estrogen-driven cancers. Utilizing robust methodologies in the QIIME 2 platform, we examined 16S rRNA vaginal and rectal microbiome data from an EC cohort: 192 women with benign gynecologic conditions, endometrial hyperplasia, or endometrial cancer. Distinct microbial compositions and community networks specific to EC were identified and related to histological grade with adjustments for EC risk factors. Vaginal health-associated Lactobacillus and Limosilactobacillus, and rectal Prevotella and Peptoniphilus, were depleted in EC, while detrimental vaginal Anaerococcus, Porphyromonas, Prevotella, Peptoniphilus, and rectal Buttiaxella were enriched. Significant bacterial features were shared between rectal and vaginal sites in EC, such as Prevotella timonensis and Peptoniphilus A. Vaginal Lactobacillus abundance contributed to less feature sharing from the rectum. Putative microbial metabolic analysis identified dysregulation of amino acid, complex carbohydrate, and hormone metabolism amongst patients with EC.

IGLoo enables comprehensive analysis and assembly of immunoglobulin heavy-chain loci in lymphoblastoid cell lines using PacBio high-fidelity reads

High-quality human genome assemblies derived from lymphoblastoid cell lines (LCLs) provide reference genomes and pangenomes for genomics studies. However, LCLs pose technical challenges for profiling immunoglobulin (IG) genes, as their IG loci contain a mixture of germline and somatically recombined haplotypes, making genotyping and assembly difficult with widely used frameworks. To address this, we introduce IGLoo, a software tool that analyzes sequence data and assemblies derived from LCLs, characterizing somatic V(D)J recombination events and identifying breakpoints and missing IG genes in the assemblies. Furthermore, IGLoo implements a reassembly framework to improve germline assembly quality by integrating information on somatic events and population structural variations in IG loci. Applying IGLoo to the assemblies from the Human Pangenome Reference Consortium, we gained valuable insights into the mechanisms, gene usage, and patterns of V(D)J recombination and the causes of assembly artifacts in the IG heavy-chain (IGH) locus, and we improved the representation of IGH assemblies.

Human mitochondrial DNA in public metagenomes: Opportunity or privacy threat?

Human DNA is unavoidably present in metagenomic analyses of human microbiomes. While current protocols remove human DNA before submission to public repositories, mitochondrial DNA (mtDNA) has been overlooked and frequently persists. We discuss the privacy risks and research opportunities associated with mtDNA, urging consideration by the scientific, ethics, and legal communities.

Combining chromosome conformation capture and exome sequencing for simultaneous detection of structural and single-nucleotide variants

BACKGROUND

Effective molecular diagnosis of congenital diseases hinges on comprehensive genomic analysis, traditionally reliant on various methodologies specific to each variant type-whole exome or genome sequencing for single nucleotide variants (SNVs), array CGH for copy-number variants (CNVs), and microscopy for structural variants (SVs).

METHODS

We introduce a novel, integrative approach combining exome sequencing with chromosome conformation capture, termed Exo-C. This method enables the concurrent identification of SNVs in clinically relevant genes and SVs across the genome and allows analysis of heterozygous and mosaic carriers. Enhanced with targeted long-read sequencing, Exo-C evolves into a cost-efficient solution capable of resolving complex SVs at base-pair accuracy.

RESULTS

Applied to 66 human samples Exo-C achieved 100% recall and 73% precision in detecting chromosomal translocations and SNVs. We further benchmarked its performance for inversions and CNVs and demonstrated its utility in detecting mosaic SVs and resolving diagnostically challenging cases.

CONCLUSIONS

Through several case studies, we demonstrate how Exo-C's multifaceted application can effectively uncover diverse causative variants and elucidate disease mechanisms in patients with rare disorders.

Co-variation of Host Gene Expression and Gut Microbiome in Intestine-Specific Spp1 Conditional Knockout Mice

Osteopontin, which is a highly phosphorylated and glycosylated acidic secreted protein encoded by the secreted phosphoprotein 1 (Spp1) gene, plays a crucial role in immune regulation, inflammatory responses, and cell adhesion. However, its impact on intestinal gene expression and gut microbiota remains underexplored. In this study, we developed an Spp1 conditional knockout mouse model to investigate alterations in the intestinal transcriptome and microbiome, with particular emphasis on changes in gene expression and predicted metabolic pathways. Our findings demonstrated that Spp1 gene conditional knockout significantly modified the expression of genes involved in immune regulation and lipid metabolism. Moreover, metagenomic analysis revealed marked shifts in gut microbial diversity and predicted the metabolic pathways associated with digestion, absorption, and lipid metabolism. These results suggest that Spp1 is instrumental in maintaining gut microbial equilibrium and in regulating host lipid metabolism and immune responses. This study offers new insights into the role of Spp1 in host-microbiota interactions and the potential foundations for developing related therapeutic strategies.

Rare predicted loss-of-function and damaging missense variants in CFHR5 associate with protection from age-related macular degeneration

Age-related macular degeneration (AMD) is a leading cause of blindness among older adults worldwide, but treatment options are limited. Genetics studies have implicated the CFH locus, containing CFH and five CFHR genes, CFHR1-5, in AMD. While CFH has been robustly linked with AMD risk, potential additional roles for the CFHR genes remain unclear, obscured by strong linkage disequilibrium across the locus. Investigating rare coding variants can help to identify causal genes in such regions. We used whole-exome sequencing data from 406,952 UK Biobank participants to examine AMD associations with genes at the CFH locus. For each gene, we used burden testing to examine associations of rare (minor-allele frequency [MAF] < 1%) predicted loss-of-function (pLoF) and predicted damaging missense variants with AMD. We considered "broadly defined AMD" (ICD-10 35.3; ncases = 10,700) and "strictly defined AMD" (dry or wet AMD; ncases = 346). Adjusting for CFH-region variants known to independently associate with AMD, we find that CFHR5 rare variant burden significantly associates with a decreased risk of broadly defined AMD (odds ratio [OR] = 0.75, p = 7 × 10-4), with this association primarily driven by pLoF variants. Furthermore, the association of CFHR5 rare variants with AMD protection is estimated to be stronger for individuals with the CFH rs1061170 AMD risk allele (p.Tyr402His [p.Y402H]; interaction p = 0.04). Corresponding analyses of strict AMD were underpowered. However, we observe that thinning of the photoreceptor layer outer segment strongly predicts strict AMD and find that CFHR5 rare variant burden is significantly associated with increased thickness of this retinal layer (+0.34 SD, p = 4 × 10-4, n = 45,365). These findings suggest CFHR5 inhibition as a potential therapeutic approach for AMD.

The genetic architecture of and evolutionary constraints on the human pelvic form

Human pelvic evolution following the human-chimpanzee divergence is thought to result in an obstetrical dilemma, a mismatch between large infant brains and narrowed female birth canals, but empirical evidence has been equivocal. By using deep learning on 31,115 dual-energy x-ray absorptiometry scans from UK Biobank, we identified 180 loci associated with seven highly heritable pelvic phenotypes. Birth canal phenotypes showed sex-specific genetic architecture, aligning with reproductive function. Larger birth canals were linked to slower walking pace and reduced back pain but increased hip osteoarthritis risk, whereas narrower birth canals were associated with reduced pelvic floor disorder risk but increased obstructed labor risk. Lastly, genetic correlation between birth canal and head widths provides evidence of coevolution between the human pelvis and brain, partially mitigating the dilemma.

Contribution of Rare and Potentially Functionally Relevant Sequence Variants in Schizophrenia Risk-Locus Xq28,distal

Duplications of the Xq28,distal locus have been described in male and female patients with schizophrenia (SCZ) or intellectual disability. The Xq28,distal locus spans eight protein-coding genes (F8, CMC4, MTCP1, BRCC3, VBP1, FUNDC2, CLIC2, and RAB39B) and is flanked by recurrent genomic breakpoints. Thus, the issue of which gene/s at this locus is/are relevant in terms of SCZ pathogenesis remains unclear. The aim of this study was to investigate the contribution of rare and potentially functionally relevant sequence variants within the Xq28,distal locus to SCZ risk using the single-molecule molecular inversion probes (smMIP) method. Targeted sequencing was performed in a cohort of 1935 patients with SCZ and 1905 controls of European ancestry. The consecutive statistical analysis addressed two main areas. On the level of the individual variants, allele counts in the patient and control cohort were systematically compared with a Fisher's exact test: (i) for the entire present study cohort; (ii) for patients and controls separated by sex; and (iii) in combination with data published by the Schizophrenia Exome Meta-Analysis (SCHEMA) consortium. On the gene-wise level, a burden analysis was performed using the X-chromosomal model of the Optimal Unified Sequence Kernel Association Test (SKAT-O), with adjustment for possible sex-specific effects. Targeted sequencing identified a total of 13 rare and potentially functional variants in four patients and 11 controls. However, neither at the level of individual rare and potentially functional variants nor at the level of the eight protein-coding genes at the Xq28,distal locus was a statistically significant enrichment in patients compared to controls observed. Although inconclusive, the present findings represent a step toward improved understanding of the contribution of X-chromosomal risk factors in neuropsychiatric disorder development, which is an underrepresented aspect of genetic studies in this field.

CUT&Tag recovers up to half of ENCODE ChIP-seq histone acetylation peaks

DNA-protein interactions have traditionally been profiled via chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq). Cleavage Under Targets & Tagmentation (CUT&Tag) is a rapidly expanding technique that enables the profiling of such interactions in situ at high sensitivity. However, thorough evaluation and benchmarking against established ChIP-seq datasets are lacking. Here, we comprehensively benchmarked CUT&Tag for H3K27ac and H3K27me3 against published ChIP-seq profiles from ENCODE in K562 cells. Combining multiple new and published CUT&Tag datasets, there was an average recall of 54% known ENCODE peaks for both histone modifications. We tested peak callers MACS2 and SEACR and identified optimal peak calling parameters. Overall, peaks identified by CUT&Tag represent the strongest ENCODE peaks and show the same functional and biological enrichments as ChIP-seq peaks identified by ENCODE. Our workflow systematically evaluates the merits of methodological adjustments, providing a benchmarking framework for the experimental design and analysis of CUT&Tag studies.

Inferring active mutational processes in cancer using single cell sequencing and evolutionary constraints

Ongoing mutagenesis in cancer drives genetic diversity throughout the natural history of cancers. As the activities of mutational processes are dynamic throughout evolution, distinguishing the mutational signatures of 'active' and 'historical' processes has important implications for studying how tumors evolve. This can aid in understanding mutagenic states at the time of presentation, and in associating active mutational process with therapeutic resistance. As bulk sequencing primarily captures historical mutational processes, we studied whether ultra-low-coverage single-cell whole-genome sequencing (scWGS), which measures the distribution of mutations across hundreds or thousands of individual cells, could enable the distinction between historical and active mutational processes. While technical challenges and data sparsity have limited mutation analysis in scWGS, we show that these data contain valuable information about dynamic mutational processes. To robustly interpret single nucleotide variants (SNVs) in scWGS, we introduce ArtiCull, a method to identify and remove SNV artifacts by leveraging evolutionary constraints, enabling reliable detection of mutations for signature analysis. Applying this approach to scWGS data from pancreatic ductal adenocarcinoma (PDAC), triple-negative breast cancer (TNBC), and high-grade serous ovarian cancer (HGSOC), we uncover temporal and spatial patterns in mutational processes. In PDAC, we observe a temporal increase in mismatch repair deficiency (MMRd). In cisplatin-treated TNBC patient-derived xenografts, we identify therapy-induced mutagenesis and inactivation of APOBEC3 activity. In HGSOC, we show distinct patterns of APOBEC3 mutagenesis, including late tumor-wide activation in one case and clade-specific enrichment in another. Additionally, we detect a clone-specific increase in SBS17 activity, in a clone previously linked to recurrence. Our findings establish ultra-low-coverage scWGS as a powerful approach for studying active mutational processes that may influence ongoing clonal evolution and therapeutic resistance.

Facilitating high throughput collections-based genomics: a comparison of DNA extraction and library building methods

While DNA barcoding methods are an increasingly important tool in biological conservation, the resource requirements of constructing reference libraries frequently reduce their efficacy. One efficient way of sourcing taxonomically validated DNA for reference libraries is to use museum collections. However, DNA degradation intrinsic to historical museum specimens can, if not addressed in the wet lab, lead to low quality data generation and severely limit scientific output. Several DNA extraction and library build methods that are designed to work with degraded DNA have been developed, although the ability to implement these methods at scale and at low cost has yet to be formally addressed. Here, the performance of widely used DNA extraction and library build methods are compared using museum specimens. We find that while our selected DNA extraction methods do not significantly differ in DNA yield, the Santa Cruz Reaction (SCR) library build method is not only the most effective at retrieving degraded DNA from museum specimens but also easily implemented at high throughput for low cost. Results highlight the importance of lab protocol on data yield. An optimised "sample to sequencing" high-throughput protocol which incorporates SCR is included to allow for easy uptake by the wider scientific community.

Combining DNA and protein alignments to improve genome annotation with LiftOn

As the number and variety of assembled genomes continue to grow, the number of annotated genomes is falling behind, particularly for eukaryotes. DNA-based mapping tools help to address this challenge, but they are only able to transfer annotation between closely related species. Here we introduce LiftOn, a homology-based software tool that integrates DNA and protein alignments to enhance the accuracy of genome-scale annotation and to allow mapping between relatively distant species. LiftOn's protein-centric algorithm considers both types of alignments, chooses optimal open reading frames, resolves overlapping gene loci, and finds additional gene copies when they exist. LiftOn can reliably transfer annotation between genomes representing members of the same species, as we demonstrate on human, mouse, honeybee, rice, and Arabidopsis thaliana It can further map annotation effectively across species pairs as far apart as mouse and rat or Drosophila melanogaster and Drosophila erecta.

Single, but not mixed dietary fibers suppress body weight gain and adiposity in high fat-fed mice

Dietary fiber can suppress excess adipose tissue and weight gain in rodents and humans when fed high fat diets. The gut microbiome is thought to have a key role, although exactly how remains unclear. In a tightly controlled murine study, we explored how different types of dietary fiber and doses affect the gut microbiota and gut epithelial gene expression. We show that 10% pectin and 10% FOS suppress high fat diet (HFD)-induced weight gain, effects not seen at 2% doses. Furthermore, 2 and 10% mixtures of dietary fiber were also without effect. Each fiber treatment stimulated a distinct gut microbiota profile at the family and operational taxonomic unit (OTU) level. Mechanistically it is likely that the single 10% fiber dose shifted selected bacteria above some threshold abundance, required to suppress body weight, which was not achieved by the 10% Mix, composed of 4 fibers each at 2.5%. Plasma levels of the gut hormone PYY were elevated by 10% pectin and FOS, but not 10% mixed fibers, and similarly RNA seq revealed some distinct effects of the 10% single fibers on gut epithelial gene expression. These data show how the ability of dietary fiber to suppress HFD-induced weight gain is dependent upon both fiber type and dose. It also shows that the microbial response to dietary fiber is distinct and that there is not a single microbial response associated with the inhibition of adiposity and weight gain. PYY seems key to the latter response, although the role of other factors such as Reg3γ and CCK needs to be explored.

Genome-wide gene expression profiles throughout human malaria parasite liver stage development in humanized mice

Gene expression of Plasmodium falciparum (Pf) liver-stage (LS) parasites has remained poorly characterized, although they are major vaccine and drug targets. Using a human liver-chimaeric mouse model and a fluorescent parasite line (PfNF54CSPGFP), we isolated PfLS and performed transcriptomics on key LS developmental phases. We linked clustered gene expression to ApiAP2, a major family of transcription factors that regulate the parasite life cycle. This provided insights into transcriptional regulation of LS infection and expression of essential LS metabolic and biosynthetic pathways. We observed expression of antigenically variant PfEMP1 proteins and the major Pf protein export machine PTEX and identified protein candidates that might be exported by LS parasites. Comparing Pf and P. vivax LS transcriptomes, we uncovered differences in their expression of sexual commitment factors. This data will aid LS research and vaccine and drug target identification for prevention of malaria infection.

The maternal X chromosome affects cognition and brain ageing in female mice

Female mammalian cells have two X chromosomes, one of maternal origin and one of paternal origin. During development, one X chromosome randomly becomes inactivated1-4. This renders either the maternal X (Xm) chromosome or the paternal X (Xp) chromosome inactive, causing X mosaicism that varies between female individuals, with some showing considerable or complete skew of the X chromosome that remains active5-7. Parent-of-X origin can modify epigenetics through DNA methylation8,9 and possibly gene expression; thus, mosaicism could buffer dysregulated processes in ageing and disease. However, whether X skew or its mosaicism alters functions in female individuals is largely unknown. Here we tested whether skew towards an active Xm chromosome influences the brain and body-and then delineated unique features of Xm neurons and Xp neurons. An active Xm chromosome impaired cognition in female mice throughout the lifespan and led to worsened cognition with age. Cognitive deficits were accompanied by Xm-mediated acceleration of biological or epigenetic ageing of the hippocampus, a key centre for learning and memory, in female mice. Several genes were imprinted on the Xm chromosome of hippocampal neurons, suggesting silenced cognitive loci. CRISPR-mediated activation of Xm-imprinted genes improved cognition in ageing female mice. Thus, the Xm chromosome impaired cognition, accelerated brain ageing and silenced genes that contribute to cognition in ageing. Understanding how Xm impairs brain function could lead to an improved understanding of heterogeneity in cognitive health in female individuals and to X-chromosome-derived pathways that protect against cognitive deficits and brain ageing.

The shared genetic architecture and evolution of human language and musical rhythm

This study aimed to test theoretical predictions over biological underpinnings of previously documented phenotypic correlations between human language-related and musical rhythm traits. Here, after identifying significant genetic correlations between rhythm, dyslexia and various language-related traits, we adapted multivariate methods to capture genetic signals common to genome-wide association studies of rhythm (N = 606,825) and dyslexia (N = 1,138,870). The results revealed 16 pleiotropic loci (P < 5 × 10-8) jointly associated with rhythm impairment and dyslexia, and intricate shared genetic and neurobiological architectures. The joint genetic signal was enriched for foetal and adult brain cell-specific regulatory regions, highlighting complex cellular composition in their shared underpinnings. Local genetic correlation with a key white matter tract (the left superior longitudinal fasciculus-I) substantiated hypotheses about auditory-motor connectivity as a genetically influenced, evolutionarily relevant neural endophenotype common to rhythm and language processing. Overall, we provide empirical evidence of multiple aspects of shared biology linking language and musical rhythm, contributing novel insight into the evolutionary relationships between human musicality and linguistic communication traits.

Epigenetic trajectory predicts development of clinical rheumatoid arthritis in ACPA+ individuals: Targeting Immune Responses for Prevention of Rheumatoid Arthritis (TIP-RA)

OBJECTIVE

The presence of autoantibodies to citrullinated protein antigens (ACPAs) in the absence of clinically-apparent inflammatory arthritis (IA) identifies individuals at-risk for developing future clinical rheumatoid arthritis (RA). However, it is unclear why some ACPA+ individuals convert to clinical RA while others do not. We explored the possibility in the Targeting Immune Responses for Prevention of Rheumatoid Arthritis (TIP-RA) study that epigenetic remodeling is part of the trajectory from an at-risk state to clinical disease and identifies novel biomarkers associated with conversion to clinical RA.

METHODS

ACPA- Controls, ACPA+ At-Risk, and Early RA individuals were followed for up to 5 years, including obtaining blood samples annually and at RA diagnosis. Peripheral blood mononuclear cells (PBMCs) were separated into CD19+ B cells, memory CD4+ T cells, and naive CD4+ T cells using antibodies and magnetic beads. Genome-wide methylation within each cell lineage was assayed using the Illumina MethylationEPIC v1.0 beadchip. ACPA+ At-Risk participants who did or did not develop RA were designated Pre-RA or Non-converters, respectively. Differentially methylated loci (DML) were selected using the Limma software package. Using the Caret package, we constructed machine learning models in test and validation cohorts and identified the most predictive loci of clinical RA conversion.

RESULTS

Cross-sectional differential methylation analysis at baseline revealed DMLs that distinguish the Pre-RA methylome from ACPA+ Non-converters, the latter which closely resembled ACPA- Controls. Genes overlapping these DMLs correspond to aberrant NOTCH signaling and DNA repair pathways in B cells. Longitudinal analysis showed that ACPA- Control and ACPA+ Non-converter methylomes are relatively constant. In contrast, the Pre-RA methylome remodeled along a dynamic RA methylome trajectory characterized by epigenetic changes in active regulatory elements. Clinical conversion to RA, defined based on diagnosis, marked an epigenetic inflection point for cell cycle pathways in B cells and adaptive immunity pathways in naive T cells. Machine learning revealed individual loci associated with RA conversion. This model significantly outperformed autoantibodies plus acute phase reactants as predictors of RA conversion.

CONCLUSION

DNA methylation is a dynamic process in ACPA+ individuals at-risk for developing RA that eventually transition to clinical disease. In contrast, non-converters and controls have stable methylomes. The accumulation of epigenetic marks over time prior to conversion to clinical RA conforms to pathways that are associated with immunity and can be used to identify potential pathogenic pathways for therapeutic targeting and/or use as prognostic biomarkers.

Homo Sapiens Chromosomal Location Ontology: A Framework for Genomic Data in Biomedical Knowledge Graphs

The Homo sapiens Chromosomal Location Ontology (HSCLO) is designed to facilitate the integration of human genomic features into biomedical knowledge graphs from releases GRCh37 and GRCh38 at multiple resolutions. HSCLO comprises two distinct versions, HSCLO37 and HSCLO38, each tailored to its respective human genome release. This ontology supports the efficient integration and analysis of human genomic data across scales ranging from entire chromosomes to individual base pairs, thereby enhancing data retrieval and interoperability within large-scale biomedical datasets. Unlike existing ontologies that primarily focus on genomic feature identification or annotation, HSCLO is specifically engineered to optimize the interoperability and scalability of genomic data within biomedical knowledge graphs. The utility and performance of HSCLO are demonstrated through a case study involving the integration of high-resolution chromatin interaction data, which reveals significant improvements in query efficiency and data linkage. HSCLO represents a valuable resource for advancing research in disease genetics, personalized medicine, and other domains that require complex genomic data integration.

Highly accurate Korean draft genomes reveal structural variation highlighting human telomere evolution

Given the presence of highly repetitive genomic regions such as subtelomeric regions, understanding human genomic evolution remains challenging. Recently, long-read sequencing technology has facilitated the identification of complex genetic variants, including structural variants (SVs), at the single-nucleotide level. Here, we resolved SVs and their underlying DNA damage-repair mechanisms in subtelomeric regions, which are among the most uncharted genomic regions. We generated ∼20 × high-fidelity long-read sequencing data from three Korean individuals and their partially phased high-quality de novo genome assemblies (contig N50: 6.3-58.2 Mb). We identified 131 138 deletion and 121 461 insertion SVs, 41.6% of which were prevalent in the East Asian population. The commonality of the SVs identified among the Korean population was examined by short-read sequencing data from 103 Korean individuals, providing the first comprehensive SV set representing the population based on the long-read assemblies. Manual investigation of 19 large subtelomeric SVs (≥5 kb) and their associated repair signatures revealed the potential repair mechanisms leading to the formation of these SVs. Our study provides mechanistic insight into human telomere evolution and can facilitate our understanding of human SV formation.

VDGE: a data repository of variation database for gene-edited animals across multiple species

Gene-edited animals are crucial for addressing fundamental questions in biology and medicine and hold promise for practical applications. In light of the rapid advancement of gene editing technologies over the past decade, a dramatically increased number of gene-edited animals have been generated. Genome editing at off-target sites can, however, introduce genomic variations, potentially leading to unintended functional consequences in these animals. So, there is an urgent need to systematically collect and collate these variations in gene-edited animals to aid data mining and integrative in-depth analyses. However, existing databases are currently insufficient to meet this need. Here, we present the Variation Database of Gene-Edited animals (VDGE, https://ngdc.cncb.ac.cn/vdge), the first open-access repository to present genomic variations and annotations in gene-edited animals, with a particular focus on larger animals such as monkeys. At present, VDGE houses 151 on-target mutations from 210 samples, and 115,710 variations identified from 107 gene-edited and wild-type animal trios through unified and standardized analysis and concurrently provides comprehensive annotation details for each variation, thus facilitating the assessment of their functional consequences and promoting mechanistic studies and practical applications for gene-edited animals.

Next-generation sequencing analysis with a population-specific human reference genome

Next-generation sequencing (NGS) has become widely available and is routinely used in basic research and clinical practice. The reference genome sequence is an essential resource for NGS analysis, and several population-specific reference genomes have recently been constructed to provide a choice to deal with the vast genetic diversity of human samples. However, resources supporting population-specific references are insufficient, and it is burdensome to perform analysis using these reference genomes. Here, we constructed a set of resources to support NGS analysis using the Japanese reference genome, JG. We created resources for variant calling, variant effect prediction, gene and repeat element annotations, read mappability and RNA-seq analysis. We also provide a resource for reference coordinate conversion for further annotation enrichment. We then provide a variant calling protocol with JG. Our resources provide a guide to prepare sufficient resources for the use of population-specific reference genomes and can facilitate the migration of reference genomes.

Normalization of gene counts affects principal components-based exploratory analysis of RNA-sequencing data

Normalization of gene expression count data is an essential step of in the analysis of RNA-sequencing data. Its statistical analysis has been mostly addressed in the context of differential expression analysis, that is in the univariate setting. However, relationships among genes and samples are better explored and quantified using multivariate exploratory data analysis tools like Principal Component Analysis (PCA). In this study we investigate how normalization impacts the PCA model and its interpretation, considering twelve different widely used normalization methods that were applied on simulated and experimental data. Correlation patterns in the normalized data were explored using both summary statistics and Covariance Simultaneous Component Analysis. The impact of normalization on the PCA solution was assessed by exploring the model complexity, the quality of sample clustering in the low-dimensional PCA space and gene ranking in the model fit to normalized data. PCA models upon normalization were interpreted in the context gene enrichment pathway analysis. We found that although PCA score plots are often similar independently form the normalization used, biological interpretation of the models can depend heavily on the normalization method applied.

Characterization of a novel heterozygous variant in the histidyl-tRNA synthetase gene associated with Charcot-Marie-Tooth disease type 2W

Heterozygous pathogenic variants in the histidyl-tRNA synthetase (HARS) gene are associated with Charcot-Marie-Tooth (CMT) type 2W disease, classified as an axonal peripheral neuropathy. To date, at least 60 variants causing CMT symptoms have been identified in seven different aminoacyl-tRNA synthetases, with eight being found in the catalytic domain of HARS. The genetic data clearly show a causative role of aminoacyl-tRNA synthetases in CMT; however, the cellular mechanisms leading to pathology can vary widely and are unknown in the case of most identified variants. Here we describe a novel HARS variant, c.412T>C; p.Y138H, identified through a CMT gene panel in a patient with peripheral neuropathy. To determine the effect of p.Y138H we employed a humanized HARS yeast model and recombinant protein biochemistry, which identified a deficiency in protein dimerization and a growth defect which shows mild but significant improvement with histidine supplementation. This raises the potential for a clinical trial of histidine.

Activation of recombinases at specific DNA loci by zinc-finger domain insertions

Recombinases have several potential advantages as genome editing tools compared to nucleases and other editing enzymes, but the process of engineering them to efficiently recombine predetermined DNA targets demands considerable investment of time and labor. Here we sought to harness zinc-finger DNA-binding domains (ZFDs) to program recombinase binding by developing fusions, in which ZFDs are inserted into recombinase coding sequences. By screening libraries of hybrid proteins, we optimized the insertion site, linker length, spacing and ZFD orientation and generated Cre-type recombinases that remain dormant unless the insertionally fused ZFD binds its target site placed in the vicinity of the recombinase binding site. The developed fusion improved targeted editing efficiencies of recombinases by four-fold and abolished measurable off-target activity in mammalian cells. The ZFD-dependent activity is transferable to a recombinase with relaxed specificity, providing the means for developing fully programmable recombinases. Our engineered recombinases provide improved genome editing tools with increased precision and efficiency.

Generation and analysis of a mouse multitissue genome annotation atlas

Generating an accurate and complete genome annotation for an organism is complex because the cells within each tissue can express a unique set of transcript isoforms from a unique set of genes. A comprehensive genome annotation should contain information on what tissues express what transcript isoforms at what level. This tissue-level isoform information can then inform a wide range of research questions as well as experiment designs. Long-read sequencing technology combined with advanced full-length cDNA library preparation methods has now achieved throughput and accuracy where generating these types of annotations is achievable. Here, we show this by generating a genome annotation of the mouse (Mus musculus). We used the nanopore-based R2C2 long-read sequencing method to generate 64 million highly accurate full-length cDNA consensus reads-averaging 5.4 million reads per tissue for a dozen tissues. Using the Mandalorion tool, we processed these reads to generate the Tissue-level Atlas of Mouse Isoforms which is available as a trackhub for the UCSC Genome Browser and contains at least one full-length isoform for the vast majority of expressed genes in each tissue.

Leveraging the T2T assembly to resolve rare and pathogenic inversions in reference genome gaps

Chromosomal inversions (INVs) are particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage-sensitive genes in cis Short-read genome sequencing (srGS) can only resolve ∼70% of cytogenetically visible inversions referred to clinical diagnostic laboratories, likely due to breakpoints in repetitive regions. Here, we study 12 inversions by long-read genome sequencing (lrGS) (n = 9) or srGS (n = 3) and resolve nine of them. In four cases, the inversion breakpoint region was missing from at least one of the human reference genomes (GRCh37, GRCh38, T2T-CHM13) and a reference agnostic analysis was needed. One of these cases, an INV9 mappable only in de novo assembled lrGS data using T2T-CHM13 disrupts EHMT1 consistent with a Mendelian diagnosis (Kleefstra syndrome 1; MIM#610253). Next, by pairwise comparison between T2T-CHM13, GRCh37, and GRCh38, as well as the chimpanzee and bonobo, we show that hundreds of megabases of sequence are missing from at least one human reference, highlighting that primate genomes contribute to genomic diversity. Aligning population genomic data to these regions indicated that these regions are variable between individuals. Our analysis emphasizes that T2T-CHM13 is necessary to maximize the value of lrGS for optimal inversion detection in clinical diagnostics. These results highlight the importance of leveraging diverse and comprehensive reference genomes to resolve unsolved molecular cases in rare diseases.

Duchenne Muscular Dystrophy in Two Half-Brothers Due to Inherited 306 Kb Inverted Insertion of 10p15.1 into Intron 44 of the Dp427m Transcript of the DMD Gene

Duchenne muscular dystrophy (DMD) is a rare genetic disorder caused by the absence of a fully functional dystrophin protein in myocytes. In skeletal muscle, the lack of dystrophin ultimately results in muscle wasting and the replacement of myocytes with fatty or fibrous tissues. In the heart, cardiomyocytes eventually fail and cause fatal cardiomyopathy. We present a case of a male patient and his younger brother with a maternally inherited inverted insertion of approximately 306 kb of chromosome 10 in the deep intronic region between exons 44 and 45 of the DMD gene, leading to Duchenne muscular dystrophy. Chromosomal microarray, comprehensive muscular dystrophy genetic testing, and whole exome sequencing were negative. Targeted transcriptome RNA sequencing at an external lab showed no aberrant splicing. Research whole genome sequencing identified the copy number gain and insertion. Subsequent reanalysis of the RNA sequencing data showed possible aberrant splicing involving DMD exons 44-45, and research RNA sequencing revealed a fusion between the DMD gene on the minus strand of chromosome X and the PFKFB3 gene on the plus strand of chromosome 10. We demonstrate that whole genome sequencing can be valuable for identifying intronic events in the DMD gene previously undetected or not reported by traditional clinical testing.

Prior Fc receptor activation primes macrophages for increased sensitivity to IgG via long-term and short-term mechanisms

Macrophages measure the "eat-me" signal immunoglobulin G (IgG) to identify targets for phagocytosis. We tested whether prior encounters with IgG influence macrophage appetite. IgG is recognized by the Fc receptor. To temporally control Fc receptor activation, we engineered an Fc receptor that is activated by the light-induced oligomerization of Cry2, triggering phagocytosis. Using this tool, we demonstrate that subthreshold Fc receptor activation primes mouse bone-marrow-derived macrophages to be more sensitive to IgG in future encounters. Macrophages that have previously experienced subthreshold Fc receptor activation eat more IgG-bound human cancer cells. Increased phagocytosis occurs by two discrete mechanisms-a short- and long-term priming. Long-term priming requires new protein synthesis and Erk activity. Short-term priming does not require new protein synthesis and correlates with an increase in Fc receptor mobility. Our work demonstrates that IgG primes macrophages for increased phagocytosis, suggesting that therapeutic antibodies may become more effective after initial priming doses.

StableLift: Optimized Germline and Somatic Variant Detection Across Genome Builds

Reference genomes are foundational to modern genomics. Our growing understanding of genome structure leads to continual improvements in reference genomes and new genome "builds" with incompatible coordinate systems. We quantified the impact of genome build on germline and somatic variant calling by analyzing tumour-normal whole-genome pairs against the two most widely used human genome builds. The average individual had a build-discordance of 3.8% for germline SNPs, 8.6% for germline SVs, 25.9% for somatic SNVs and 49.6% for somatic SVs. Build-discordant variants are not simply false-positives: 47% were verified by targeted resequencing. Build-discordant variants were associated with specific genomic and technical features in variant- and algorithm-specific patterns. We leveraged these patterns to create StableLift, an algorithm that predicts cross-build stability with AUROCs of 0.934 ± 0.029. These results call for significant caution in cross-build analyses and for use of StableLift as a computationally efficient solution to mitigate inter-build artifacts.

Cell-free DNA from ascites identifies clinically relevant variants and tumour evolution in patients with advanced ovarian cancer

The emergence of targeted therapies has transformed ovarian cancer treatment. However, biomarker profiling for precision medicine is limited by access to quality, tumour-enriched tissue samples. The use of cell-free DNA (cfDNA) in ascites presents a potential solution to this challenge. In this study, next-generation sequencing was performed on ascites-derived cfDNA samples (26 samples from 15 human participants with ovarian cancer), with matched DNA from ascites-derived tumour cells (n = 5) and archived formalin-fixed paraffin-embedded (FFPE) tissue (n = 5). Similar tumour purity and variant detection were achieved with cfDNA compared to FFPE and ascites cell DNA. Analysis of large-scale genomic alterations, loss of heterozygosity and tumour mutation burden identified six cases of high genomic instability (including four with pathogenic BRCA1 and BRCA2 mutations). Copy number profiles and subclone prevalence changed between sequential ascites samples, particularly in a case where deletions and chromothripsis in Chr17p13.1 and Chr8q resulted in changes in clinically relevant TP53 and MYC variants over time. Ascites cfDNA identified clinically actionable information, concordant to tissue biopsies, enabling opportunistic molecular profiling. This advocates for analysis of ascites cfDNA in lieu of accessing tumour tissue via biopsy.

Repeat-Rich Regions Cause False-Positive Detection of NUMTs: A Case Study in Amphibians Using an Improved Cane Toad Reference Genome

Mitochondrial DNA (mtDNA) has been widely used in genetics research for decades. Contamination from nuclear DNA of mitochondrial origin (NUMTs) can confound studies of phylogenetic relationships and mtDNA heteroplasmy. Homology searches with mtDNA are widely used to detect NUMTs in the nuclear genome. Nevertheless, false-positive detection of NUMTs is common when handling repeat-rich sequences, while fragmented genomes might result in missing true NUMTs. In this study, we investigated different NUMT detection methods and how the quality of the genome assembly affects them. We presented an improved nuclear genome assembly (aRhiMar1.3) of the invasive cane toad (Rhinella marina) with additional long-read Nanopore and 10× linked-read sequencing. The final assembly was 3.47 Gb in length with 91.3% of tetrapod universal single-copy orthologs (n = 5,310), indicating the gene-containing regions were well assembled. We used 3 complementary methods (NUMTFinder, dinumt, and PALMER) to study the NUMT landscape of the cane toad genome. All 3 methods yielded consistent results, showing very few NUMTs in the cane toad genome. Furthermore, we expanded NUMT detection analyses to other amphibians and confirmed a weak relationship between genome size and the number of NUMTs present in the nuclear genome. Amphibians are repeat-rich, and we show that the number of NUMTs found in highly repetitive genomes is prone to inflation when using homology-based detection without filters. Together, this study provides an exemplar of how to robustly identify NUMTs in complex genomes when confounding effects on mtDNA analyses are a concern.

Droplet Hi-C enables scalable, single-cell profiling of chromatin architecture in heterogeneous tissues

Current methods for analyzing chromatin architecture are not readily scalable to heterogeneous tissues. Here we introduce Droplet Hi-C, which uses a commercial microfluidic device for high-throughput, single-cell chromatin conformation profiling in droplets. Using Droplet Hi-C, we mapped the chromatin architecture of the mouse cortex and analyzed gene regulatory programs in major cortical cell types. In addition, we used this technique to detect copy number variations, structural variations and extrachromosomal DNA in human glioblastoma, colorectal and blood cancer cells, revealing clonal dynamics and other oncogenic events during treatment. We refined the technique to allow joint profiling of chromatin architecture and transcriptome in single cells, facilitating exploration of the links between chromatin architecture and gene expression in both normal tissues and tumors. Thus, Droplet Hi-C both addresses critical gaps in chromatin analysis of heterogeneous tissues and enhances understanding of gene regulation.

LINC01432 binds to CELF2 in newly diagnosed multiple myeloma promoting short progression-free survival to standard therapy

Multiple Myeloma (MM) is an incurable form of cancer that arises from malignant plasma cells, with over 35,000 new cases diagnosed annually in the United States. While there are a growing number of approved therapies, MM remains incurable and nearly all patients will relapse and exhaust available treatments. Mechanisms for disease progression are unclear and little is known regarding the role of long non-coding RNAs (lncRNA) in mediating disease progression and response to treatment. Here, we used transcriptome sequencing to compare newly diagnosed MM (NDMM) patients who had short progression-free survival (PFS) to standard first-line treatment (PFS < 24 months) to patients who had prolonged PFS (PFS > 24 months). We identified 157 differentially upregulated lncRNAs with short PFS and focused our efforts on characterizing the most upregulated lncRNA, LINC01432. We investigated LINC01432 to show that its overexpression significantly increases cell viability and reduces apoptosis, while knockdown significantly reduces viability and increases apoptosis. Next, we show that LINC01432 directly interacts with the RNA binding protein, CELF2. Lastly, we showed that LINC01432-targeted locked nucleic acid antisense oligonucleotides reduce viability and increases apoptosis. In summary, this fundamental study identified lncRNAs associated with short PFS to standard NDMM treatment and further characterized LINC01432.

Graphasing: phasing diploid genome assembly graphs with single-cell strand sequencing

Haplotype information is crucial for biomedical and population genetics research. However, current strategies to produce de novo haplotype-resolved assemblies often require either difficult-to-acquire parental data or an intermediate haplotype-collapsed assembly. Here, we present Graphasing, a workflow which synthesizes the global phase signal of Strand-seq with assembly graph topology to produce chromosome-scale de novo haplotypes for diploid genomes. Graphasing readily integrates with any assembly workflow that both outputs an assembly graph and has a haplotype assembly mode. Graphasing performs comparably to trio phasing in contiguity, phasing accuracy, and assembly quality, outperforms Hi-C in phasing accuracy, and generates human assemblies with over 18 chromosome-spanning haplotypes.

Toward building a comprehensive human pan-genome: The SEN-GENOME project

The human reference genome (GRCh38), primarily sourced from individuals of European descent, falls short in capturing the vast genetic diversity across global populations. Efforts to diversify the reference genome face challenges in accessibility and representation, exacerbating the scarcity of African genomic data crucial for studying diseases prevalent in these populations. Sherman et al. proposed constructing reference genomes tailored to distinct human sub-populations. Their African Pan-Genome initiative highlighted substantial genetic variation missing from the GRCh38 human reference genome, emphasizing the necessity for population-specific genomes. In response, local initiatives like the Senegalese Genome project (SEN-GENOME) have emerged to document the genomes of historically overlooked populations. SEN-GENOME embodies community-driven decentralized research. With meticulous recruitment criteria and ethical practices, it aims to sequence 1,000 genomes from 31 ethnolinguistic groups, in the fourteen administrative regions of Senegal, fostering local genomic research tailored to the region. The key to SEN-GENOME's success is its commitment to local governance of data, capacity building, and integration with broader pan-genome projects in Africa. Despite the complexities of data harmonization and sharing, our collaborative efforts are aligned with common goals, ensuring steady progress toward a comprehensive human pan-genome. We invite and welcome collaboration with other research entities to achieve this shared vision. In summary, local initiatives such as SEN-GENOME are pivotal in bridging genomic disparities, offering pathways to equitable and inclusive genomic research. Collaborative endeavors guided by a collective vision for human health will propel us toward a more encompassing understanding of the human genome and better health through genomic medicine.

Genetic complexity of killer-cell immunoglobulin-like receptor genes in human pangenome assemblies

The killer-cell immunoglobulin-like receptor (KIR) gene complex, a highly polymorphic region of the human genome that encodes proteins involved in immune responses, poses strong challenges in genotyping owing to its remarkable genetic diversity and structural intricacy. Accurate analysis of KIR alleles, including their structural variations, is crucial for understanding their roles in various immune responses. Leveraging the high-quality genome assemblies from the Human Pangenome Reference Consortium (HPRC), we present a novel bioinformatic tool, the structural KIR annoTator (SKIRT), to investigate gene diversity and facilitate precise KIR allele analysis. In 47 HPRC-phased assemblies, SKIRT identifies a recurrent novel KIR2DS4/3DL1 fusion gene in the paternal haplotype of HG02630 and maternal haplotype of NA19240. Additionally, SKIRT accurately identifies eight structural variants and 15 novel nonsynonymous alleles, all of which are independently validated using short-read data or quantitative polymerase chain reaction. Our study has discovered a total of 570 novel alleles, among which eight haplotypes harbor at least one KIR gene duplication, six haplotypes have lost at least one framework gene, and 75 out of 94 haplotypes (79.8%) carry at least five novel alleles, thus confirming KIR genetic diversity. These findings are pivotal in providing insights into KIR gene diversity and serve as a solid foundation for understanding the functional consequences of KIR structural variations. High-resolution genome assemblies offer unprecedented opportunities to explore polymorphic regions that are challenging to investigate using short-read sequencing methods. The SKIRT pipeline emerges as a highly efficient tool, enabling the comprehensive detection of the complete spectrum of KIR alleles within human genome assemblies.

Stress increases sperm respiration and motility in mice and men

Semen quality and fertility has declined over the last 50 years, corresponding to ever-increasing environmental stressors. However, the cellular mechanisms involved and their impact on sperm functions remain unknown. In a repeated sampling human cohort study, we identify a significant effect of prior perceived stress to increase sperm motility 2-3 months following stress, timing that expands upon our previous studies revealing significant stress-associated changes in sperm RNA important for fertility. We mechanistically examine this post-stress timing in mice using an in vitro stress model in the epididymal epithelial cells responsible for sperm maturation and find 7282 differentially H3K27me3 bound DNA regions involving genes critical for mitochondrial and metabolic pathways. Further, prior stress exposure significantly changes the composition and size of epithelial cell-secreted extracellular vesicles that when incubated with mouse sperm, increase mitochondrial respiration and sperm motility, adding to our prior work showing impacts on embryo development. Together, these studies identify a time-dependent, translational signaling pathway that communicates stress experience to sperm, ultimately affecting reproductive functions.

Rare genomic copy number variants implicate new candidate genes for bicuspid aortic valve

Bicuspid aortic valve (BAV), the most common congenital heart defect, is a major cause of aortic valve disease requiring valve interventions and thoracic aortic aneurysms predisposing to acute aortic dissections. The spectrum of BAV ranges from early onset valve and aortic complications (EBAV) to sporadic late onset disease. Rare genomic copy number variants (CNVs) have previously been implicated in the development of BAV and thoracic aortic aneurysms. We determined the frequency and gene content of rare CNVs in EBAV probands (n = 272) using genome-wide SNP microarray analysis and three complementary CNV detection algorithms (cnvPartition, PennCNV, and QuantiSNP). Unselected control genotypes from the Database of Genotypes and Phenotypes were analyzed using identical methods. We filtered the data to select large genic CNVs that were detected by multiple algorithms. Findings were replicated in a BAV cohort with late onset sporadic disease (n = 5040). We identified 3 large and rare (< 1,1000 in controls) CNVs in EBAV probands. The burden of CNVs intersecting with genes known to cause BAV when mutated was increased in case-control analysis. CNVs intersecting with GATA4 and DSCAM were enriched in cases, recurrent in other datasets, and segregated with disease in families. In total, we identified potentially pathogenic CNVs in 9% of EBAV cases, implicating alterations of candidate genes at these loci in the pathogenesis of BAV.

Comprehensive network modeling approaches unravel dynamic enhancer-promoter interactions across neural differentiation

BACKGROUND

Increasing evidence suggests that a substantial proportion of disease-associated mutations occur in enhancers, regions of non-coding DNA essential to gene regulation. Understanding the structures and mechanisms of the regulatory programs this variation affects can shed light on the apparatuses of human diseases.

RESULTS

We collect epigenetic and gene expression datasets from seven early time points during neural differentiation. Focusing on this model system, we construct networks of enhancer-promoter interactions, each at an individual stage of neural induction. These networks serve as the base for a rich series of analyses, through which we demonstrate their temporal dynamics and enrichment for various disease-associated variants. We apply the Girvan-Newman clustering algorithm to these networks to reveal biologically relevant substructures of regulation. Additionally, we demonstrate methods to validate predicted enhancer-promoter interactions using transcription factor overexpression and massively parallel reporter assays.

CONCLUSIONS

Our findings suggest a generalizable framework for exploring gene regulatory programs and their dynamics across developmental processes; this includes a comprehensive approach to studying the effects of disease-associated variation on transcriptional networks. The techniques applied to our networks have been published alongside our findings as a computational tool, E-P-INAnalyzer. Our procedure can be utilized across different cellular contexts and disorders.

Beyond the Human Genome Project: The Age of Complete Human Genome Sequences and Pangenome References

The Human Genome Project was an enormous accomplishment, providing a foundation for countless explorations into the genetics and genomics of the human species. Yet for many years, the human genome reference sequence remained incomplete and lacked representation of human genetic diversity. Recently, two major advances have emerged to address these shortcomings: complete gap-free human genome sequences, such as the one developed by the Telomere-to-Telomere Consortium, and high-quality pangenomes, such as the one developed by the Human Pangenome Reference Consortium. Facilitated by advances in long-read DNA sequencing and genome assembly algorithms, complete human genome sequences resolve regions that have been historically difficult to sequence, including centromeres, telomeres, and segmental duplications. In parallel, pangenomes capture the extensive genetic diversity across populations worldwide. Together, these advances usher in a new era of genomics research, enhancing the accuracy of genomic analysis, paving the path for precision medicine, and contributing to deeper insights into human biology.

Metagenome-wide characterization of shared antimicrobial resistance genes in sympatric people and lemurs in rural Madagascar

Background

Tracking the spread of antibiotic resistant bacteria is critical to reduce global morbidity and mortality associated with human and animal infections. There is a need to understand the role that wild animals in maintenance and transfer of antibiotic resistance genes (ARGs).

Methods

This study used metagenomics to identify and compare the abundance of bacterial species and ARGs detected in the gut microbiomes from sympatric humans and wild mouse lemurs in a forest-dominated, roadless region of Madagascar near Ranomafana National Park. We examined the contribution of human geographic location toward differences in ARG abundance and compared the genomic similarity of ARGs between host source microbiomes.

Results

Alpha and beta diversity of species and ARGs between host sources were distinct but maintained a similar number of detectable ARG alleles. Humans were differentially more abundant for four distinct tetracycline resistance-associated genes compared to lemurs. There was no significant difference in human ARG diversity from different locations. Human and lemur microbiomes shared 14 distinct ARGs with highly conserved in nucleotide identity. Synteny of ARG-associated assemblies revealed a distinct multidrug-resistant gene cassette carrying dfrA1 and aadA1 present in human and lemur microbiomes without evidence of geographic overlap, suggesting that these resistance genes could be widespread in this ecosystem. Further investigation into intermediary processes that maintain drug-resistant bacteria in wildlife settings is needed.

IGLoo: Profiling the Immunoglobulin Heavy chain locus in Lymphoblastoid Cell Lines with PacBio High-Fidelity Sequencing reads

New high-quality human genome assemblies derived from lymphoblastoid cell lines (LCLs) provide reference genomes and pangenomes for genomics studies. However, the characteristics of LCLs pose technical challenges to profiling immunoglobulin (IG) genes. IG loci in LCLs contain a mixture of germline and somatically recombined haplotypes, making them difficult to genotype or assemble accurately. To address these challenges, we introduce IGLoo, a software tool that implements novel methods for analyzing sequence data and genome assemblies derived from LCLs. IGLoo characterizes somatic V(D)J recombination events in the sequence data and identifies the breakpoints and missing IG genes in the LCL-based assemblies. Furthermore, IGLoo implements a novel reassembly framework to improve germline assembly quality by integrating information about somatic events and population structural variantions in the IG loci. We applied IGLoo to study the assemblies from the Human Pangenome Reference Consortium, providing new insights into the mechanisms, gene usage, and patterns of V(D)J recombination, causes of assembly fragmentation in the IG heavy chain (IGH) locus, and improved representation of the IGH assemblies.

Impact of genome build on RNA-seq interpretation and diagnostics

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network and Genomics Research to Elucidate the Genetics of Rare Disease Consortium. Across six routinely collected biospecimens, 61% of quantified genes were not influenced by genome build. However, we identified 1,492 genes with build-dependent quantification, 3,377 genes with build-exclusive expression, and 9,077 genes with annotation-specific expression across six routinely collected biospecimens, including 566 clinically relevant and 512 known OMIM genes. Further, we demonstrate that between builds for a given gene, a larger difference in quantification is well correlated with a larger change in expression outlier calling. Combined, we provide a database of genes impacted by build choice and recommend that transcriptomics-guided analyses and diagnoses are cross referenced with these data for robustness.

Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample

Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.

Vaginal and rectal microbiome contribute to genital inflammation in chronic pelvic pain

BACKGROUND

Chronic pelvic pain (CPP) is a multifactorial syndrome that can substantially affect a patient's quality of life. Endometriosis is one cause of CPP, and alterations of the immune and microbiome profiles have been observed in patients with endometriosis. The objective of this pilot study was to investigate differences in the vaginal and gastrointestinal microbiomes and cervicovaginal immune microenvironment in patients with CPP and endometriosis diagnosis compared to those with CPP without endometriosis and no CPP.

METHODS

Vaginal swabs, rectal swabs, and cervicovaginal lavages (CVL) were collected among individuals undergoing gynecologic laparoscopy. Participants were grouped based on patients seeking care for chronic pain and/or pathology results: CPP and endometriosis (CPP-Endo) (n = 35), CPP without endometriosis (n = 23), or patients without CPP or endometriosis (controls) (n = 15). Sensitivity analyses were performed on CPP with endometriosis location, stage, and co-occurring gynecologic conditions (abnormal uterine bleeding, fibroids). 16S rRNA sequencing was performed to profile the microbiome, and a panel of soluble immune mediators was quantified using a multiplex assay. Statistical analysis was conducted with SAS, R, MicrobiomeAnalyst, MetaboAnalyst, and QIIME 2.

RESULTS

Significant differences were observed between participants with CPP alone, CPP-Endo, and surgical controls for body mass index, ethnicity, diagnosis of ovarian cysts, and diagnosis of fibroids. In rectal microbiome analysis, both CPP alone and CPP-Endo exhibited lower alpha diversity than controls, and both CPP groups revealed enrichment of irritable bowel syndrome-associated bacteria. CPP-Endo exhibited an increased abundance of vaginal Streptococcus anginosus and rectal Ruminococcus. Patients with CPP and endometrioma (s) demonstrated increased vaginal Streptococcus, Lactobacillus, and Prevotella compared to other endometriosis sites. Further, abnormal uterine bleeding was associated with an increased abundance of bacterial vaginosis-associated bacteria. Immunoproteomic profiles were distinctly clustered by CPP alone and CPP-Endo compared to controls. CPP-Endo was enriched in TNF⍺, MDC, and IL-1⍺.

CONCLUSIONS

Vaginal and rectal microbiomes were observed to differ between patients with CPP alone and CPP with endometriosis, which may be useful in personalized treatment for individuals with CPP and endometriosis from those with other causes of CPP. Further investigation is warranted in patients with additional co-occurring conditions, such as AUB/fibroids, which add additional complexity to these conditions and reveal the enrichment of distinct pathogenic bacteria in both mucosal sites. This study provides foundational microbiome-immunoproteomic knowledge related to chronic pelvic pain, endometriosis, and co-occurring gynecologic conditions that can help improve the treatment of patients seeking care for pain.

Insights into non-crossover recombination from long-read sperm sequencing

Meiotic recombination is a fundamental process that generates genetic diversity by creating new combinations of existing alleles. Although human crossovers have been studied at the pedigree, population and single-cell level, the more frequent non-crossover events that lead to gene conversion are harder to study, particularly at the individual level. Here we show that single high-fidelity long sequencing reads from sperm can capture both crossovers and non-crossovers, allowing effectively arbitrary sample sizes for analysis from one male. Using fifteen sperm samples from thirteen donors we demonstrate variation between and within donors for the rates of different types of recombination. Intriguingly, we observe a tendency for non-crossover gene conversions to occur upstream of nearby PRDM9 binding sites, whereas crossover locations have a slight downstream bias. We further provide evidence for two distinct non-crossover processes. One gives rise to the vast majority of non-crossovers with mean conversion tract length under 50bp, which we suggest is an outcome of standard PRDM9-induced meiotic recombination. In contrast ~2% of non-crossovers have much longer mean tract length, and potentially originate from the same process as complex events with more than two haplotype switches, which is not associated with PRDM9 binding sites and is also seen in somatic cells.

Every CFTR variant counts - Target-capture based next-generation-sequencing for molecular diagnosis in the German CF Registry

BACKGROUND

In times of genotype guided therapy options, a total of 3.2 % of people with CF (pwCF) in the German CF Registry[1] only have one or no CFTR-variant detected after genetic analysis. Additionally, genetic data in the Registry can be documented as free text and can therefore be prone to error. In order to allow the greatest possible amount of pwCF access to modern therapies, we conducted a re-evaluation of free text entries and established a custom-whole-CFTR-locus NGS-approach for all pwCF who remained without genetic confirmation afterwards.

METHODS

To this end, we assembled 731 free text variants of 655 pwCF in the German CF Registry. All variants were evaluated using ClinVar, HGMD and CFTR1/2, corrected in the Registries' database and uploaded to ClinVar. PwCF whose diagnosis remained uncertain as well as additional pwCF or pwCFTR-RD that were assembled through a nationwide call for testing of unclear cases were offered genetic analysis. Samples were analysed using a target-capture based NGS-custom-design-panel covering the entire CFTR-locus.

RESULTS

Evaluation of free text variants led to the discovery of 43 variants not formerly reported in the context of CF. The Registries' dropdown list was extended by 497 variants and over 500 pwCF were provided with their most up-to-date genotype. Samples of 47 pwCF/pwCFTR-RD were sequenced via NGS with an overall success rate of 61.7 %, resulting in implementation of entire CFTR-genotyping into routine diagnostics.

CONCLUSION

Entire CFTR-genotyping can greatly increase the genetic diagnostic rate of pwCF/pwCFTR-RD and should be considered after inconspicuous CFTR screening panels in CFTR-diagnostics.

LINC01432 binds to CELF2 in newly diagnosed multiple myeloma promoting short progression-free survival to standard therapy

Multiple Myeloma (MM) is a highly prevalent and incurable form of cancer that arises from malignant plasma cells, with over 35,000 new cases diagnosed annually in the United States. While there are a growing number of approved therapies, MM remains incurable and nearly all patients will relapse and exhaust all available treatment options. Mechanisms for disease progression are unclear and in particular, little is known regarding the role of long non-coding RNAs (lncRNA) in mediating disease progression and response to treatment. In this study, we used transcriptome sequencing to compare newly diagnosed MM patients who had short progression-free survival (PFS) to standard first-line treatment (PFS < 24 months) to patients who had prolonged PFS (PFS > 24 months). We identified 157 differentially upregulated lncRNAs with short PFS and focused our efforts on characterizing the most upregulated lncRNA, LINC01432. We investigated LINC01432 overexpression and CRISPR/Cas9 knockdown in MM cell lines to show that LINC01432 overexpression significantly increases cell viability and reduces apoptosis, while knockdown significantly reduces viability and increases apoptosis, supporting the clinical relevance of this lncRNA. Next, we used individual-nucleotide resolution cross-linking immunoprecipitation with RT-qPCR to show that LINC01432 directly interacts with the RNA binding protein, CELF2. Lastly, we showed that LINC01432-targeted locked nucleic acid antisense oligonucleotides reduce viability and increases apoptosis. In summary, this fundamental study identified lncRNAs associated with short PFS to standard NDMM treatment and further characterized LINC01432, which inhibits apoptosis.