A Parent-of-Origin Effect Impacts the Phenotype in Low Penetrance Retinoblastoma Families Segregating the c.1981C>T/p.Arg661Trp Mutation of RB1

Retinoblastoma (Rb), the most common pediatric intraocular neoplasm, results from inactivation of both alleles of the RB1 tumor suppressor gene. The second allele is most commonly lost, as demonstrated by loss of heterozygosity studies. RB1 germline carriers usually develop bilateral tumors, but some Rb families display low penetrance and variable expressivity. In order to decipher the underlying mechanisms, 23 unrelated low penetrance pedigrees segregating the common c.1981C>T/p.Arg661Trp mutation and other low penetrance mutations were studied. In families segregating the c.1981C>T mutation, we demonstrated, for the first time, a correlation between the gender of the transmitting carrier and penetrance, as evidenced by Fisher’s exact test: the probability of being unaffected is 90.3% and 32.5% when the mutation is inherited from the mother and the father, respectively (p-value = 7.10⁻⁷). Interestingly, a similar correlation was observed in families segregating other low penetrance alleles. Consequently, we investigated the putative involvement of an imprinted, modifier gene in low penetrance Rb. We first ruled out a MED4-driven mechanism by MED4 methylation and expression analyses. We then focused on the differentially methylated CpG85 island located in intron 2 of RB1 and showing parent-of-origin-specific DNA methylation. This differential methylation promotes expression of the maternal c.1981C>T allele. We propose that the maternally inherited c.1981C>T/p.Arg661Trp allele retains sufficient tumor suppressor activity to prevent retinoblastoma development. In contrast, when the mutation is paternally transmitted, the low residual activity would mimic a null mutation and subsequently lead to retinoblastoma. This implies that the c.1981C>T mutation is not deleterious per se but needs to be destabilized in order to reach pRb haploinsufficiency and initiate tumorigenesis. We suggest that this phenomenon might be a general mechanism to explain phenotypic differences in low penetrance Rb families.     

Rate of Gene Silencing at Duplicate Loci: A Theoretical Study and Interpretation of Data from Tetraploid Fishes

A large-scale simulation has been conducted on the rate of gene loss at duplicate loci under irreversible mutation. It is found that tight linkage does not provide a strong sheltering effect, as thought by previous authors; indeed, the mean loss time for the case of tight linkage is of the same order of magnitude as that for no linkage, as long as Nu is not much larger than 1, where N is the effective population size and u the mutation rate. When Nu is 0.01 or less, the two loci behave almost as neutral loci, regardless of linkage, and the mean loss time is about only half the mean extinction time for a neutral allele under irreversible mutation. However, the former becomes two or more times larger than the latter when Nu ≥ 1.——In the simulation, the sojourn times in the frequency intervals (0, 0.01) and (0.99, 1) and the time for the frequency of the null allele to reach 0.99 at one of the two loci have also been recorded. The results show that the population is monomorphic for the normal allele most of the time if Nu ≤ 0.01, but polymorphic for the null and the normal alleles most of the time if Nu ≥ 0.1.——The distribution of the frequency of the null allele in an equilibrium tetraploid population has been studied analytically. The present results have been applied to interpret data from some fish groups that are of tetraploid origin, and a model for explaining the slow rate of gene loss in these fishes is proposed.     

Population dynamics of sex-determining alleles in honey bees and self-incompatibility alleles in plants

Mathematical theories of the population dynamics of sex-determining alleles in honey bees are developed. It is shown that in an infinitely large population the equilibrium frequency of a sex allele is 1/n, where n is the number of alleles in the population, and the asymptotic rate of approach to this equilibrium is 2/(3n) per generation. Formulae for the distribution of allele frequencies and the effective and actual numbers of alleles that can be maintained in a finite population are derived by taking into account the population size and mutation rate. It is shown that the allele frequencies in a finite population may deviate considerably from 1/n. Using these results, available data on the number of sex alleles in honey bee populations are discussed. It is also shown that the number of self-incompatibility alleles in plants can be studied in a much simpler way by the method used in this paper. A brief discussion about general overdominant selection is presented.     

Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits

The obesity epidemic is responsible for a substantial economic burden in developed countries and is a major risk factor for type 2 diabetes and cardiovascular disease. The disease is the result not only of several environmental risk factors, but also of genetic predisposition. To take advantage of recent advances in gene-mapping technology, we executed a genome-wide association scan to identify genetic variants associated with obesity-related quantitative traits in the genetically isolated population of Sardinia. Initial analysis suggested that several SNPs in the FTO and PFKP genes were associated with increased BMI, hip circumference, and weight. Within the FTO gene, rs9930506 showed the strongest association with BMI (p = 8.6 ×10⁻⁷), hip circumference (p = 3.4 × 10⁻⁸), and weight (p = 9.1 × 10⁻⁷). In Sardinia, homozygotes for the rare “G” allele of this SNP (minor allele frequency = 0.46) were 1.3 BMI units heavier than homozygotes for the common “A” allele. Within the PFKP gene, rs6602024 showed very strong association with BMI (p = 4.9 × 10⁻⁶). Homozygotes for the rare “A” allele of this SNP (minor allele frequency = 0.12) were 1.8 BMI units heavier than homozygotes for the common “G” allele. To replicate our findings, we genotyped these two SNPs in the GenNet study. In European Americans (N = 1,496) and in Hispanic Americans (N = 839), we replicated significant association between rs9930506 in the FTO gene and BMI (p-value for meta-analysis of European American and Hispanic American follow-up samples, p = 0.001), weight (p = 0.001), and hip circumference (p = 0.0005). We did not replicate association between rs6602024 and obesity-related traits in the GenNet sample, although we found that in European Americans, Hispanic Americans, and African Americans, homozygotes for the rare “A” allele were, on average, 1.0–3.0 BMI units heavier than homozygotes for the more common “G” allele. In summary, we have completed a whole genome–association scan for three obesity-related quantitative traits and report that common genetic variants in the FTO gene are associated with substantial changes in BMI, hip circumference, and body weight. These changes could have a significant impact on the risk of obesity-related morbidity in the general population.     

The Characteristic Trajectory of a Fixing Allele: A Consequence of Fictitious Selection That Arises from Conditioning

This work is concerned with the historical progression, to fixation, of an allele in a finite population. This progression is characterized by the average frequency trajectory of alleles that achieve fixation before a given time, T. Under a diffusion analysis, the average trajectory, conditional on fixation by time T, is shown to be equivalent to the average trajectory in an unconditioned problem involving additional selection. We call this additional selection “fictitious selection”; it plays the role of a selective force in the unconditioned problem but does not exist in reality. It is a consequence of conditioning on fixation. The fictitious selection is frequency dependent and can be very large compared with any real selection that is acting. We derive an approximation for the characteristic trajectory of a fixing allele, when subject to real additive selection, from an unconditioned problem, where the total selection is a combination of real and fictitious selection. Trying to reproduce the characteristic trajectory from the action of additive selection, in an infinite population, can lead to estimates of the strength of the selection that deviate from the real selection by >1000% or have the opposite sign. Strong evolutionary forces may be invoked in problems where conditioning has been carried out, but these forces may largely be an outcome of the conditioning and hence may not have a real existence. The work presented here clarifies these issues and provides two useful tools for future analyses: the characteristic trajectory of a fixing allele and the force that primarily drives this, namely fictitious selection. These should prove useful in a number of areas of interest including coalescence with selection, experimental evolution, time series analyses of ancient DNA, game theory in finite populations, and the historical dynamics of selected alleles in wild populations.     

Evidence of Natural Selection Acting on a Polymorphic Hybrid Incompatibility Locus in Mimulus

As a common cause of reproductive isolation in diverse taxa, hybrid incompatibilities are fundamentally important to speciation. A key question is which evolutionary forces drive the initial substitutions within species that lead to hybrid dysfunction. Previously, we discovered a simple genetic incompatibility that causes nearly complete male sterility and partial female sterility in hybrids between the two closely related yellow monkeyflower species Mimulus guttatus and M. nasutus. In this report, we fine map the two major incompatibility loci—hybrid male sterility 1 (hms1) and hybrid male sterility 2 (hms2)—to small nuclear genomic regions (each <70 kb) that include strong candidate genes. With this improved genetic resolution, we also investigate the evolutionary dynamics of hms1 in a natural population of M. guttatus known to be polymorphic at this locus. Using classical genetic crosses and population genomics, we show that a 320-kb region containing the hms1 incompatibility allele has risen to intermediate frequency in this population by strong natural selection. This finding provides direct evidence that natural selection within plant species can lead to hybrid dysfunction between species.     

Mutation Spectrum of Common Deafness-Causing Genes in Patients with Non-Syndromic Deafness in the Xiamen Area,China

In China, approximately 30,000 babies are born with hearing impairment each year. However, the molecular factors causing congenital hearing impairment in the Xiamen area of Fujian province have not been evaluated. To provide accurate genetic testing and counseling in the Xiamen area, we investigated the molecular etiology of non-syndromic deafness in a deaf population from Xiamen. Unrelated students with hearing impairment (n = 155) who attended Xiamen Special Education School in Fujian Province were recruited for this study. Three common deafness-related genes, GJB2, SLC26A4, and mtDNA12SrRNA, were analyzed using all-exon sequencing. GJB2 mutations were detected in 27.1% (42/155) of the entire cohort. The non-syndromic hearing loss (NSHL) hotspot mutations c.109G>A (p.V37I) and c.235delC were found in this population, whereas the Caucasian hotspot mutation c.35delG was not. The allelic frequency of the c.109G>A mutation was 9.03% (28/310), slightly higher than that of c.235delC (8.39%, 26/310), which is the most common GJB2 mutation in most areas of China. The allelic frequency of the c.109G>A mutation was significantly higher in this Xiamen’s deaf population than that in previously reported cohorts (P = 0.00). The SLC26A4 mutations were found in 16.77% (26/155) of this cohort. The most common pathogenic allele was c.IVS7-2A>G (6.13%, 19/310), and the second most common was the c.1079C>T (p.A360V) mutation (1.94%, 6/310) which has rarely been reported as a hotspot mutation in other studies. The mutation rate of mtDNA12SrRNA in this group was 3.87% (6/155), all being the m.A1555G mutation. These findings show the specificity of the common deaf gene-mutation spectrum in this area. According to this study, there were specific hotspot mutations in Xiamen deaf patients. Comprehensive sequencing analysis of the three common deaf genes can help portray the mutation spectrum and develop optimal testing strategies for deaf patients in this area.     

Deficient and Null Variants of SERPINA1 Are Proteotoxic in a Caenorhabditis elegans Model of α1-Antitrypsin Deficiency

α1-antitrypsin deficiency (ATD) predisposes patients to both loss-of-function (emphysema) and gain-of-function (liver cirrhosis) phenotypes depending on the type of mutation. Although the Z mutation (ATZ) is the most prevalent cause of ATD, >120 mutant alleles have been identified. In general, these mutations are classified as deficient (<20% normal plasma levels) or null (<1% normal levels) alleles. The deficient alleles, like ATZ, misfold in the ER where they accumulate as toxic monomers, oligomers and aggregates. Thus, deficient alleles may predispose to both gain- and loss-of-function phenotypes. Null variants, if translated, typically yield truncated proteins that are efficiently degraded after being transiently retained in the ER. Clinically, null alleles are only associated with the loss-of-function phenotype. We recently developed a C. elegans model of ATD in order to further elucidate the mechanisms of proteotoxicity (gain-of-function phenotype) induced by the aggregation-prone deficient allele, ATZ. The goal of this study was to use this C. elegans model to determine whether different types of deficient and null alleles, which differentially affect polymerization and secretion rates, correlated to any extent with proteotoxicity. Animals expressing the deficient alleles, Mmalton, Siiyama and S (ATS), showed overall toxicity comparable to that observed in patients. Interestingly, Siiyama expressing animals had smaller intracellular inclusions than ATZ yet appeared to have a greater negative effect on animal fitness. Surprisingly, the null mutants, although efficiently degraded, showed a relatively mild gain-of-function proteotoxic phenotype. However, since null variant proteins are degraded differently and do not appear to accumulate, their mechanism of proteotoxicity is likely to be different to that of polymerizing, deficient mutants. Taken together, these studies showed that C. elegans is an inexpensive tool to assess the proteotoxicity of different AT variants using a transgenic approach.     

Soft Selective Sweeps in Complex Demographic Scenarios

Adaptation from de novo mutation can produce so-called soft selective sweeps, where adaptive alleles of independent mutational origin sweep through the population at the same time. Population genetic theory predicts that such soft sweeps should be likely if the product of the population size and the mutation rate toward the adaptive allele is sufficiently large, such that multiple adaptive mutations can establish before one has reached fixation; however, it remains unclear how demographic processes affect the probability of observing soft sweeps. Here we extend the theory of soft selective sweeps to realistic demographic scenarios that allow for changes in population size over time. We first show that population bottlenecks can lead to the removal of all but one adaptive lineage from an initially soft selective sweep. The parameter regime under which such “hardening” of soft selective sweeps is likely is determined by a simple heuristic condition. We further develop a generalized analytical framework, based on an extension of the coalescent process, for calculating the probability of soft sweeps under arbitrary demographic scenarios. Two important limits emerge within this analytical framework: In the limit where population-size fluctuations are fast compared to the duration of the sweep, the likelihood of soft sweeps is determined by the harmonic mean of the variance effective population size estimated over the duration of the sweep; in the opposing slow fluctuation limit, the likelihood of soft sweeps is determined by the instantaneous variance effective population size at the onset of the sweep. We show that as a consequence of this finding the probability of observing soft sweeps becomes a function of the strength of selection. Specifically, in species with sharply fluctuating population size, strong selection is more likely to produce soft sweeps than weak selection. Our results highlight the importance of accurate demographic estimates over short evolutionary timescales for understanding the population genetics of adaptation from de novo mutation.     

Ages,extinction times,and first passage probabilities for a multiallele diffusion model with irreversible mutation

We derive some new results for diffusion models in population genetics of the “infinite-alleles” type. Assuming present allelic frequencies known, we find expected values for times which may be interpreted either as the age of the kth oldest allele or the time to extinction of the “kth-to-last” allele to be lost. We also find some conditional expectations and probabilities related to the order of extinction or creation of the alleles. Computations illustrate the dependence of the expected times on the mutation rate.     

Worldwide Distribution of the MYH9 Kidney Disease Susceptibility Alleles and Haplotypes: Evidence of Historical Selection in Africa

MYH9 was recently identified as renal susceptibility gene (OR 3–8, p<10⁻⁸) for major forms of kidney disease disproportionately affecting individuals of African descent. The risk haplotype (E-1) occurs at much higher frequencies in African Americans (≥60%) than in European Americans (<4%), revealing a genetic basis for a major health disparity. The population distributions of MYH9 risk alleles and the E-1 risk haplotype and the demographic and selective forces acting on the MYH9 region are not well explored. We reconstructed MYH9 haplotypes from 4 tagging single nucleotide polymorphisms (SNPs) spanning introns 12–23 using available data from HapMap Phase II, and by genotyping 938 DNAs from the Human Genome Diversity Panel (HGDP). The E-1 risk haplotype followed a cline, being most frequent within sub-Saharan African populations (range 50–80%), less frequent in populations from the Middle East (9–27%) and Europe (0–9%), and rare or absent in Asia, the Americas, and Oceania. The fixation indexes (F_ST) for pairwise comparisons between the risk haplotypes for continental populations were calculated for MYH9 haplotypes; F_ST ranged from 0.27–0.40 for Africa compared to other continental populations, possibly due to selection. Uniquely in Africa, the Yoruba population showed high frequency extended haplotype length around the core risk allele (C) compared to the alternative allele (T) at the same locus (rs4821481, iHs = 2.67), as well as high population differentiation (F_{ST(CEU vs. YRI)} = 0.51) in HapMap Phase II data, also observable only in the Yoruba population from HGDP (F_ST = 0.49), pointing to an instance of recent selection in the genomic region. The population-specific divergence in MYH9 risk allele frequencies among the world''s populations may prove important in risk assessment and public health policies to mitigate the burden of kidney disease in vulnerable populations.     

Reversibility and the age of an allele. I. Moran's infinitely many neutral alleles model

A population of constant size is subjected to mutation, such that each mutant is of a new allelic type. For the particular population model studied in this paper, the age of an allele, whose present frequency is known, is a random variable with distribution independent of the frequencies of other alleles. As a consequence of reversibility of the population process, the age of an allele, from the past to the present, has the same distribution as its time to extinction, from the present into the future. This verifies, and re-interprets, certain diffusion approximations found by Kimura and Ohta [Genetics 75, 199–212 (1973)] and Maruyama [Genet. Res. Cambridge 23, 137–143 (1974)].     

Mosaic and Intronic Mutations in TSC1/TSC2 Explain the Majority of TSC Patients with No Mutation Identified by Conventional Testing

Tuberous sclerosis complex (TSC) is an autosomal dominant tumor suppressor gene syndrome due to germline mutations in either TSC1 or TSC2. 10–15% of TSC individuals have no mutation identified (NMI) after thorough conventional molecular diagnostic assessment. 53 TSC subjects who were NMI were studied using next generation sequencing to search for mutations in these genes. Blood/saliva DNA including parental samples were available from all subjects, and skin tumor biopsy DNA was available from six subjects. We identified mutations in 45 of 53 subjects (85%). Mosaicism was observed in the majority (26 of 45, 58%), and intronic mutations were also unusually common, seen in 18 of 45 subjects (40%). Seventeen (38%) mutations were seen at an allele frequency < 5%, five at an allele frequency < 1%, and two were identified in skin tumor biopsies only, and were not seen at appreciable frequency in blood or saliva DNA. These findings illuminate the extent of mosaicism in TSC, indicate the importance of full gene coverage and next generation sequencing for mutation detection, show that analysis of TSC-related tumors can increase the mutation detection rate, indicate that it is not likely that a third TSC gene exists, and enable provision of genetic counseling to the substantial population of TSC individuals who are currently NMI.     

Statistical Mechanics and the Evolution of Polygenic Quantitative Traits

The evolution of quantitative characters depends on the frequencies of the alleles involved, yet these frequencies cannot usually be measured. Previous groups have proposed an approximation to the dynamics of quantitative traits, based on an analogy with statistical mechanics. We present a modified version of that approach, which makes the analogy more precise and applies quite generally to describe the evolution of allele frequencies. We calculate explicitly how the macroscopic quantities (i.e., quantities that depend on the quantitative trait) depend on evolutionary forces, in a way that is independent of the microscopic details. We first show that the stationary distribution of allele frequencies under drift, selection, and mutation maximizes a certain measure of entropy, subject to constraints on the expectation of observable quantities. We then approximate the dynamical changes in these expectations, assuming that the distribution of allele frequencies always maximizes entropy, conditional on the expected values. When applied to directional selection on an additive trait, this gives a very good approximation to the evolution of the trait mean and the genetic variance, when the number of mutations per generation is sufficiently high (4Nμ > 1). We show how the method can be modified for small mutation rates (4Nμ → 0). We outline how this method describes epistatic interactions as, for example, with stabilizing selection.     

Development and characterization of microsatellite markers in sawih tree (Duabanga moluccana Blume) using ISSR-suppression PCR techniques

Duabanga moluccana (or locally known as sawih) is an indigenous fast growing tropical tree species that confers various advantages for the timber industry and for planted forests development. In this paper, we isolated and characterized 8 polymorphic microsatellite markers from the D. moluccana genome using ISSR-suppression PCR techniques. The number of alleles and PIC values ranged from 3 to 8 alleles per locus and from 0.488 to 0.792, respectively. Three microsatellite loci were deviated from Hardy-Weinberg equilibrium (P < 0.05). The transferability rate ranged from 24 to 100 % among the three indigenous tree species tested. This indicates that the newly developed microsatellite markers would be useful tools for population genetic studies on D. moluccana and other indigenous tree species.     

Identification of a Novel Homozygous Mutation,TMPRSS3: c.535G>A,in a Tibetan Family with Autosomal Recessive Non-Syndromic Hearing Loss

Different ethnic groups have distinct mutation spectrums associated with inheritable deafness. In order to identify the mutations responsible for congenital hearing loss in the Tibetan population, mutation screening for 98 deafness-related genes by microarray and massively parallel sequencing of captured target exons was conducted in one Tibetan family with familiar hearing loss. A homozygous mutation, TMPRSS3: c.535G>A, was identified in two affected brothers. Both parents are heterozygotes and an unaffected sister carries wild type alleles. The same mutation was not detected in 101 control Tibetan individuals. This missense mutation results in an amino acid change (p.Ala179Thr) at a highly conserved site in the scavenger receptor cysteine rich (SRCR) domain of the TMPRSS3 protein, which is essential for protein-protein interactions. Thus, this mutation likely affects the interactions of this transmembrane protein with extracellular molecules. According to our bioinformatic analyses, the TMPRSS3: c.535G>A mutation might damage protein function and lead to hearing loss. These data suggest that the homozygous mutation TMPRSS3: c.535G>A causes prelingual hearing loss in this Tibetan family. This is the first TMPRSS3 mutation found in the Chinese Tibetan population.