Estimating allele age

Population genomics unravels a lag phase during the global fall armyworm invasion

The time that elapsed between the initial introduction and the proliferation of an invasive species is referred to as the lag phase. The identification of the lag phase is critical for generating plans for pest management and for the prevention of biosecurity failure. However, lag phases have been identified mostly through retrospective searches of historical records. The agricultural pest fall armyworm (FAW; Spodoptera frugiperda) is native to the New World. FAW invasion was first reported from West Africa in 2016, then it spread quickly through Africa, Asia, and Oceania. Here, using population genomics approaches, we demonstrate that the FAW invasion involved an undocumented lag phase. Invasive FAW populations have negative signs of genomic Tajima's D, and invasive population-specific genetic variations have particularly decreased Tajima's D, supporting a substantial amount of time for the generation of new mutations in introduced FAW populations. Model-based diffusion approximations support the existence of a period with a cessation of gene flow between native and invasive FAW populations. Taken together, these results provide strong support for the presence of a lag phase during the FAW invasion. These results show the usefulness of using population genomics analyses to identify lag phases in biological invasions.

Allele age estimators designed for whole genome datasets show only a modest decrease in accuracy when applied to whole exome datasets

Personalized genomics in the healthcare system is becoming increasingly accessible as the costs of sequencing decreases. With the increase in number of genomes, larger numbers of rare variants are being discovered and much work is being done to identify their functional impacts in relation to disease phenotypes. One way to characterize these variants is to estimate the time the mutation entered the population. However, allele age estimators such as Relate, Genealogical Estimator of Variant Age, and time of coalescence, were developed based on the assumption that datasets include the entire genome. We examined the performance of each of these estimators on simulated exome data under a neutral constant population size model and found that each provides usable estimates of allele age from whole-exome datasets. To test the robustness of these methods, analyses were undertaken to simulate data under a population expansion model and background selection. Relate performs the best amongst all three estimators with Pearson coefficients of 0.64 and 0.68 (neutral constant and expansion population model) with a 17 percent and 15 percent drop in accuracy between whole genome and whole exome estimations. Of the three estimators, Relate is best able to parallelize to yield quick results with little resources, however even Relate is only able to scale to thousands of samples making it unable to match the hundreds of thousands of samples being currently released. While more work is needed to expand the capabilities of current methods of estimating allele age, these methods estimate the age of mutations with a modest decrease in performance.

A novel splice-site FHOD3 founder variant is a common cause of hypertrophic cardiomyopathy in the population of the Balkans-A cohort study

Founder variants in sarcomere protein genes account for a significant proportion of disease-causing variants in patients with hypertrophic cardiomyopathy (HCM). However, information on founder variants in non-sarcomeric protein genes, such as FHOD3, which have only recently been associated with HCM, remains scarce. In this study, we conducted a retrospective analysis of exome sequencing data of 134 probands with HCM for recurrent pathogenic variants. We discovered a novel likely pathogenic variant c.1646+2T>C in FHOD3 in heterozygous state in eight probands with HCM and confirmed its presence in seven additional relatives. Individuals with this variant had a wide range of ages at onset of the disease (4-63 years). No adverse cardiac events were observed. Haplotype analysis revealed that the individuals with this variant shared a genomic region of approximately 5 Mbp surrounding the variant, confirming the founder effect of the variant. FHOD3 c.1646+2T>C is estimated to have arisen 58 generations ago (95% CI: 45-81) in a common ancestor living on the Balkans. A founder FHOD3 c.1646+2T>C variant is the second most common genetic variant in our cohort of patients with HCM, occurring in 16% of probands with a known genetic cause of HCM, which represents a substantially higher proportion than the currently estimated 0.5-2% for causal FHOD3 variants. Our study broadens the understanding of the genetic causes of HCM and may improve the diagnosis of this condition, particularly in patients from the Balkans.

Linking New Alleles at the Oscillator Loci to Flowering and Expansion of Asian Rice

The central oscillator is believed to be the key mechanism by which plants adapt to new environments. However, impacts from hybridization, the natural environment, and human selection have rarely been assessed on the oscillator of a crop. Here, from clearly identified alleles at oscillator loci (OsCCA1/LHY, OsPRR95, OsPRR37, OsPRR59, and OsPRR1) in ten diverse genomes of Oryza sativa, additional accessions, and functional analysis, we show that rice's oscillator was rebuilt primarily by new alleles from recombining parental sequences and subsequent 5' or/and coding mutations. New alleles may exhibit altered transcript levels from that of a parental allele and are transcribed variably among genetic backgrounds and natural environments in RIL lines. Plants carrying more expressed OsCCA1_a and less transcribed OsPRR1_e flower early in the paddy field. 5' mutations are instrumental in varied transcription, as shown by EMSA tests on one deletion at the 5' region of highly transcribed OsPRR1_a. Compared to relatively balanced mutations at oscillator loci of Arabidopsis thaliana, 5' mutations of OsPRR37 (and OsCCA1 to a less degree) were under negative selection while those of OsPRR1 alleles were under strong positive selection. Together, range expansion of Asian rice can be elucidated by human selection on OsPRR1 alleles via local flowering time-yield relationships.

The BRCA1 c.4096+1G>A Is a Founder Variant Which Originated in Ancient Times

Approximately 30-50% of hereditary breast and ovarian cancer (HBOC) is due to the presence of germline pathogenic variants in the BRCA1 (OMIM 113705) and BRCA2 (OMIM 600185) onco-suppressor genes, which are involved in DNA damage response. Women who carry pathogenic BRCA1 variants are particularly likely to develop breast cancer (BC) and ovarian cancer (OC), with a 45-79 percent and 39-48 percent chance, respectively. The BRCA1 c.4096+1G>A variant has been frequently ascertained in Tuscany, Italy, and it has also been detected in other Italian regions and other countries. Its pathogenetic status has been repeatedly changed from a variant of uncertain significance, to pathogenic, to likely pathogenic. In our study, 48 subjects (38 of whom are carriers) from 27 families were genotyped with the Illumina OncoArray Infinium platform (533,531 SNPs); a 20 Mb region (24.6 cM) around BRCA1, including 4130 SNPs (21 inside BRCA1) was selected for haplotype analysis. We used a phylogenetic method to estimate the time to the most recent common ancestor (MRCA) of BRCA1 c.4096+1G>A founder pathogenic variant. This analysis suggests that the MRCA lived about 155 generations ago-around 3000 years ago.

Analyses of allele age and fitness impact reveal human beneficial alleles to be older than neutral controls

A classic population genetic prediction is that alleles experiencing directional selection should swiftly traverse allele frequency space, leaving detectable reductions in genetic variation in linked regions. However, despite this expectation, identifying clear footprints of beneficial allele passage has proven to be surprisingly challenging. We addressed the basic premise underlying this expectation by estimating the ages of large numbers of beneficial and deleterious alleles in a human population genomic data set. Deleterious alleles were found to be young, on average, given their allele frequency. However, beneficial alleles were older on average than non-coding, non-regulatory alleles of the same frequency. This finding is not consistent with directional selection and instead indicates some type of balancing selection. Among derived beneficial alleles, those fixed in the population show higher local recombination rates than those still segregating, consistent with a model in which new beneficial alleles experience an initial period of balancing selection due to linkage disequilibrium with deleterious recessive alleles. Alleles that ultimately fix following a period of balancing selection will leave a modest 'soft' sweep impact on the local variation, consistent with the overall paucity of species-wide 'hard' sweeps in human genomes.

Recurrent mutation in the ancestry of a rare variant

Recurrent mutation produces multiple copies of the same allele which may be co-segregating in a population. Yet, most analyses of allele-frequency or site-frequency spectra assume that all observed copies of an allele trace back to a single mutation. We develop a sampling theory for the number of latent mutations in the ancestry of a rare variant, specifically a variant observed in relatively small count in a large sample. Our results follow from the statistical independence of low-count mutations, which we show to hold for the standard neutral coalescent or diffusion model of population genetics as well as for more general coalescent trees. For populations of constant size, these counts are distributed like the number of alleles in the Ewens sampling formula. We develop a Poisson sampling model for populations of varying size and illustrate it using new results for site-frequency spectra in an exponentially growing population. We apply our model to a large data set of human SNPs and use it to explain dramatic differences in site-frequency spectra across the range of mutation rates in the human genome.

Insight into the Natural History of Pathogenic Variant c.919-2A>G in the SLC26A4 Gene Involved in Hearing Loss: The Evidence for Its Common Origin in Southern Siberia (Russia)

Pathogenic variants in the SLC26A4 gene leading to nonsyndromic recessive deafness (DFNB4), or Pendred syndrome, are some of the most common causes of hearing loss worldwide. Earlier, we found a high proportion of SLC26A4-related hearing loss with prevailing pathogenic variant c.919-2A>G (69.3% among all mutated SLC26A4 alleles that have been identified) in Tuvinian patients belonging to the indigenous Turkic-speaking Siberian people living in the Tyva Republic (Southern Siberia, Russia), which implies a founder effect in the accumulation of c.919-2A>G in Tuvinians. To evaluate a possible common origin of c.919-2A>G, we genotyped polymorphic STR and SNP markers, intragenic and flanking SLC26A4, in patients homozygous for c.919-2A>G and in healthy controls. The common STR and SNP haplotypes carrying c.919-2A>G were revealed, which convincingly indicates the origin of c.919-2A>G from a single ancestor, supporting a crucial role of the founder effect in the c.919-2A>G prevalence in Tuvinians. Comparison analysis with previously published data revealed the identity of the small SNP haplotype (~4.5 kb) in Tuvinian and Han Chinese carriers of c.919-2A>G, which suggests their common origin from founder chromosomes. We assume that c.919-2A>G could have originated in the geographically close territories of China or Tuva and subsequently spread to other regions of Asia. In addition, the time intervals of the c.919-2A>G occurrence in Tuvinians were roughly estimated.

Global Distribution of Founder Variants Associated with Non-Syndromic Hearing Impairment

The genetic etiology of non-syndromic hearing impairment (NSHI) is highly heterogeneous with over 124 distinct genes identified. The wide spectrum of implicated genes has challenged the implementation of molecular diagnosis with equal clinical validity in all settings. Differential frequencies of allelic variants in the most common NSHI causal gene, gap junction beta 2 (GJB2), has been described as stemming from the segregation of a founder variant and/or spontaneous germline variant hot spots. We aimed to systematically review the global distribution and provenance of founder variants associated with NSHI. The study protocol was registered on PROSPERO, the International Prospective Register of Systematic Reviews, with the registration number "CRD42020198573". Data from 52 reports, involving 27,959 study participants from 24 countries, reporting 56 founder pathogenic or likely pathogenic (P/LP) variants in 14 genes (GJB2, GJB6, GSDME, TMC1, TMIE, TMPRSS3, KCNQ4, PJVK, OTOF, EYA4, MYO15A, PDZD7, CLDN14, and CDH23), were reviewed. Varied number short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs) were used for haplotype analysis to identify the shared ancestral informative markers in a linkage disequilibrium and variants' origins, age estimates, and common ancestry computations in the reviewed reports. Asia recorded the highest number of NSHI founder variants (85.7%; 48/56), with variants in all 14 genes, followed by Europe (16.1%; 9/56). GJB2 had the highest number of ethnic-specific P/LP founder variants. This review reports on the global distribution of NSHI founder variants and relates their evolution to population migration history, bottleneck events, and demographic changes in populations linked with the early evolution of deleterious founder alleles. International migration and regional and cultural intermarriage, coupled to rapid population growth, may have contributed to re-shaping the genetic architecture and structural dynamics of populations segregating these pathogenic founder variants. We have highlighted and showed the paucity of data on hearing impairment (HI) variants in Africa, establishing unexplored opportunities in genetic traits.

Museum Genomics of an Agricultural Super-Pest, the Colorado Potato Beetle, Leptinotarsa decemlineata (Chrysomelidae), Provides Evidence of Adaptation from Standing Variation

Despite extensive research on agricultural pests, our knowledge about their evolutionary history is often limited. A mechanistic understanding of the demographic changes and modes of adaptation remains an important goal, as it improves our understanding of organismal responses to environmental change and our ability to sustainably manage pest populations. Emerging genomic datasets now allow for characterization of demographic and adaptive processes, but face limits when they are drawn from contemporary samples, especially in the context of strong demographic change, repeated selection, or adaptation involving modest shifts in allele frequency at many loci. Temporal sampling, however, can improve our ability to reconstruct evolutionary events. Here, we leverage museum samples to examine whether population genomic diversity and structure has changed over time, and to identify genomic regions that appear to be under selection. We focus on the Colorado potato beetle (CPB), Leptinotarsa decemlineata (Say 1824; Coleoptera: Chrysomelidae), which is widely regarded as a super-pest due to its rapid, and repeated, evolution to insecticides. By combining whole genome resequencing data from 78 museum samples with modern sampling, we demonstrate that CPB expanded rapidly in the 19th century, leading to a reduction in diversity and limited genetic structure from the Midwest to Northeast United States. Temporal genome scans provide extensive evidence for selection acting in resistant field populations in Wisconsin and New York, including numerous known insecticide resistance genes. We also validate these results by showing that known selective sweeps in modern populations are identified by our genome scan. Perhaps most importantly, temporal analysis indicates selection on standing genetic variation, as we find evidence for parallel evolution in the two geographical regions. Parallel evolution involves a range of phenotypic traits not previously identified as under selection in CPB, such as reproductive and morphological functional pathways that might be important for adaptation to agricultural habitats.

Modeling the spatiotemporal spread of beneficial alleles using ancient genomes

Ancient genome sequencing technologies now provide the opportunity to study natural selection in unprecedented detail. Rather than making inferences from indirect footprints left by selection in present-day genomes, we can directly observe whether a given allele was present or absent in a particular region of the world at almost any period of human history within the last 10,000 years. Methods for studying selection using ancient genomes often rely on partitioning individuals into discrete time periods or regions of the world. However, a complete understanding of natural selection requires more nuanced statistical methods which can explicitly model allele frequency changes in a continuum across space and time. Here we introduce a method for inferring the spread of a beneficial allele across a landscape using two-dimensional partial differential equations. Unlike previous approaches, our framework can handle time-stamped ancient samples, as well as genotype likelihoods and pseudohaploid sequences from low-coverage genomes. We apply the method to a panel of published ancient West Eurasian genomes to produce dynamic maps showcasing the inferred spread of candidate beneficial alleles over time and space. We also provide estimates for the strength of selection and diffusion rate for each of these alleles. Finally, we highlight possible avenues of improvement for accurately tracing the spread of beneficial alleles in more complex scenarios.

MPYS Modulates Fatty Acid Metabolism and Immune Tolerance at Homeostasis Independent of Type I IFNs

MPYS/STING (stimulator of IFN genes) senses cyclic dinucleotides (CDNs), generates type I IFNs, and plays a critical role in infection, inflammation, and cancer. In this study, analyzing genotype and haplotype data from the 1000 Genomes Project, we found that the R71H-G230A-R293Q (HAQ) MPYS allele frequency increased 57-fold in East Asians compared with sub-Saharan Africans. Meanwhile, the G230A-R293Q (AQ) allele frequency decreased by 98% in East Asians compared with sub-Saharan Africans. We propose that the HAQ and AQ alleles underwent a natural selection during the out-of-Africa migration. We used mouse models of HAQ and AQ to investigate the underlying mechanism. We found that the mice carrying the AQ allele, which disappeared in East Asians, had normal CDN-type I IFN responses. Adult AQ mice, however, had less fat mass than did HAQ or wild-type mice on a chow diet. AQ epididymal adipose tissue had increased regulatory T cells and M2 macrophages with protein expression associated with enhanced fatty acid oxidation. Conditional knockout mice and adoptive cell transfer indicate a macrophage and regulatory T cell-intrinsic role of MPYS in fatty acid metabolism. Mechanistically, AQ/IFNAR1-/- mice had a similar lean phenotype as for the AQ mice. MPYS intrinsic tryptophan fluorescence revealed that the R71H change increased MPYS hydrophilicity. Lastly, we found that the second transmembrane (TM) and the TM2-TM3 linker region of MPYS interact with activated fatty acid, fatty acyl-CoA. In summary, studying the evolution of the human MPYS gene revealed an MPYS function in modulating fatty acid metabolism that may be critical during the out-of-Africa migration.

Genetic basis of speciation and adaptation: from loci to causative mutations

Does evolution proceed in small steps or large leaps? How repeatable is evolution? How constrained is the evolutionary process? Answering these long-standing questions in evolutionary biology is indispensable for both understanding how extant biodiversity has evolved and predicting how organisms and ecosystems will respond to changing environments in the future. Understanding the genetic basis of phenotypic diversification and speciation in natural populations is key to properly answering these questions. The leap forward in genome sequencing technologies has made it increasingly easier to not only investigate the genetic architecture but also identify the variant sites underlying adaptation and speciation in natural populations. Furthermore, recent advances in genome editing technologies are making it possible to investigate the functions of each candidate gene in organisms from natural populations. In this article, we discuss how these recent technological advances enable the analysis of causative genes and mutations and how such analysis can help answer long-standing evolutionary biology questions.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

A recurrent ZP1 variant is responsible for Oocyte Maturation Defect with degenerated oocytes in infertile females

A female factor is present in approximately 70% of couple infertility, often due to ovulatory disorders. In oocyte maturation defect (OMD), affected patients have a primary infertility with normal menstrual cycles but produce no oocyte, degenerated (atretic) or abnormal oocytes blocked at different stages of maturation. Four genes have so far been associated with OMD: PATL2, TUBB8, WEE2, and ZP1. In our initial study, 6 out of 23 OMD subjects were shown to carry the same PATL2 homozygous loss of function variant and one patient had a TUBB8 truncating variant. Here, we included four additional OMD patients and reanalyzed all 27 subjects. In addition to the seven patients with a previously identified defect, five carried the same deleterious homozygous ZP1 variant (c.1097G>A; p.Arg366Gln). All the oocytes from ZP1‐associated patients appeared shriveled and dark indicating that the abnormal ZP1 protein induced oocyte death and degeneration. Overall ZP1‐associated patients had degenerated or absent oocytes contrary to PATL2‐associated subjects who had immature oocytes blocked mainly at the germinal vesicle stage. In this cohort of North African OMD patients, whole exome sequencing permitted to diagnose 44% of the patients studied and to identify a new frequent ZP1 variant.

CLCC1 c. 75C>A Mutation in Pakistani Derived Retinitis Pigmentosa Families Likely Originated With a Single Founder Mutation 2,000-5,000 Years Ago

Background: A CLCC1 c. 75C > A (p.D25E) mutation has been associated with autosomal recessive pigmentosa in patients in and from Pakistan. CLCC1 is ubiquitously expressed, and knockout models of this gene in zebrafish and mice are lethal in the embryonic period, suggesting that possible retinitis pigmentosa mutations in this gene might be limited to those leaving partial activity. In agreement with this hypothesis, the mutation is the only CLCC1 mutation associated with retinitis pigmentosa to date, and all identified patients with this mutation share a common SNP haplotype surrounding the mutation, suggesting a common founder.

Methods: SNPs were genotyped by a combination of WGS and Sanger sequencing. The original founder haplotype, and recombination pathways were delineated by examination to minimize recombination events. Mutation age was estimated by four methods including an explicit solution, an iterative approach, a Bayesian approach and an approach based solely on ancestral segment lengths using high density SNP data.

Results: All members of each of the nine families studied shared a single autozygous SNP haplotype for the CLCC1 region ranging from approximately 1–3.5 Mb in size. The haplotypes shared by the families could be derived from a single putative ancestral haplotype with at most two recombination events. Based on the haplotype and Gamma analysis, the estimated age of the founding mutation varied from 79 to 196 generations, or approximately 2,000–5,000 years, depending on the markers used in the estimate. The DMLE (Bayesian) estimates ranged from 2,160 generations assuming a population growth rate of 0–309 generations assuming a population growth rate of 2% with broad 95% confidence intervals.

Conclusion: These results provide insight into the origin of the CLCC1 mutation in the Pakistan population. This mutation is estimated to have occurred 2000–5,000 years ago and has been transmitted to affected families of Pakistani origin in geographically dispersed locations around the world. This is the only mutation in CLCC1 identified to date, suggesting that the CLCC1 gene is under a high degree of constraint, probably imposed by functional requirements for this gene during embryonic development.

Adaptive eQTLs reveal the evolutionary impacts of pleiotropy and tissue-specificity while contributing to health and disease

Large numbers of expression quantitative trait loci (eQTLs) have recently been identified in humans, and many of these regulatory variants have large allele frequency differences between populations. Here, we conducted genome-wide scans of selection to identify adaptive eQTLs (i.e., eQTLs with large population branch statistics). We then tested if tissue pleiotropy affects whether eQTLs are more or less likely to be adaptive and identified tissues that have been key targets of positive selection during the last 100,000 years. Top adaptive eQTL outliers include rs1043809, rs66899053, and rs2814778 (a SNP that is associated with malaria resistance). We found that effect sizes of eQTLs were negatively correlated with population branch statistics and that adaptive eQTLs affect two-thirds as many tissues as do non-adaptive eQTLs. Because the tissue breadth of an eQTL can be viewed as a measure of pleiotropy, these results imply that pleiotropy inhibits adaptation. The proportion of eQTLs that are adaptive varies by tissue, and we found that eQTLs that regulate expression in testis, thyroid, blood, or sun-exposed skin are enriched for signatures of positive selection. By contrast, eQTLs that regulate expression in the cerebrum or female-specific tissues have a relative lack of adaptive outliers. Scans of selections also reveal that many adaptive eQTLs are closely linked to disease-associated loci. Taken together, our results indicate that eQTLs have played an important role in recent human evolution.

An ancient positively selected BMPRIB missense variant increases litter size of Mongolian sheep populations following latitudinal gradient

New gene mutation origination is a driving force for the evolution of organisms. The effect of FecB mutation in BMPRIB gene on the litter size of sheep has been well known for a long time, each copy of the mutant allele increases litter size by 0.4-0.5. However, the origin and adaptive evolution mechanism of FecB mutation are still unclear. Here we carried on the thorough analysis on evolutionary features of BMPRIB gene and found that 150 species as a whole is under purifying selection while sheep lineage shows evidence of positive selection. The results of allele age estimation revealed that the FecB mutation in Mongolian sheep of China originated in Mongolian Plateau at about 5000 years ago. Due the two shape drops in temperature subsequently, Mongolian sheep migrated from north to south following the northern nomadic people. Accordingly, the FecB mutant allele frequency increased, with the lowest in sheep locating at Mongolian plateau (0.01) and the highest in sheep locating at Yangtze River valley (0.96). In conclusion, the FecB mutation in Mongolian sheep of China originated in Mongolian Plateau at about 5000 years ago, and the differentiated litter size of Mongolian sheep might be the result of adaptation to various environments during the migration following latitudinal gradient. This study may well exemplify selection on an ancient variation triggered by drastic ecological shifts, and is also helpful to analyze the adaptive evolution mechanism of economic traits of domestic animals and identify major genes and molecular markers.

Current Developments in Detection of Identity-by-Descent Methods and Applications

Identity-by-descent (IBD), the detection of shared segments inherited from a common ancestor, is a fundamental concept in genomics with broad applications in the characterization and analysis of genomes. While historically the concept of IBD was extensively utilized through linkage analyses and in studies of founder populations, applications of IBD-based methods subsided during the genome-wide association study era. This was primarily due to the computational expense of IBD detection, which becomes increasingly relevant as the field moves toward the analysis of biobank-scale datasets that encompass individuals from highly diverse backgrounds. To address these computational barriers, the past several years have seen new methodological advances enabling IBD detection for datasets in the hundreds of thousands to millions of individuals, enabling novel analyses at an unprecedented scale. Here, we describe the latest innovations in IBD detection and describe opportunities for the application of IBD-based methods across a broad range of questions in the field of genomics.

The Aldehyde Dehydrogenase ALDH2*2 Allele, Associated with Alcohol Drinking Behavior, Dates Back to Prehistoric Times

Human alcohol-consumption behavior is partly genetically encoded. The alcohol consumption of 987 residents in Keelung, Taiwan, was evaluated by using the Alcohol Use Disorder Identification Test (AUDIT). We assessed ~750,000 genomic variants of 71 residents who drank hazardously (AUDIT score ≥ 8) and 126 residents who did not drink in their daily lives (AUDIT score = 0), using high-density single nucleotide polymorphism (SNP) arrays. The rs671 G > A manifests the highest significance of the association with drinking behavior (Fisher’s exact P = 8.75 × 10⁻⁹). It is a pleiotropic, non-synonymous variant in the aldehyde dehydrogenase 2 (ALDH2) gene. The minor allele “A”, commonly known as ALDH2*2, is associated with non-drinkers. Intriguingly, identity-by-descent haplotypes encompassing genomic regions with a median length of 1.6 (0.6–2.0) million nucleotide bases were found in all study participants with either heterozygous or homozygous ALDH2*2 (n = 81 and 13, respectively). We also analyzed a public-domain dataset with genome-wide genotypes of 2000 participants in Guangzhou, a coastal city in Southern China. Among them, 175 participants have homozygous ALDH2*2 genotype, and again, long ALDH2*2-carrying haplotypes were found in all 175 participants without exceptions. The median length of the ALDH2*2-carrying haplotype is 1.7 (0.5–2.8) million nucleotide bases. The haplotype lengths in the Keelung and Guangzhou cohorts combined indicate that the origin of the ALDH2*2 allele dates back to 7935 (7014–9381) years ago. In conclusion, the rs671 G > A is the leading genomic variant associated with the long-term drinking behavior among residents of Keelung, Taiwan. The ALDH2*2 allele has been in Asian populations since prehistoric times.

Contribution of the mutation T865G in TPH1 gene to the genetic potentiality of housed Mongolian sheep to year-round breeding

Seasonal breeding is widespread in sheep and significantly affects the development of the housed sheep industry. To improve and balance the reproduction performance of sheep, year-round breeding has the goal of modern sheep farming. The tryptophan hydroxylase (TPH), which initiates and regulates biosynthesis of melatonin, is an important player in the formation of mammalian year-round breeding. However, little is known about its role in regulation of sheep seasonal breeding. In this study, a missense mutation, T865G in TPH1 gene was detected in 328 individuals of six Mongolian sheep groups. It was positively selected among Mongolian sheep. This mutation may appear between 13,683 and 350,973 years ago and only exist in Hu sheep now. In Hu sheep, the frequency of allele T was 89.66%, and that of allele G was 10.34%. The TPH1 protein structure and property analysis suggested that this mutation from T to G affect the three-dimensional structure and reduce the hydropathicity of catalytic core. When the allele is T, the protein activity is twice that of the allele G, and their difference was significant (p < .05). In conclusion, T865G is an ancient mutation of TPH1 gene and affects the function of TPH protein, which may contribute to the genetic potentiality of Mongolian sheep to year-round breeding.

Evolution and genomic signatures of spontaneous somatic mutation in Drosophila intestinal stem cells

Spontaneous mutations can alter tissue dynamics and lead to cancer initiation. Although large-scale sequencing projects have illuminated processes that influence somatic mutation and subsequent tumor evolution, the mutational dynamics operating in the very early stages of cancer development are currently not well understood. To explore mutational processes in the early stages of cancer evolution, we exploited neoplasia arising spontaneously in the Drosophila intestine. Analysing whole-genome sequencing data with a dedicated bioinformatic pipeline, we found neoplasia formation to be driven largely through the inactivation of Notch by structural variants, many of which involve highly complex genomic rearrangements. The genome-wide mutational burden in neoplasia was found to be similar to that of several human cancers. Finally, we identified genomic features associated with spontaneous mutation, and defined the evolutionary dynamics and mutational landscape operating within intestinal neoplasia over the short lifespan of the adult fly. Our findings provide unique insight into mutational dynamics operating over a short timescale in the genetic model system, Drosophila melanogaster.

Identification and Characterization of the Most Common Genetic Variant Responsible for Acephalic Spermatozoa Syndrome in Men Originating from North Africa

Acephalic spermatozoa syndrome (ASS) is a rare but extremely severe type of teratozoospermia, defined by the presence of a majority of headless flagella and a minority of tail-less sperm heads in the ejaculate. Like the other severe monomorphic teratozoospermias, ASS has a strong genetic basis and is most often caused by bi-allelic variants in SUN5 (Sad1 and UNC84 domain-containing 5). Using whole exome sequencing (WES), we investigated a cohort of nine infertile subjects displaying ASS. These subjects were recruited in three centers located in France and Tunisia, but all originated from North Africa. Sperm from subjects carrying candidate genetic variants were subjected to immunofluorescence analysis and transmission electron microscopy. Moreover, fluorescent in situ hybridization (FISH) was performed on sperm nuclei to assess their chromosomal content. Variant filtering permitted us to identify the same SUN5 homozygous frameshift variant (c.211+1_211+2dup) in 7/9 individuals (78%). SUN5 encodes a protein localized on the posterior part of the nuclear envelope that is necessary for the attachment of the tail to the sperm head. Immunofluorescence assays performed on sperm cells from three mutated subjects revealed a total absence of SUN5, thus demonstrating the deleterious impact of the identified variant on protein expression. Transmission electron microscopy showed a conserved flagellar structure and a slightly decondensed chromatin. FISH did not highlight a higher rate of chromosome aneuploidy in spermatozoa from SUN5 patients compared to controls, indicating that intra-cytoplasmic sperm injection (ICSI) can be proposed for patients carrying the c.211+1_211+2dup variant. These results suggest that the identified SUN5 variant is the main cause of ASS in the North African population. Consequently, a simple and inexpensive genotyping of the 211+1_211+2dup variant could be beneficial for affected men of North African origin before resorting to more exhaustive genetic analyses.

Age and Origin of the Founder Antithrombin Budapest 3 (p.Leu131Phe) Mutation; Its High Prevalence in the Roma Population and Its Association With Cardiovascular Diseases

Background: Antithrombin (AT) is one of the most important regulator of hemostasis. AT Budapest 3 (ATBp3) is a prevalent type II heparin-binding site (IIHBS) deficiency due to founder effect. Thrombosis is a complex disease including arterial (ATE) and venous thrombotic events (VTE) and the Roma population, the largest ethnic minority in Europe has increased susceptibility to these diseases partly due to their unfavorable genetic load. We aimed to calculate the age and origin of ATBp3 and to explore whether the frequency of it is higher in the Roma population as compared with the general population from the corresponding geographical area. We investigated the association of ATBp3 with thrombotic events in well-defined patients' populations in order to refine the recommendation when testing for ATBp3 is useful.

Methods and Results: Prevalence of ATBp3, investigated in large samples (n = 1,000 and 1,185 for general Hungarian and Roma populations, respectively) was considerably high, almost 3%, among Roma and the founder effect was confirmed in their samples, while it was absent in the Hungarian general population. Age of ATBp3—as calculated by analysis of 8 short tandem repeat sequences surrounding SERPINC1—was dated back to XVII Century, when Roma migration in Central and Eastern Europe occurred. In our IIHBS cohort (n = 230), VTE was registered in almost all ATBp3 homozygotes (93%) and in 44% of heterozygotes. ATE occurred with lower frequency in ATBp3 (around 6%); it was rather associated with AT Basel (44%). All patients with ATE were young at the time of diagnosis. Upon investigating consecutive young (<40 years) patients with ATE (n = 92) and VTE (n = 110), the presence of ATBp3 was remarkable.

Conclusions: ATBp3, a 400-year-old founder mutation is prevalent in Roma population and its Roma origin can reasonably be assumed. By the demonstration of the presence of ATBp3 in ATE patients, we draw the attention to consider type IIHBS AT deficiency in the background of not only VTE but also ATE, especially in selected populations as young patients without advanced atherosclerosis. We recommend including the investigation of ATBp3 as part of thrombosis risk assessment and stratification in Roma individuals.

Estimation of Natural Selection and Allele Age from Time Series Allele Frequency Data Using a Novel Likelihood-Based Approach

Temporally spaced genetic data allow for more accurate inference of population genetic parameters and hypothesis testing on the recent action of natural selection. In this work, we develop a novel likelihood-based method for jointly estimating selection coefficient and allele age from time series data of allele frequencies. Our approach is based on a hidden Markov model where the underlying process is a Wright-Fisher diffusion conditioned to survive until the time of the most recent sample. This formulation circumvents the assumption required in existing methods that the allele is created by mutation at a certain low frequency. We calculate the likelihood by numerically solving the resulting Kolmogorov backward equation backward in time while reweighting the solution with the emission probabilities of the observation at each sampling time point. This procedure reduces the two-dimensional numerical search for the maximum of the likelihood surface, for both the selection coefficient and the allele age, to a one-dimensional search over the selection coefficient only. We illustrate through extensive simulations that our method can produce accurate estimates of the selection coefficient and the allele age under both constant and nonconstant demographic histories. We apply our approach to reanalyze ancient DNA data associated with horse base coat colors. We find that ignoring demographic histories or grouping raw samples can significantly bias the inference results.

High Rates of Three Common GJB2 Mutations c.516G>C, c.-23+1G>A, c.235delC in Deaf Patients from Southern Siberia Are Due to the Founder Effect

The mutations in the GJB2 gene (13q12.11, MIM 121011) encoding transmembrane protein connexin 26 (Cx26) account for a significant portion of hereditary hearing loss worldwide. Earlier we found a high prevalence of recessive GJB2 mutations c.516G>C, c.-23+1G>A, c.235delC in indigenous Turkic-speaking Siberian peoples (Tuvinians and Altaians) from the Tyva Republic and Altai Republic (Southern Siberia, Russia) and proposed the founder effect as a cause for their high rates in these populations. To reconstruct the haplotypes associated with each of these mutations, the genotyping of polymorphic genetic markers both within and flanking the GJB2 gene was performed in 28 unrelated individuals homozygous for c.516G>C (n = 18), c.-23+1G>A (n = 6), or c.235delC (n = 4) as well as in the ethnically matched controls (62 Tuvinians and 55 Altaians) without these mutations. The common haplotypes specific for mutations c.516G>C, c.-23+1G>A, or c.235delC were revealed implying a single origin of each of these mutations. The age of mutations estimated by the DMLE+ v2.3 software and the single marker method is discussed in relation to ethnic history of Tuvinians and Altaians. The data obtained in this study support a crucial role of the founder effect in the high prevalence of GJB2 mutations c.516G>C, c.-23+1G>A, c.235delC in indigenous populations of Southern Siberia.

Do Genes Associated with Dyslexia of Chinese Characters Evolve Neutrally?

Dyslexia, or reading disability, is found to have a genetic basis, and several related genes have been reported. We investigated whether natural selection has acted on single nucleotide polymorphisms (SNPs) that were reported to be associated with risk/non-risk for the reading disability of Chinese characters. We applied recently developed 2D SFS-based statistics to SNP data of East Asian populations to examine whether there is any sign of selective sweep. While neutrality was not rejected for most SNPs, significant signs of selection were detected for two linkage disequilibrium (LD) regions containing the reported SNPs of GNPTAB and DCDC2. Furthermore, we searched for a selection target site among the SNPs in these LD regions, because a causal site is not necessarily a reported SNP but could instead be a tightly linked site. In both LD regions, we found candidate target sites, which may have an effect on expression regulation and have been selected, although which genes these SNPs affect remains unknown. Because most people were not engaged in reading until recently, it is unlikely that there has been selective pressure on reading ability itself. Consistent with this, our results suggest a possibility of genetic hitchhiking, whereby alleles of the reported SNPs may have increased in frequency together with the selected target, which could have functions for other genes and traits apart from reading ability.

Prevalence and founder effect of the BRCA1 p.(Val1833Met) variant in the Greek population, with further evidence for pathogenicity and risk modification

PURPOSE

Multiple lines of evidence have suggested a likely causative role in breast/ovarian cancer (BrCa/OvCa) predisposition for the BRCA1 p.(Val1833Met) variant, predominantly found among Greek patients. Our aim was to study the variant's prevalence and founder effect on the Greek population, while providing additional data for its pathogenicity.

METHODS

We genotyped 3531 BrCa/OvCa patients using Sanger and next generation sequencing, as well as 1558 healthy, age-matched females with real-time PCR. Carriers underwent haplotype analysis to determine a founder effect. A co-segregation analysis was applied to estimate the likelihood ratio for pathogenicity.

RESULTS

In total, 27 BrCa/OvCa patients (0.77%; 27/3531) were found to carry the p.(Val1833Met) variant. No carriers were identified in the control group diagnosis. A common shared haplotype, spanning 2.76 Mb on chromosome 17 was demonstrated among carriers, establishing the founder effect. BRCA1, p.(Val1833Met) is possibly a disease-associated variant, supported by a likelihood ratio of 1.88, while a correlation to ovarian cancer is suspected.

CONCLUSIONS

Altogether, BRCA1, p.(Val1833Met) variant is a Greek founder and is very likely to predispose for BrCa/OvCa. Therefore, such carriers should be counselled accordingly, with clinical recommendations supporting surveillance and risk-reduction strategies, while providing the option for targeted therapeutic interventions.

An estimator of first coalescent time reveals selection on young variants and large heterogeneity in rare allele ages among human populations

Allele age has long been a focus of population genetic research, primarily because it can be an important clue to the fitness effects of an allele. By virtue of their effects on fitness, alleles under directional selection are expected to be younger than neutral alleles of the same frequency. We developed a new coalescent-based estimator of a close proxy for allele age, the time when a copy of an allele first shares common ancestry with other chromosomes in a sample not carrying that allele. The estimator performs well, including for the very rarest of alleles that occur just once in a sample, with a bias that is typically negative. The estimator is mostly insensitive to population demography and to factors that can arise in population genomic pipelines, including the statistical phasing of chromosomes. Applications to 1000 Genomes Data and UK10K genome data confirm predictions that singleton alleles that alter proteins are significantly younger than those that do not, with a greater difference in the larger UK10K dataset, as expected. The 1000 Genomes populations varied markedly in their distributions for singleton allele ages, suggesting that these distributions can be used to inform models of demographic history, including recent events that are only revealed by their impacts on the ages of very rare alleles.

Chloroplast DNA variation in a hyperdiverse tropical tree community

We investigate chloroplast DNA variation in a hyperdiverse community of tropical rainforest trees in French Guiana, focusing on patterns of intraspecific and interspecific variation. We test whether a species genetic diversity is higher when it has congeners in the community with which it can exchange genes and if shared haplotypes are more frequent in genetically diverse species, as expected in the presence of introgression.We sampled a total of 1,681 individual trees from 472 species corresponding to 198 genera and sequenced them at a noncoding chloroplast DNA fragment.Polymorphism was more frequent in species that have congeneric species in the study site than in those without congeners (30% vs. 12%). Moreover, more chloroplast haplotypes were shared with congeners in polymorphic species than in monomorphic ones (44% vs. 28%).Despite large heterogeneities caused by genus-specific behaviors in patterns of hybridization, these results suggest that the higher polymorphism in the presence of congeners is caused by local introgression rather than by incomplete lineage sorting. Our findings suggest that introgression has the potential to drive intraspecific genetic diversity in species-rich tropical forests.

Genetic disease risks can be misestimated across global populations

BACKGROUND

Accurate assessment of health disparities requires unbiased knowledge of genetic risks in different populations. Unfortunately, most genome-wide association studies use genotyping arrays and European samples. Here, we integrate whole genome sequence data from global populations, results from thousands of genome-wide association studies (GWAS), and extensive computer simulations to identify how genetic disease risks can be misestimated.

RESULTS

In contrast to null expectations, we find that risk allele frequencies at known disease loci are significantly different for African populations compared to other continents. Strikingly, ancestral risk alleles are found at 9.51% higher frequency in Africa, and derived risk alleles are found at 5.40% lower frequency in Africa. By simulating GWAS with different study populations, we find that non-African cohorts yield disease associations that have biased allele frequencies and that African cohorts yield disease associations that are relatively free of bias. We also find empirical evidence that genotyping arrays and SNP ascertainment bias contribute to continental differences in risk allele frequencies. Because of these causes, polygenic risk scores can be grossly misestimated for individuals of African descent. Importantly, continental differences in risk allele frequencies are only moderately reduced if GWAS use whole genome sequences and hundreds of thousands of cases and controls. Finally, comparisons between uncorrected and corrected genetic risk scores reveal the benefits of considering whether risk alleles are ancestral or derived.

CONCLUSIONS

Our results imply that caution must be taken when extrapolating GWAS results from one population to predict disease risks in another population.

Recent Activity in Expanding Populations and Purifying Selection Have Shaped Transposable Element Landscapes across Natural Accessions of the Mediterranean Grass Brachypodium distachyon

Transposable element (TE) activity has emerged as a major cause of variation in genome size and structure among species. To what extent TEs contribute to genetic variation and divergence within species, however, is much less clear, mainly because population genomic data have so far only been available for the classical model organisms. In this study, we use the annual Mediterranean grass Brachypodium distachyon to investigate TE dynamics in natural populations. Using whole-genome sequencing data for 53 natural accessions, we identified more than 5,400 TE polymorphisms across the studied genomes. We found, first, that while population bottlenecks and expansions have shaped genetic diversity in B. distachyon, these events did not lead to lineage-specific activations of TE families, as observed in other species. Instead, the same families have been active across the species range and TE activity is homogeneous across populations, indicating the presence of conserved regulatory mechanisms. Second, almost half of the TE insertion polymorphisms are accession-specific, most likely because of recent activity in expanding populations and the action of purifying selection. And finally, although TE insertion polymorphisms are underrepresented in and around genes, more than 1,000 of them occur in genic regions and could thus contribute to functional divergence. Our study shows that while TEs in B. distachyon are “well-behaved” compared with TEs in other species with larger genomes, they are an abundant source of lineage-specific genetic variation and may play an important role in population divergence and adaptation.

A Homozygous Ancestral SVA-Insertion-Mediated Deletion in WDR66 Induces Multiple Morphological Abnormalities of the Sperm Flagellum and Male Infertility

Multiple morphological abnormalities of the sperm flagellum (MMAF) is a severe form of male infertility defined by the presence of a mosaic of anomalies, including short, bent, curled, thick, or absent flagella, resulting from a severe disorganization of the axoneme and of the peri-axonemal structures. Mutations in DNAH1, CFAP43, and CFAP44, three genes encoding axoneme-related proteins, have been described to account for approximately 30% of the MMAF cases reported so far. Here, we searched for pathological copy-number variants in whole-exome sequencing data from a cohort of 78 MMAF-affected subjects to identify additional genes associated with MMAF. In 7 of 78 affected individuals, we identified a homozygous deletion that removes the two penultimate exons of WDR66 (also named CFAP251), a gene coding for an axonemal protein preferentially localized in the testis and described to localize to the calmodulin- and spoke-associated complex at the base of radial spoke 3. Sequence analysis of the breakpoint region revealed in all deleted subjects the presence of a single chimeric SVA (SINE-VNTR-Alu) at the breakpoint site, suggesting that the initial deletion event was potentially mediated by an SVA insertion-recombination mechanism. Study of Trypanosoma WDR66's ortholog (TbWDR66) highlighted high sequence and structural analogy with the human protein and confirmed axonemal localization of the protein. Reproduction of the human deletion in TbWDR66 impaired flagellar movement, thus confirming WDR66 as a gene associated with the MMAF phenotype and highlighting the importance of the WDR66 C-terminal region.

Fine-tuning of Approximate Bayesian Computation for human population genomics

Approximate Bayesian Computation (ABC) is a flexible statistical tool widely applied to addressing a variety of questions regarding the origin and evolution of humans. The significant growth of genomic scale data from diverse geographic populations has facilitated the use of ABC in modelling the complex processes that underlie human demography and local adaptation. However, a fundamental issue still remains in how to efficiently capture patterns of genetic variation with a set of summary statistics in order to achieve better approximation of Bayesian inference. Here, we review recent advances in ABC methodology and its applications for human population genomics, with a particular focus on optimal tuning of ABC approaches for different types of genetic data and different sets of evolutionary parameters.

Nonequilibrium Neutral Theory for Hitchhikers

Selective sweep is a phenomenon of reduced variation at presumably neutrally evolving sites (hitchhikers) in the genome that is caused by the spread of a selected allele at a linked focal site, and is widely used to test for action of positive selection. Nonetheless, selective sweep may also provide an unprecedented opportunity for studying nonequilibrium properties of the neutral variation itself. We have demonstrated this possibility in relation to ancient selective sweep for modern human-specific changes and ongoing selective sweep for local population-specific changes.

Estimating Time to the Common Ancestor for a Beneficial Allele

The haplotypes of a beneficial allele carry information about its history that can shed light on its age and the putative cause for its increase in frequency. Specifically, the signature of an allele's age is contained in the pattern of variation that mutation and recombination impose on its haplotypic background. We provide a method to exploit this pattern and infer the time to the common ancestor of a positively selected allele following a rapid increase in frequency. We do so using a hidden Markov model which leverages the length distribution of the shared ancestral haplotype, the accumulation of derived mutations on the ancestral background, and the surrounding background haplotype diversity. Using simulations, we demonstrate how the inclusion of information from both mutation and recombination events increases accuracy relative to approaches that only consider a single type of event. We also show the behavior of the estimator in cases where data do not conform to model assumptions, and provide some diagnostics for assessing and improving inference. Using the method, we analyze population-specific patterns in the 1000 Genomes Project data to estimate the timing of adaptation for several variants which show evidence of recent selection and functional relevance to diet, skin pigmentation, and morphology in humans.

Inverse European Latitudinal Cline at the timeless Locus of Drosophila melanogaster Reveals Selection on a Clock Gene: Population Genetics of ls-tim

The spread of adaptive genetic variants in populations is a cornerstone of evolutionary theory but with relatively few biologically well-understood examples. Previous work on the ls-tim variant of timeless, which encodes the light-sensitive circadian regulator in Drosophila melanogaster, suggests that it may have originated in southeastern Italy. Flies characterized by the new allele show photoperiod-related phenotypes likely to be adaptive in seasonal environments. ls-tim may be spreading from its point of origin in Italy by directional selection, but there are alternative explanations for its observed clinal geographical distribution, including balancing selection and demography. From population analyses of ls-tim frequencies collected on the eastern side of the Iberian Peninsula, we show that ls-tim frequencies are inverted compared with those in Italy. This pattern is consistent with a scenario of directional selection rather than latitude-associated balancing selection. Neutrality tests further reveal the signature of directional selection at the ls-tim site, which is reduced a few kb pairs either side of ls-tim. A reanalysis of allele frequencies from a large number of microsatellite loci do not demonstrate any frequent ls-tim-like spatial patterns, so a general demographic effect or population expansion from southeastern Italy cannot readily explain current ls-tim frequencies. Finally, a revised estimate of the age of ls-tim allele using linkage disequilibrium and coalescent-based approaches reveals that it may be only 300 to 3000 years old, perhaps explaining why it has not yet gone to fixation. ls-tim thus provides a rare temporal snapshot of a new allele that has come under selection before it reaches equilibrium.

Allele Age Under Non-Classical Assumptions is Clarified by an Exact Computational Markov Chain Approach

Determination of the age of an allele based on its population frequency is a well-studied problem in population genetics, for which a variety of approximations have been proposed. We present a new result that, surprisingly, allows the expectation and variance of allele age to be computed exactly (within machine precision) for any finite absorbing Markov chain model in a matter of seconds. This approach makes none of the classical assumptions (e.g., weak selection, reversibility, infinite sites), exploits modern sparse linear algebra techniques, integrates over all sample paths, and is rapidly computable for Wright-Fisher populations up to N_e = 100,000. With this approach, we study the joint effect of recurrent mutation, dominance, and selection, and demonstrate new examples of “selective strolls” where the classical symmetry of allele age with respect to selection is violated by weakly selected alleles that are older than neutral alleles at the same frequency. We also show evidence for a strong age imbalance, where rare deleterious alleles are expected to be substantially older than advantageous alleles observed at the same frequency when population-scaled mutation rates are large. These results highlight the under-appreciated utility of computational methods for the direct analysis of Markov chain models in population genetics.

RNA expression and disease tolerance are associated with a "keystone mutation" in the ochre sea star Pisaster ochraceus

An overdominant mutation in an intron of the elongation factor 1-α (EF1A) gene in the sea star Pisaster ochraceus has shown itself to mediate tolerance to "sea star wasting disease", a pandemic that has significantly reduced sea star populations on the Pacific coast of North America. Here we use RNA sequencing of healthy individuals to identify differences in constitutive expression of gene regions that may help explain this tolerance phenotype. Our results show that individuals carrying this mutation have lower expression at a large contingent of gene regions. Individuals without this mutation also appear to have a greater cellular response to temperature stress, which has been implicated in the outbreak of sea star wasting disease. Given the ecological significance of P. ochraceus, these results may be useful in predicting the evolutionary and demographic future for Pacific intertidal communities.

A novel RLBP1 gene geographical area-related mutation present in a young patient with retinitis punctata albescens

Background

Autosomal recessive forms of retinitis punctata albescens (RPA) have been described. RPA is characterized by progressive retinal degeneration due to alteration in visual cycle and consequent deposit of photopigments in retinal pigment epithelium. Five loci have been linked to RPA onset. Among these, the retinaldehyde-binding protein 1 gene, RLBP1, is the most frequently involved and several founder mutations were reported. We report results of a genetic molecular investigation performed on a large Sicilian family in which appears a young woman with RPA.

Results

The proband is in homozygous condition for a novel RLBP1 single-pair deletion, and her healthy parents, both heterozygous, are not consanguineous. Thenovelc.398delC (p.P133Qfs*258) involves the exon 6 and leads to a premature stop codon, resulting in a truncated protein entirely missing of CRAL-TRIO lipid-binding domain.

Pedigree analysis showed other non-consanguineous relatives heterozygous for the same mutation in the family. Extension of mutation research in the native town of the proband revealed its presence also in healthy subjects, in a heterozygous condition.

Conclusions

A novel RLBP1 truncating mutation was detected in a young girl affected by RPA. Although her parents are not consanguineous, the mutation was observed in a homozygous condition. Being them native of the same small Sicilian town of Fiumedinisi, the hypothesis of a geographical area-related mutation was assessed and confirmed.

A novel recessive TTN founder variant is a common cause of distal myopathy in the Serbian population

Variants in the TTN gene have been associated with distal myopathies and other distinctive phenotypes involving skeletal and cardiac muscle. Through whole-exome sequencing we identified a novel stop-gain variant (c.107635C>T, p.(Gln35879Ter)) in the TTN gene, coding a part of the M-line of titin, in 14 patients with autosomal recessive distal myopathy and Serbian ancestry. All patients share a common 1 Mb core haplotype associated with c.107635C>T, suggesting a founder variant. In compound heterozygotes, nine other TTN variants were identified: four stop-gain, three frameshift, one missense and one splice donor variant. Patients homozygous for the common variant did not show significant clinical differences to the compound heterozygous patients. The clinical presentation of all patients was an adult onset distal myopathy with predominant lower limb involvement. In addition, most patients had normal to mildly elevated serum creatine kinase levels, myopathic electromyograms, normal cardiologic and respiratory tests and muscle pathology consistent with a dystrophic process. In this study, we describe a distinct phenotype for patients with distal myopathy associated with novel recessive TTN variants including a Serbian founder variant. Our results expand the phenotypic and genetic spectrum of titinopathies and will facilitate the diagnosis of this condition in patients of Serbian origin.

Identification and population history of CYP4V2 mutations in patients with Bietti crystalline corneoretinal dystrophy

To identify known and novel CYP4V2 mutations in patients with Bietti crystalline cornea (BCD), expand the spectrum of CYP4V2 mutations, and characterize the population history of the c.802-8_810del17insGC mutation common in Asian populations, genomic DNA was isolated from peripheral blood samples from 58 unrelated patients with clinical diagnoses of BCD. Exons and flanking intronic regions of the CYP4V2 gene were dideoxy DNA sequenced. Nonpathogenic polymorphisms were excluded and known mutations were identified by sequencing 192 unaffected individuals from similar ethnic backgrounds and examination of online databases. The age of the c.802-8_810del17insGC mutation was estimated using three independent approaches. A total of 28 CYP4V2 mutations, 9 of which were novel, were detected in the 58 patients with BCD. These included 19 missense, 4 nonsense, 2 deletion, 2 splice site, and 1 insertion-deletion mutations. Two missense variants of uncertain significance were also detected. The age of the c.802-8_810del17insGC mutation was estimated to be 1040-8200 generations in the Chinese and 300-1100 generations in the Japanese populations. These results expand the mutation spectrum of CYP4V2, and provide insight into the origin of the c.802-8_810del17insGC mutation in the Chinese population and its transmission to the Japanese population.

A frameshift deletion in the sarcomere gene MYL4 causes early-onset familial atrial fibrillation

AIMS

Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia in man, causing substantial morbidity and mortality with a major worldwide public health impact. It is increasingly recognized as a highly heritable condition. This study aimed to determine genetic risk factors for early-onset AF.

METHODS AND RESULTS

We sequenced the whole genomes of 8453 Icelanders and imputed genotypes of the 25.5 million sequence variants we discovered into 1799 Icelanders with early-onset AF (diagnosed before 60 years of age) and 337 453 controls. Each sequence variant was tested for association based on multiplicative and recessive inheritance models. We discovered a rare frameshift deletion in the myosin MYL4 gene (c.234delC) that associates with early-onset AF under a recessive mode of inheritance (allelic frequency = 0.58%). We found eight homozygous carriers of the mutation, all of whom had early-onset AF. Six of the homozygotes were diagnosed by the age of 30 and the remaining two in their 50s. Three of the homozygotes had received pacemaker implantations due to sick sinus syndrome, three had suffered an ischemic stroke, and one suffered sudden cardiac death.

CONCLUSIONS

Through a population approach we found a loss of function mutation in the myosin gene MYL4 that, in the homozygous state, is completely penetrant for early-onset AF. The finding may provide novel mechanistic insight into the pathophysiology of this complex arrhythmia.

Haplotype analysis reveals that the recurrent BRCA1 deletion of exons 23 and 24 is a Greek founder mutation

A recurrent large genomic rearrangement (LGR) encompassing exons 23 and 24 of the BRCA1 gene has been identified in breast-ovarian cancer families of Greek origin. Its breakpoints have been determined as c.5406 + 664_*8273del11052 (RefSeq: NM_007294.3) and a diagnostic polymerase chain reaction (PCR) has been set up for rapid screening. In a series of 2,092 high-risk families completely screened for BRCA1 and BRCA2 germline mutations, we have found the deletion in 35 families (1.68%), representing 7.83% of the mutations identified in both genes and 10.3% of the total BRCA1 mutations. In order to characterize this deletion as a founder mutation, haplotype analysis was conducted in 60 carriers from 35 families, using three BRCA1 intragenic microsatellite markers and four markers surrounding the BRCA1 locus. Our results demonstrate a common shared core disease-associated haplotype of 2.89Mb. Our calculations estimate that the deletion has originated from a common ancestor 1450 years ago, which most probably inhabited the Asia Minor area. The particular (LGR) is the third mutation of such type that is proven to have a Greek founder effect in the Greek population, illustrating the necessity for LGRs testing in individuals of Greek descent.

Founder mutation in KCNJ10 in Pakistani patients with EAST syndrome

Background

EAST syndrome is an autosomal recessive disorder caused by loss‐of‐function mutations in the gene KCNJ10. Among the 14 pathogenic mutations described so far, the p.R65P mutation stands out as the most frequent one and is particularly associated with patients of Pakistani origin. As a result we aimed to establish the existence of a potential founder effect in the Pakistani population.

Methods

To this end, we genotyped 12 patients from seven families and we compared disease haplotypes with ethnically matched control chromosomes. This haplotype was used together with demographic data for Pakistan to estimate the age of this founder mutation.

Results

We identified a small homozygous 0.694 Mb region around the KCNJ10 p.R65P mutation that had identical haplotypes in all of the patients which were completely absent in the control sample. Based on current demographic data and knowledge about disease frequency, we estimate that this particular p.R65P mutation arose 20 generations (about 500 years) ago.

Conclusion

By knowing the prevalent mutation in a given population more efficient diagnostics can be performed and the families can benefit from specific counseling.

A Founder Effect for the HGD G360R Mutation in Italy: Implications for a Regional Screening of Alkaptonuria

We sought to establish rapid and specific genotyping methods for G360R mutation and for seven tightly linked markers in the homogentisate dioxygenase gene to address the question of whether G360R is a mutational hot spot or the result of a founder effect, as it has been repeatedly found in alkaptonuric patients from a geographic isolate in Italy.For G360R and single nucleotide polymorphism genotyping, high-resolution melting analysis was performed. Microsatellites were analysed by multiplex PCR and capillary electrophoresis. To investigate the natural history of the G360R mutation, we genotyped markers in 52 controls and in 8 unrelated patients from the UK and USA, who also segregated the G360R mutation, and calculated its age using DMLE+2.3 software.A distinct G360R-bearing haplotype was identified in all patients of Caucasian descent. Estimated mutation age was 545 generations (95% credible set, 402-854), suggesting that G360R arose in an ancestor who lived 8,000-10,000 years BC. Archaeological, historical and demographic data support that a G360R carrier has settled the remote valley where present-day population might have a heterozygote frequency of at least 6%.Given the late health-threatening complications of alkaptonuria and a cure within reach, inhabitants of this isolate would benefit from screening and genetic counselling.

A study of common Mendelian disease carriers across ageing British cohorts: meta-analyses reveal heterozygosity for alpha 1-antitrypsin deficiency increases respiratory capacity and height

BACKGROUND

Several recessive Mendelian disorders are common in Europeans, including cystic fibrosis (CFTR), medium-chain-acyl-Co-A-dehydrogenase deficiency (ACADM), phenylketonuria (PAH) and alpha 1-antitrypsin deficiency (SERPINA1).

METHODS

In a multicohort study of >19,000 older individuals, we investigated the relevant phenotypes in heterozygotes for these genes: lung function (forced expiratory volume in 1 second (FEV1), forced vital capacity (FVC)) for CFTR and SERPINA1; cognitive measures for ACADM and PAH; and physical capability for ACADM, PAH and SERPINA1.

RESULTS

Findings were mostly negative but lung function in SERPINA1 (protease inhibitor (PI) Z allele, rs28929474) showed enhanced FEV1 and FVC (0.13 z-score increase in FEV1 (p=1.7 × 10(-5)) and 0.16 z-score increase in FVC (p=5.2 × 10(-8))) in PI-MZ individuals. Height adjustment (a known, strong correlate of FEV1 and FVC) revealed strong positive height associations of the Z allele (1.50 cm increase in height (p=3.6 × 10(-10))).

CONCLUSIONS

The PI-MZ rare (2%) SNP effect is nearly four times greater than the 'top' common height SNP in HMGA2. However, height only partially attenuates the SERPINA1-FEV1 or FVC association (around 50%) and vice versa. Height SNP variants have recently been shown to be positively selected collectively in North versus South Europeans, while the Z allele high frequency is localised to North Europe. Although PI-ZZ is clinically disadvantageous to lung function, PI-MZ increases both height and respiratory function; potentially a balanced polymorphism. Partial blockade of PI could conceivably form part of a future poly-therapeutic approach in very short children. The notion that elastase inhibition should benefit patients with chronic obstructive pulmonary disease may also merit re-evaluation. PI is already a therapeutic target: our findings invite a reconsideration of the optimum level in respiratory care and novel pathway potential for development of agents for the management of growth disorders.

An Ashkenazi founder mutation in the PKHD1 gene

Autosomal recessive polycystic kidney disease (ARPKD) is usually detected late in pregnancies in embryos with large echogenic kidneys accompanied by oligohydramnios. Hundreds of private pathogenic variants have been identified in the large PKHD1 gene in various populations. Yet, because of the large size of the gene, segregation analysis of microsatellite polymorphic markers residing in the PKDH1 locus has commonly been utilized for prenatal diagnosis. Keeping in mind the limitations of this strategy, we utilized it for testing 7 families with affected fetuses or newborns, of which in 5 at least one parent was Ashkenazi, and identified that the same haplotype was shared by the majority of the Ashkenazi parents (7/9). This led us to suspect that they carry the same founder mutation. Whole Exome analysis of DNA from a fetus of one of the families detected an already known pathogenic variant c.3761_3762delCCinsG, an indel variant resulting in frameshift (p.Ala1254GlyfsX49). This variant was detected in 9 parents (5 families), of them 7 individuals were Ashkenazi and one Moroccan Jew who shared the same haplotype, and one Ashkenazi, who carried the same variant on a recombinant haplotype. Screening for this variant in 364 Ashkenazi individuals detected 2 carriers. These findings suggest that although c.3761_3762delCCinsG is considered one of the frequent variants detected in unrelated individuals, and was thought to have occurred independently on various haplotypes, it is in fact a founder mutation in the Ashkenazi population.

Estimating the Ages of Selection Signals from Different Epochs in Human History

Genetic variation harbors signatures of natural selection driven by selective pressures that are often unknown. Estimating the ages of selection signals may allow reconstructing the history of environmental changes that shaped human phenotypes and diseases. We have developed an approximate Bayesian computation (ABC) approach to estimate allele ages under a model of selection on new mutations and under demographic models appropriate for human populations. We have applied it to two resequencing data sets: An ultra-high depth data set from a relatively small sample of unrelated individuals and a lower depth data set in a larger sample with transmission information. In addition to evaluating the accuracy of our method based on simulations, for each SNP, we assessed the consistency between the posterior probabilities estimated by the ABC approach and the ancient DNA record, finding good agreement between the two types of data and methods. Applying this ABC approach to data for eight single nucleotide polymorphisms (SNPs), we were able to rule out an onset of selection prior to the dispersal out-of-Africa for three of them and more recent than the spread of agriculture for an additional three SNPs.

Expanding the Spectrum of Founder Mutations Causing Isolated Gonadotropin-Releasing Hormone Deficiency

CONTEXT

Loss of function (LoF) mutations in more than 20 genes are now known to cause isolated GnRH deficiency (IGD) in humans. Most causal IGD mutations are typically private, ie, limited to a single individual/pedigree. However, somewhat paradoxically, four IGD genes (GNRH1, TAC3, PROKR2, and GNRHR) have been shown to harbor LoF founder mutations that are shared by multiple unrelated individuals. It is not known whether similar founder mutations occur in other IGD genes.

OBJECTIVE

The objective of the study was to determine whether shared deleterious mutations in IGD-associated genes represent founder alleles.

SETTING

This study was an international collaboration among academic medical centers.

METHODS

IGD patients with shared mutations, defined as those documented in three or more unrelated probands in 14 IGD-associated genes, were identified from various academic institutions, the Human Gene Mutation Database, and literature reports by other international investigators. Haplotypes of single-nucleotide polymorphisms and short tandem repeats surrounding the mutations were constructed to assess genetic ancestry.

RESULTS

A total of eight founder mutations in five genes, GNRHR (Q106R, R262Q, R139H), TACR3 (W275X), PROKR2 (R85H), FGFR1 (R250Q, G687R), and HS6ST1 (R382W) were identified. Most founder alleles were present at low frequency in the general population. The estimated age of these mutant alleles ranged from 1925 to 5600 years and corresponded to the time of rapid human population expansion.

CONCLUSIONS

We have expanded the spectrum of founder alleles associated with IGD to a total of eight founder mutations. In contrast to the approximately 9000-year-old PROKR2 founder allele that may confer a heterozygote advantage, the rest of the founder alleles are relatively more recent in origin, in keeping with the timing of recent human population expansion and any selective heterozygote advantage of these alleles requires further evaluation.