Search

The genomes of present-day humans outside Africa originated almost entirely from a single out-migration ~ 50,000-70,000 years ago, followed by mixture with Neanderthals contributing ~ 2% to all non-Africans. However, the details of this initial migration remain poorly understood because no ancient DNA analyses are available from this key time period, and interpretation of present-day autosomal data is complicated due to subsequent population movements/reshaping. One locus, however, does retain male-specific information from this early period: the Y chromosome, where a detailed calibrated phylogeny has been constructed. Three present-day Y lineages were carried by the initial migration: the rare haplogroup D, the moderately rare C, and the very common FT lineage which now dominates most non-African populations. Here, we show that phylogenetic analyses of haplogroup C, D and FT sequences, including very rare deep-rooting lineages, together with phylogeographic analyses of ancient and present-day non-African Y chromosomes, all point to East/Southeast Asia as the origin 50,000-55,000 years ago of all known surviving non-African male lineages (apart from recent migrants). This observation contrasts with the expectation of a West Eurasian origin predicted by a simple model of expansion from a source near Africa, and can be interpreted as resulting from extensive genetic drift in the initial population or replacement of early western Y lineages from the east, thus informing and constraining models of the initial expansion.

We describe the variation in copy number of a ~ 10 kb region overlapping the long intergenic noncoding RNA (lincRNA) gene, TTTY22, within the IR3 inverted repeat on the short arm of the human Y chromosome, leading to individuals with 0-3 copies of this region in the general population. Variation of this CNV is common, with 266 individuals having 0 copies, 943 (including the reference sequence) having 1, 23 having 2 copies, and two having 3 copies, and was validated by breakpoint PCR, fibre-FISH, and 10× Genomics Chromium linked-read sequencing in subsets of 1234 individuals from the 1000 Genomes Project. Mapping the changes in copy number to the phylogeny of these Y chromosomes previously established by the Project identified at least 20 mutational events, and investigation of flanking paralogous sequence variants showed that the mutations involved flanking sequences in 18 of these, and could extend over > 30 kb of DNA. While either gene conversion or double crossover between misaligned sister chromatids could formally explain the 0-2 copy events, gene conversion is the more likely mechanism, and these events include the longest non-allelic gene conversion reported thus far. Chromosomes with three copies of this CNV have arisen just once in our data set via another mechanism: duplication of 420 kb that places the third copy 230 kb proximal to the existing proximal copy. Our results establish gene conversion as a previously under-appreciated mechanism of generating copy number changes in humans and reveal the exceptionally large size of the conversion events that can occur.

We present 42 new Y-chromosomal sequences from diverse Indian tribal and non-tribal populations, including the Jarawa and Onge from the Andaman Islands, which are analysed within a calibrated Y-chromosomal phylogeny incorporating South Asian (in total 305 individuals) and worldwide (in total 1286 individuals) data from the 1000 Genomes Project. In contrast to the more ancient ancestry in the South than in the North that has been claimed, we detected very similar coalescence times within Northern and Southern non-tribal Indian populations. A closest neighbour analysis in the phylogeny showed that Indian populations have an affinity towards Southern European populations and that the time of divergence from these populations substantially predated the Indo-European migration into India, probably reflecting ancient shared ancestry rather than the Indo-European migration, which had little effect on Indian male lineages. Among the tribal populations, the Birhor (Austro-Asiatic-speaking) and Irula (Dravidian-speaking) are the nearest neighbours of South Asian non-tribal populations, with a common origin in the last few millennia. In contrast, the Riang (Tibeto-Burman-speaking) and Andamanese have their nearest neighbour lineages in East Asia. The Jarawa and Onge shared haplogroup D lineages with each other within the last ~7000 years, but had diverged from Japanese haplogroup D Y-chromosomes ~53000 years ago, most likely by a split from a shared ancestral population. This analysis suggests that Indian populations have complex ancestry which cannot be explained by a single expansion model.

Y-chromosomal variation in West Asian populations has so far been studied in less detail than in the neighboring Europe. Here, we analyzed 598 Y-chromosomes from two West Asian subregions-Transcaucasia and the Armenian plateau-using 40 Y-SNPs and 17 Y-STRs and combined them with previously published data from the region. The West Asian populations fell into two clusters: upland populations from the Anatolian, Armenian and Iranian plateaus, and lowland populations from the Levant, Mesopotamia and the Arabian Peninsula. This geographic subdivision corresponds with the linguistic difference between Indo-European and Turkic speakers, on the one hand, and Semitic speakers, on the other. This subdivision could be traced back to the Neolithic epoch, when upland populations from the Anatolian and Iranian plateaus carried similar haplogroup spectra but did not overlap with lowland populations from the Levant. We also found that the initial gene pool of the Armenian motherland population has been well preserved in most groups of the Armenian Diaspora. In view of the contribution of West Asians to the autosomal gene pool of the steppe Yamnaya archaeological culture, we sequenced a large portion of the Y-chromosome in haplogroup R1b samples from present-day East European steppe populations. The ancient Yamnaya samples are located on the "eastern" R-GG400 branch of haplogroup R1b-L23, showing that the paternal descendants of the Yamnaya still live in the Pontic steppe and that the ancient Yamnaya population was not an important source of paternal lineages in present-day West Europeans.

High-altitude adaptation in Tibetans is influenced by introgression of a 32.7-kb haplotype from the Denisovans, an extinct branch of archaic humans, lying within the endothelial PAS domain protein 1 (EPAS1), and has also been reported in Sherpa. We genotyped 19 variants in this genomic region in 1507 Eurasian individuals, including 1188 from Bhutan and Nepal residing at altitudes between 86 and 4550 m above sea level. Derived alleles for five SNPs characterizing the core Denisovan haplotype (AGGAA) were present at high frequency not only in Tibetans and Sherpa, but also among many populations from the Himalayas, showing a significant correlation with altitude (Spearman's correlation coefficient = 0.75, p value 3.9 × 10(-11)). Seven East- and South-Asian 1000 Genomes Project individuals shared the Denisovan haplotype extending beyond the 32-kb region, enabling us to refine the haplotype structure and identify a candidate regulatory variant (rs370299814) that might be interacting in an additive manner with the derived G allele of rs150877473, the variant previously associated with high-altitude adaptation in Tibetans. Denisovan-derived alleles were also observed at frequencies of 3-14% in the 1000 Genomes Project African samples. The closest African haplotype is, however, separated from the Asian high-altitude haplotype by 22 mutations whereas only three mutations, including rs150877473, separate the Asians from the Denisovan, consistent with distant shared ancestry for African and Asian haplotypes and Denisovan adaptive introgression.

We have assessed copy number variation (CNV) in the male-specific part of the human Y chromosome discovered by array comparative genomic hybridization (array-CGH) in 411 apparently healthy UK males, and validated the findings using SNP genotype intensity data available for 149 of them. After manual curation taking account of the complex duplicated structure of Y-chromosomal sequences, we discovered 22 curated CNV events considered validated or likely, mean 0.93 (range 0-4) per individual. 16 of these were novel. Curated CNV events ranged in size from <1 kb to >3 Mb, and in frequency from 1/411 to 107/411. Of the 24 protein-coding genes or gene families tested, nine showed CNV. These included a large duplication encompassing the AMELY and TBL1Y genes that probably has no phenotypic effect, partial deletions of the TSPY cluster and AZFc region that may influence spermatogenesis, and other variants with unknown functional implications, including abundant variation in the number of RBMY genes and/or pseudogenes, and a novel complex duplication of two segments overlapping the AZFa region and including the 3' end of the UTY gene.

Some individuals with a particular disease-causing mutation or genotype fail to express most if not all features of the disease in question, a phenomenon that is known as 'reduced (or incomplete) penetrance'. Reduced penetrance is not uncommon; indeed, there are many known examples of 'disease-causing mutations' that fail to cause disease in at least a proportion of the individuals who carry them. Reduced penetrance may therefore explain not only why genetic diseases are occasionally transmitted through unaffected parents, but also why healthy individuals can harbour quite large numbers of potentially disadvantageous variants in their genomes without suffering any obvious ill effects. Reduced penetrance can be a function of the specific mutation(s) involved or of allele dosage. It may also result from differential allelic expression, copy number variation or the modulating influence of additional genetic variants in cis or in trans. The penetrance of some pathogenic genotypes is known to be age- and/or sex-dependent. Variable penetrance may also reflect the action of unlinked modifier genes, epigenetic changes or environmental factors. At least in some cases, complete penetrance appears to require the presence of one or more genetic variants at other loci. In this review, we summarize the evidence for reduced penetrance being a widespread phenomenon in human genetics and explore some of the molecular mechanisms that may help to explain this enigmatic characteristic of human inherited disease.

The evolutionary history of variation in the human Rh blood group system, determined by variants in the RHD and RHCE genes, has long been an unresolved puzzle in human genetics. Prior to medical treatments and interventions developed in the last century, the D-positive (RhD positive) children of D-negative (RhD negative) women were at risk for hemolytic disease of the newborn, if the mother produced anti-D antibodies following sensitization to the blood of a previous D-positive child. Given the deleterious fitness consequences of this disease, the appreciable frequencies in European populations of the responsible RHD gene deletion variant (for example, 0.43 in our study) seem surprising. In this study, we used new molecular and genomic data generated from four HapMap population samples to test the idea that positive selection for an as-of-yet unknown fitness benefit of the RHD deletion may have offset the otherwise negative fitness effects of hemolytic disease of the newborn. We found no evidence that positive natural selection affected the frequency of the RHD deletion. Thus, the initial rise to intermediate frequency of the RHD deletion in European populations may simply be explained by genetic drift/founder effect, or by an older or more complex sweep that we are insufficiently powered to detect. However, our simulations recapitulate previous findings that selection on the RHD deletion is frequency dependent and weak or absent near 0.5. Therefore, once such a frequency was achieved, it could have been maintained by a relatively small amount of genetic drift. We unexpectedly observed evidence for positive selection on the C allele of RHCE in non-African populations (on chromosomes with intact copies of the RHD gene) in the form of an unusually high F( ST ) value and the high frequency of a single haplotype carrying the C allele. RhCE function is not well understood, but the C/c antigenic variant is clinically relevant and can result in hemolytic disease of the newborn, albeit much less commonly and severely than that related to the D-negative blood type. Therefore, the potential fitness benefits of the RHCE C allele are currently unknown but merit further exploration.

We have investigated whether regions of the genome showing signs of positive selection in scans based on haplotype structure also show evidence of positive selection when sequence-based tests are applied, whether the target of selection can be localized more precisely, and whether such extra evidence can lead to increased biological insights. We used two tools: simulations under neutrality or selection, and experimental investigation of two regions identified by the HapMap2 project as putatively selected in human populations. Simulations suggested that neutral and selected regions should be readily distinguished and that it should be possible to localize the selected variant to within 40 kb at least half of the time. Re-sequencing of two ~300 kb regions (chr4:158Mb and chr10:22Mb) lacking known targets of selection in HapMap CHB individuals provided strong evidence for positive selection within each and suggested the micro-RNA gene hsa-miR-548c as the best candidate target in one region, and changes in regulation of the sperm protein gene SPAG6 in the other.

We have surveyed 15 high-altitude adaptation candidate genes for signals of positive selection in North Caucasian highlanders using targeted re-sequencing. A total of 49 unrelated Daghestani from three ethnic groups (Avars, Kubachians, and Laks) living in ancient villages located at around 2,000 m above sea level were chosen as the study population. Caucasian (Adygei living at sea level, N = 20) and CEU (CEPH Utah residents with ancestry from northern and western Europe; N = 20) were used as controls. Candidate genes were compared with 20 putatively neutral control regions resequenced in the same individuals. The regions of interest were amplified by long-PCR, pooled according to individual, indexed by adding an eight-nucleotide tag, and sequenced using the Illumina GAII platform. 1,066 SNPs were called using false discovery and false negative thresholds of ~6%. The neutral regions provided an empirical null distribution to compare with the candidate genes for signals of selection. Two genes stood out. In Laks, a non-synonymous variant within HIF1A already known to be associated with improvement in oxygen metabolism was rediscovered, and in Kubachians a cluster of 13 SNPs located in a conserved intronic region within EGLN1 showing high population differentiation was found. These variants illustrate both the common pathways of adaptation to high altitude in different populations and features specific to the Daghestani populations, showing how even a mildly hypoxic environment can lead to genetic adaptation.

The caste system has persisted in Indian Hindu society for around 3,500 years. Like the Y chromosome, caste is defined at birth, and males cannot change their caste. In order to investigate the genetic consequences of this system, we have analysed male-lineage variation in a sample of 227 Indian men of known caste, 141 from the Jaunpur district of Uttar Pradesh and 86 from the rest of India. We typed 131 Y-chromosomal binary markers and 16 microsatellites. We find striking evidence for male substructure: in particular, Brahmins and Kshatriyas (but not other castes) from Jaunpur each show low diversity and the predominance of a single distinct cluster of haplotypes. These findings confirm the genetic isolation and drift within the Jaunpur upper castes, which are likely to result from founder effects and social factors. In the other castes, there may be either larger effective population sizes, or less strict isolation, or both.

Arab forces conquered the Indus Delta region in 711 AD: and, although a Muslim state was established there, their influence was barely felt in the rest of South Asia at that time. By the end of the tenth century, Central Asian Muslims moved into India from the northwest and expanded throughout the subcontinent. Muslim communities are now the largest minority religion in India, comprising more than 138 million people in a predominantly Hindu population of over one billion. It is unclear whether the Muslim expansion in India was a purely cultural phenomenon or had a genetic impact on the local population. To address this question from a male perspective, we typed eight microsatellite loci and 16 binary markers from the Y chromosome in 246 Muslims from Andhra Pradesh, and compared them to published data on 4,204 males from East Asia, Central Asia, other parts of India, Sri Lanka, Pakistan, Iran, the Middle East, Turkey, Egypt and Morocco. We find that the Muslim populations in general are genetically closer to their non-Muslim geographical neighbors than to other Muslims in India, and that there is a highly significant correlation between genetics and geography (but not religion). Our findings indicate that, despite the documented practice of marriage between Muslim men and Hindu women, Islamization in India did not involve large-scale replacement of Hindu Y chromosomes. The Muslim expansion in India was predominantly a cultural change and was not accompanied by significant gene flow, as seen in other places, such as China and Central Asia.

Previous studies of human Y-chromosomal single-nucleotide polymorphisms (Y-SNPs) established a link between the extant Y-SNP haplogroup distribution and the prehistoric demography of Europe. By contrast, our analysis of seven rapidly evolving Y-chromosomal short tandem repeat loci (Y-STRs) in over 12,700 samples from 91 different locations in Europe reveals a signature of more recent historic events, not previously detected by other genetic markers. Cluster analysis based upon molecular variance yields two clearly identifiable sub-clusters of Western and Eastern European Y-STR haplotypes, and a diverse transition zone in central Europe, where haplotype spectra change more rapidly with longitude than with latitude. This and other observed patterns of Y-STR similarity may plausibly be related to particular historical incidents, including, for example, the expansion of the Franconian and Ottoman Empires. We conclude that Y-STRs may be capable of resolving male genealogies to an unparalleled degree and could therefore provide a useful means to study local population structure and recent demographic history.

We have analyzed eight Y-chromosomal binary markers (YAP, RPS4Y(711), M9, M175, LINE1, SRY(+465), 47z, and M95) and three Y-STR markers (DYS390, DYS391, and DYS393) in 738 males from 11 ethnic groups in east Asia in order to study the male lineage history of Korea. Haplogroup DE-YAP was found at a high frequency only in Japan but was also present at low frequencies in northeast Asia, including 2.5% in Korea, suggesting a northern origin for these chromosomes. Haplogroup C-RPS4Y(711) was present in Korea and Manchuria at moderate frequencies: higher than in populations from southeast Asia, but lower than those in the northeast, which may imply a northern Asian expansion of these lineages, perhaps from Mongolia or Siberia. The major Y-chromosomal expansions in east Asia were those of haplogroup O-M175 (and its sublineages). This haplogroup is likely to have originated in southern east Asia and subsequently expanded to all of east Asia. The moderate frequency of one sublineage in the Koreans, haplogroup O-LINE1 (12.5%), could be a result of interaction with Chinese populations. The age of another sublineage, haplogroup O-SRY(+465), and Y-STR haplotype diversity provide evidence for relatively recent male migration, originally from China, through Korea into Japan. In conclusion, the distribution pattern of Y-chromosomal haplogroups reveals the complex origin of the Koreans, resulting from genetic contributions involving the northern Asian settlement and range expansions mostly from southern-to-northern China.

Migration and admixture history of populations have always been curious and an interesting theme. The West Coast of India harbours a rich diversity, bestowing various ethno-linguistic groups, with many of them having well-documented history of migrations. The Roman Catholic is one such distinct group, whose origin was much debated. While some historians and anthropologists relating them to ancient group of Gaud Saraswat Brahmins, others relating them for being members of the Jews Lost Tribes in the first Century migration to India. Historical records suggests that this community was later forcibly converted to Christianity by the Portuguese in Goa during the Sixteenth Century. Till date, no genetic study was done on this group to infer their origin and genetic affinity. Hence, we analysed 110 Roman Catholics from three different locations of West Coast of India including Goa, Kumta and Mangalore using both uniparental and autosomal markers to understand their genetic history. We found that the Roman Catholics have close affinity with the Indo-European linguistic groups, particularly Brahmins. Additionally, we detected genetic signal of Jews in the linkage disequilibrium-based admixture analysis, which was absent in other Indo-European populations, who are inhabited in the same geographical regions. Haplotype-based analysis suggests that the Roman Catholics consist of South Asian-specific ancestry and showed high drift. Ancestry-specific historical population size estimation points to a possible bottleneck around the time of Goan inquisition (fifteenth century). Analysis of the Roman Catholics data along with ancient DNA data of Neolithic and bronze age revealed that the Roman Catholics fits well in a basic model of ancient ancestral composition, typical of most of the Indo-European caste groups of India. Mitochondrial DNA (mtDNA) analysis suggests that most of the Roman Catholics have aboriginal Indian maternal genetic ancestry; while the Y chromosomal DNA analysis indicates high frequency of R1a lineage, which is predominant in groups with higher ancestral North Indian (ANI) component. Therefore, we conclude that the Roman Catholics of Goa, Kumta and Mangalore regions are the remnants of very early lineages of Brahmin community of India, having Indo-Europeans genetic affinity along with cryptic Jewish admixture, which needs to be explored further., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Primary immunodeficiencies (PID) are a diverse group of genetic disorders caused by inadequate development and function of immune system. Identifying genetic etiology is important for genetic counselling and treatment decisions. Clinical relevance of genetic variants is a complex problem depending on gene-specific and variant specific genotype-phenotype interactions. To address this challenge, we aimed to characterize the pathogenic landscape of PID genes by combining the analysis of germline variations reported in ClinVar and HGMD ®  and identification of damaging variations available in dbSNP. We generated a joint ClinVar/HGMD database, which included 111,940 variants, among them 32,452 were classified as pathogenic/likely pathogenic. From a total of 5,415,794 bi- or multiallelic variants in PID genes recorded in dbSNP, we retrieved 38,291 high impact (HI) biallelic variants with presumably disruptive impact in the protein, of them 25,500 variants were not present in ClinVar/HGMD. Using a functional prediction algorithm, we additionally identified 28,507 deleterious and 56,016 neutral missense variants among dbSNP variants and created a collection of damaging and neutral variations in PID genes, not currently present in ClinVar/HGMD, with their allele frequencies and mappings to protein domains. The distribution of pathogenic variants from ClinVar/HGMD, HI variants and deleterious missense variants from dbSNP was analyzed in the context of hereditary pattern and gene specific metrics, such as pLI and haploinsufficiency. Our report summarized data on complex gene-specific variability in PID genes and might be useful for the identification of the most promising variants and gene regions for further study., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

The suite of phenotypic diversity across geographically distributed human populations is the outcome of genetic drift, gene flow, and natural selection throughout human evolution. Human genetic variation underlying local biological adaptations to selective pressures is incompletely characterized. With the emergence of population genetics modeling of large-scale genomic data derived from diverse populations, scientists are able to map signatures of natural selection in the genome in a process known as selection mapping. Inferred selection signals further can be used to identify candidate functional alleles that underlie putative adaptive phenotypes. Phenotypic association, fine mapping, and functional experiments facilitate the identification of candidate adaptive alleles. Functional investigation of candidate adaptive variation using novel techniques in molecular biology is slowly beginning to unravel how selection signals translate to changes in biology that underlie the phenotypic spectrum of our species. In addition to informing evolutionary hypotheses of adaptation, the discovery and functional annotation of adaptive alleles also may be of clinical significance. While selection mapping efforts in non-European populations are growing, there remains a stark under-representation of diverse human populations in current public genomic databases, of both clinical and non-clinical cohorts. This lack of inclusion limits the study of human biological variation. Identifying and functionally validating candidate adaptive alleles in more global populations is necessary for understanding basic human biology and human disease.

Aims/hypothesis: Glucagon-expressing pancreatic alpha cells have attracted much attention for their plasticity to transdifferentiate into insulin-producing beta cells; however, it remains unclear precisely when, and from where, alpha cells emerge and what regulates alpha cell fate. We therefore explored the spatial and transcriptional heterogeneity of alpha cell differentiation using a novel time-resolved reporter system., Methods: We established the mouse model, 'Gcg-Timer', in which newly generated alpha cells can be distinguished from more-differentiated cells by their fluorescence. Fluorescence imaging and transcriptome analysis were performed with Gcg-Timer mice during the embryonic and postnatal stages., Results: Fluorescence imaging and flow cytometry demonstrated that green fluorescence-dominant cells were present in Gcg-Timer mice at the embryonic and neonatal stages but not after 1 week of age, suggesting that alpha cell neogenesis occurs during embryogenesis and early neonatal stages under physiological conditions. Transcriptome analysis of Gcg-Timer embryos revealed that the mRNAs related to angiogenesis were enriched in newly generated alpha cells. Histological analysis revealed that some alpha cells arise close to the pancreatic ducts, whereas the others arise away from the ducts and adjacent to the blood vessels. Notably, when the glucagon signal was suppressed by genetic ablation or by chemicals, such as neutralising glucagon antibody, green-dominant cells emerged again in adult mice., Conclusions/interpretation: Novel time-resolved analysis with Gcg-Timer reporter mice uncovered spatiotemporal features of alpha cell neogenesis that will enhance our understanding of cellular identity and plasticity within the islets., Data Availability: Raw and processed RNA sequencing data for this study has been deposited in the Gene Expression Omnibus under accession number GSE229090., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Diabetes mellitus (DM) has become a significant health concern with an increasing rate of morbidity and mortality worldwide. India ranks second in the number of diabetes cases in the world. The increasing burden of DM can be explained by genetic predisposition of Indians to type 2 diabetes mellitus (T2DM) coupled with rapid urbanization and socio-economic development in the last 3 decades leading to drastic changes in lifestyle. Environment and lifestyle changes contribute to T2DM development by altering epigenetic processes such as DNA methylation, histone post-translational modifications, and long non-coding RNAs, all of which regulate chromatin structure and gene expression. Although the genetic predisposition of Indians to T2DM is well established, how environmental and genetic factors interact and lead to T2DM is not well understood. In this review, we discuss the prevalence of diabetes and its complications across different states in India and how various risk factors contribute to its pathogenesis. The review also highlights the role of genetic predisposition among the Indian population and epigenetic factors involved in the etiology of diabetes. Lastly, we review current treatments and emphasize the knowledge gap with respect to genetic and epigenetic factors in the Indian context. Further understanding of the genetic and epigenetic determinants will help in risk prediction and prevention as well as therapeutic interventions, which will improve the clinical management of diabetes and associated macro- and micro-vascular complications., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Brain arteriovenous malformation (BAVM) is a rare but serious cerebrovascular disease whose pathogenesis has not been fully elucidated. Studies have found that epigenetic regulation, genetic variation and their signaling pathways, immune inflammation, may be the cause of BAVM the main reason. This review comprehensively analyzes the key pathways and inflammatory factors related to BAVMs, and explores their interplay with epigenetic regulation and genetics. Studies have found that epigenetic regulation such as DNA methylation, non-coding RNAs and m6A RNA modification can regulate endothelial cell proliferation, apoptosis, migration and damage repair of vascular malformations through different target gene pathways. Gene defects such as KRAS, ACVRL1 and EPHB4 lead to a disordered vascular environment, which may promote abnormal proliferation of blood vessels through ERK, NOTCH, mTOR, Wnt and other pathways. PDGF-B and PDGFR-β were responsible for the recruitment of vascular adventitial cells and smooth muscle cells in the extracellular matrix environment of blood vessels, and played an important role in the pathological process of BAVM. Recent single-cell sequencing data revealed the diversity of various cell types within BAVM, as well as the heterogeneous expression of vascular-associated antigens, while neutrophils, macrophages and cytokines such as IL-6, IL-1, TNF-α, and IL-17A in BAVM tissue were significantly increased. Currently, there are no specific drugs targeting BAVMs, and biomarkers for BAVM formation, bleeding, and recurrence are lacking clinically. Therefore, further studies on molecular biological mechanisms will help to gain insight into the pathogenesis of BAVM and develop potential therapeutic strategies., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Trypsinogen (PRSS1, PRSS2) copy number gains and regulatory variants have both been proposed to elevate pancreatitis risk through a gene dosage effect (i.e., by increasing the expression of wild-type protein). However, to date, their impact on pancreatitis risk has not been thoroughly evaluated whilst the underlying pathogenic mechanisms remain to be explicitly investigated in mouse models. Genetic studies of the rare trypsinogen duplication and triplication copy number variants (CNVs), and the common rs10273639C variant, were collated from PubMed and/or ClinVar. Mouse studies that analyzed the influence of a transgenically expressed wild-type human PRSS1 or PRSS2 gene on the development of pancreatitis were identified from PubMed. The genetic effects of the different risk genotypes, in terms of odds ratios, were calculated wherever appropriate. The genetic effects of the rare trypsinogen duplication and triplication CNVs were also evaluated by reference to their associated disease subtypes. We demonstrate a positive correlation between increased trypsinogen gene dosage and pancreatitis risk in the context of the rare duplication and triplication CNVs, and between the level of trypsinogen expression and disease risk in the context of the heterozygous and homozygous rs10273639C-tagged genotypes. We retrospectively identify three mouse transgenic studies that are informative in relation to the pathogenic mechanism underlying the trypsinogen gene dosage effect in pancreatitis. Trypsinogen gene dosage correlates with pancreatitis risk across genetic and transgenic studies, highlighting the fundamental role of dysregulated expression of wild-type trypsinogen in the etiology of pancreatitis. Specifically downregulating trypsinogen expression in the pancreas may serve as a potential therapeutic and/or prevention strategy for pancreatitis., (© 2022. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Epistatic interactions complicate the identification of variants involved in phenotypic effect. In-depth knowledge in modifiers and in pathogenic variants would benefit the mechanistic studies on the genetic basis of complex traits. We systematically compared the modifier variants which have evidence of modifier effect with the pathogenic variants from multiple angles. Our study found that genomic loci of modifier variations differ from pathogenic loci in many aspects, such as population genetics statistics, epigenetic features, evolutionary characteristics and functional properties of the variations. Genes containing modifier variation(s) exhibit higher probability of being haploinsufficient and higher probability of recessive disease causation, and they are relatively more important in network communication. Furthermore, we reinforced that co-expression analysis is an effective methodology to predict functional associations between modifier genes and their potential target genes. In many aspects, we detected statistically significant differences between modifier variants/genes and pathogenic variants/genes, and investigated relationships between modifiers and their potential targets. Our results offer some actionable insights that may provide appropriate guidelines to clinical genetics and researchers to elucidate the molecular mechanism underlying the human phenotypic variation., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Structural variation includes a change in copy number, orientation, or location of a part of the genome. Copy number variants (CNVs) are a common cause of genetic hearing loss, comprising nearly 20% of diagnosed cases. While large deletions involving the gene STRC are the most common pathogenic CNVs, a significant proportion of known hearing loss genes also contain pathogenic CNVs. In this review, we provide an overview of currently used methods for detection of CNVs in genes known to cause hearing loss including molecular techniques such as multiplex ligation probe amplification (MLPA) and digital droplet polymerase chain reaction (ddPCR), array-CGH and single-nucleotide polymorphism (SNP) arrays, as well as techniques for detection of CNVs using next-generation sequencing data analysis including targeted gene panel, exome, and genome sequencing data. In addition, in this review, we compile published data on pathogenic hearing loss CNVs to provide an up-to-date overview. We show that CNVs have been identified in 29 different non-syndromic hearing loss genes. An understanding of the contribution of CNVs to genetic hearing loss is critical to the current diagnosis of hearing loss and is crucial for future gene therapies. Thus, evaluation for CNVs is required in any modern pipeline for genetic diagnosis of hearing loss., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

A non-negligible proportion of human pathogenic variants are known to be present as wild type in at least some non-human mammalian species. The standard explanation for this finding is that molecular mechanisms of compensatory epistasis can alleviate the mutations' otherwise pathogenic effects. Examples of compensated variants have been described in the literature but the interacting residue(s) postulated to play a compensatory role have rarely been ascertained. In this study, the examination of five human X-chromosomally encoded proteins (FIX, GLA, HPRT1, NDP and OTC) allowed us to identify several candidate compensated variants. Strong evidence for a compensated/compensatory pair of amino acids in the coagulation FIXa protein (involving residues 270 and 271) was found in a variety of mammalian species. Both amino acid residues are located within the 60-loop, spatially close to the 39-loop that performs a key role in coagulation serine proteases. To understand the nature of the underlying interactions, molecular dynamics simulations were performed. The predicted conformational change in the 39-loop consequent to the Glu270Lys substitution (associated with hemophilia B) appears to impair the protein's interaction with its substrate but, importantly, such steric hindrance is largely mitigated in those proteins that carry the compensatory residue (Pro271) at the neighboring amino acid position., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Fetal growth is an important determinant of cardiometabolic disease risk during childhood and adulthood. The genetic architecture of fetal growth remains largely understudied in ancestrally diverse populations. We conducted genome-wide admixture mapping scan and analysis of genetic ancestry among Hispanic American, African American, European American, and Asian American pregnant women to identify genetic loci associated with fetal growth measures across 13-40 weeks gestation. Fetal growth measures were associated with genome-wide average African, European, Amerindigenous and East Asian ancestry proportions (P ranged from10 -3 to 4.8 × 10 -2 ). Admixture mapping analysis identified ten African ancestry loci and three Amerindigenous ancestry loci significantly associated with fetal growth measures at Bonferroni-corrected levels of significance (P ranged from 2.18 × 10 -8 to 3.71 × 10 -6 ). At the chr2q23.3-24.2 locus in which higher African ancestry was associated with long bone (femur and humerus) lengths, the T allele of rs13030825 (GALNT13) was associated with longer humerus length in African Americans (β = 0.44, P = 6.25 × 10 -6 at week 27; β = 0.39, P = 7.72 × 10 -5 at week 40). The rs13030825 SNP accounted for most of the admixture association at the chr2q23.3-24.2 locus and has substantial allele frequency difference between African and European reference samples (F ST  = 0.55, P = 0.03). Regulatory annotation shows that rs13030825 overlaps with the serum response factor (SRF) transcription factor previously implicated in postnatal bone development of mice. Overall, we identified ancestry-related maternal genetic loci that influence fetal growth, shedding light on molecular pathways that regulate fetal growth and potential effects on health across the lifespan.Clinical trials registration ClinicalTrials.gov, NCT00912132.

Copy number variation (CNV) related disorders tend to show complex phenotypic profiles that do not match known diseases. This makes it difficult to ascertain their underlying molecular basis. A potential solution is to compare the affected genomic regions for multiple patients that share a pathological phenotype, looking for commonalities. Here, we present a novel approach to associate phenotypes with functional systems, in terms of GO categories and KEGG and Reactome pathways, based on patient data. The approach uses genomic and phenomic data from the same patients, finding shared genomic regions between patients with similar phenotypes. These regions are mapped to genes to find associated functional systems. We applied the approach to analyse patients in the DECIPHER database with de novo CNVs, finding functional systems associated with most phenotypes, often due to mutations affecting related genes in the same genomic region. Manual inspection of the ten top-scoring phenotypes found multiple FunSys connections supported by the previous studies for seven of them. The workflow also produces reports focussed on the genes and FunSys connected to the different phenotypes, alongside patient-specific reports, which give details of the associated genes and FunSys for each individual in the cohort. These can be run in "confidential" mode, preserving patient confidentiality. The workflow presented here can be used to associate phenotypes with functional systems using data at the level of a whole cohort of patients, identifying important connections that could not be found when considering them individually. The full workflow is available for download, enabling it to be run on any patient cohort for which phenotypic and CNV data are available.

We have analyzed human genetic diversity in 33 Old World populations including 23 populations obtained through Genographic Project studies. A set of 1,536 SNPs in five X chromosome regions were genotyped in 1,288 individuals (mostly males). We use a novel analysis employing subARG network construction with recombining chromosomal segments. Here, a subARG is constructed independently for each of five gene-free regions across the X chromosome, and the results are aggregated across them. For PCA, MDS and ancestry inference with STRUCTURE, the subARG is processed to obtain feature vectors of samples and pairwise distances between samples. The observed population structure, estimated from the five short X chromosomal segments, supports genome-wide frequency-based analyses: African populations show higher genetic diversity, and the general trend of shared variation is seen across the globe from Africa through Middle East, Europe, Central Asia, Southeast Asia, and East Asia in broad patterns. The recombinational analysis was also compared with established methods based on SNPs and haplotypes. For haplotypes, we also employed a fixed-length approach based on information-content optimization. Our recombinational analysis suggested a southern migration route out of Africa, and it also supports a single, rapid human expansion from Africa to East Asia through South Asia.