Recombination between heterologous human acrocentric chromosomes.

The short arms of the human acrocentric chromosomes 13, 14, 15, 21 and 22 (SAACs) share large homologous regions, including ribosomal DNA repeats and extended segmental duplications1,2. Although the resolution of these regions in the first complete assembly of a human genome-the Telomere-to-Telomere Consortium's CHM13 assembly (T2T-CHM13)-provided a model of their homology3, it remained unclear whether these patterns were ancestral or maintained by ongoing recombination exchange. Here we show that acrocentric chromosomes contain pseudo-homologous regions (PHRs) indicative of recombination between non-homologous sequences. Utilizing an all-to-all comparison of the human pangenome from the Human Pangenome Reference Consortium4 (HPRC), we find that contigs from all of the SAACs form a community. A variation graph5 constructed from centromere-spanning acrocentric contigs indicates the presence of regions in which most contigs appear nearly identical between heterologous acrocentric chromosomes in T2T-CHM13. Except on chromosome 15, we observe faster decay of linkage disequilibrium in the pseudo-homologous regions than in the corresponding short and long arms, indicating higher rates of recombination6,7. The pseudo-homologous regions include sequences that have previously been shown to lie at the breakpoint of Robertsonian translocations8, and their arrangement is compatible with crossover in inverted duplications on chromosomes 13, 14 and 21. The ubiquity of signals of recombination between heterologous acrocentric chromosomes seen in the HPRC draft pangenome suggests that these shared sequences form the basis for recurrent Robertsonian translocations, providing sequence and population-based confirmation of hypotheses first developed from cytogenetic studies 50 years ago9.

A draft human pangenome reference.

Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals1. These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample.

Building pangenome graphs.

Pangenome graphs can represent all variation between multiple reference genomes, but current approaches to build them exclude complex sequences or are based upon a single reference. In response, we developed the PanGenome Graph Builder, a pipeline for constructing pangenome graphs without bias or exclusion. The PanGenome Graph Builder uses all-to-all alignments to build a variation graph in which we can identify variation, measure conservation, detect recombination events and infer phylogenetic relationships.

A novel quantitative trait locus implicates Msh3 in the propensity for genome-wide short tandem repeat expansions in mice.

Short tandem repeats (STRs) are a class of rapidly mutating genetic elements typically characterized by repeated units of 1-6 bp. We leveraged whole-genome sequencing data for 152 recombinant inbred (RI) strains from the BXD family of mice to map loci that modulate genome-wide patterns of new mutations arising during parent-to-offspring transmission at STRs. We defined quantitative phenotypes describing the numbers and types of germline STR mutations in each strain and performed quantitative trait locus (QTL) analyses for each of these phenotypes. We identified a locus on Chromosome 13 at which strains inheriting the C57BL/6J (B) haplotype have a higher rate of STR expansions than those inheriting the DBA/2J (D) haplotype. The strongest candidate gene in this locus is Msh3, a known modifier of STR stability in cancer and at pathogenic repeat expansions in mice and humans, as well as a current drug target against Huntington's disease. The D haplotype at this locus harbors a cluster of variants near the 5' end of Msh3, including multiple missense variants near the DNA mismatch recognition domain. In contrast, the B haplotype contains a unique retrotransposon insertion. The rate of expansion covaries positively with Msh3 expression-with higher expression from the B haplotype. Finally, detailed analysis of mutation patterns showed that strains carrying the B allele have higher expansion rates, but slightly lower overall total mutation rates, compared with those with the D allele, particularly at tetranucleotide repeats. Our results suggest an important role for inherited variants in Msh3 in modulating genome-wide patterns of germline mutations at STRs.

Efficient dynamic variation graphs.

MOTIVATION: Pangenomics is a growing field within computational genomics. Many pangenomic analyses use bidirected sequence graphs as their core data model. However, implementing and correctly using this data model can be difficult, and the scale of pangenomic datasets can be challenging to work at. These challenges have impeded progress in this field. RESULTS: Here, we present a stack of two C++ libraries, libbdsg and libhandlegraph, which use a simple, field-proven interface, designed to expose elementary features of these graphs while preventing common graph manipulation mistakes. The libraries also provide a Python binding. Using a diverse collection of pangenome graphs, we demonstrate that these tools allow for efficient construction and manipulation of large genome graphs with dense variation. For instance, the speed and memory usage are up to an order of magnitude better than the prior graph implementation in the VG toolkit, which has now transitioned to using libbdsg's implementations. AVAILABILITY AND IMPLEMENTATION: libhandlegraph and libbdsg are available under an MIT License from https://github.com/vgteam/libhandlegraph and https://github.com/vgteam/libbdsg.

Maternal exome analysis for the diagnosis of oocyte maturation defects and early embryonic developmental arrest.

RESEARCH QUESTION: Can a methodology be developed for case selection and whole-exome sequencing (WES) analysis of women who are infertile owing to recurrent oocyte maturation defects (OOMD) and/or preimplantation embryo lethality (PREMBL)? DESIGN: Data were collected from IVF patients attending the Istanbul Memorial Hospital (2015-2021). A statistical methodology to identify infertile endophenotypes (recurrent low oocyte maturation rate, low fertilization rate and preimplantation developmental arrest) was developed using a large IVF dataset (11,221 couples). Twenty-eight infertile women with OOMD/PREMBL were subsequently enrolled for WES on their genomic DNA. Pathogenic variants were prioritized using a custom-made bioinformatic pipeline set to minimize false-positive discoveries through resampling in control cohorts (the Human Genome Diversity Project and 1343 whole-exome sequences from oocyte donors). Individual single-cell RNA sequencing data from 18 human metaphase II (MII) oocytes and antral granulosa cells was used for genome-wide validation. WES and bioinformatics were performed at Igenomix and the National Research Council, Italy. RESULTS: Variant prioritization analysis identified 265 unique variants in 248 genes (average 22.4 per sample). Of the genes harbouring high-impact variants 78% were expressed by MII oocytes and/or antral granulosa cells, significantly higher than for random sample of controls (odds ratio = 5, Fisher's exact P = 0.0004). Seven of the 28 women (25%) were homozygous carriers of missense pathogenic variants in known candidate genes for OOMD/PREMBL, including PATL2, NLRP5 (n = 2),TLE6, PADI6, TUBB8 and TRIP13. Furthermore, novel gene-disease associations were identified. In fact, one woman with a low oocyte maturation rate was a homozygous carrier of high-impact variants in ENSA, an essential gene for prophase I meiotic transition in mice. CONCLUSIONS: This analytical framework could reveal known and new genes associated with isolated recurrent OOMD/PREMBL, providing essential indications for scaling this strategy to larger studies.

Identification of sex determination genes and their evolution in Phlebotominae sand flies (Diptera, Nematocera).

BACKGROUND: Phlebotomine sand flies (Diptera, Nematocera) are important vectors of several pathogens, including Leishmania parasites, causing serious diseases of humans and dogs. Despite their importance as disease vectors, most aspects of sand fly biology remain unknown including the molecular basis of their reproduction and sex determination, aspects also relevant for the development of novel vector control strategies. RESULTS: Using comparative genomics/transcriptomics data mining and transcriptional profiling, we identified the sex determining genes in phlebotomine sand flies and proposed the first model for the sex determination cascade of these insects. For all the genes identified, we produced manually curated gene models, developmental gene expression profile and performed evolutionary molecular analysis. We identified and characterized, for the first time in a Nematocera species, the transformer (tra) homolog which exhibits both conserved and novel features. The analysis of the tra locus in sand flies and its expression pattern suggest that this gene is able to autoregulate its own splicing, as observed in the fruit fly Ceratitis capitata and several other insect species. CONCLUSIONS: Our results permit to fill the gap about sex determination in sand flies, contribute to a better understanding of this developmental pathway in Nematocera and open the way for the identification of sex determining orthologs in other species of this important Diptera sub-order. Furthermore, the sex determination genes identified in our work also provide the opportunity of future biotechnological applications to control natural population of sand flies, reducing their impact on public health.

Vitamin B12 ameliorates the phenotype of a mouse model of DiGeorge syndrome.

Abstract: Pathological conditions caused by reduced dosage of a gene, such as gene haploinsufficiency, can potentially be reverted by enhancing the expression of the functional allele. In practice, low specificity of therapeutic agents, or their toxicity reduces their clinical applicability. Here, we have used a high throughput screening (HTS) approach to identify molecules capable of increasing the expression of the gene Tbx1, which is involved in one of the most common gene haploinsufficiency syndromes, the 22q11.2 deletion syndrome. Surprisingly, we found that one of the two compounds identified by the HTS is the vitamin B12. Validation in a mouse model demonstrated that vitamin B12 treatment enhances Tbx1 gene expression and partially rescues the haploinsufficiency phenotype. These results lay the basis for preclinical and clinical studies to establish the effectiveness of this drug in the human syndrome.

Development and validation of a comprehensive genomic diagnostic tool for myeloid malignancies.

The diagnosis of hematologic malignancies relies on multidisciplinary workflows involving morphology, flow cytometry, cytogenetic, and molecular genetic analyses. Advances in cancer genomics have identified numerous recurrent mutations with clear prognostic and/or therapeutic significance to different cancers. In myeloid malignancies, there is a clinical imperative to test for such mutations in mainstream diagnosis; however, progress toward this has been slow and piecemeal. Here we describe Karyogene, an integrated targeted resequencing/analytical platform that detects nucleotide substitutions, insertions/deletions, chromosomal translocations, copy number abnormalities, and zygosity changes in a single assay. We validate the approach against 62 acute myeloid leukemia, 50 myelodysplastic syndrome, and 40 blood DNA samples from individuals without evidence of clonal blood disorders. We demonstrate robust detection of sequence changes in 49 genes, including difficult-to-detect mutations such as FLT3 internal-tandem and mixed-lineage leukemia (MLL) partial-tandem duplications, and clinically significant chromosomal rearrangements including MLL translocations to known and unknown partners, identifying the novel fusion gene MLL-DIAPH2 in the process. Additionally, we identify most significant chromosomal gains and losses, and several copy neutral loss-of-heterozygosity mutations at a genome-wide level, including previously unreported changes such as homozygosity for DNMT3A R882 mutations. Karyogene represents a dependable genomic diagnosis platform for translational research and for the clinical management of myeloid malignancies, which can be readily adapted for use in other cancers.

IFITM3 restricts the morbidity and mortality associated with influenza.

The 2009 H1N1 influenza pandemic showed the speed with which a novel respiratory virus can spread and the ability of a generally mild infection to induce severe morbidity and mortality in a subset of the population. Recent in vitro studies show that the interferon-inducible transmembrane (IFITM) protein family members potently restrict the replication of multiple pathogenic viruses. Both the magnitude and breadth of the IFITM proteins' in vitro effects suggest that they are critical for intrinsic resistance to such viruses, including influenza viruses. Using a knockout mouse model, we now test this hypothesis directly and find that IFITM3 is essential for defending the host against influenza A virus in vivo. Mice lacking Ifitm3 display fulminant viral pneumonia when challenged with a normally low-pathogenicity influenza virus, mirroring the destruction inflicted by the highly pathogenic 1918 'Spanish' influenza. Similar increased viral replication is seen in vitro, with protection rescued by the re-introduction of Ifitm3. To test the role of IFITM3 in human influenza virus infection, we assessed the IFITM3 alleles of individuals hospitalized with seasonal or pandemic influenza H1N1/09 viruses. We find that a statistically significant number of hospitalized subjects show enrichment for a minor IFITM3 allele (SNP rs12252-C) that alters a splice acceptor site, and functional assays show the minor CC genotype IFITM3 has reduced influenza virus restriction in vitro. Together these data reveal that the action of a single intrinsic immune effector, IFITM3, profoundly alters the course of influenza virus infection in mouse and humans.

Revisiting the thrifty gene hypothesis via 65 loci associated with susceptibility to type 2 diabetes.

We have investigated the evidence for positive selection in samples of African, European, and East Asian ancestry at 65 loci associated with susceptibility to type 2 diabetes (T2D) previously identified through genome-wide association studies. Selection early in human evolutionary history is predicted to lead to ancestral risk alleles shared between populations, whereas late selection would result in population-specific signals at derived risk alleles. By using a wide variety of tests based on the site frequency spectrum, haplotype structure, and population differentiation, we found no global signal of enrichment for positive selection when we considered all T2D risk loci collectively. However, in a locus-by-locus analysis, we found nominal evidence for positive selection at 14 of the loci. Selection favored the protective and risk alleles in similar proportions, rather than the risk alleles specifically as predicted by the thrifty gene hypothesis, and may not be related to influence on diabetes. Overall, we conclude that past positive selection has not been a powerful influence driving the prevalence of T2D risk alleles.

An Ethnolinguistic and Genetic Perspective on the Origins of the Dravidian-Speaking Brahui in Pakistan.

Pakistan is a part of South Asia that modern humans encountered soon after they left Africa ~50 - 70,000 years ago. Approximately 9,000 years ago they began establishing cities that eventually expanded to represent the Harappan culture, rivalling the early city states of Mesopotamia. The modern state constitutes the north western land mass of the Indian sub-continent and is now the abode of almost 200 million humans representing many ethnicities and linguistic groups. Studies utilising autosomal, Y chromosomal and mitochondrial DNA markers in selected Pakistani populations revealed a mixture of Western Eurasian-, South- and East Asian-specific lineages, some of which were unequivocally associated with past migrations. Overall in Pakistan, genetic relationships are generally predicted more accurately by geographic proximity than linguistic origin. The Dravidian-speaking Brahui population are a prime example of this. They currently reside in south-western Pakistan, surrounded by Indo-Europeans speakers with whom they share a common genetic origin. In contrast, the Hazara share the highest affinity with East Asians, despite their Indo-European linguistic affiliation. In this report we reexamine the genetic origins of the Brahuis, and compare them with diverse populations from India, including several Dravidian-speaking groups, and present a genetic perspective on ethnolinguistic groups in present-day Pakistan. Given the high affinity of Brahui to the other Indo-European Pakistani populations and the absence of population admixture with any of the examined Indian Dravidian groups, we conclude that Brahui are an example of cultural (linguistic) retention following a major population replacement.

Monoamine Oxidase A gene polymorphisms and self reported aggressive behaviour in a Pakistani ethnic group.

OBJECTIVE: To investigate the association of monoamine oxidase Agene polymorphisms with aggression. METHODS: The study was conducted in an ethnic community in Lahore, Pakistan, from August 2008 to December 2009 on the basis of data that was collected through a questionnaire between August 2004 and September 2005. It analysed 10 single nucleotide polymorphisms of monoamine oxidase A in unrelated males from the same ethnic background who were administered a Punjabi translation of the Buss and Perry aggression questionnaire. SPSS 13 was used for statistical analysis. RESULTS: Of the total 133 haplotypes studied, 52(39%) were Haplotype A, 58(43.6%) B, 8(6%) C, 3(2.3%) D, 9(6.8%) E and 3(2.3%) F. The six haplotypes were analysed for association with scores of the four subscales of the aggression questionnaire and multivariate analysis of variance showed no significant differences (p>0.05 each) in the error variances of the total scores and scores for three of the sub-scales across the haplotypes. The variance was significantly different only for the anger sub-scale (p<0.05). CONCLUSIONS: The association of an extended haplotype with low levels of self-reported aggression in this study should assist in characterisation of functional variants responsible for non-aggressive behaviour in male subjects.

Genome-wide meta-analysis of common variant differences between men and women.

The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency >0.05) single-nucleotide polymorphisms (SNPs). Allele frequencies were compared between men and women for directly-typed and imputed variants within each study. Forward-time simulations for unlinked, neutral, autosomal, common loci were performed under the demographic model for European populations with a fixed sex ratio and a random mating scheme to assess the probability of detecting significant allele frequency differences. We do not detect any genome-wide significant (P < 5 × 10(-8)) common SNP differences between men and women in this well-powered meta-analysis. The simulated data provided results entirely consistent with these findings. This large-scale investigation across ~115 000 individuals shows no detectable contribution from common genetic variants to the observed skew in the sex ratio. The absence of sex-specific differences is useful in guiding genetic association study design, for example when using mixed controls for sex-biased traits.

Genetic landscape of populations along the Silk Road: admixture and migration patterns.

BACKGROUND: The ancient Silk Road has been a trading route between Europe and Central Asia from the 2(nd) century BCE to the 15(th) century CE. While most populations on this route have been characterized, the genetic background of others remains poorly understood, and little is known about past migration patterns. The scientific expedition "Marco Polo" has recently collected genetic and phenotypic data in six regions (Georgia, Armenia, Azerbaijan, Uzbekistan, Kazakhstan, Tajikistan) along the Silk Road to study the genetics of a number of phenotypes. RESULTS: We characterized the genetic structure of these populations within a worldwide context. We observed a West-East subdivision albeit the existence of a genetic component shared within Central Asia and nearby populations from Europe and Near East. We observed a contribution of up to 50% from Europe and Asia to most of the populations that have been analyzed. The contribution from Asia dates back to ~25 generations and is limited to the Eastern Silk Road. Time and direction of this contribution are consistent with the Mongolian expansion era. CONCLUSIONS: We clarified the genetic structure of six populations from Central Asia and suggested a complex pattern of gene flow among them. We provided a map of migration events in time and space and we quantified exchanges among populations. Altogether these novel findings will support the future studies aimed at understanding the genetics of the phenotypes that have been collected during the Marco Polo campaign, they will provide insights into the history of these populations, and they will be useful to reconstruct the developments and events that have shaped modern Eurasians genomes.

A comprehensive analysis of SARS-CoV-2 missense mutations indicates that all possible amino acid replacements in the viral proteins occurred within the first two-and-a-half years of the pandemic.

The surveillance of COVID-19 pandemic has led to the determination of millions of genome sequences of the SARS-CoV-2 virus, with the accumulation of a wealth of information never collected before for an infectious disease. Exploring the information retrieved from the GISAID database reporting at that time >13 million genome sequences, we classified the 141,639 unique missense mutations detected in the first two-and-a-half years (up to October 2022) of the pandemic. Notably, our analysis indicates that 98.2 % of all possible conservative amino acid replacements occurred. Even non-conservative mutations were highly represented (73.9 %). For a significant number of residues (3 %), all possible replacements with the other nineteen amino acids have been observed. These observations strongly indicate that, in this time interval, the virus explored all possible alternatives in terms of missense mutations for all sites of its polypeptide chain and that those that are not observed severely affect SARS-CoV-2 integrity. The implications of the present findings go well beyond the structural biology of SARS-CoV-2 as the huge amount of information here collected and classified may be valuable for the elucidation of the sequence-structure-function relationships in proteins.

Pangenome reconstruction in rats enhances genotype-phenotype mapping and novel variant discovery.

The HXB/BXH family of recombinant inbred rat strains is a unique genetic resource that has been extensively phenotyped over 25 years, resulting in a vast dataset of quantitative molecular and physiological phenotypes. We built a pangenome graph from 10x Genomics Linked-Read data for 31 recombinant inbred rats to study genetic variation and association mapping. The pangenome includes 0.2Gb of sequence that is not present the reference mRatBN7.2, confirming the capture of substantial additional variation. We validated variants in challenging regions, including complex structural variants resolving into multiple haplotypes. Phenome-wide association analysis of validated SNPs uncovered variants associated with glucose/insulin levels and hippocampal gene expression. We propose an interaction between Pirl1l1, chromogranin expression, TNF-α levels, and insulin regulation. This study demonstrates the utility of linked-read pangenomes for comprehensive variant detection and mapping phenotypic diversity in a widely used rat genetic reference panel.

Three-dimensional environment sensitizes pancreatic cancer cells to the anti-proliferative effect of budesonide by reprogramming energy metabolism.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is the most lethal cancer with an aggressive metastatic phenotype and very poor clinical prognosis. Interestingly, a lower occurrence of PDAC has been described in individuals with severe and long-standing asthma. Here we explored the potential link between PDAC and the glucocorticoid (GC) budesonide, a first-line therapy to treat asthma. METHODS: We tested the effect of budesonide and the classical GCs on the morphology, proliferation, migration and invasiveness of patient-derived PDAC cells and pancreatic cancer cell lines, using 2D and 3D cultures in vitro. Furthermore, a xenograft model was used to investigate the effect of budesonide on PDAC tumor growth in vivo. Finally, we combined genome-wide transcriptome analysis with genetic and pharmacological approaches to explore the mechanisms underlying budesonide activities in the different environmental conditions. RESULTS: We found that in 2D culture settings, high micromolar concentrations of budesonide reduced the mesenchymal invasive/migrating features of PDAC cells, without affecting proliferation or survival. This activity was specific and independent of the Glucocorticoid Receptor (GR). Conversely, in a more physiological 3D environment, low nanomolar concentrations of budesonide strongly reduced PDAC cell proliferation in a GR-dependent manner. Accordingly, we found that budesonide reduced PDAC tumor growth in vivo. Mechanistically, we demonstrated that the 3D environment drives the cells towards a general metabolic reprogramming involving protein, lipid, and energy metabolism (e.g., increased glycolysis dependency). This metabolic change sensitizes PDAC cells to the anti-proliferative effect of budesonide, which instead induces opposite changes (e.g., increased mitochondrial oxidative phosphorylation). Finally, we provide evidence that budesonide inhibits PDAC growth, at least in part, through the tumor suppressor CDKN1C/p57Kip2. CONCLUSIONS: Collectively, our study reveals that the microenvironment influences the susceptibility of PDAC cells to GCs and provides unprecedented evidence for the anti-proliferative activity of budesonide on PDAC cells in 3D conditions, in vitro and in vivo. Our findings may explain, at least in part, the reason for the lower occurrence of pancreatic cancer in asthmatic patients and suggest a potential suitability of budesonide for clinical trials as a therapeutic approach to fight pancreatic cancer.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.