Axonally synthesized ATF4 transmits a neurodegenerative signal across brain regions.

In Alzheimer's disease (AD) brain, exposure of axons to Aß causes pathogenic changes that spread retrogradely by unknown mechanisms, affecting the entire neuron. We found that locally applied Aß1-42 initiates axonal synthesis of a defined set of proteins including the transcription factor ATF4. Inhibition of local translation and retrograde transport or knockdown of axonal Atf4 mRNA abolished Aß-induced ATF4 transcriptional activity and cell loss. Aß1-42 injection into the dentate gyrus (DG) of mice caused loss of forebrain neurons whose axons project to the DG. Protein synthesis and Atf4 mRNA were upregulated in these axons, and coinjection of Atf4 siRNA into the DG reduced the effects of Aß1-42 in the forebrain. ATF4 protein and transcripts were found with greater frequency in axons in the brain of AD patients. These results reveal an active role for intra-axonal translation in neurodegeneration and identify ATF4 as a mediator for the spread of AD pathology.

Optical genome mapping with genome sequencing identifies subtelomeric Xq28 deletion and inserted 7p22.3 duplication in a male with multisystem developmental disorder.

We report a 17-year-old male with supravalvular stenosis, initial failure to thrive and delayed early development, short stature, acromelia, dysmorphic facial features, hypertelorism, macrocephaly, syringomyelia, hypertension, and anxiety disorder. Fluorescent in situ hybridization (FISH), chromosomal microarray analysis (CMA), and exome sequencing (ES) were nondiagnostic. Combined optical genome mapping (OGM) and genome sequencing (GS) showed a complex rearrangement including an X chromosome with a 22.5 kb deletion in band Xq28 replaced by a 61.4 kb insertion of duplicated chromosome 7p22.3 material. The deletion removes the distal 3' untranslated region (UTR) of FUNDC2, the entire CMC4 and MTCP1, and the first five exons of BRCC3. Transcriptome analysis revealed absent expression of CMC4 and MTCP1 and BRCC3 with normal transcript level of FUNDC2. The inserted duplication includes only one known gene: UNCX. Similar overlapping Xq28 deletions have been reported to be associated with Moyamoya disease (MMD), short stature, hypergonadotropic hypogonadism (HH), and facial dysmorphism. Although he has short stature, our patient does not have signs of Moyamoya arteriopathy or hypogonadism. The structurally abnormal X chromosome was present in his mother, but not in his unaffected brother, maternal uncle, or maternal grandparents. We propose that the combination of his absent Xq28 and duplicated 7p22.3 genomic material is responsible for his phenotype. This case highlights the potential of combined OGM and GS for detecting complex structural variants compared with standard of care genetic testing such as CMA and ES.

The genotypic spectrum of ALDH7A1 mutations resulting in pyridoxine dependent epilepsy: A common epileptic encephalopathy.

Pyridoxine dependent epilepsy (PDE) is a treatable epileptic encephalopathy characterized by a positive response to pharmacologic doses of pyridoxine. Despite seizure control, at least 75% of individuals have intellectual disability and developmental delay. Current treatment paradigms have resulted in improved cognitive outcomes emphasizing the importance of an early diagnosis. As genetic testing is increasingly accepted as first tier testing for epileptic encephalopathies, we aimed to provide a comprehensive overview of ALDH7A1 mutations that cause PDE. The genotypes, ethnic origin and reported gender was collected from 185 subjects with a diagnosis of PDE. The population frequency for the variants in this report and the existing literature were reviewed in the Genome Aggregation Database (gnomAD). Novel variants identified in population databases were also evaluated through in silico prediction software and select variants were over-expressed in an E.coli-based expression system to measure α-aminoadipic semialdehyde dehydrogenase activity and production of α-aminoadipic acid. This study adds 47 novel variants to the literature resulting in a total of 165 reported pathogenic variants. Based on this report, in silico predictions, and general population data, we estimate an incidence of approximately 1:64,352 live births. This report provides a comprehensive overview of known ALDH7A1 mutations that cause PDE, and suggests that PDE may be more common than initially estimated. Due to the relative high frequency of the disease, the likelihood of under-diagnosis given the wide clinical spectrum and limited awareness among clinicians as well as the cognitive improvement noted with early treatment, newborn screening for PDE may be warranted.

A mutation abolishing the ZMPSTE24 cleavage site in prelamin A causes a progeroid disorder.

In 1994 in the Journal of Cell Science, Hennekes and Nigg reported that changing valine to arginine at the endoproteolytic cleavage site in chicken prelamin A abolishes its conversion to lamin A. The consequences of this mutation in an organism have remained unknown. We now report that the corresponding mutation in a human subject leads to accumulation of prelamin A and causes a progeroid disorder. Next generation sequencing of the subject and her parents' exomes identified a de novo mutation in the lamin A/C gene (LMNA) that resulted in a leucine to arginine amino acid substitution at residue 647 in prelamin A. The subject's fibroblasts accumulated prelamin A, a farnesylated protein, which led to an increased percentage of cultured cells with morphologically abnormal nuclei. Treatment with a protein farnesyltransferase inhibitor improved abnormal nuclear morphology. This case demonstrates that accumulation of prelamin A, independent of the loss of function of ZMPSTE24 metallopeptidase that catalyzes processing of prelamin A, can cause a progeroid disorder and that a cell biology assay could be used in precision medicine to identify a potential therapy.

Inherited selective cobalamin malabsorption in Komondor dogs associated with a CUBN splice site variant.

BACKGROUND: Three Komondor dogs in a small family and 3 sporadic cases exhibited a constellation of signs that included juvenile-onset of failure-to-thrive, inappetence, vomiting and/or diarrhea, and weakness. In each we documented dyshematopoiesis, increased anion gap, methylmalonic acidemia/-uria, and serum cobalamin deficiency. Urine protein electrophoresis demonstrated excretion of cubam ligands. All clinical signs and metabolic abnormalities, except proteinuria, were reversed by regular parenteral cobalamin administration. The pattern of occurrence and findings in the disorder suggested an autosomal recessive inheritance of cobalamin malabsorption with proteinuria, a condition in humans called Imerslund-Gräsbeck syndrome. The purpose of this study was to determine the molecular cause of this disorder in Komondors. RESULTS: Whole genome sequencing of two affected Komondor dogs of unknown relatedness and one parent and a clinically-normal littermate of an affected dog revealed a pathogenic single-base change in the CUBN intron 55 splice donor consensus sequence (NM_001003148.1: c.8746 + 1G > A) that was homozygous in affected dogs and heterozygous in the unaffected parents. Alleles of the variant co-segregated with alleles of the disease locus in the entire family and all more distantly-related sporadic cases. A population study using a simple allele-specific DNA test indicated mutant allele frequencies of 8.3 and 4.5% among North American and Hungarian Komondors, respectively. CONCLUSIONS: DNA testing can be used diagnostically in Komondors when clinical signs are suggestive of cobalamin deficiency or to inform Komondor breeders prospectively and prevent occurrence of future affected dogs. This represents the third cubilin variant causing inherited selective cobalamin malabsorption in a large animal ortholog of human Imerslund-Gräsbeck syndrome.

The proper splicing of RNAi factors is critical for pericentric heterochromatin assembly in fission yeast.

Heterochromatin preferentially assembles at repetitive DNA elements, playing roles in transcriptional silencing, recombination suppression, and chromosome segregation. The RNAi machinery is required for heterochromatin assembly in a diverse range of organisms. In fission yeast, RNA splicing factors are also required for pericentric heterochromatin assembly, and a prevailing model is that splicing factors provide a platform for siRNA generation independently of their splicing activity. Here, by screening the fission yeast deletion library, we discovered four novel splicing factors that are required for pericentric heterochromatin assembly. Sequencing total cellular RNAs from the strongest of these mutants, cwf14Δ, showed intron retention in mRNAs of several RNAi factors. Moreover, introducing cDNA versions of RNAi factors significantly restored pericentric heterochromatin in splicing mutants. We also found that mutations of splicing factors resulted in defective telomeric heterochromatin assembly and mis-splicing the mRNA of shelterin component Tpz1, and that replacement of tpz1+ with its cDNA partially rescued heterochromatin defects at telomeres in splicing mutants. Thus, proper splicing of RNAi and shelterin factors contributes to heterochromatin assembly at pericentric regions and telomeres.

Analysis of the ABCA4 genomic locus in Stargardt disease.

Autosomal recessive Stargardt disease (STGD1, MIM 248200) is caused by mutations in the ABCA4 gene. Complete sequencing of ABCA4 in STGD patients identifies compound heterozygous or homozygous disease-associated alleles in 65-70% of patients and only one mutation in 15-20% of patients. This study was designed to find the missing disease-causing ABCA4 variation by a combination of next-generation sequencing (NGS), array-Comparative Genome Hybridization (aCGH) screening, familial segregation and in silico analyses. The entire 140 kb ABCA4 genomic locus was sequenced in 114 STGD patients with one known ABCA4 exonic mutation revealing, on average, 200 intronic variants per sample. Filtering of these data resulted in 141 candidates for new mutations. Two variants were detected in four samples, two in three samples, and 20 variants in two samples, the remaining 117 new variants were detected only once. Multimodal analysis suggested 12 new likely pathogenic intronic ABCA4 variants, some of which were specific to (isolated) ethnic groups. No copy number variation (large deletions and insertions) was detected in any patient suggesting that it is a very rare event in the ABCA4 locus. Many variants were excluded since they were not conserved in non-human primates, were frequent in African populations and, therefore, represented ancestral, and not disease-associated, variants. The sequence variability in the ABCA4 locus is extensive and the non-coding sequences do not harbor frequent mutations in STGD patients of European-American descent. Defining disease-associated alleles in the ABCA4 locus requires exceptionally well characterized large cohorts and extensive analyses by a combination of various approaches.

The Glc7 phosphatase subunit of the cleavage and polyadenylation factor is essential for transcription termination on snoRNA genes.

Glc7, the yeast protein phosphatase 1, is a component of the cleavage and polyadenylation factor (CPF). Here we show that downregulation of Glc7, or its dissociation from CPF in the absence of CPF subunits Ref2 or Swd2, results in similar snoRNA termination defects. Overexpressing a C-terminal fragment of Sen1, a superfamily I helicase required for snoRNA termination, suppresses the growth and termination defects associated with loss of Swd2 or Ref2, but not Glc7. Suppression by Sen1 requires nuclear localization and direct interaction with Glc7, which can dephosphorylate Sen1 in vitro. The suppressing fragment, and in a similar manner full-length Sen1, copurifies with the snoRNA termination factors Nrd1 and Nab3, suggesting loss of Glc7 from CPF can be compensated by recruiting Glc7 to Nrd1-Nab3 through Sen1. Swd2 is also a subunit of the Set1c histone H3K4 methyltransferase complex and is required for its stability and optimal methyltransferase activity.

Tls1 regulates splicing of shelterin components to control telomeric heterochromatin assembly and telomere length.

Heterochromatin preferentially forms at repetitive DNA elements through RNAi-mediated targeting of histone-modifying enzymes. It was proposed that splicing factors interact with the RNAi machinery or regulate the splicing of repeat transcripts to directly participate in heterochromatin assembly. Here, by screening the fission yeast deletion library, we comprehensively identified factors required for telomeric heterochromatin assembly, including a novel gene tls1+. Purification of Tls1 and mass spectrometry analysis of its interacting proteins show that Tls1 associates with the spliceosome subunit Brr2. RNA sequencing analysis shows that the splicing of a subset of mRNAs are affected in tls1Δ cells, including mRNAs of shelterin components rap1+ and poz1+. Importantly, replacing rap1+ and poz1+ with their cDNAs significantly alleviated heterochromatin defects of tls1Δ cells, suggesting that the missplicing of shelterin components is the cause of such defects, and that splicing factors regulate telomeric heterochromatin through the proper splicing of heterochromatin factors. In addition to its role in telomeric heterochromatin assembly, Tls1-mediated splicing of shelterin mRNAs also regulates telomere length. Given that its human homologue C9ORF78 also associates with the spliceosome and is overexpressed in multiple cancer cell lines, our results suggest that C9ORF78 overexpression might alter the proper splicing of genes during cancer progression.

Comparative anatomy of chromosomal domains with imprinted and non-imprinted allele-specific DNA methylation.

Allele-specific DNA methylation (ASM) is well studied in imprinted domains, but this type of epigenetic asymmetry is actually found more commonly at non-imprinted loci, where the ASM is dictated not by parent-of-origin but instead by the local haplotype. We identified loci with strong ASM in human tissues from methylation-sensitive SNP array data. Two index regions (bisulfite PCR amplicons), one between the C3orf27 and RPN1 genes in chromosome band 3q21 and the other near the VTRNA2-1 vault RNA in band 5q31, proved to be new examples of imprinted DMRs (maternal alleles methylated) while a third, between STEAP3 and C2orf76 in chromosome band 2q14, showed non-imprinted haplotype-dependent ASM. Using long-read bisulfite sequencing (bis-seq) in 8 human tissues we found that in all 3 domains the ASM is restricted to single differentially methylated regions (DMRs), each less than 2kb. The ASM in the C3orf27-RPN1 intergenic region was placenta-specific and associated with allele-specific expression of a long non-coding RNA. Strikingly, the discrete DMRs in all 3 regions overlap with binding sites for the insulator protein CTCF, which we found selectively bound to the unmethylated allele of the STEAP3-C2orf76 DMR. Methylation mapping in two additional genes with non-imprinted haplotype-dependent ASM, ELK3 and CYP2A7, showed that the CYP2A7 DMR also overlaps a CTCF site. Thus, two features of imprinted domains, highly localized DMRs and allele-specific insulator occupancy by CTCF, can also be found in chromosomal domains with non-imprinted ASM. Arguing for biological importance, our analysis of published whole genome bis-seq data from hES cells revealed multiple genome-wide association study (GWAS) peaks near CTCF binding sites with ASM.

Isoform-level brain expression profiling of the spermidine/spermine N1-Acetyltransferase1 (SAT1) gene in major depression and suicide.

Low brain expression of the spermidine/spermine N-1 acetyltransferase (SAT1) gene, the rate-limiting enzyme involved in catabolism of polyamines that mediate the polyamine stress response (PSR), has been reported in depressed suicides. However, it is unknown whether this effect is associated with depression or with suicide and whether all or only specific isoforms expressed by SAT1, such as the primary 171 amino acid protein-encoding transcript (SSAT), or an alternative splice variant (SSATX) that is involved in SAT1 regulated unproductive splicing and transcription (RUST), are involved. We applied next generation sequencing (RNA-seq) to assess gene-level, isoform-level, and exon-level SAT1 expression differences between healthy controls (HC, N = 29), DSM-IV major depressive disorder suicides (MDD-S, N = 21) and MDD non-suicides (MDD, N = 9) in the dorsal lateral prefrontal cortex (Brodmann Area 9, BA9) of medication-free individuals postmortem. Using small RNA-seq, we also examined miRNA species putatively involved in SAT1 post-transcriptional regulation. A DSM-IV diagnosis was made by structured interview. Toxicology and history ruled out recent psychotropic medication. At the gene-level, we found low SAT1 expression in both MDD-S (vs. HC, p = 0.002) and MDD (vs. HC, p = 0.002). At the isoform-level, reductions in MDD-S (vs. HC) were most pronounced in four transcripts including SSAT and SSATX, while reductions in MDD (vs. HC) were pronounced in three transcripts, one of which was reduced in MDD relative to MDD-S (all p < 0.1 FDR corrected). We did not observe evidence for differential exon-usage (i.e. splicing) nor differences in miRNA expression. Results replicate the finding of low SAT1 brain expression in depressed suicides in an independent sample and implicate low SAT1 brain expression in MDD independent of suicide. Low expressions of both SSAT and SATX isoforms suggest that shared transcriptional mechanisms involved in RUST may account for low SAT1 brain expression in depressed suicides. Future studies are required to understand the functions and regulation of SAT1 isoforms, and how they relate to the pathogenesis of MDD and suicide.

LMNA variants cause cytoplasmic distribution of nuclear pore proteins in Drosophila and human muscle.

Mutations in the human LMNA gene, encoding A-type lamins, give rise to laminopathies, which include several types of muscular dystrophy. Here, heterozygous sequence variants in LMNA, which result in single amino-acid substitutions, were identified in patients exhibiting muscle weakness. To assess whether the substitutions altered lamin function, we performed in vivo analyses using a Drosophila model. Stocks were generated that expressed mutant forms of the Drosophila A-type lamin modeled after each variant. Larvae were used for motility assays and histochemical staining of the body-wall muscle. In parallel, immunohistochemical analyses were performed on human muscle biopsy samples from the patients. In control flies, muscle-specific expression of the wild-type A-type lamin had no apparent affect. In contrast, expression of the mutant A-type lamins caused dominant larval muscle defects and semi-lethality at the pupal stage. Histochemical staining of larval body wall muscle revealed that the mutant A-type lamin, B-type lamins, the Sad1p, UNC-84 domain protein Klaroid and nuclear pore complex proteins were mislocalized to the cytoplasm. In addition, cytoplasmic actin filaments were disorganized, suggesting links between the nuclear lamina and the cytoskeleton were disrupted. Muscle biopsies from the patients showed dystrophic histopathology and architectural abnormalities similar to the Drosophila larvae, including cytoplasmic distribution of nuclear envelope proteins. These data provide evidence that the Drosophila model can be used to assess the function of novel LMNA mutations and support the idea that loss of cellular compartmentalization of nuclear proteins contributes to muscle disease pathogenesis.

Tetratricopeptide repeat domain 7A (TTC7A) mutation in a newborn with multiple intestinal atresia and combined immunodeficiency.

In the past year, two centers reported autosomal recessive mutations in tetratricopeptide repeat domain 7A (TTC7A) gene in patients with multiple intestinal atresia and immunodeficiency. Here, we present clinical progress of an infant with multiple intestinal atresia and combined immunodeficiency who carries novel compound heterozygote mutations in TTC7A gene.

Detection of Constitutional Structural Variants by Optical Genome Mapping: A Multisite Study of Postnatal Samples.

Optical genome mapping is a high-resolution technology that can detect all types of structural variations in the genome. This second phase of a multisite study compares the performance of optical genome mapping and current standard-of-care methods for diagnostic testing of individuals with constitutional disorders, including neurodevelopmental impairments and congenital anomalies. Among the 627 analyses in phase 2, 405 were of retrospective samples supplied by five diagnostic centers in the United States and 94 were prospective samples collected over 18 months by two diagnostic centers (June 2021 to October 2022). Additional samples represented a family cohort to determine inheritance (n = 119) and controls (n = 9). Full concordance of results between optical genome mapping and one or more standard-of-care diagnostic tests was 98.6% (618/627), with partial concordance in an additional 1.1% (7/627).

Multisite Assessment of Optical Genome Mapping for Analysis of Structural Variants in Constitutional Postnatal Cases.

This study compares optical genome mapping (OGM) performed at multiple sites with current standard-of-care (SOC) methods used in clinical cytogenetics. This study included 50 negative controls and 359 samples from individuals (patients) with suspected genetic conditions referred for cytogenetic testing. OGM was performed using the Saphyr system and Bionano Access software version 1.7. Structural variants, including copy number variants, aneuploidy, and regions of homozygosity, were detected and classified according to American College of Medical Genetics and Genomics guidelines. Repeated expansions in FMR1 and contractions in facioscapulohumeral dystrophy 1 were also analyzed. OGM results were compared with SOC for technical concordance, clinical classification concordance, intrasite and intersite reproducibility, and ability to provide additional, clinically relevant information. Across five testing sites, 98.8% (404/409) of samples yielded successful OGM data for analysis and interpretation. Overall, technical concordance for OGM to detect previously reported SOC results was 99.5% (399/401). The blinded analysis and variant classification agreement between SOC and OGM was 97.6% (364/373). Replicate analysis of 130 structural variations was 100% concordant. On the basis of this demonstration of the analytic validity and clinical utility of OGM by this multisite assessment, the authors recommend this technology as an alternative to existing SOC tests for rapid detection and diagnosis in postnatal constitutional disorders.

Discovery of the BMPR1A promoter and germline mutations that cause juvenile polyposis.

Juvenile polyposis (JP) is an autosomal dominant hamartomatous polyposis syndrome where affected individuals are predisposed to colorectal and upper gastrointestinal cancer. Forty-five percent of JP patients have mutations or deletions involving the coding regions of SMAD4 and BMPR1A, but the genetic basis of other cases is unknown. We set out to identify the JP gene in a large kindred having 10 affected members without SMAD4 or BMPR1A coding sequence mutations or deletions. We found a germline deletion segregating in all affected members, mapping 119 kb upstream of the coding region of BMPR1A by multiplex ligation-dependent probe amplification and comparative genomic hybridization. To further understand the genomic structure of BMPR1A, we performed 5' RACE from lymphoblastoid cell lines and normal colon tissue, which revealed four non-coding (NC) exons and two putative promoters. Further analysis of this deletion showed that it encompassed 12 433 bp, including one promoter and NC exon. The activities of each promoter and deletion constructs were evaluated by luciferase assays, and the stronger promoter sequence analyzed for changes in JP patients without SMAD4 or BMPR1A alterations. A total of 6 of 65 JP probands were found to have mutations affecting this promoter. All probands examined had diminished BMPR1A protein by ELISA, and all promoter mutations but one led to significantly reduced luciferase activity relative to the wild-type promoter reporter. We conclude that we have identified the promoter for BMPR1A, in which mutations may be responsible for as many as 10% of JP cases with unknown mutations.

2q24 deletions: further characterization of clinical findings and their relation to the SCN cluster.

As the resolution of molecular cytogenetic methods continues to improve, it has become increasingly possible to refine genotype-phenotype correlations based upon gene involvement. We report three new patients with nonrecurrent deletions involving subbands of 2q24. These patients were referred for evaluation of developmental delay, but were found to have unique, nonoverlapping clinical features. Patient 1 presented with infantile seizures, microcephaly, and brain anomalies, along with facial dysmorphism, growth retardation, neuromuscular scoliosis, and later with developmental regression. Array comparative genomic hybridization (aCGH) detected an 8 Mb interstitial deletion encompassing the neuronal sodium channel (SCN) gene cluster. Patient 2 presented with growth retardation, congenital heart defect, and hypotonia. Patient 3 presented with developmental delay and behavioral problems. Patients 2 and 3 had no history of seizures, microcephaly, or brain anomalies and were found to have deletions of 2q24, ∼8 Mb and <500 kb respectively, centromeric to and outside the SCN cluster. It has been demonstrated that mutations and copy number variants (CNVs) affecting the SCN gene cluster result in severe, early-onset seizures. It is however, less clear whether haploinsufficiency of regions outside the SCN cluster may result in phenotypically recognizable and clinically significant features. We discuss additional dosage sensitive genes that may exist outside the SCN cluster. Our and published data indicate that 2q24 deletions not involving the SCN cluster are associated with fewer neurobehavioral problems, but may predispose to congenital malformations.

Aerobic exercise elicits clinical adaptations in myotonic dystrophy type 1 patients independently of pathophysiological changes.

BackgroundMyotonic dystrophy type 1 (DM1) is a complex life-limiting neuromuscular disorder characterized by severe skeletal muscle atrophy, weakness, and cardiorespiratory defects. Exercised DM1 mice exhibit numerous physiological benefits that are underpinned by reduced CUG foci and improved alternative splicing. However, the efficacy of physical activity in patients is unknown.MethodsEleven genetically diagnosed DM1 patients were recruited to examine the extent to which 12 weeks of cycling can recuperate clinical and physiological metrics. Furthermore, we studied the underlying molecular mechanisms through which exercise elicits benefits in skeletal muscle of DM1 patients.RESULTSDM1 was associated with impaired muscle function, fitness, and lung capacity. Cycling evoked several clinical, physical, and metabolic advantages in DM1 patients. We highlight that exercise-induced molecular and cellular alterations in patients do not conform with previously published data in murine models and propose a significant role of mitochondrial function in DM1 pathology. Finally, we discovered a subset of small nucleolar RNAs (snoRNAs) that correlated to indicators of disease severity.ConclusionWith no available cures, our data support the efficacy of exercise as a primary intervention to partially mitigate the clinical progression of DM1. Additionally, we provide evidence for the involvement of snoRNAs and other noncoding RNAs in DM1 pathophysiology.Trial registrationThis trial was approved by the HiREB committee (no. 7901) and registered under ClinicalTrials.gov (NCT04187482).FundingNeil and Leanne Petroff. Canadian Institutes of Health Research Foundation (no. 143325).

The genetic legacy of the Hunyadi descendants.

The Hunyadi family is one of the most influential families in the history of Central Europe in the 14th-16th centuries. The family's prestige was established by Johannes Hunyadi, a Turk-beater who rose to the position of governor of the Kingdom of Hungary. His second son, Matthias Hunyadi, became the elected ruler of the Kingdom of Hungary in 1458. The Hunyadi family had unknown origin. Moreover, Matthias failed to found a dynasty because of lacking a legitimate heir and his illegitimate son Johannes Corvinus was unable to obtain the crown. His grandson, Christophorus Corvinus, died in childhood, thus the direct male line of the family ended. In the framework of on interdisciplinary research, we have determined the whole genome sequences of Johannes Corvinus and Christophorus Corvinus by next-generation sequencing technology. Both of them carried the Y-chromosome haplogroup is E1b1b1a1b1a6a1c ∼, which is widespread in Eurasia. The father-son relationship was verified using the classical STR method and whole genome data. Christophorus Corvinus belongs to the rare, sporadically occurring T2c1+146 mitochondrial haplogroup, most frequent around the Mediterranean, while his father belongs to the T2b mitochondrial haplogroup, widespread in Eurasia, both are consistent with the known origin of the mothers. Archaeogenomic analysis indicated that the Corvinus had an ancient European genome composition. Based on the reported genetic data, it will be possible to identify all the other Hunyadi family member, whose only known grave site is known, but who are resting assorted with several other skeletons.

The genetic origin of Huns, Avars, and conquering Hungarians.

Huns, Avars, and conquering Hungarians were migration-period nomadic tribal confederations that arrived in three successive waves in the Carpathian Basin between the 5th and 9th centuries. Based on the historical data, each of these groups are thought to have arrived from Asia, although their exact origin and relation to other ancient and modern populations have been debated. Recently, hundreds of ancient genomes were analyzed from Central Asia, Mongolia, and China, from which we aimed to identify putative source populations for the above-mentioned groups. In this study, we have sequenced 9 Hun, 143 Avar, and 113 Hungarian conquest period samples and identified three core populations, representing immigrants from each period with no recent European ancestry. Our results reveal that this "immigrant core" of both Huns and Avars likely originated in present day Mongolia, and their origin can be traced back to Xiongnus (Asian Huns), as suggested by several historians. On the other hand, the "immigrant core" of the conquering Hungarians derived from an earlier admixture of Mansis, early Sarmatians, and descendants of late Xiongnus. We have also shown that a common "proto-Ugric" gene pool appeared in the Bronze Age from the admixture of Mezhovskaya and Nganasan people, supporting genetic and linguistic data. In addition, we detected shared Hun-related ancestry in numerous Avar and Hungarian conquest period genetic outliers, indicating a genetic link between these successive nomadic groups. Aside from the immigrant core groups, we identified that the majority of the individuals from each period were local residents harboring "native European" ancestry.