A common sequence motif associated with recombination hot spots and genome instability in humans

The formation and propagation of human Robertsonian chromosomes

Robertsonian chromosomes are a type of variant chromosome found commonly in nature. Present in one in 800 humans, these chromosomes can underlie infertility, trisomies, and increased cancer incidence. Recognized cytogenetically for more than a century, their origins have remained mysterious. Recent advances in genomics allowed us to assemble three human Robertsonian chromosomes completely. We identify a common breakpoint and epigenetic changes in centromeres that provide insight into the formation and propagation of common Robertsonian translocations. Further investigation of the assembled genomes of chimpanzee and bonobo highlights the structural features of the human genome that uniquely enable the specific crossover event that creates these chromosomes. Resolving the structure and epigenetic features of human Robertsonian chromosomes at a molecular level paves the way to understanding how chromosomal structural variation occurs more generally, and how chromosomes evolve.

Male-biased recombination at chromosome ends in a songbird revealed by precisely mapping crossover positions

Recombination plays a crucial role in evolution by generating novel haplotypes and disrupting linkage between genes, thereby enhancing the efficiency of selection. Here, we analyze the genomes of 12 great reed warblers (Acrocephalus arundinaceus) in a 3-generation pedigree to identify precise crossover positions along the chromosomes. We located more than 200 crossovers and found that these were highly concentrated toward the telomeric ends of the chromosomes. Apart from this major pattern in the recombination landscape, we found significantly higher frequencies of crossovers in genic compared with intergenic regions, and in exons compared with introns. Moreover, while the number of recombination events was similar between the sexes, the crossovers were located significantly closer to the ends of paternal compared with maternal chromosomes. In conclusion, our study of the great reed warbler revealed substantial variation in crossover frequencies within chromosomes, with a distinct bias toward the sub-telomeric regions, particularly on the paternal side. These findings emphasize the importance of thoroughly screening the entire length of chromosomes to characterize the recombination landscape and uncover potential sex-biases in recombination.

A fine-scale genetic map of the Japanese population

Genetic maps are fundamental resources for linkage and association studies. A fine-scale genetic map can be constructed by inferring historical recombination events from the genome-wide structure of linkage disequilibrium-a non-random association of alleles among loci-by using population-scale sequencing data. We constructed a fine-scale genetic map and identified recombination hotspots from 10 092 551 bi-allelic high-quality autosomal markers segregating among 150 unrelated Japanese individuals whose genotypes were determined by high-coverage (30×) whole-genome sequencing, and the genotype quality was carefully controlled by using their parents' and offspring's genotypes. The pedigree information was also utilized for haplotype phasing. The resulting genome-wide recombination rate profiles were concordant with those of the worldwide population on a broad scale, and the resolution was much improved. We identified 9487 recombination hotspots and confirmed the enrichment of previously known motifs in the hotspots. Moreover, we demonstrated that the Japanese genetic map improved the haplotype phasing and genotype imputation accuracy for the Japanese population. The construction of a population-specific genetic map will help make genetics research more accurate.

Polymorphism Identification in the Coding Sequences (ORFs) of the Porcine Pregnancy-Associated Glycoprotein 2-like Gene Subfamily in Pigs

Pregnancy-associated glycoproteins (PAGs) are a polygenic family with many scattered genes and pseudogenes resulting from the duplication or fusion of a pseudogene with expression beginning in the trophoblast during the peri-implantation period and continuing in the trophectoderm. In this study, single-nucleotide polymorphism (SNP) and insertion/deletion (InDels) in the open reading frame (nine exons) of crossbreed pigs are reported for the first time. Novel SNPs/InDels were researched using genomic DNA templates isolated from the leukocytes of crossbreed pigs (N = 25), which were amplified, gel-out-purified, and sequenced. Sixteen SNPs and one InDel (g.6961_6966 Ins TGCCAA) were identified in the crossbreed pigs. In silico analysis revealed that among 16 SNPs, only 10 SNPs cause amino acid (aa) substitutions, and InDel codes asparagine (N298) and alanine (A299). The results provide a novel broad-based database (main pattern) that will be critical for future research into the possible correlations between the SNP genotypes of the pPAG2-L subfamily in pigs of various breeds whose reproductive traits are known.

Novel crossover and recombination hotspots massively spread across primate genomes

BACKGROUND

The recombination landscape and subsequent natural selection have vast consequences forevolution and speciation. However, most of the crossover and recombination hotspots are yet to be discovered. We previously reported the relevance of C and G trinucleotide two-repeat units (CG-TTUs) in crossovers and recombination.

METHODS

On a genome-wide scale, here we mapped all combinations of A and T trinucleotide two-repeat units (AT-TTUs) in human, consisting of AATAAT, ATAATA, ATTATT, TTATTA, TATTAT, and TAATAA. We also compared a number of the colonies formed by the AT-TTUs (distance between consecutive AT-TTUs < 500 bp) in several other primates and mouse.

RESULTS

We found that the majority of the AT-TTUs (> 96%) resided in approximately 1.4 million colonies, spread throughout the human genome. In comparison to the CG-TTU colonies, the AT-TTU colonies were significantly more abundant and larger in size. Pure units and overlapping units of the pure units were readily detectable in the same colonies, signifying that the units were the sites of unequal crossover. We discovered dynamic sharedness of several of the colonies across the primate species studied, which mainly reached maximum complexity and size in human.

CONCLUSIONS

We report novel crossover and recombination hotspots of the finest molecular resolution, massively spread and shared across the genomes of human and several other primates. With respect to crossovers and recombination, these genomes are far more dynamic than previously envisioned.

Understanding the Genetic Basis of Variation in Meiotic Recombination: Past, Present, and Future

Meiotic recombination is a fundamental feature of sexually reproducing species. It is often required for proper chromosome segregation and plays important role in adaptation and the maintenance of genetic diversity. The molecular mechanisms of recombination are remarkably conserved across eukaryotes, yet meiotic genes and proteins show substantial variation in their sequence and function, even between closely related species. Furthermore, the rate and distribution of recombination shows a huge diversity within and between chromosomes, individuals, sexes, populations, and species. This variation has implications for many molecular and evolutionary processes, yet how and why this diversity has evolved is not well understood. A key step in understanding trait evolution is to determine its genetic basis-that is, the number, effect sizes, and distribution of loci underpinning variation. In this perspective, I discuss past and current knowledge on the genetic basis of variation in recombination rate and distribution, explore its evolutionary implications, and present open questions for future research.

Effects of demographic history on recombination hotspots in soybean

Recombination is the primary mechanism underlying genetic improvement in populations and allows plant breeders to create new allelic combinations for agronomic improvement. Soybean [Glycine max (L.) Merr.] has gone through multiple genetic bottlenecks that have significantly affected its genetic diversity, linkage disequilibrium, and altered allele frequencies. To investigate the impact of genetic bottlenecks on recombination hotspots in soybeans, historical recombination was studied in three soybean populations. The populations were wild soybean [Glycine soja (Sieb. and Zucc.)], landraces, and North American elite soybean cultivars that have been genotyped with the SoySNP50K BeadChip. While each population after a genetic bottleneck had an increased average haplotype block size, they did not have a significant difference in the number of hotspots between each population. Instead, the increase in observed haplotype block size is likely due to an elimination of individuals that contained historical recombination at hotspots which decreased the observed rate of recombination for the hotspot after each genetic bottleneck. Conversely, heterochromatic DNA which has an increased haplotype block size compared to euchromatic DNA had a significantly different number of hotspots but not a significant difference in the average hotspot recombination rate. Previously identified genomic motifs associated with hotspots were also associated with hotspots found in the historical populations suggesting a common mechanism. This characterization of historical recombination hotspots in soybeans provides further insights into the effect genetic bottlenecks and selection have on recombination hotspots.

Recombination hotspots in Neurospora crassa controlled by idiomorphic sequences and meiotic silencing

Genes regulating recombination in specific chromosomal intervals of Neurospora crassa were described in the 1960s, but the mechanism is still unknown. For each of the rec-1, rec-2, and rec-3 genes, a single copy of the putative dominant allele, for example, rec-2SL found in St Lawrence OR74 A wild type, reduces recombination in chromosomal regions specific to that gene. However, when we sequenced the recessive allele, rec-2LG (derived from the Lindegren 1A wild type), we found that a 10 kb region in rec-2SL strains was replaced by a 2.7 kb unrelated sequence, making the "alleles" idiomorphs. When we introduced sad-1, a mutant lacking the RNA-dependent RNA polymerase that silences unpaired coding regions during meiosis into crosses heterozygous rec-2SL/rec-2LG, it increased recombination, indicating that meiotic silencing of a gene promoting recombination is responsible for dominant suppression of recombination. Consistent with this, mutation of rec-2LG by Repeat-Induced Point mutation generated an allele with multiple stop codons in the predicted rec-2 gene, which does not promote recombination and is recessive to rec-2LG. Sad-1 also relieves suppression of recombination in relevant target regions, in crosses heterozygous for rec-1 alleles and in crosses heterozygous for rec-3 alleles. We conclude that for all 3 known rec genes, 1 allele appears dominant only because meiotic silencing prevents the product of the active, "recessive," allele from stimulating recombination during meiosis. In addition, the proposed amino acid sequence of REC-2 suggests that regulation of recombination in Neurospora differs from any currently known mechanism.

SPO-Seq: An Accessible Method for Efficient Evaluation of Spo11 Catalytic Activity and Profiling Meiotic DSB Hotspots in Saccharomyces cerevisiae

Meiotic recombination is a key process facilitating the formation of crossovers and the exchange of genetic material between homologous chromosomes in early meiosis. This involves controlled double-strand breaks (DSBs) formation catalyzed by Spo11. DSBs exhibit a preferential location in specific genomic regions referred to as hotspots, and their variability is tied to varying Spo11 activity levels. We have refined a ChIP-Seq technique, called SPO-Seq, to map Spo11-specific DSB formation in Saccharomyces cerevisiae. The chapter describes our streamlined approach and the developed bioinformatic tools for processing data and comparing with existing DSB hotspot maps. Through this combined experimental and computational approach, we aim to enhance our understanding of meiotic recombination and genetic exchange processes in budding yeast, with the potential to expand this methodology to other organisms by applying a few modifications.

Meiotic DNA breaks drive multifaceted mutagenesis in the human germ line

Meiotic recombination commences with hundreds of programmed DNA breaks; however, the degree to which they are accurately repaired remains poorly understood. We report that meiotic break repair is eightfold more mutagenic for single-base substitutions than was previously understood, leading to de novo mutation in one in four sperm and one in 12 eggs. Its impact on indels and structural variants is even higher, with 100- to 1300-fold increases in rates per break. We uncovered new mutational signatures and footprints relative to break sites, which implicate unexpected biochemical processes and error-prone DNA repair mechanisms, including translesion synthesis and end joining in meiotic break repair. We provide evidence that these mechanisms drive mutagenesis in human germ lines and lead to disruption of hundreds of genes genome wide.

Cancer Associated PRDM9: Implications for Linking Genomic Instability and Meiotic Recombination

The PR domain-containing 9 or PRDM9 is a gene recognized for its fundamental role in meiosis, a process essential for forming reproductive cells. Recent findings have implicated alterations in the PRDM9, particularly its zinc finger motifs, in the onset and progression of cancer. This association is manifested through genomic instability and the misregulation of genes critical to cell growth, proliferation, and differentiation. In our comprehensive study, we harnessed advanced bioinformatic mining tools to delve deep into the intricate relationship between PRDM9F and cancer. We analyzed 136,752 breakpoints and found an undeniable association between specific PRDM9 motifs and the occurrence of double-strand breaks, a phenomenon evidenced in every cancer profile examined. Utilizing R statistical querying and the Regioner package, 55 unique sequence variations of PRDM9 were statistically correlated with cancer, from a pool of 1024 variations. A robust analysis using the Enrichr tool revealed prominent associations with various cancer types. Moreover, connections were noted with specific phenotypic conditions and molecular functions, underlining the pervasive influence of PRDM9 variations in the biological spectrum. The Reactome tool identified 25 significant pathways associated with cancer, offering insights into the mechanistic underpinnings linking PRDM9 to cancer progression. This detailed analysis not only confirms the pivotal role of PRDM9 in cancer development, but also unveils a complex network of biological processes influenced by its variations. The insights gained lay a solid foundation for future research aimed at deciphering the mechanistic pathways of PRDM9, offering prospects for targeted interventions and innovative therapeutic approaches in cancer management.

High population frequencies of MICA copy number variations originate from independent recombination events

MICA is a stress-induced ligand of the NKG2D receptor that stimulates NK and T cell responses and was identified as a key determinant of anti-tumor immunity. The MICA gene is located inside the MHC complex and is in strong linkage disequilibrium with HLA-B. While an HLA-B*48-linked MICA deletion-haplotype was previously described in Asian populations, little is known about other MICA copy number variations. Here, we report the genotyping of more than two million individuals revealing high frequencies of MICA duplications (1%) and MICA deletions (0.4%). Their prevalence differs between ethnic groups and can rise to 2.8% (Croatia) and 9.2% (Mexico), respectively. Targeted sequencing of more than 70 samples indicates that these copy number variations originate from independent nonallelic homologous recombination events between segmental duplications upstream of MICA and MICB. Overall, our data warrant further investigation of disease associations and consideration of MICA copy number data in oncological study protocols.

Take a walk on the KRAB side

Canonical Krüppel-associated box (KRAB)-containing zinc finger proteins (KZFPs) act as major repressors of transposable elements (TEs) via the KRAB-mediated recruitment of the heterochromatin scaffold KRAB-associated protein (KAP)1. KZFP genes emerged some 420 million years ago in the last common ancestor of coelacanth, lungfish, and tetrapods, and dramatically expanded to give rise to lineage-specific repertoires in contemporary species paralleling their TE load and turnover. However, the KRAB domain displays sequence and function variations that reveal repeated diversions from a linear TE-KZFP trajectory. This Review summarizes current knowledge on the evolution of KZFPs and discusses how ancestral noncanonical KZFPs endowed with variant KRAB, SCAN or DUF3669 domains have been utilized to achieve KAP1-independent functions.

The chromatin determinants and Ph1 gene effect at wheat sites with contrasting recombination frequency

INTRODUCTION

Meiotic recombination is one of the most important processes of evolution and adaptation to environmental conditions. Even though there is substantial knowledge about proteins involved in the process, targeting specific DNA loci by the recombination machinery is not well understood.

OBJECTIVES

This study aims to investigate a wheat recombination hotspot (H1) in comparison with a "regular" recombination site (Rec7) on the sequence and epigenetic level in conditions with functional and non-functional Ph1 locus.

METHODS

The DNA sequence, methylation pattern, and recombination frequency were analyzed for the H1 and Rec7 in three mapping populations derived by crossing introgressive wheat line 8.1 with cv. Chinese Spring (with Ph1 and ph1 alleles) and cv. Tähti.

RESULTS

The H1 and Rec7 loci are 1.586 kb and 2.538 kb long, respectively. High-density mapping allowed to delimit the Rec7 and H1 to 19 and 574 bp and 593 and 571 bp CO sites, respectively. A new method (ddPing) allowed screening recombination frequency in almost 66 thousand gametes. The screening revealed a 5.94-fold higher recombination frequency at the H1 compared to the Rec7. The H1 was also found out of the Ph1 control, similarly as gamete distortion. The recombination was strongly affected by larger genomic rearrangements but not by the SNP proximity. Moreover, chromatin markers for open chromatin and DNA hypomethylation were found associated with crossover occurrence except for the CHH methylation.

CONCLUSION

Our results, for the first time, allowed study of wheat recombination directly on sequence, shed new light on chromatin landmarks associated with particular recombination sites, and deepened knowledge about role of the Ph1 locus in control of wheat recombination processes. The results are suggesting more than one recombination control pathway. Understanding this phenomenon may become a base for more efficient wheat genome manipulation, gene pool enrichment, breeding, and study processes of recombination itself.

Single-cell multi-omics sequencing of human spermatogenesis reveals a DNA demethylation event associated with male meiotic recombination

Human spermatogenesis is a highly ordered process; however, the roles of DNA methylation and chromatin accessibility in this process remain largely unknown. Here by simultaneously investigating the chromatin accessibility, DNA methylome and transcriptome landscapes using the modified single-cell chromatin overall omic-scale landscape sequencing approach, we revealed that the transcriptional changes throughout human spermatogenesis were correlated with chromatin accessibility changes. In particular, we identified a set of transcription factors and cis elements with potential functions. A round of DNA demethylation was uncovered upon meiosis initiation in human spermatogenesis, which was associated with male meiotic recombination and conserved between human and mouse. Aberrant DNA hypermethylation could be detected in leptotene spermatocytes of certain nonobstructive azoospermia patients. Functionally, the intervention of DNA demethylation affected male meiotic recombination and fertility. Our work provides multi-omics landscapes of human spermatogenesis at single-cell resolution and offers insights into the association between DNA demethylation and male meiotic recombination.

The Recombination Hotspot Paradox: Co-evolution between PRDM9 and its target sites

Recombination often concentrates in small regions called recombination hotspots where recombination is much higher than the genome's average. In many vertebrates, including humans, gene PRDM9 specifies which DNA motifs will be the target for breaks that initiate recombination, ultimately determining the location of recombination hotspots. Because the sequence that breaks (allowing recombination) is converted into the sequence that does not break (preventing recombination), the latter sequence is over-transmitted to future generations and recombination hotspots are self-destructive. Given their self-destructive nature, recombination hotspots should eventually become extinct in genomes where they are found. While empirical evidence shows that individual hotspots do become inactive over time (die), hotspots are abundant in many vertebrates: a contradiction called the Recombination Hotspot Paradox. What saves recombination hotspots from their foretold extinction? Here we formulate a co-evolutionary model of the interaction among sequence-specific gene conversion, fertility selection, and recurrent mutation. We find that allelic frequencies oscillate leading to stable limit cycles. From a biological perspective this means that when fertility selection is weaker than gene conversion, it cannot stop individual hotspots from dying but can save them from extinction by driving their re-activation (resuscitation). In our model, mutation balances death and resuscitation of hotspots, thus maintaining their number over evolutionary time. Interestingly, we find that multiple alleles result in oscillations that are chaotic and multiple targets in oscillations that are asynchronous between targets thus helping to maintain the average genomic recombination probability constant. Furthermore, we find that the level of expression of PRDM9 should control for the fraction of targets that are hotspots and the overall temperature of the genome. Therefore, our co-evolutionary model improves our understanding of how hotspots may be replaced, thus contributing to solve the Recombination Hotspot Paradox. From a more applied perspective our work provides testable predictions regarding the relation between mutation probability and fertility selection with life expectancy of hotspots.

Conserved γδ T cell selection by BTNL proteins limits progression of human inflammatory bowel disease

Murine intraepithelial γδ T cells include distinct tissue-protective cells selected by epithelial butyrophilin-like (BTNL) heteromers. To determine whether this biology is conserved in humans, we characterized the colonic γδ T cell compartment, identifying a diverse repertoire that includes a phenotypically distinct subset coexpressing T cell receptor Vγ4 and the epithelium-binding integrin CD103. This subset was disproportionately diminished and dysregulated in inflammatory bowel disease, whereas on-treatment CD103+γδ T cell restoration was associated with sustained inflammatory bowel disease remission. Moreover, CD103+Vγ4+cell dysregulation and loss were also displayed by humans with germline BTNL3/BTNL8 hypomorphism, which we identified as a risk factor for penetrating Crohn's disease (CD). Thus, BTNL-dependent selection and/or maintenance of distinct tissue-intrinsic γδ T cells appears to be an evolutionarily conserved axis limiting the progression of a complex, multifactorial, tissue-damaging disease of increasing global incidence.

Chromosomal Instability in Genome Evolution: From Cancer to Macroevolution

The integrity of the genome is crucial for the survival of all living organisms. However, genomes need to adapt to survive certain pressures, and for this purpose use several mechanisms to diversify. Chromosomal instability (CIN) is one of the main mechanisms leading to the creation of genomic heterogeneity by altering the number of chromosomes and changing their structures. In this review, we will discuss the different chromosomal patterns and changes observed in speciation, in evolutional biology as well as during tumor progression. By nature, the human genome shows an induction of diversity during gametogenesis but as well during tumorigenesis that can conclude in drastic changes such as the whole genome doubling to more discrete changes as the complex chromosomal rearrangement chromothripsis. More importantly, changes observed during speciation are strikingly similar to the genomic evolution observed during tumor progression and resistance to therapy. The different origins of CIN will be treated as the importance of double-strand breaks (DSBs) or the consequences of micronuclei. We will also explain the mechanisms behind the controlled DSBs, and recombination of homologous chromosomes observed during meiosis, to explain how errors lead to similar patterns observed during tumorigenesis. Then, we will also list several diseases associated with CIN, resulting in fertility issues, miscarriage, rare genetic diseases, and cancer. Understanding better chromosomal instability as a whole is primordial for the understanding of mechanisms leading to tumor progression.

2022 William Allan Award

This article is based on the address given by the author at the 2022 meeting of The American Society of Human Genetics (ASHG) in Los Angeles, CA. The video of the original address can be found at the ASHG website.

Genetic Redundancy in Rye Shows in a Variety of Ways

Fifty years ago Susumu Ohno formulated the famous C-value paradox, which states that there is no correlation between the physical sizes of the genome, i.e., the amount of DNA, and the complexity of the organism, and highlighted the problem of genome redundancy. DNA that does not have a positive effect on the fitness of organisms has been characterized as "junk or selfish DNA". The controversial concept of junk DNA remains viable. Rye is a convenient subject for yet another test of the correctness and scientific significance of this concept. The genome of cultivated rye, Secale cereale L., is considered one of the largest among species of the tribe Triticeae and thus it tops the average angiosperm genome and the genomes of its closest evolutionary neighbors, such as species of barley, Hordeum (by approximately 30-35%), and diploid wheat species, Triticum (approximately 25%). The review provides an analysis of the structural organization of various regions of rye chromosomes with a description of the molecular mechanisms contributing to their size increase during evolution and the classes of DNA sequences involved in these processes. The history of the development of the concept of eukaryotic genome redundancy is traced and the current state of this problem is discussed.

Recombination and mutation shape variations in the major histocompatibility complex

The major histocompatibility complex (MHC) is closely associated with numerous diseases, but its high degree of polymorphism complicates the discovery of disease-associated variants. In principle, recombination and de novo mutations are two critical factors responsible for MHC polymorphisms. However, direct evidence for this hypothesis is lacking. Here, we report the generation of fine-scale MHC recombination and de novo mutation maps of ∼5 Mb by deep sequencing (> 100×) of the MHC genome for 17 MHC recombination and 30 non-recombination Han Chinese families (a total of 190 individuals). Recombination hotspots and Han-specific breakpoints are located in close proximity at haplotype block boundaries. The average MHC de novo mutation rate is higher than the genome-wide de novo mutation rate, particularly in MHC recombinant individuals. Notably, mutation and recombination generated polymorphisms are located within and outside linkage disequilibrium regions of the MHC, respectively, and evolution of the MHC locus was mainly controlled by positive selection. These findings provide insights on the evolutionary causes of the MHC diversity and may facilitate the identification of disease-associated genetic variants.

Further understanding of paternal uniparental disomy in Beckwith-Wiedemann syndrome

INTRODUCTION

Paternal uniparental disomy of chromosome 11 (upd(11)pat) accounts for up to 20% of molecularly confirmed Beckwith-Wiedemann spectrum (BWSp) cases. It belongs to the BWSp subgroup with the second highest tumor risk, and therefore needs particular awareness in research, diagnostics and clinical management.

AREAS COVERED

We overview the contribution of paternal (mosaic) uniparental disomy of chromosome 11 (UPD, upd(11)pat) and mosaic paternal uniparental diploidy in patients with Beckwith-Wiedemann features. The review comprises the current knowledge on their formation and their molecular and clinical consequences. Accordingly, the consequences for diagnostic testing and clinical monitoring are compiled.

EXPERT OPINION

The necessity to diagnostically identify and thus discriminate genome-wide paternal uniparental disomy, and upd(11)pat becomes obvious, due to the differences in the clinical course, disease prognosis, and treatment. In particular, monitoring of tumor development by liquid biopsy might be a promising option in the future. From the research point of view, it should be addressed why 11p is prone to mitotic recombination and thus also provide to the role of upd(11) as second hit in tumorigenesis.

Epigenetic targeting of transposon relics: beating the dead horses of the genome?

Transposable elements (TEs) have been seen as selfish genetic elements that can propagate in a host genome. Their propagation success is however hindered by a combination of mechanisms such as mutations, selection, and their epigenetic silencing by the host genome. As a result, most copies of TEs in a given genome are dead relics: their sequence is too degenerated to allow any transposition. Nevertheless, these TE relics often, but not always, remain epigenetically silenced, and if not to prevent transposition anymore, one can wonder the reason for this phenomenon. The mere self-perpetuating loop inherent to epigenetic silencing could alone explain that even when inactive, TE copies remain silenced. Beyond this process, nevertheless, antagonistic selective forces are likely to act on TE relic silencing. Especially, without the benefit of preventing transposition, TE relic silencing may prove deleterious to the host fitness, suggesting that the maintenance of TE relic silencing is the result of a fine, and perhaps case-by-case, evolutionary trade-off between beneficial and deleterious effects. Ultimately, the release of TE relics silencing may provide a 'safe' ground for adaptive epimutations to arise. In this review, we provide an overview of these questions in both plants and animals.

Variation in fine-scale recombination rate in temperature-evolved Drosophila melanogaster populations in response to selection

Meiotic recombination plays a critical evolutionary role in maintaining fitness in response to selective pressures due to changing environments. Variation in recombination rate has been observed amongst and between species and populations and within genomes across numerous taxa. Studies have demonstrated a link between changes in recombination rate and selection, but the extent to which fine-scale recombination rate varies between evolved populations during the evolutionary period in response to selection is under active research. Here, we utilize a set of 3 temperature-evolved Drosophila melanogaster populations that were shown to have diverged in several phenotypes, including recombination rate, based on the temperature regime in which they evolved. Using whole-genome sequencing data from these populations, we generated linkage disequilibrium-based fine-scale recombination maps for each population. With these maps, we compare recombination rates and patterns among the 3 populations and show that they have diverged at fine scales but are conserved at broader scales. We further demonstrate a correlation between recombination rates and genomic variation in the 3 populations. Lastly, we show variation in localized regions of enhanced recombination rates, termed warm spots, between the populations with these warm spots and associated genes overlapping areas previously shown to have diverged in the 3 populations due to selection. These data support the existence of recombination modifiers in these populations which are subject to selection during evolutionary change.

Analysis of archaic human haplotypes suggests that 5hmC acts as an epigenetic guide for NCO recombination

Background

Non-crossover (NCO) refers to a mechanism of homologous recombination in which short tracks of DNA are copied between homologue chromatids. The allelic changes are typically restricted to one or few SNPs, which potentially allow for the gradual adaptation and maturation of haplotypes. It is assumed to be a stochastic process but the analysis of archaic and modern human haplotypes revealed a striking variability in local NCO recombination rates.

Methods

NCO recombination rates of 1.9 million archaic SNPs shared with Denisovan hominids were defined by a linkage study and correlated with functional and genomic annotations as well as ChIP-Seq data from modern humans.

Results

We detected a strong correlation between NCO recombination rates and the function of the respective region: low NCO rates were evident in introns and quiescent intergenic regions but high rates in splice sites, exons, 5′- and 3′-UTRs, as well as CpG islands. Correlations with ChIP-Seq data from ENCODE and other public sources further identified epigenetic modifications that associated directly with these recombination events. A particularly strong association was observed for 5-hydroxymethylcytosine marks (5hmC), which were enriched in virtually all of the functional regions associated with elevated NCO rates, including CpG islands and ‘poised’ bivalent regions.

Conclusion

Our results suggest that 5hmC marks may guide the NCO machinery specifically towards functionally relevant regions and, as an intermediate of oxidative demethylation, may open a pathway for environmental influence by specifically targeting recently opened gene loci.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01353-9.

Chromatin accessibility shapes meiotic recombination in mouse primordial germ cells through assisting double-strand breaks and loop formation

Meiotic recombination is a driver of evolution, and aberrant recombination is a major contributor to aneuploidy in mammals. Mechanism of recombination remains elusive yet. Here, we present a computational analysis to explore recombination-related dynamics of chromatin accessibility in mouse primordial germ cells (PGCs). Our data reveals that: (1) recombination hotspots which get accessible at meiosis-specific DNase I-hypersensitive sites (DHSs) only when PGCs enter meiosis are located preferentially in intronic and distal intergenic regions; (2) stable DHSs maintained stably across PGC differentiation are enriched by CTCF motifs and CTCF binding and mediate chromatin loop formation; (3) compared with the specific DHSs aroused at meiotic stage, stable DHSs are largely encoded in DNA sequence and also enriched by epigenetic marks; (4) PRDM9 is likely to target nucleosome-occupied hotspot regions and remodels local chromatin structure to make them accessible for recombination machinery; and (5) cells undergoing meiotic recombination are deficient in TAD structure and chromatin loop arrays are organized regularly along the axis formed between homologous chromosomes. Taken together, by analyzing DHS-related DNA features, epigenetic marks and 3D genome structure, we revealed some specific roles of chromatin accessibility in recombination, which would expand our understanding of recombination mechanism.

Adaptive Control of the Meiotic Recombination Landscape by DNA Site-dependent Hotspots With Implications for Evolution

Meiosis is an essential component of the sexual life cycle in eukaryotes. The independent assortment of chromosomes in meiosis increases genetic diversity at the level of whole chromosomes and meiotic recombination increases genetic diversity within chromosomes. The resulting variability fuels evolution. Interestingly, global mapping of recombination in diverse taxa revealed dramatic changes in its frequency distribution between closely related species, subspecies, and even isolated populations of the same species. New insight into mechanisms for these evolutionarily rapid changes has come from analyses of environmentally induced plasticity of recombination in fission yeast. Many different DNA sites, and where identified their binding/activator proteins, control the positioning of recombination at hotspots. Each different class of hotspots functions as an independently controlled rheostat that modulates rates of recombination over a broad dynamic range in response to changing conditions. Together, this independent modulation can rapidly and dramatically alter the global frequency distribution of recombination. This process likely contributes substantially to (i.e., can largely explain) evolutionarily rapid, Prdm9-independent changes in the recombination landscape. Moreover, the precise control mechanisms allow cells to dynamically favor or disfavor newly arising combinations of linked alleles in response to changing extracellular and intracellular conditions, which has striking implications for the impacts of meiotic recombination on evolution.

Deep learning identifies and quantifies recombination hotspot determinants

MOTIVATION

Recombination is one of the essential genetic processes for sexually reproducing organisms, which can happen more frequently in some regions, called recombination hotspots. Although several factors, such as PRDM9 binding motifs, are known to be related to the hotspots, their contributions to the recombination hotspots have not been quantified, and other determinants are yet to be elucidated. Here, we propose a computational method, RHSNet, based on deep learning and signal processing, to identify and quantify the hotspot determinants in a purely data-driven manner, utilizing datasets from various studies, populations, sexes and species.

RESULTS

RHSNet can significantly outperform other sequence-based methods on multiple datasets across different species, sexes and studies. In addition to being able to identify hotspot regions and the well-known determinants accurately, more importantly, RHSNet can quantify the determinants that contribute significantly to the recombination hotspot formation in the relation between PRDM9 binding motif, histone modification and GC content. Further cross-sex, cross-population and cross-species studies suggest that the proposed method has the generalization power and potential to identify and quantify the evolutionary determinant motifs.

AVAILABILITY AND IMPLEMENTATION

https://github.com/frankchen121212/RHSNet.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

On the Base Composition of Transposable Elements

Transposable elements exhibit a base composition that is often different from the genomic average and from hosts’ genes. The most common compositional bias is towards Adenosine and Thymine, although this bias is not universal, and elements with drastically different base composition can coexist within the same genome. The AT-richness of transposable elements is apparently maladaptive because it results in poor transcription and sub-optimal translation of proteins encoded by the elements. The cause(s) of this unusual base composition remain unclear and have yet to be investigated. Here, I review what is known about the nucleotide content of transposable elements and how this content can affect the genome of their host as well as their own replication. The compositional bias of transposable elements could result from several non-exclusive processes including horizontal transfer, mutational bias, and selection. It appears that mutation alone cannot explain the high AT-content of transposons and that selection plays a major role in the evolution of the compositional bias. The reason why selection would favor a maladaptive nucleotide content remains however unexplained and is an area of investigation that clearly deserves attention.

High-speed rail model reveals the gene tandem amplification mediated by short repeated sequence in eukaryote

The occurrence of gene duplication/amplification (GDA) provide potential material for adaptive evolution with environmental stress. Several molecular models have been proposed to explain GDA, recombination via short stretches of sequence similarity plays a crucial role. By screening genomes for such events, we propose a “SRS (short repeated sequence) *N + unit + SRS*N” amplified unit under USCE (unequal sister-chromatid exchange) for tandem amplification mediated by SRS with different repeat numbers in eukaryotes. The amplified units identified from 2131 well-organized amplification events that generate multi gene/element copy amplified with subsequent adaptive evolution in the respective species. Genomic data we analyzed showed dynamic changes among related species or subspecies or plants from different ecotypes/strains. This study clarifies the characteristics of variable copy number SRS on both sides of amplified unit under USCE mechanism, to explain well-organized gene tandem amplification under environmental stress mediated by SRS in all eukaryotes.

PRDM9-directed recombination hotspots depleted near meiotically transcribed genes

PRDM9 drives recombination hotspots in some mammals, including mice and apes. Non-functional orthologs of PRDM9 are present in a wide variety of vertebrates, but why it is functionally maintained in some lineages is not clear. One possible explanation is that PRDM9 plays a role in ensuring that meiosis is successful. During meiosis, available DNA may be a limiting resource given the tight packaging of chromosomes and could lead to competition between two key processes: meiotic transcription and recombination. Here we explore this potential competition and the role that PRDM9 might play in their interaction. Leveraging existing mouse genomic data, we use resampling schemes that simulate shuffled features along the genome and models that account for the rarity of features in the genome, to test if PRDM9 influences interactions between recombination hotspots and meiotic transcription in a whole genome framework. We also explored possible DNA sequence motifs associated to clusters of hotspots not tied to transcription or PRDM9. We find evidence of competition between meiotic transcription and recombination, with PRDM9 appearing to relocate recombination to avoid said conflict. We also find that retrotransposons may be playing a role in directing hotspots in the absence of other factors.

Mutational bias in spermatogonia impacts the anatomy of regulatory sites in the human genome

Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9 binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF1, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (≥5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

A Chinese pedigree with glucocorticoid remediable aldosteronism

Glucocorticoid-remediable aldosteronism (GRA) is an autosomal-dominant inherited aldosteronism that is often accompanied by early-onset hypertension. GRA is caused by the unequal crossover of the 11β-hydroxylase (CYP11B1) and aldosterone synthase (CYP11B2) genes. As a result of chimeric gene duplication, aldosterone is ectopically synthesized in the adrenal zona fasciculata under the control of adrenocorticotropic hormone (ACTH). Here, we describe a Chinese pedigree with three affected subjects. Both the uncle and nephew were hospitalized in our hospital due to early-onset hypertension (onset <20 years old) and were diagnosed with primary aldosteronism (PA). Their laboratory test results revealed hyperaldosteronism, hyporeninemia, a high plasma aldosterone to renin (ARR) ratio, and normal serum potassium (K⁺). Captopril failed to suppress aldosterone secretion. This family had a strong paternal history of hypertension. Thirteen members underwent gene testing, and three of them were found to be GRA positive. Through long-extension PCR (XL-PCR) and direct sequencing, we identified the CYP11B1/CYP11B2 chimeric gene, and with unequal crossover breakpoints located between intron 2 of CYP11B1 and exon 3 of CYP11B2 in the three patients. Low-dose dexamethasone was effective. This is the first family report of GRA in northern China. Moreover, a case of GRA combined with a CACNA1H gene mutation is reported for the first time. We found that dihydropyridine calcium channel blockers (CCBs) combined with aldosterone receptor antagonists exerted good therapeutic effects in controlling blood pressure in GRA patients for whom glucocorticoid therapy was not an option.

Natural variation identifies SNI1, the SMC5/6 component, as a modifier of meiotic crossover in Arabidopsis

Meiotic recombination plays a fundamental role in shaping genetic diversity in eukaryotes. Extensive variation in crossover rate exists between populations and species. The identity of modifier loci and their roles in genome evolution remain incompletely understood. We explored natural variation in Arabidopsis crossover and identified SNI1 as the causal gene underlying a major modifier locus. To date, SNI1 had no known role in crossover. SNI1 is a component of the SMC5/6 complex that is closely related to cohesin and condensin. Arabidopsis sni1 and other SMC5/6 mutants show similar effects on the interference-independent crossover pathway. Hence, our findings demonstrate that the SMC5/6 complex, which is known for its role in DNA damage repair, is also important for control of meiotic crossover.

The frequency and distribution of meiotic crossovers are tightly controlled; however, variation in this process can be observed both within and between species. Using crosses of two natural Arabidopsis thaliana accessions, Col and Ler, we mapped a crossover modifier locus to semidominant polymorphisms in SUPPRESSOR OF NPR1-1 INDUCIBLE 1 (SNI1), which encodes a component of the SMC5/6 complex. The sni1 mutant exhibits a modified pattern of recombination across the genome with crossovers elevated in chromosome distal regions but reduced in pericentromeres. Mutations in SNI1 result in reduced crossover interference and can partially restore the fertility of a Class I crossover pathway mutant, which suggests that the protein affects noninterfering crossover repair. Therefore, we tested genetic interactions between SNI1 and both RECQ4 and FANCM DNA helicases, which showed that additional Class II crossovers observed in the sni1 mutant are FANCM independent. Furthermore, genetic analysis of other SMC5/6 mutants confirms the observations of crossover redistribution made for SNI1. The study reveals the importance of the SMC5/6 complex in ensuring the proper progress of meiotic recombination in plants.

Epigenetic Marks and Variation of Sequence-Based Information Along Genomic Regions Are Predictive of Recombination Hot/Cold Spots in Saccharomyces cerevisiae

Characterization and identification of recombination hotspots provide important insights into the mechanism of recombination and genome evolution. In contrast with existing sequence-based models for predicting recombination hotspots which were defined in a ORF-based manner, here, we first defined recombination hot/cold spots based on public high-resolution Spo11-oligo-seq data, then characterized them in terms of DNA sequence and epigenetic marks, and finally presented classifiers to identify hotspots. We found that, in addition to some previously discovered DNA-based features like GC-skew, recombination hotspots in yeast can also be characterized by some remarkable features associated with DNA physical properties and shape. More importantly, by using DNA-based features and several epigenetic marks, we built several classifiers to discriminate hotspots from coldspots, and found that SVM classifier performs the best with an accuracy of ∼92%, which is also the highest among the models in comparison. Feature importance analysis combined with prediction results show that epigenetic marks and variation of sequence-based features along the hotspots contribute dominantly to hotspot identification. By using incremental feature selection method, an optimal feature subset that consists of much less features was obtained without sacrificing prediction accuracy.

Genetic variation in recombination rate in the pig

Background

Meiotic recombination results in the exchange of genetic material between homologous chromosomes. Recombination rate varies between different parts of the genome, between individuals, and is influenced by genetics. In this paper, we assessed the genetic variation in recombination rate along the genome and between individuals in the pig using multilocus iterative peeling on 150,000 individuals across nine genotyped pedigrees. We used these data to estimate the heritability of recombination and perform a genome-wide association study of recombination in the pig.

Results

Our results confirmed known features of the recombination landscape of the pig genome, including differences in genetic length of chromosomes and marked sex differences. The recombination landscape was repeatable between lines, but at the same time, there were differences in average autosome-wide recombination rate between lines. The heritability of autosome-wide recombination rate was low but not zero (on average 0.07 for females and 0.05 for males). We found six genomic regions that are associated with recombination rate, among which five harbour known candidate genes involved in recombination: RNF212, SHOC1, SYCP2, MSH4 and HFM1.

Conclusions

Our results on the variation in recombination rate in the pig genome agree with those reported for other vertebrates, with a low but nonzero heritability, and the identification of a major quantitative trait locus for recombination rate that is homologous to that detected in several other species. This work also highlights the utility of using large-scale livestock data to understand biological processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00643-0.

Identification and characterization of inheritable structural variations induced by ion beam radiations in rice

High energy ion beams are effective physical mutagens for mutation induction in plants. Due to their high linear energy transfer (LET) property, they are known to generate single nucleotide variations (SNVs) and insertion/deletions (InDels, <50 bp) as well as structural variations (SVs). However, due to the technical difficulties to identify SVs, studies on ion beam induced SVs by genome sequencing have so far been limited in numbers and inadequate in nature, and knowledge of SVs is scarce with regards to their characteristics. In the present study, we identified and validated SVs in six M₄ plants (designated as Ar_50, Ar_100, C_150, C_200, Ne_50 and Ne_100 according to ion beam types and irradiation doses), two each induced by argon (⁴⁰Ar¹⁸⁺), carbon (¹²C⁶⁺) and neon (²⁰Ne¹⁰⁺) ion beams and performed in depth analyses of their characteristics. In total, 22 SVs were identified and validated, consisting of 11 deletions, 1 duplication, and 4 intra-chromosomal and 6 inter-chromosomal translocations. There were several SVs larger than 1 kbp. The SVs were distributed across the whole genome with an aggregation with SNVs and InDels only in the Ne_50 mutants. An enrichment of a 11-bp wide G-rich DNA motif 'GAAGGWGGRGG' was identified around the SV breakpoints. Three mechanisms might be involved in the SV formation, i.e., the expansion of tandem repeats, transposable element insertion, and non-allelic homologous recombination. Put together, the present study provides a preliminary view of SVs induced by Ar, C and Ne ion beam radiations, and as a pilot study, it contributes to our understanding of how SVs might form after ion beam irradiation in rice.

The Structural, Functional and Evolutionary Impact of Transposable Elements in Eukaryotes

Transposable elements (TEs) are nearly ubiquitous in eukaryotes. The increase in genomic data, as well as progress in genome annotation and molecular biology techniques, have revealed the vast number of ways mobile elements have impacted the evolution of eukaryotes. In addition to being the main cause of difference in haploid genome size, TEs have affected the overall organization of genomes by accumulating preferentially in some genomic regions, by causing structural rearrangements or by modifying the recombination rate. Although the vast majority of insertions is neutral or deleterious, TEs have been an important source of evolutionary novelties and have played a determinant role in the evolution of fundamental biological processes. TEs have been recruited in the regulation of host genes and are implicated in the evolution of regulatory networks. They have also served as a source of protein-coding sequences or even entire genes. The impact of TEs on eukaryotic evolution is only now being fully appreciated and the role they may play in a number of biological processes, such as speciation and adaptation, remains to be deciphered.

Human L1 Transposition Dynamics Unraveled with Functional Data Analysis

Long INterspersed Elements-1 (L1s) constitute >17% of the human genome and still actively transpose in it. Characterizing L1 transposition across the genome is critical for understanding genome evolution and somatic mutations. However, to date, L1 insertion and fixation patterns have not been studied comprehensively. To fill this gap, we investigated three genome-wide data sets of L1s that integrated at different evolutionary times: 17,037 de novo L1s (from an L1 insertion cell-line experiment conducted in-house), and 1,212 polymorphic and 1,205 human-specific L1s (from public databases). We characterized 49 genomic features-proxying chromatin accessibility, transcriptional activity, replication, recombination, etc.-in the ±50 kb flanks of these elements. These features were contrasted between the three L1 data sets and L1-free regions using state-of-the-art Functional Data Analysis statistical methods, which treat high-resolution data as mathematical functions. Our results indicate that de novo, polymorphic, and human-specific L1s are surrounded by different genomic features acting at specific locations and scales. This led to an integrative model of L1 transposition, according to which L1s preferentially integrate into open-chromatin regions enriched in non-B DNA motifs, whereas they are fixed in regions largely free of purifying selection-depleted of genes and noncoding most conserved elements. Intriguingly, our results suggest that L1 insertions modify local genomic landscape by extending CpG methylation and increasing mononucleotide microsatellite density. Altogether, our findings substantially facilitate understanding of L1 integration and fixation preferences, pave the way for uncovering their role in aging and cancer, and inform their use as mutagenesis tools in genetic studies.

MendelVar: gene prioritization at GWAS loci using phenotypic enrichment of Mendelian disease genes

Motivation

Gene prioritization at human GWAS loci is challenging due to linkage-disequilibrium and long-range gene regulatory mechanisms. However, identifying the causal gene is crucial to enable identification of potential drug targets and better understanding of molecular mechanisms. Mapping GWAS traits to known phenotypically relevant Mendelian disease genes near a locus is a promising approach to gene prioritization.

Results

We present MendelVar, a comprehensive tool that integrates knowledge from four databases on Mendelian disease genes with enrichment testing for a range of associated functional annotations such as Human Phenotype Ontology, Disease Ontology and variants from ClinVar. This open web-based platform enables users to strengthen the case for causal importance of phenotypically matched candidate genes at GWAS loci. We demonstrate the use of MendelVar in post-GWAS gene annotation for type 1 diabetes, type 2 diabetes, blood lipids and atopic dermatitis.

Availability and implementation

MendelVar is freely available at https://mendelvar.mrcieu.ac.uk

Supplementary information

Supplementary data are available at Bioinformatics online.

The Transposable Element Environment of Human Genes Differs According to Their Duplication Status and Essentiality

Transposable elements (TEs) are major components of eukaryotic genomes and represent approximately 45% of the human genome. TEs can be important sources of novelty in genomes and there is increasing evidence that TEs contribute to the evolution of gene regulation in mammals. Gene duplication is an evolutionary mechanism that also provides new genetic material and opportunities to acquire new functions. To investigate how duplicated genes are maintained in genomes, here, we explored the TE environment of duplicated and singleton genes. We found that singleton genes have more short-interspersed nuclear elements and DNA transposons in their vicinity than duplicated genes, whereas long-interspersed nuclear elements and long-terminal repeat retrotransposons have accumulated more near duplicated genes. We also discovered that this result is highly associated with the degree of essentiality of the genes with an unexpected accumulation of short-interspersed nuclear elements and DNA transposons around the more-essential genes. Our results underline the importance of taking into account the TE environment of genes to better understand how duplicated genes are maintained in genomes.

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

Length Polymorphism and Methylation Status of UPS29 Minisatellite of the ACAP3 Gene as Molecular Biomarker of Epilepsy. Sex Differences in Seizure Types and Symptoms

Epilepsy is a neurological disease with different clinical forms and inter-individuals heterogeneity, which may be associated with genetic and/or epigenetic polymorphisms of tandem-repeated noncoding DNA. These polymorphisms may serve as predictive biomarkers of various forms of epilepsy. ACAP3 is the protein regulating morphogenesis of neurons and neuronal migration and is an integral component of important signaling pathways. This study aimed to carry out an association analysis of the length polymorphism and DNA methylation of the UPS29 minisatellite of the ACAP3 gene in patients with epilepsy. We revealed an association of short UPS29 alleles with increased risk of development of symptomatic and cryptogenic epilepsy in women, and also with cerebrovascular pathologies, structural changes in the brain, neurological status, and the clinical pattern of seizures in both women and men. The increase of frequency of hypomethylated UPS29 alleles in men with symptomatic epilepsy, and in women with both symptomatic and cryptogenic epilepsy was observed. For patients with hypomethylated UPS29 alleles, we also observed structural changes in the brain, neurological status, and the clinical pattern of seizures. These associations had sex-specific nature similar to a genetic association. In contrast with length polymorphism epigenetic changes affected predominantly the long UPS29 allele. We suppose that genetic and epigenetic alterations UPS29 can modify ACAP3 expression and thereby affect the development and clinical course of epilepsy.

Extreme enrichment of VNTR-associated polymorphicity in human subtelomeres: genes with most VNTRs are predominantly expressed in the brain

The human genome harbors numerous structural variants (SVs) which, due to their repetitive nature, are currently underexplored in short-read whole-genome sequencing approaches. Using single-molecule, real-time (SMRT) long-read sequencing technology in combination with FALCON-Unzip, we generated a de novo assembly of the diploid genome of a 115-year-old Dutch cognitively healthy woman. We combined this assembly with two previously published haploid assemblies (CHM1 and CHM13) and the GRCh38 reference genome to create a compendium of SVs that occur across five independent human haplotypes using the graph-based multi-genome aligner REVEAL. Across these five haplotypes, we detected 31,680 euchromatic SVs (>50 bp). Of these, ~62% were comprised of repetitive sequences with 'variable number tandem repeats' (VNTRs), ~10% were mobile elements (Alu, L1, and SVA), while the remaining variants were inversions and indels. We observed that VNTRs with GC-content >60% and repeat patterns longer than 15 bp were 21-fold enriched in the subtelomeric regions (within 5 Mb of the ends of chromosome arms). VNTR lengths can expand to exceed a critical length which is associated with impaired gene transcription. The genes that contained most VNTRs, of which PTPRN2 and DLGAP2 are the most prominent examples, were found to be predominantly expressed in the brain and associated with a wide variety of neurological disorders. Repeat-induced variation represents a sizeable fraction of the genetic variation in human genomes and should be included in investigations of genetic factors associated with phenotypic traits, specifically those associated with neurological disorders. We make available the long and short-read sequence data of the supercentenarian genome, and a compendium of SVs as identified across 5 human haplotypes.

Disentangling the determinants of transposable elements dynamics in vertebrate genomes using empirical evidences and simulations

The interactions between transposable elements (TEs) and their hosts constitute one of the most profound co-evolutionary processes found in nature. The population dynamics of TEs depends on factors specific to each TE families, such as the rate of transposition and insertional preference, the demographic history of the host and the genomic landscape. How these factors interact has yet to be investigated holistically. Here we are addressing this question in the green anole (Anolis carolinensis) whose genome contains an extraordinary diversity of TEs (including non-LTR retrotransposons, SINEs, LTR-retrotransposons and DNA transposons). We observed a positive correlation between recombination rate and frequency of TEs and densities for LINEs, SINEs and DNA transposons. For these elements, there was a clear impact of demography on TE frequency and abundance, with a loss of polymorphic elements and skewed frequency spectra in recently expanded populations. On the other hand, some LTR-retrotransposons displayed patterns consistent with a very recent phase of intense amplification. To determine how demography, genomic features and intrinsic properties of TEs interact we ran simulations using SLiM3. We determined that i) short TE insertions are not strongly counter-selected, but long ones are, ii) neutral demographic processes, linked selection and preferential insertion may explain positive correlations between average TE frequency and recombination, iii) TE insertions are unlikely to have been massively recruited in recent adaptation. We demonstrate that deterministic and stochastic processes have different effects on categories of TEs and that a combination of empirical analyses and simulations can disentangle these mechanisms.

MRLR: unraveling high-resolution meiotic recombination by linked reads

MOTIVATION

Meiotic recombination facilitates the transmission of exchanged genetic material between homologous chromosomes and plays a crucial role in increasing the genetic variations in eukaryotic organisms. In humans, thousands of crossover events have been identified by genotyping related family members. However, most of these crossover regions span tens to hundreds of kb, which is not sufficient resolution to accurately identify the crossover breakpoints in a typical trio family.

RESULTS

We have developed MRLR, a software using 10X linked reads to identify crossover events at a high resolution. By reconstructing the gamete genome, MRLR only requires a trio family dataset and can efficiently discover the crossover events. Using MRLR, we revealed a fine-scale pattern of crossover regions in six human families. From the two closest heterozygous alleles around the crossovers, we determined that MRLR achieved a median resolution 4.5 kb. This method can delineate a genome-wide landscape of crossover events at a precise scale, which is important for both functional and genomic features analysis of meiotic recombination.

AVAILABILITY AND IMPLEMENTATION

MRLR is freely available at https://github.com/ChongLab/MRLR, implemented in Perl.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Insights into variation in meiosis from 31,228 human sperm genomes

Meiosis, although essential for reproduction, is also variable and error-prone: rates of chromosome crossover vary among gametes, between the sexes, and among humans of the same sex, and chromosome missegregation leads to abnormal chromosome numbers (aneuploidy)1-8. To study diverse meiotic outcomes and how they covary across chromosomes, gametes and humans, we developed Sperm-seq, a way of simultaneously analysing the genomes of thousands of individual sperm. Here we analyse the genomes of 31,228 human gametes from 20 sperm donors, identifying 813,122 crossovers and 787 aneuploid chromosomes. Sperm donors had aneuploidy rates ranging from 0.01 to 0.05 aneuploidies per gamete; crossovers partially protected chromosomes from nondisjunction at the meiosis I cell division. Some chromosomes and donors underwent more-frequent nondisjunction during meiosis I, and others showed more meiosis II segregation failures. Sperm genomes also manifested many genomic anomalies that could not be explained by simple nondisjunction. Diverse recombination phenotypes-from crossover rates to crossover location and separation, a measure of crossover interference-covaried strongly across individuals and cells. Our results can be incorporated with earlier observations into a unified model in which a core mechanism, the variable physical compaction of meiotic chromosomes, generates interindividual and cell-to-cell variation in diverse meiotic phenotypes.

An equivariant Bayesian convolutional network predicts recombination hotspots and accurately resolves binding motifs

MOTIVATION

Convolutional neural networks (CNNs) have been tremendously successful in many contexts, particularly where training data are abundant and signal-to-noise ratios are large. However, when predicting noisily observed phenotypes from DNA sequence, each training instance is only weakly informative, and the amount of training data is often fundamentally limited, emphasizing the need for methods that make optimal use of training data and any structure inherent in the process.

RESULTS

Here we show how to combine equivariant networks, a general mathematical framework for handling exact symmetries in CNNs, with Bayesian dropout, a version of Monte Carlo dropout suggested by a reinterpretation of dropout as a variational Bayesian approximation, to develop a model that exhibits exact reverse-complement symmetry and is more resistant to overtraining. We find that this model combines improved prediction consistency with better predictive accuracy compared to standard CNN implementations and state-of-art motif finders. We use our network to predict recombination hotspots from sequence, and identify binding motifs for the recombination-initiation protein PRDM9 previously unobserved in this data, which were recently validated by high-resolution assays. The network achieves a predictive accuracy comparable to that attainable by a direct assay of the H3K4me3 histone mark, a proxy for PRDM9 binding.

AVAILABILITY AND IMPLEMENTATION

https://github.com/luntergroup/EquivariantNetworks.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Determining the impact of uncharacterized inversions in the human genome by droplet digital PCR

Despite the interest in characterizing genomic variation, the presence of large repeats at the breakpoints hinders the analysis of many structural variants. This is especially problematic for inversions, since there is typically no gain or loss of DNA. Here, we tested novel linkage-based droplet digital PCR (ddPCR) assays to study 20 inversions ranging from 3.1 to 742 kb flanked by inverted repeats (IRs) up to 134 kb long. Of those, we validated 13 inversions predicted by different genome-wide techniques. In addition, we obtained new experimental human population information across 95 African, European, and East Asian individuals for 16 inversions, including four already validated variants without high-throughput genotyping methods. Through comparison with previous data, independent replicates and both inversion breakpoints, we demonstrate that the technique is highly accurate and reproducible. Most studied inversions are widespread across continents, and their frequency is negatively correlated with genetic length. Moreover, all except two show clear signs of being recurrent, and we could better define the factors affecting recurrence levels and estimate the inversion rate across the genome. Finally, the generated genotypes have allowed us to check inversion functional effects, validating gene expression differences reported before for two inversions and finding new candidate associations. Therefore, the developed methodology makes it possible to screen these and other complex genomic variants quickly in a large number of samples for the first time, highlighting the importance of direct genotyping to assess their potential consequences and clinical implications.