Gene conversion: mechanisms, evolution and human disease

Evolutionary variation in gene conversion at the avian MHC is explained by fluctuating selection, gene copy numbers and life history

The major histocompatibility complex (MHC) multigene family encodes key pathogen-recognition molecules of the vertebrate adaptive immune system. Hyper-polymorphism of MHC genes is de novo generated by point mutations, but new haplotypes may also arise by re-shuffling of existing variation through intra- and inter-locus gene conversion. Although the occurrence of gene conversion at the MHC has been known for decades, we still have limited understanding of its functional importance. Here, I took advantage of extensive genetic resources (~9000 sequences) to investigate broad scale macroevolutionary patterns in gene conversion processes at the MHC across nearly 200 avian species. Gene conversion was found to constitute a universal mechanism in birds, as 83% of species showed footprints of gene conversion at either MHC class and 25% of all allelic variants were attributed to gene conversion. Gene conversion processes were stronger at MHC-II than MHC-I, but inter-specific variation at both MHC classes was explained by similar evolutionary scenarios, reflecting fluctuating selection towards different optima and drift. Gene conversion showed uneven phylogenetic distribution across birds and was driven by gene copy number variation, supporting significant role of inter-locus gene conversion processes in the evolution of the avian MHC. Finally, MHC gene conversion was stronger in species with fast life histories (high fecundity) and in long-distance migrants, likely reflecting variation in population sizes and host-pathogen coevolutionary dynamics. The results provide a robust comparative framework for understanding macroevolutionary variation in gene conversion at the avian MHC and reinforce important contribution of this mechanism to functional MHC diversity.

Selected Lark Mitochondrial Genomes Provide Insights into the Evolution of Second Control Region with Tandem Repeats in Alaudidae (Aves, Passeriformes)

The control region (CR) regulates the replication and transcription of the mitochondrial genome (mitogenome). Some avian mitogenomes possess two CRs, and the second control region (CR2) may enhance replication and transcription; however, the CR2 in lark mitogenome appears to be undergoing loss and is accompanied by tandem repeats. Here, we characterized six lark mitogenomes from Alaudala cheleensis, Eremophila alpestris, Alauda razae, and Calandrella cinerea and reconstructed the phylogeny of Passerida. Through further comparative analysis among larks, we traced the evolutionary process of CR2. The mitochondrial gene orders were conserved in all published lark mitogenomes, with Cytb-trnT-CR1-trnP-ND6-trnE-remnant CR2 with tandem repeat-trnF-rrnS. Phylogenetic analysis revealed Alaudidae and Panuridae are sister groups at the base of Sylvioidea, and sporadic losses of CR2 may occur in their common ancestor. CR sequence and phylogeny analysis indicated CR2 tandem repeats were generated within CR2, originating in the ancestor of all larks, rather than inherited from CR1. The secondary structure comparison of tandem repeat units within and between species suggested slipped-strand mispairing and DNA turnover as suitable models for explaining the origin and evolution of these repeats. This study reveals the evolutionary process of the CR2 containing tandem repeat in Alaudidae, providing reference for understanding the evolutionary characteristics and dynamics of tandem repeats.

Low-input PacBio sequencing generates high-quality individual fly genomes and characterizes mutational processes

Long-read sequencing, exemplified by PacBio, revolutionizes genomics, overcoming challenges like repetitive sequences. However, the high DNA requirement ( > 1 µg) is prohibitive for small organisms. We develop a low-input (100 ng), low-cost, and amplification-free library-generation method for PacBio sequencing (LILAP) using Tn5-based tagmentation and DNA circularization within one tube. We test LILAP with two Drosophila melanogaster individuals, and generate near-complete genomes, surpassing preexisting single-fly genomes. By analyzing variations in these two genomes, we characterize mutational processes: complex transpositions (transposon insertions together with extra duplications and/or deletions) prefer regions characterized by non-B DNA structures, and gene conversion of transposons occurs on both DNA and RNA levels. Concurrently, we generate two complete assemblies for the endosymbiotic bacterium Wolbachia in these flies and similarly detect transposon conversion. Thus, LILAP promises a broad PacBio sequencing adoption for not only mutational studies of flies and their symbionts but also explorations of other small organisms or precious samples.

Enigmatic Nodal and Lefty gene repertoire discrepancy: Latent evolutionary history revealed by vertebrate-wide phylogeny

Homology in vertebrate body plans is traditionally ascribed to the high-level conservation of regulatory components within the genetic programs governing them, particularly during the "phylotypic stage." However, advancements in embryology and molecular phylogeny have unveiled the dynamic nature of gene repertoires responsible for early development. Notably, the Nodal and Lefty genes, members of the transforming growth factor-beta superfamily producing intercellular signaling molecules and crucial for left-right (L-R) symmetry breaking, exhibit distinctive features within their gene repertoires. These features encompass among-species gene repertoire variations resulting from gene gain and loss, as well as gene conversion. Despite their significance, these features have been largely unexplored in a phylogenetic context, but accumulating genome-wide sequence information is allowing the scrutiny of these features. It has exposed hidden paralogy between Nodal1 and Nodal2 genes resulting from differential gene loss in amniotes. In parallel, the tandem cluster of Lefty1 and Lefty2 genes, which was thought to be confined to mammals, is observed in sharks and rays, with an unexpected phylogenetic pattern. This article provides a comprehensive review of the current understanding of the origins of these vertebrate gene repertoires and proposes a revised nomenclature based on the elucidated history of vertebrate genome evolution.

A Unifying Hypothesis for the Genome Dynamics Proposed to Underlie Neuropsychiatric Phenotypes

The sheer number of gene variants and the extent of the observed clinical and molecular heterogeneity recorded in neuropsychiatric disorders (NPDs) could be due to the magnified downstream effects initiated by a smaller group of genomic higher-order alterations in response to endogenous or environmental stress. Chromosomal common fragile sites (CFS) are functionally linked with microRNAs, gene copy number variants (CNVs), sub-microscopic deletions and duplications of DNA, rare single-nucleotide variants (SNVs/SNPs), and small insertions/deletions (indels), as well as chromosomal translocations, gene duplications, altered methylation, microRNA and L1 transposon activity, and 3-D chromosomal topology characteristics. These genomic structural features have been linked with various NPDs in mostly isolated reports and have usually only been viewed as areas harboring potential candidate genes of interest. The suggestion to use a higher level entry point (the 'fragilome' and associated features) activated by a central mechanism ('stress') for studying NPD genetics has the potential to unify the existing vast number of different observations in this field. This approach may explain the continuum of gene findings distributed between affected and unaffected individuals, the clustering of NPD phenotypes and overlapping comorbidities, the extensive clinical and molecular heterogeneity, and the association with certain other medical disorders.

Shared evolutionary processes shape landscapes of genomic variation in the great apes

For at least the past 5 decades, population genetics, as a field, has worked to describe the precise balance of forces that shape patterns of variation in genomes. The problem is challenging because modeling the interactions between evolutionary processes is difficult, and different processes can impact genetic variation in similar ways. In this paper, we describe how diversity and divergence between closely related species change with time, using correlations between landscapes of genetic variation as a tool to understand the interplay between evolutionary processes. We find strong correlations between landscapes of diversity and divergence in a well-sampled set of great ape genomes, and explore how various processes such as incomplete lineage sorting, mutation rate variation, GC-biased gene conversion and selection contribute to these correlations. Through highly realistic, chromosome-scale, forward-in-time simulations, we show that the landscapes of diversity and divergence in the great apes are too well correlated to be explained via strictly neutral processes alone. Our best fitting simulation includes both deleterious and beneficial mutations in functional portions of the genome, in which 9% of fixations within those regions is driven by positive selection. This study provides a framework for modeling genetic variation in closely related species, an approach which can shed light on the complex balance of forces that have shaped genetic variation.

Role for gene conversion in the evolution of cell-surface antigens of the malaria parasite Plasmodium falciparum

While the malaria parasite Plasmodium falciparum has low average genome-wide diversity levels, likely due to its recent introduction from a gorilla-infecting ancestor (approximately 10,000 to 50,000 years ago), some genes display extremely high diversity levels. In particular, certain proteins expressed on the surface of human red blood cell-infecting merozoites (merozoite surface proteins (MSPs)) possess exactly 2 deeply diverged lineages that have seemingly not recombined. While of considerable interest, the evolutionary origin of this phenomenon remains unknown. In this study, we analysed the genetic diversity of 2 of the most variable MSPs, DBLMSP and DBLMSP2, which are paralogs (descended from an ancestral duplication). Despite thousands of available Illumina WGS datasets from malaria-endemic countries, diversity in these genes has been hard to characterise as reads containing highly diverged alleles completely fail to align to the reference genome. To solve this, we developed a pipeline leveraging genome graphs, enabling us to genotype them at high accuracy and completeness. Using our newly- resolved sequences, we found that both genes exhibit 2 deeply diverged lineages in a specific protein domain (DBL) and that one of the 2 lineages is shared across the genes. We identified clear evidence of nonallelic gene conversion between the 2 genes as the likely mechanism behind sharing, leading us to propose that gene conversion between diverged paralogs, and not recombination suppression, can generate this surprising genealogy; a model that is furthermore consistent with high diversity levels in these 2 genes despite the strong historical P. falciparum transmission bottleneck.

A 25-year odyssey of genomic technology advances and structural variant discovery

This perspective focuses on advances in genome technology over the last 25 years and their impact on germline variant discovery within the field of human genetics. The field has witnessed tremendous technological advances from microarrays to short-read sequencing and now long-read sequencing. Each technology has provided genome-wide access to different classes of human genetic variation. We are now on the verge of comprehensive variant detection of all forms of variation for the first time with a single assay. We predict that this transition will further transform our understanding of human health and biology and, more importantly, provide novel insights into the dynamic mutational processes shaping our genomes.

Complex genotypes in family with metachromatic leukodystrophy: Effect of trans and cis mutations distribution on the phenotype variability

Metachromatic leukodystrophy (MLD) is a severe metabolic disorder caused by the deficient activity of arylsulfatase A due to ARSA gene mutations. According to the age of onset, MLD is classified into three forms: infantile, juvenile, and adult. In our study, we aimed to perform a genetic analysis for two siblings with juvenile MLD for a better characterization of the molecular mechanisms behind the disease. A consanguineous family including two MLD patients (PII.1 and PII.2) was enrolled in our study. The diagnosis was made based on the clinical and neuroimaging investigations. The sequencing of ARSA gene was performed followed by in silico analysis. Besides, the cis/trans distribution of the variants was verified through a PCR-RFLP. The ARSA gene sequencing revealed three known variants, two exonic c.1055A > G and c.1178C > G and an intronic one (c.1524 + 95A > G) in the 3'UTR region. All variants were present at heterozygous state in the two siblings and their mother. The assessment of the cis/trans distribution showed the presence of these variants in cis within the mother, while PII.2 and PII.2 present the c.1055A > G/c.1524 + 95A > G and the c.1178C > G in trans. Additionally, PII.1 harbored a de novo novel missense variant c.1119G > T, whose pathogenicity was supported by our predictive results. Our genetic findings, supported by a clinical examination, confirmed the affection of the mother by the adult MLD. Our results proved the implication of the variable distribution of the found variants in the age of MLD onset. Besides, we described a variable severity between the two siblings due to the de novo pathogenic variant. In conclusion, we identified a complex genotype of ARSA variants within two MLD siblings with a variable severity due to a de novo variant present in one of them. Our results allowed the establishment of an adult MLD diagnosis and highlighted the importance of an assessment of the trans/cis distribution in the cases of complex genotypes.

Rad53 regulates the lifetime of Rdh54 at homologous recombination intermediates

Rdh54 is a conserved DNA translocase that participates in homologous recombination (HR), DNA checkpoint adaptation, and chromosome segregation. Saccharomyces cerevisiae Rdh54 is a known target of the Mec1/Rad53 signaling axis, which globally protects genome integrity during DNA metabolism. While phosphorylation of DNA repair proteins by Mec1/Rad53 is critical for HR progression little is known about how specific post translational modifications alter HR reactions. Phosphorylation of Rdh54 is linked to protection of genomic integrity but the consequences of modification remain poorly understood. Here, we demonstrate that phosphorylation of the Rdh54 C-terminus by the effector kinase Rad53 regulates Rdh54 clustering activity as revealed by single molecule imaging. This stems from phosphorylation dependent and independent interactions between Rdh54 and Rad53. Genetic assays reveal that loss of phosphorylation leads to phenotypic changes resulting in loss-of-heterozygosity (LOH) outcomes. Our data highlight Rad53 as a key regulator of HR intermediates through activation and attenuation of Rdh54 motor function.

Enhancement of Repeat-Mediated Deletion Rearrangement Induced by Particle Irradiation in a RecA-Dependent Manner in Escherichia coli

Repeat-mediated deletion (RMD) rearrangement is a major source of genome instability and can be deleterious to the organism, whereby the intervening sequence between two repeats is deleted along with one of the repeats. RMD rearrangement is likely induced by DNA double-strand breaks (DSBs); however, it is unclear how the complexity of DSBs influences RMD rearrangement. Here, a transgenic Escherichia coli strain K12 MG1655 with a lacI repeat-controlled amp activation was used while taking advantage of particle irradiation, such as proton and carbon irradiation, to generate different complexities of DSBs. Our research confirmed the enhancement of RMD under proton and carbon irradiation and revealed a positive correlation between RMD enhancement and LET. In addition, RMD enhancement could be suppressed by an intermolecular homologous sequence, which was regulated by its composition and length. Meanwhile, RMD enhancement was significantly stimulated by exogenous λ-Red recombinase. Further results investigating its mechanisms showed that the enhancement of RMD, induced by particle irradiation, occurred in a RecA-dependent manner. Our finding has a significant impact on the understanding of RMD rearrangement and provides some clues for elucidating the repair process and possible outcomes of complex DNA damage.

Unzipped chromosome-level genomes reveal allopolyploid nematode origin pattern as unreduced gamete hybridization

The formation and consequences of polyploidization in animals with clonal reproduction remain largely unknown. Clade I root-knot nematodes (RKNs), characterized by parthenogenesis and allopolyploidy, show a widespread geographical distribution and extensive agricultural destruction. Here, we generated 4 unzipped polyploid RKN genomes and identified a putative novel alternative telomeric element. Then we reconstructed 4 chromosome-level assemblies and resolved their genome structures as AAB for triploid and AABB for tetraploid. The phylogeny of subgenomes revealed polyploid RKN origin patterns as hybridization between haploid and unreduced gametes. We also observed extensive chromosomal fusions and homologous gene expression decrease after polyploidization, which might offset the disadvantages of clonal reproduction and increase fitness in polyploid RKNs. Our results reveal a rare pathway of polyploidization in parthenogenic polyploid animals and provide a large number of high-precision genetic resources that could be used for RKN prevention and control.

Fungal antigenic variation using mosaicism and reassortment of subtelomeric genes' repertoires

Surface antigenic variation is crucial for major pathogens that infect humans. To escape the immune system, they exploit various mechanisms. Understanding these mechanisms is important to better prevent and fight the deadly diseases caused. Those used by the fungus Pneumocystis jirovecii that causes life-threatening pneumonia in immunocompromised individuals remain poorly understood. Here, though this fungus is currently not cultivable, our detailed analysis of the subtelomeric sequence motifs and genes encoding surface proteins suggests that the system involves the reassortment of the repertoire of ca. 80 non-expressed genes present in each strain, from which single genes are retrieved for mutually exclusive expression. Dispersion of the new repertoires, supposedly by healthy carrier individuals, appears very efficient because identical alleles are observed in patients from different countries. Our observations reveal a unique strategy of antigenic variation. They also highlight the possible role in genome rearrangements of small imperfect mirror sequences forming DNA triplexes.

Unexpectedly high mutation rate of cyp11b1 compared to cyp21a2 in randomly-selected turkish women: a large screening study

PURPOSE

Congenital adrenal hyperplasia (CAH) is a group of autosomal recessive disorders resulting from enzyme deficiencies associated with steroidogenesis. The clinical presentation of non-classic CAH (NCAH) in females is often indistinguishable from other hyperandrogenic disorders like polycystic ovary syndrome (PCOS). The data on the prevalence of NCAH in unselected women in the literature is scanty. The research aimed to evaluate the prevalence of NCAH, carrier frequencies, and the correlation between clinical symptoms and genotype in Turkish women.

METHODS

The study group comprised two hundred and seventy randomly-selected unrelated asymptomatic women of reproductive age (18-45). Subjects were recruited from female blood donors. All volunteers underwent clinical examination and hormone measurements. The protein-encoding exons and exon-intron boundaries of the CYP21A2, CYP11B1, HSD3β2 and CYP21A2 promoter were sequenced by direct DNA sequencing.

RESULTS

After genotyping, seven (2.2%) individuals were diagnosed with NCAH. The heterozygous carrier frequencies of CYP21A2, CYP21A2 promoter, CYP11B1, and HSD3β2 genes with 34, 34, 41, and 1 pathologic mutation were determined at 12.6%, 12.6%, 15.2%, and 0.37% of volunteers, respectively. Gene-conversion (GC) frequencies between CYP21A2/CYP21A1P and CYP11B1/CYP11B2 were determined as 10.4% and 14.8%, respectively.

CONCLUSION

Despite GC-derived higher mutation frequency determined in the CYP11B1 gene, the reason for the low frequency of NCAH due to 11OHD compared to 21OHD might be that gene-conversion arises with active CYP11B2 rather than an inactive pseudogene. HSD3β1 exhibits high homology with HSD3β2 located on the same chromosome; remarkably, it demonstrates low heterozygosity and no GC, most probably the outcome of a tissue-specific expression pattern.

Systematic analysis of paralogous regions in 41,755 exomes uncovers clinically relevant variation

The short lengths of short-read sequencing reads challenge the analysis of paralogous genomic regions in exome and genome sequencing data. Most genetic variants within these homologous regions therefore remain unidentified in standard analyses. Here, we present a method (Chameleolyser) that accurately identifies single nucleotide variants and small insertions/deletions (SNVs/Indels), copy number variants and ectopic gene conversion events in duplicated genomic regions using whole-exome sequencing data. Application to a cohort of 41,755 exome samples yields 20,432 rare homozygous deletions and 2,529,791 rare SNVs/Indels, of which we show that 338,084 are due to gene conversion events. None of the SNVs/Indels are detectable using regular analysis techniques. Validation by high-fidelity long-read sequencing in 20 samples confirms >88% of called variants. Focusing on variation in known disease genes leads to a direct molecular diagnosis in 25 previously undiagnosed patients. Our method can readily be applied to existing exome data.

Shared evolutionary processes shape landscapes of genomic variation in the great apes

For at least the past five decades population genetics, as a field, has worked to describe the precise balance of forces that shape patterns of variation in genomes. The problem is challenging because modelling the interactions between evolutionary processes is difficult, and different processes can impact genetic variation in similar ways. In this paper, we describe how diversity and divergence between closely related species change with time, using correlations between landscapes of genetic variation as a tool to understand the interplay between evolutionary processes. We find strong correlations between landscapes of diversity and divergence in a well sampled set of great ape genomes, and explore how various processes such as incomplete lineage sorting, mutation rate variation, GC-biased gene conversion and selection contribute to these correlations. Through highly realistic, chromosome-scale, forward-in-time simulations we show that the landscapes of diversity and divergence in the great apes are too well correlated to be explained via strictly neutral processes alone. Our best fitting simulation includes both deleterious and beneficial mutations in functional portions of the genome, in which 9% of fixations within those regions is driven by positive selection. This study provides a framework for modelling genetic variation in closely related species, an approach which can shed light on the complex balance of forces that have shaped genetic variation.

Inference of Coalescence Times and Variant Ages Using Convolutional Neural Networks

Accurate inference of the time to the most recent common ancestor (TMRCA) between pairs of individuals and of the age of genomic variants is key in several population genetic analyses. We developed a likelihood-free approach, called CoalNN, which uses a convolutional neural network to predict pairwise TMRCAs and allele ages from sequencing or SNP array data. CoalNN is trained through simulation and can be adapted to varying parameters, such as demographic history, using transfer learning. Across several simulated scenarios, CoalNN matched or outperformed the accuracy of model-based approaches for pairwise TMRCA and allele age prediction. We applied CoalNN to settings for which model-based approaches are under-developed and performed analyses to gain insights into the set of features it uses to perform TMRCA prediction. We next used CoalNN to analyze 2,504 samples from 26 populations in the 1,000 Genome Project data set, inferring the age of ∼80 million variants. We observed substantial variation across populations and for variants predicted to be pathogenic, reflecting heterogeneous demographic histories and the action of negative selection. We used CoalNN's predicted allele ages to construct genome-wide annotations capturing the signature of past negative selection. We performed LD-score regression analysis of heritability using summary association statistics from 63 independent complex traits and diseases (average N=314k), observing increased annotation-specific effects on heritability compared to a previous allele age annotation. These results highlight the effectiveness of using likelihood-free, simulation-trained models to infer properties of gene genealogies in large genomic data sets.

Phenotype reversion as "natural gene therapy" in Fanconi anemia by a gene conversion event

Somatic mosaicism appears as a recurrent phenomenon among patients suffering from Fanconi anemia (FA), but its direct prognostic significance mostly remains an open question. The clinical picture of FA mosaic subjects could indeed vary from just mild features to severe hematologic failure. Here, we illustrate the case of a proband whose FA familiarity, modest signs (absence of hematological anomalies and fertility issues), and chromosome fragility test transition to negative overtime were suggestive of somatic mosaicism. In line with this hypothesis, genetic testing on patient's peripheral blood and buccal swab reported the presence of the only FANCA paternal variant (FANCA:c.2638C>T, p. Arg880*) and of both parental alleles (the additional FANCA:c.3164G>A, p. Arg1055Gln), respectively. Moreover, the SNP analysis performed on the same biological specimens allowed us to attribute the proband's mosaicism status to a possible gene conversion mechanism. Our case clearly depicts the positive association between somatic mosaicism and the proband's favorable clinical course due to the occurrence of the reversion event at the hematopoietic stem cell level. Since this condition concerns only a limited subgroup of FA individuals, the accurate evaluation of the origin and extent of clonality would be key to steer clinicians toward the most appropriate therapeutic decision for their FA mosaic patients.

Detection of Mosaic Absence of Heterozygosity (AOH) Using Low-Pass Whole Genome Sequencing in Prenatal Diagnosis: A Preliminary Report

Objective: Mosaicism is a common biological phenomenon in organisms and has been reported in many types of chromosome abnormalities, including the absence of heterozygosity (AOH). Due to the detection limitations of the sequencing approach, mosaic AOH events are rarely assessed in clinical cases. Herein, we report the performance of mosaic AOH identification using a low-pass (5~8-fold) WGS method (termed 'CMA-seq', an abbreviation for 'Chromosome Analysis by Sequencing') in fetal genetic diagnosis. Methods: Thirty AOH-negative, eleven constitutional AOH, and three mosaic AOH samples were collected as training data sets to develop the algorithm and evaluate the suitable thresholds for distinguishing mosaic AOH. Twenty-four new chromosomal aberrant cases, along with sixteen constitutional AOH samples, which were previously ascertained via the SNP-array-based method, were used as a validation data set to measure the performance in terms of sensitivity and specificity of this algorithm. Results: A new statistic, 'D-value', was implemented to identify and distinguish constitutional and mosaic AOH events. The reporting thresholds for constitutional and mosaic AOH were also established. In the validation set consisting of 24 new cases, seven constitutional AOH cases and 1 mosaic AOH case were successfully identified, indicating that the results were consistent with those of the SNP-array-based method. The results of all sixteen constitutional AOH validation samples also met the threshold requirements. Conclusions: In this study, we developed a new bioinformatic algorithm to accurately distinguish mosaic AOH from constitutional AOH by low-pass WGS. However, due to the small sample size of the training data set, the algorithm proposed in this manuscript still needs further refinements.

Interlocus Gene Conversion, Natural Selection, and Paralog Homogenization

Following a duplication, the resulting paralogs tend to diverge. While mutation and natural selection can accelerate this process, they can also slow it. Here, we quantify the paralog homogenization that is caused by point mutations and interlocus gene conversion (IGC). Among 164 duplicated teleost genes, the median percentage of postduplication codon substitutions that arise from IGC rather than point mutation is estimated to be between 7% and 8%. By differentiating between the nonsynonymous codon substitutions that homogenize the protein sequences of paralogs and the nonhomogenizing nonsynonymous substitutions, we estimate the homogenizing nonsynonymous rates to be higher for 163 of the 164 teleost data sets as well as for all 14 data sets of duplicated yeast ribosomal protein-coding genes that we consider. For all 14 yeast data sets, the estimated homogenizing nonsynonymous rates exceed the synonymous rates.

The human fungal pathogen Aspergillus fumigatus can produce the highest known number of meiotic crossovers

Sexual reproduction involving meiosis is essential in most eukaryotes. This produces offspring with novel genotypes, both by segregation of parental chromosomes as well as crossovers between homologous chromosomes. A sexual cycle for the opportunistic human pathogenic fungus Aspergillus fumigatus is known, but the genetic consequences of meiosis have remained unknown. Among other Aspergilli, it is known that A. flavus has a moderately high recombination rate with an average of 4.2 crossovers per chromosome pair, whereas A. nidulans has in contrast a higher rate with 9.3 crossovers per chromosome pair. Here, we show in a cross between A. fumigatus strains that they produce an average of 29.9 crossovers per chromosome pair and large variation in total map length across additional strain crosses. This rate of crossovers per chromosome is more than twice that seen for any known organism, which we discuss in relation to other genetic model systems. We validate this high rate of crossovers through mapping of resistance to the laboratory antifungal acriflavine by using standing variation in an undescribed ABC efflux transporter. We then demonstrate that this rate of crossovers is sufficient to produce one of the common multidrug resistant haplotypes found in the cyp51A gene (TR34/L98H) in crosses among parents harboring either of 2 nearby genetic variants, possibly explaining the early spread of such haplotypes. Our results suggest that genomic studies in this species should reassess common assumptions about linkage between genetic regions. The finding of an unparalleled crossover rate in A. fumigatus provides opportunities to understand why these rates are not generally higher in other eukaryotes.

Technological and computational approaches to detect somatic mosaicism in epilepsy

Lesional epilepsy is a common and severe disease commonly associated with malformations of cortical development, including focal cortical dysplasia and hemimegalencephaly. Recent advances in sequencing and variant calling technologies have identified several genetic causes, including both short/single nucleotide and structural somatic variation. In this review, we aim to provide a comprehensive overview of the methodological advancements in this field while highlighting the unresolved technological and computational challenges that persist, including ultra-low variant allele fractions in bulk tissue, low availability of paired control samples, spatial variability of mutational burden within the lesion, and the issue of false-positive calls and validation procedures. Information from genetic testing in focal epilepsy may be integrated into clinical care to inform histopathological diagnosis, postoperative prognosis, and candidate precision therapies.

Repeat-Induced Point Mutation and Gene Conversion Coinciding with Heterochromatin Shape the Genome of a Plant-Pathogenic Fungus

Meiosis is associated with genetic changes in the genome-via recombination, gene conversion, and mutations. The occurrence of gene conversion and mutations during meiosis may further be influenced by the chromatin conformation, similar to the effect of the chromatin conformation on the mitotic mutation rate. To date, however, the exact distribution and type of meiosis-associated changes and the role of the chromatin conformation in this context are largely unexplored. Here, we determine recombination, gene conversion, and de novo mutations using whole-genome sequencing of all meiotic products of 23 individual meioses in Zymoseptoria tritici, an important pathogen of wheat. We confirm a high genome-wide recombination rate of 65 centimorgan (cM)/Mb and see higher recombination rates on the accessory compared to core chromosomes. A substantial fraction of 0.16% of all polymorphic markers was affected by gene conversions, showing a weak GC-bias and occurring at higher frequency in regions of constitutive heterochromatin, indicated by the histone modification H3K9me3. The de novo mutation rate associated with meiosis was approximately three orders of magnitude higher than the corresponding mitotic mutation rate. Importantly, repeat-induced point mutation (RIP), a fungal defense mechanism against duplicated sequences, is active in Z. tritici and responsible for the majority of these de novo meiotic mutations. Our results indicate that the genetic changes associated with meiosis are a major source of variability in the genome of an important plant pathogen and shape its evolutionary trajectory. IMPORTANCE The impact of meiosis on the genome composition via gene conversion and mutations is mostly poorly understood, in particular, for non-model species. Here, we sequenced all four meiotic products for 23 individual meioses and determined the genetic changes caused by meiosis for the important fungal wheat pathogen Zymoseptoria tritici. We found a high rate of gene conversions and an effect of the chromatin conformation on gene conversion rates. Higher conversion rates were found in regions enriched with the H3K9me3-a mark for constitutive heterochromatin. Most importantly, meiosis was associated with a much higher frequency of de novo mutations than mitosis; 78% of the meiotic mutations were caused by repeat-induced point mutations-a fungal defense mechanism against duplicated sequences. In conclusion, the genetic changes associated with meiosis are therefore a major factor shaping the genome of this fungal pathogen.

Limitations of gene editing assessments in human preimplantation embryos

Range of DNA repair in response to double-strand breaks induced in human preimplantation embryos remains uncertain due to the complexity of analyzing single- or few-cell samples. Sequencing of such minute DNA input requires a whole genome amplification that can introduce artifacts, including coverage nonuniformity, amplification biases, and allelic dropouts at the target site. We show here that, on average, 26.6% of preexisting heterozygous loci in control single blastomere samples appear as homozygous after whole genome amplification indicative of allelic dropouts. To overcome these limitations, we validate on-target modifications seen in gene edited human embryos in embryonic stem cells. We show that, in addition to frequent indel mutations, biallelic double-strand breaks can also produce large deletions at the target site. Moreover, some embryonic stem cells show copy-neutral loss of heterozygosity at the cleavage site which is likely caused by interallelic gene conversion. However, the frequency of loss of heterozygosity in embryonic stem cells is lower than in blastomeres, suggesting that allelic dropouts is a common whole genome amplification outcome limiting genotyping accuracy in human preimplantation embryos.

Identification of miR-3182 and miR-3143 target genes involved in the cell cycle as a novel approach in TNBC treatment: A systems biology approach

Triple-negative breast cancer (TNBC) is an aggressive subtype of breast cancer with a poor prognosis, lacking therapeutic targets. miRNAs play crucial roles in TNBC through regulating various mechanisms, including cellular growth and proliferation. This study aims to identify critical target genes of two novel miRNAs (miR-3143 and miR-3182) involved in the cell cycle of TNBC as possible therapeutic targets and investigates their regulatory and therapeutic roles through a systems biology approach and in vitro experiment. Datasets related to the TNBC cell line (MDA-MB-231) were screened and retrieved, and Gene regulatory networks were constructed. Significant regulatory motifs were detected and analyzed using the FANMOD and Cytoscape analyzer, and the clusters and seeds were identified using the MCODE. Functional enrichment analysis was also performed using DAVID and STRING. The most critical genes were determined using the analysis of GRN motifs and PPI clusters. The essential genes involved in the cell cycle were selected and verified using the bc-GenExMiner v4.7. We overexpressed miR-3143 and miR-3182 in the MDA-MB-231 cell line using human umbilical cord mesenchymal stem cell (HUCMSC)-miRNA loaded exosomes, and the expression of the critical target genes was investigated using RT-qPCR. We identified eight critical genes as potential therapeutic targets. Their expression decreased by overexpression of miR-3143 and miR-3182 in RT-qPCR. The identified critical genes have probably significant roles in the pathogenesis of TNBC through the cell cycle. We suggest that the overexpression of miR-3143 and miR-3182 could be a new therapeutic candidate in TNBC and is worth more investigation.

Genomic structural variation: A complex but important driver of human evolution

Structural variants (SVs)-including duplications, deletions, and inversions of DNA-can have significant genomic and functional impacts but are technically difficult to identify and assay compared with single-nucleotide variants. With the aid of new genomic technologies, it has become clear that SVs account for significant differences across and within species. This phenomenon is particularly well-documented for humans and other primates due to the wealth of sequence data available. In great apes, SVs affect a larger number of nucleotides than single-nucleotide variants, with many identified SVs exhibiting population and species specificity. In this review, we highlight the importance of SVs in human evolution by (1) how they have shaped great ape genomes resulting in sensitized regions associated with traits and diseases, (2) their impact on gene functions and regulation, which subsequently has played a role in natural selection, and (3) the role of gene duplications in human brain evolution. We further discuss how to incorporate SVs in research, including the strengths and limitations of various genomic approaches. Finally, we propose future considerations in integrating existing data and biospecimens with the ever-expanding SV compendium propelled by biotechnology advancements.

Characterization of the immunoglobulin lambda chain locus from diverse populations reveals extensive genetic variation

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-Cap. From these 16 individuals, we identified significant allelic diversity, including 36 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

CFH and CFHR structural variants in atypical Hemolytic Uremic Syndrome: Prevalence, genomic characterization and impact on outcome

Introduction

Atypical hemolytic uremic syndrome (aHUS) is a rare disease that manifests with microangiopathic hemolytic anemia, thrombocytopenia, and acute renal failure, and is associated with dysregulation of the alternative complement pathway. The chromosomal region including CFH and CFHR1-5 is rich in repeated sequences, favoring genomic rearrangements that have been reported in several patients with aHUS. However, there are limited data on the prevalence of uncommon CFH-CFHR genomic rearrangements in aHUS and their impact on disease onset and outcomes.

Methods

In this study, we report the results of CFH-CFHR Copy Number Variation (CNV) analysis and the characterization of resulting structural variants (SVs) in a large cohort of patients, including 258 patients with primary aHUS and 92 with secondary forms.

Results

We found uncommon SVs in 8% of patients with primary aHUS: 70% carried rearrangements involving CFH alone or CFH and CFHR (group A; n=14), while 30% exhibited rearrangements including only CFHRs (group B; n=6). In group A, 6 patients presented CFH::CFHR1 hybrid genes, 7 patients carried duplications in the CFH-CFHR region that resulted either in the substitution of the last CFHR1 exon(s) with those of CFH (CFHR1::CFH reverse hybrid gene) or in an internal CFH duplication. In group A, the large majority of aHUS acute episodes not treated with eculizumab (12/13) resulted in chronic ESRD; in contrast, anti-complement therapy induced remission in 4/4 acute episodes. aHUS relapse occurred in 6/7 grafts without eculizumab prophylaxis and in 0/3 grafts with eculizumab prophylaxis. In group B, 5 subjects had the CFHR3_1-5::CFHR4₁₀ hybrid gene and one had 4 copies of CFHR1 and CFHR4. Compared with group A, patients in group B exhibited a higher prevalence of additional complement abnormalities and earlier disease onset. However, 4/6 patients in this group underwent complete remission without eculizumab treatment. In secondary forms we identified uncommon SVs in 2 out of 92 patients: the CFHR3_1-5::CFHR4₁₀ hybrid and a new internal duplication of CFH.

Discussion

In conclusion, these data highlight that uncommon CFH-CFHR SVs are frequent in primary aHUS and quite rare in secondary forms. Notably, genomic rearrangements involving the CFH are associated with a poor prognosis but carriers respond to anti-complement therapy.

High-Quality Genomes of Pangolins: Insights into the Molecular Basis of Scale Formation and Adaption to Myrmecophagous Diet

Pangolins are one of nature's most fascinating species being scales covered and myrmecophagous diet, yet relatively little is known about the molecular basis. Here, we combine the multi-omics, evolution, and fundamental proteins feature analysis of both Chinese and Malayan pangolins, highlighting the molecular mechanism of both myrmecophagous diet and scale formation, representing a fascinating evolutionary strategy to occupy the unique ecological niches. In contrast to conserved organization of epidermal differentiation complex, pangolin has undergone large scale variation and gene loss events causing expression pattern and function conversion that contribute to cornified epithelium structures on stomach to adapt myrmecophagous diet. Our assemblies also enable us to discover large copies number of high glycine-tyrosine keratin-associated proteins (HGT-KRTAPs). In addition, highly homogenized tandem array, amino content, and the specific expression pattern further validate the strong connection between the molecular mechanism of scale hardness and HGT-KRTAPs.

A Revised Perspective on the Evolution of Troponin I and Troponin T Gene Families in Vertebrates

The troponin (Tn) complex, responsible for the Ca2+ activation of striated muscle, is composed of three interacting protein subunits: TnC, TnI, and TnT, encoded by TNNC, TNNI, and TNNT genes. TNNI and TNNT are sister gene families, and in mammals the three TNNI paralogs (TNNI1, TNNI2, TNNI3), which encode proteins with tissue-specific expression, are each in close genomic proximity with one of the three TNNT paralogs (TNNT2, TNNT3, TNNT1, respectively). It has been widely presumed that all vertebrates broadly possess genes of these same three classes, although earlier work has overlooked jawless fishes (cyclostomes) and cartilaginous fishes (chimeras, rays, and sharks), which are distantly related to other jawed vertebrates. With a new phylogenetic and synteny analysis of a diverse array of vertebrates including these taxonomic groups, we define five distinct TNNI classes (TNNI1-5), with TNNI4 and TNNI5 being only present in non-amniote vertebrates and typically found in tandem, and four classes of TNNT (TNNT1-4). These genes are located in four genomic loci that were generated by the 2R whole-genome duplications. TNNI3, encoding "cardiac TnI" in tetrapods, was independently lost in cartilaginous and ray-finned fishes. Instead, ray-finned fishes predominantly express TNNI1 in the heart. TNNI5 is highly expressed in shark hearts and contains a N-terminal extension similar to that of TNNI3 found in tetrapod hearts. Given that TNNI3 and TNNI5 are distantly related, this supports the hypothesis that the N-terminal extension may be an ancestral feature of vertebrate TNNI and not an innovation unique to TNNI3, as has been commonly believed.

The genome loading model for the origin and maintenance of sex in eukaryotes

Understanding why sexual reproduction-which involves syngamy (union of gametes) and meiosis-emerged and how it has subsisted for millions of years remains a fundamental problem in biology. Considered as the essence of sex, meiotic recombination is initiated by a DNA double-strand break (DSB) that forms on one of the pairing homologous chromosomes. This DNA lesion is subsequently repaired by gene conversion, the non-reciprocal transfer of genetic information from the intact homolog. A major issue is which of the pairing homologs undergoes DSB formation. Accumulating evidence shows that chromosomal sites where the pairing homologs locally differ in size, i.e., are heterozygous for an insertion or deletion, often display disparity in gene conversion. Biased conversion tends to duplicate insertions and lose deletions. This suggests that DSB is preferentially formed on the "shorter" homologous region, which thereby acts as the recipient for DNA transfer. Thus, sex primarily functions as a genome (re)loading mechanism. It ensures the restoration of formerly lost DNA sequences (deletions) and allows the efficient copying and, mainly in eukaryotes, subsequent spreading of newly emerged sequences (insertions) arising initially in an individual genome, even if they confer no advantage to the host. In this way, sex simultaneously repairs deletions and increases genetic variability underlying adaptation. The model explains a remarkable increase in DNA content during the evolution of eukaryotic genomes.

RNA polymerase III directly participates in DNA homologous recombination

A recent study showed that RNA transcription is directly involved in DNA homologous recombination (HR). The first step in HR is end resection, which degrades a few kilobases or more from the 5'-end strand at DNA breaks, but the 3'-end strand remains strictly intact. Such protection of the 3'-end strand is achieved by the transient formation of an RNA-DNA hybrid structure. The RNA strand in the hybrid is newly synthesized by RNA polymerase III. The revelation of the existence of an RNA-DNA hybrid intermediate should further help resolve several long-standing questions of HR. In this article, we also put forward our views on some controversial issues related to RNA-DNA hybrids, RNA polymerases, and the protection of 3'-end strands.

Overcoming the Pitfalls of Next-Generation Sequencing-Based Molecular Diagnosis of Shwachman-Diamond Syndrome

Shwachman-Diamond syndrome (SDS) is the second most common cause of exocrine pancreatic insufficiency, and 90% of patients carry mutations in the SBDS gene, the most common being the c.183_184delinsCT and c.258+2T>C variants. However, precise detection of these most contributory variants by conventional short-read next-generation sequencing data analysis is limited because of the SBDS/SBDSP1 highly homologous sequences. In this study, an efficient approach was established to infer the haplotype of SBDS based on the expectation-maximization algorithm. The workflow was retrospectively applied to detect the two most common SBDS variants in a Chinese SDS high-risk cohort, and a systematic comparison of variant detection results was performed between the workflow and conventional next-generation sequencing analysis based on Sanger sequencing validation. Among the Chinese SDS high-risk cohort (n = 47) and their available parents (n = 64), the established workflow improved the diagnostic rate for these two variants by 27.7% (95% CI, 15.6%-42.6%) compared with conventional analysis. For overall variant detection, the established workflow achieved 100% (95% CI, 92.5%-100%) concordance with Sanger sequencing, whereas conventional analysis showed only 65.8% accuracy; these results included 25.2% with missed variant calls, 7.2% with diagnosed but inaccurate variant calls, and 1.8% with false-positive calls. With its favorable result in both SDS patient diagnosis and carrier detection performance, the provided workflow showed its potential in clinical application for SDS molecular diagnosis.

Structure and evolution of the squamate major histocompatibility complex as revealed by two Anolis lizard genomes

The major histocompatibility complex (MHC) is an important genomic region for adaptive immunity and has long been studied in ecological and evolutionary contexts, such as disease resistance and mate and kin selection. The MHC has been investigated extensively in mammals and birds but far less so in squamate reptiles, the third major radiation of amniotes. We localized the core MHC genomic region in two squamate species, the green anole (Anolis carolinensis) and brown anole (A. sagrei), and provide the first detailed characterization of the squamate MHC, including the presence and ordering of known MHC genes in these species and comparative assessments of genomic structure and composition in MHC regions. We find that the Anolis MHC, located on chromosome 2 in both species, contains homologs of many previously-identified mammalian MHC genes in a single core MHC region. The repetitive element composition in anole MHC regions was similar to those observed in mammals but had important distinctions, such as higher proportions of DNA transposons. Moreover, longer introns and intergenic regions result in a much larger squamate MHC region (11.7 Mb and 24.6 Mb in the green and brown anole, respectively). Evolutionary analyses of MHC homologs of anoles and other representative amniotes uncovered generally monophyletic relationships between species-specific homologs and a loss of the peptide-binding domain exon 2 in one of two mhc2β gene homologs of each anole species. Signals of diversifying selection in each anole species was evident across codons of mhc1, many of which appear functionally relevant given known structures of this protein from the green anole, chicken, and human. Altogether, our investigation fills a major gap in understanding of amniote MHC diversity and evolution and provides an important foundation for future squamate-specific or vertebrate-wide investigations of the MHC.

Hereditary Transthyretin-Related Amyloidosis: Genetic Heterogeneity and Early Personalized Gene Therapy

Point mutations of the transthyretin (TTR) gene are related with hereditary amyloidosis (hATTR). The number of people affected by this rare disease is only partially estimated. The real impact of somatic mosaicism and other genetic factors on expressivity, complexity, progression, and transmission of the disease should be better investigated. The relevance of this rare disease is increasing and many efforts have been made to improve the time to diagnosis and to estimate the real number of cases in endemic and non-endemic areas. In this context, somatic mosaicism should be better investigated to explain the complexity of the heterogeneity of the hATTR clinical features, to better estimate the number of new cases, and to focus on early and personalized gene therapy. Gene therapy can potentially improve the living conditions of affected individuals and is one of the central goals in research on amyloidosis related to the TTR gene, with the advantage of overcoming liver transplantation as the sole treatment for hATTR disease.

Individual Genetic Heterogeneity

Genetic variation has been widely covered in literature, however, not from the perspective of an individual in any species. Here, a synthesis of genetic concepts and variations relevant for individual genetic constitution is provided. All the different levels of genetic information and variation are covered, ranging from whether an organism is unmixed or hybrid, has variations in genome, chromosomes, and more locally in DNA regions, to epigenetic variants or alterations in selfish genetic elements. Genetic constitution and heterogeneity of microbiota are highly relevant for health and wellbeing of an individual. Mutation rates vary widely for variation types, e.g., due to the sequence context. Genetic information guides numerous aspects in organisms. Types of inheritance, whether Mendelian or non-Mendelian, zygosity, sexual reproduction, and sex determination are covered. Functions of DNA and functional effects of variations are introduced, along with mechanism that reduce and modulate functional effects, including TARAR countermeasures and intraindividual genetic conflict. TARAR countermeasures for tolerance, avoidance, repair, attenuation, and resistance are essential for life, integrity of genetic information, and gene expression. The genetic composition, effects of variations, and their expression are considered also in diseases and personalized medicine. The text synthesizes knowledge and insight on individual genetic heterogeneity and organizes and systematizes the central concepts.

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.

Mutation accumulation opposes polymorphism: supergenes and the curious case of balanced lethals

Supergenes offer spectacular examples of long-term balancing selection in nature, but their origin and maintenance remain a mystery. Reduced recombination between arrangements, a critical aspect of many supergenes, protects adaptive multi-trait phenotypes but can lead to mutation accumulation. Mutation accumulation can stabilize the system through the emergence of associative overdominance (AOD), destabilize the system, or lead to new evolutionary outcomes. One outcome is the formation of maladaptive balanced lethal systems, where only heterozygotes remain viable and reproduce. We investigated the conditions under which these different outcomes occur, assuming a scenario of introgression after divergence. We found that AOD aided the invasion of a new supergene arrangement and the establishment of a polymorphism. However, this polymorphism was easily destabilized by further mutation accumulation, which was often asymmetric, disrupting the quasi-equilibrium state. Mechanisms that accelerated degeneration tended to amplify asymmetric mutation accumulation between the supergene arrangements and vice-versa. As the evolution of balanced lethal systems requires symmetric degeneration of both arrangements, this leaves only restricted conditions for their evolution, namely small population sizes and low rates of gene conversion. The dichotomy between the persistence of polymorphism and degeneration of supergene arrangements likely underlies the rarity of balanced lethal systems in nature.

This article is part of the theme issue ‘Genomic architecture of supergenes: causes and evolutionary consequences’.

Genome doubling enabled the expansion of yeast vesicle traffic pathways

Vesicle budding and fusion in eukaryotes depend on a suite of protein types, such as Arfs, Rabs, coats and SNAREs. Distinct paralogs of these proteins act at distinct intracellular locations, suggesting a link between gene duplication and the expansion of vesicle traffic pathways. Genome doubling, a common source of paralogous genes in fungi, provides an ideal setting in which to explore this link. Here we trace the fates of paralog doublets derived from the 100-Ma-old hybridization event that gave rise to the whole genome duplication clade of budding yeast. We find that paralog doublets involved in specific vesicle traffic functions and pathways are convergently retained across the entire clade. Vesicle coats and adaptors involved in secretory and early-endocytic pathways are retained as doublets, at rates several-fold higher than expected by chance. Proteins involved in later endocytic steps and intra-Golgi traffic, including the entire set of multi-subunit and coiled-coil tethers, have reverted to singletons. These patterns demonstrate that selection has acted to expand and diversify the yeast vesicle traffic apparatus, across species and time.

Conversion between duplicated genes generated by polyploidization contributes to the divergence of poplar and willow

BACKGROUND

Gene conversion has an important effect on duplicate genes produced by polyploidization. Poplar (Populus trichocarpa) and willow (Salix brachista) are leading models and excellent green plants in the Salicaceae. Although much attention has been paid to the evolution of duplicated genes in poplar and willow, the role of conversion between duplicates generated from polyploidization remains poorly understood.

RESULTS

Here, through genomic synteny analyses, we identified duplicate genes generated by the Salicaceae common tetraploidization (SCT) in the poplar and willow genomes. We estimated that at least 0.58% and 0.25% of poplar and willow duplicates were affected by whole-gene conversion after the poplar-willow divergence, with more (5.73% and 2.66%) affected by partial-gene conversion. Moreover, we found that the converted duplicated genes were unevenly distributed on each chromosome in the two genomes, and the well-preserved homoeologous chromosome regions may facilitate the conversion of duplicates. Notably, we found that conversion maintained the similarity of duplicates, likely contributing to the conservation of certain sequences, but is essentially accelerated the rate of evolution and increased species divergence. In addition, we found that converted duplicates tended to have more similar expression patterns than nonconverted duplicates. We found that genes associated with multigene families were preferentially converted. We also found that the genes encoding conserved structural domains associated with specific traits exhibited a high frequency of conversion.

CONCLUSIONS

Extensive conversion between duplicate genes generated from the SCT contributes to the diversification of the family Salicaceae and has had long-lasting effects on those genes with important biological functions.

Robust and accurate estimation of paralog-specific copy number for duplicated genes using whole-genome sequencing

The human genome contains hundreds of low-copy repeats (LCRs) that are challenging to analyze using short-read sequencing technologies due to extensive copy number variation and ambiguity in read mapping. Copy number and sequence variants in more than 150 duplicated genes that overlap LCRs have been implicated in monogenic and complex human diseases. We describe a computational tool, Parascopy, for estimating the aggregate and paralog-specific copy number of duplicated genes using whole-genome sequencing (WGS). Parascopy is an efficient method that jointly analyzes reads mapped to different repeat copies without the need for global realignment. It leverages multiple samples to mitigate sequencing bias and to identify reliable paralogous sequence variants (PSVs) that differentiate repeat copies. Analysis of WGS data for 2504 individuals from diverse populations showed that Parascopy is robust to sequencing bias, has higher accuracy compared to existing methods and enables prioritization of pathogenic copy number changes in duplicated genes.

Unraveling Structural Rearrangements of the CFH Gene Cluster in Atypical Hemolytic Uremic Syndrome Patients Using Molecular Combing and Long-Fragment Targeted Sequencing

Complement factor H (CFH) and its related proteins have an essential role in regulating the alternative pathway of the complement system. Mutations and structural variants (SVs) of the CFH gene cluster, consisting of CFH and its five related genes (CFHR1-5), have been reported in renal pathologies as well as in complex immune diseases like age-related macular degeneration and systemic lupus erythematosus. SV analysis of this cluster is challenging because of its high degree of sequence homology. Following first-line next-generation sequencing gene panel sequencing, we applied Genomic Vision's Molecular Combing Technology to detect and visualize SVs within the CFH gene cluster and resolve its structural haplotypes completely. This approach was tested in three patients with atypical hemolytic uremic syndrome and known SVs and 18 patients with atypical hemolytic uremic syndrome or complement factor 3 glomerulopathy with unknown CFH gene cluster haplotypes. Three SVs, a CFH/CFHR1 hybrid gene in two patients and a rare heterozygous CFHR4/CFHR1 deletion in trans with the common CFHR3/CFHR1 deletion in a third patient, were newly identified. For the latter, the breakpoints were determined using a targeted enrichment approach for long DNA fragments (Samplix Xdrop) in combination with Oxford Nanopore sequencing. Molecular combing in addition to next-generation sequencing was able to improve the molecular genetic yield in this pilot study. This (cost-)effective approach warrants validation in larger cohorts with CFH/CFHR-associated disease.

Frequent somatic gene conversion as a mechanism for loss of heterozygosity in tumor suppressor genes

The major processes in carcinogenesis include the inactivation of tumor-suppressor genes (TSGs). Although Knudson's two-hit model requires two independent inactivating mutations, perhaps more frequently, a TSG inactivation can occur through a loss of heterozygosity (LOH) of an inactivating mutation. Deletion and uniparental disomy (UPD) have been well documented as LOH mechanisms, but the role of gene conversion is poorly understood. Here, we developed a simple algorithm to detect somatic gene conversion from short-read sequencing data. We applied it to 6285 cancer patient samples, from which 4978 somatic mutations that underwent gene conversion to achieve LOH were found. This number accounted for 14.8% of the total LOH mutations. We further showed that LOH by gene conversion was enriched in TSGs compared with non-TSG genes, showing a significant contribution of gene conversion to carcinogenesis.

Supergene origin and maintenance in Atlantic cod

Supergenes are sets of genes that are inherited as a single marker and encode complex phenotypes through their joint action. They are identified in an increasing number of organisms, yet their origins and evolution remain enigmatic. In Atlantic cod, four megabase-scale supergenes have been identified and linked to migratory lifestyle and environmental adaptations. Here we investigate the origin and maintenance of these four supergenes through analysis of whole-genome-sequencing data, including a new long-read-based genome assembly for a non-migratory Atlantic cod individual. We corroborate the finding that chromosomal inversions underlie all four supergenes, and we show that they originated at different times between 0.40 and 1.66 million years ago. We reveal gene flux between supergene haplotypes where migratory and stationary Atlantic cod co-occur and conclude that this gene flux is driven by gene conversion, on the basis of an increase in GC content in exchanged sites. Additionally, we find evidence for double crossover between supergene haplotypes, leading to the exchange of an ~275 kilobase fragment with genes potentially involved in adaptation to low salinity in the Baltic Sea. Our results suggest that supergenes can be maintained over long timescales in the same way as hybridizing species, through the selective purging of introduced genetic variation.

Atlantic cod carries four supergenes linked to migratory lifestyle and environmental adaptations. Using whole-genome sequencing, the authors show that the genome inversions that underlie the supergenes originated at different times and show gene flux between supergene haplotypes.

Episodes of Rapid Recovery of the Functional Activity of the ras85D Gene in the Evolutionary History of Phylogenetically Distant Drosophila Species

As assemblies of genomes of new species with varying degrees of relationship appear, it becomes obvious that structural rearrangements of the genome, such as inversions, translocations, and transposon movements, are an essential and often the main source of evolutionary variation. In this regard, the following questions arise. How conserved are the regulatory regions of genes? Do they have a common evolutionary origin? And how and at what rate is the functional activity of genes restored during structural changes in the promoter region? In this article, we analyze the evolutionary history of the formation of the regulatory region of the ras85D gene in different lineages of the genus Drosophila, as well as the participation of mobile elements in structural rearrangements and in the replacement of specific areas of the promoter region with those of independent evolutionary origin. In the process, we substantiate hypotheses about the selection of promoter elements from a number of frequently repeated motifs with different degrees of degeneracy in the ancestral sequence, as well as about the restoration of the minimum required set of regulatory sequences using a conversion mechanism or similar.

Evolution of host-microbe cell adherence by receptor domain shuffling

Stable adherence to epithelial surfaces is required for colonization by diverse host-associated microbes. Successful attachment of pathogenic microbes to host cells via adhesin molecules is also the first step in many devastating infections. Despite the primacy of epithelial adherence in establishing host-microbe associations, the evolutionary processes that shape this crucial interface remain enigmatic. Carcinoembryonic antigen-related cell adhesion molecules (CEACAMs) encompass a multifunctional family of vertebrate cell surface proteins which are recurrent targets of bacterial adhesins at epithelial barriers. Here, we show that multiple members of the primate CEACAM family exhibit evidence of repeated natural selection at protein surfaces targeted by bacteria, consistent with pathogen-driven evolution. Divergence of CEACAM proteins between even closely related great apes is sufficient to control molecular interactions with a range of bacterial adhesins. Phylogenetic analyses further reveal that repeated gene conversion of CEACAM extracellular domains during primate divergence plays a key role in limiting bacterial adhesin host tropism. Moreover, we demonstrate that gene conversion has continued to shape CEACAM diversity within human populations, with abundant human CEACAM1 variants mediating evasion of adhesins from pathogenic Neisseria. Together this work reveals a mechanism by which gene conversion shapes first contact between microbes and animal hosts.

Trillions of bacteria live in and on the human body. Most of them are harmless but some can cause serious infections. To grow in or on the body, bacteria often attach to proteins on the surface of cells that make up the lining of tissues like the gut or the throat. In some cases, bacteria use these proteins to invade the cells causing an infection. Genetic mutations in the genes encoding these proteins that protect against infection are more likely to be passed on to future generations. This may lead to rapid spread of these beneficial genes in a population.

A family of proteins called CEACAMs are frequent targets of infection-causing bacteria. These proteins have been shown to play a role in cancer progression. But they also play many helpful roles in the body, including helping transmit messages between cells, aiding cell growth, and helping the immune system recognize pathogens. Scientists are not sure if these multi-tasking CEACAM proteins can evolve to evade bacteria without affecting their other roles.

Baker et al. show that CEACAM proteins targeted by bacteria have undergone rapid evolution in primates. In the experiments, human genes encoding CEACAMs were compared with equivalent genes from 19 different primates. Baker et al. found the changes in human and primate CEACAMs often occur through a process called gene conversion. Gene conversion occurs when DNA sections are copied and pasted from one gene to another. Using laboratory experiments, they showed that some of these changes enabled CEACAM proteins to prevent certain harmful bacteria from binding.

The experiments suggest that some versions of CEACAM genes may protect humans or other primates against bacterial infections. Studies in natural populations are needed to test if this is the case. Learning more about how CEACAM proteins evolve and what they do may help scientists better understand the role they play in cancer and help improve cancer care. Studying CEACAM evolution may also help scientists understand how bacteria and other pathogens drive protein evolution in the body.

The underlying cause of the simple virilizing phenotype in patients with 21-hydroxylase deficiency harboring P31L variant

OBJECTIVE

To analyze the relationship between genotype and phenotype in 21-Hydroxylase deficiency patients harboring P31L variant and the underlying mechanism.

METHODS

A total of 29 Chinese patients with 21-OHD harboring P31L variant were recruited, and the detailed clinical features of the patients were extracted and analyzed retrospectively. The TA clone combined with sequencing of the region containing the promotor and exon1 of CYP21A2 was performed to determine whether the variants in promotor and P31L aligned in cis. We further compared the clinical characteristics of 21-OHD patients between the promoter variant group and no promoter variant group.

RESULTS

Among the 29 patients diagnosed with 21-OHD harboring P31L variant, the incidence of classical simple virilizing form was 62.1%. Thirteen patients owned promoter variants (1 homozygote and 12 heterozygote) and all exhibited SV form. The promoter variants and the P31L variant were located in the same mutant allele as validated by TA cloning and sequencing. There were statistically significant differences in clinical phenotype and 17-OHP level between the patients with and without promoter region variations (P<0.05).

CONCLUSION

There exists high incidence (57.4%) of SV form among the 21-OHD patients harboring P31L variant, and the underlying mechanism is partially due to both the promoter variants and P31L aligning in cis on one allele. Further sequencing of promoter region will provide important hints for the explanation of phenotype in patients harboring P31L.