Female-to-Male Breeding Ratio in Modern Humans—an Analysis Based on Historical Recombinations

Was the past genetic contribution of women and men to the current human population equal? Was polygyny (excess of breeding women) present among hominid lineages? We addressed these questions by measuring the ratio of population recombination rates between the X chromosome and the autosomes, ρ_X/ρ_A. The X chromosome recombines only in female meiosis, whereas autosomes undergo crossovers in both sexes; thus, ρ_X/ρ_A reflects the female-to-male breeding ratio, β. We estimated β from ρ_X/ρ_A inferred from genomic diversity data and calibrated with recombination rates derived from pedigree data. For the HapMap populations, we obtained β of 1.4 in the Yoruba from West Africa, 1.3 in Europeans, and 1.1 in East Asian samples. These values are consistent with a high prevalence of monogamy and limited polygyny in human populations. More mutations occur during male meiosis as compared to female meiosis at the rate ratio referred to as α. We show that at α ≠ 1, the divergence rates and genetic diversities of the X chromosome relative to the autosomes are complex functions of both α and β, making their independent estimation difficult. Because our estimator of β does not require any knowledge of the mutation rates, our approach should allow us to dissociate the effects of α and β on the genetic diversity and divergence rate ratios of the sex chromosomes to the autosomes.     

Accelerated Adaptive Evolution on a Newly Formed X Chromosome

Sex chromosomes originated from ordinary autosomes, and their evolution is characterized by continuous gene loss from the ancestral Y chromosome. Here, we document a new feature of sex chromosome evolution: bursts of adaptive fixations on a newly formed X chromosome. Taking advantage of the recently formed neo-X chromosome of Drosophila miranda, we compare patterns of DNA sequence variation at genes located on the neo-X to genes on the ancestral X chromosome. This contrast allows us to draw inferences of selection on a newly formed X chromosome relative to background levels of adaptation in the genome while controlling for demographic effects. Chromosome-wide synonymous diversity on the neo-X is reduced 2-fold relative to the ancestral X, as expected under recent and recurrent directional selection. Several statistical tests employing various features of the data consistently identify 10%–15% of neo-X genes as targets of recent adaptive evolution but only 1%–3% of genes on the ancestral X. In addition, both the rate of adaptation and the fitness effects of adaptive substitutions are estimated to be roughly an order of magnitude higher for neo-X genes relative to genes on the ancestral X. Thus, newly formed X chromosomes are not passive players in the evolutionary process of sex chromosome differentiation, but respond adaptively to both their sex-biased transmission and to Y chromosome degeneration, possibly through demasculinization of their gene content and the evolution of dosage compensation.     

High genetic diversity at the extreme range edge: nucleotide variation at nuclear loci in Scots pine (Pinus sylvestris L.) in Scotland

Nucleotide polymorphism at 12 nuclear loci was studied in Scots pine populations across an environmental gradient in Scotland, to evaluate the impacts of demographic history and selection on genetic diversity. At eight loci, diversity patterns were compared between Scottish and continental European populations. At these loci, a similar level of diversity (θ_sil=∼0.01) was found in Scottish vs mainland European populations, contrary to expectations for recent colonization, however, less rapid decay of linkage disequilibrium was observed in the former (ρ=0.0086±0.0009, ρ=0.0245±0.0022, respectively). Scottish populations also showed a deficit of rare nucleotide variants (multi-locus Tajima''s D=0.316 vs D=−0.379) and differed significantly from mainland populations in allelic frequency and/or haplotype structure at several loci. Within Scotland, western populations showed slightly reduced nucleotide diversity (π_tot=0.0068) compared with those from the south and east (0.0079 and 0.0083, respectively) and about three times higher recombination to diversity ratio (ρ/θ=0.71 vs 0.15 and 0.18, respectively). By comparison with results from coalescent simulations, the observed allelic frequency spectrum in the western populations was compatible with a relatively recent bottleneck (0.00175 × 4N_e generations) that reduced the population to about 2% of the present size. However, heterogeneity in the allelic frequency distribution among geographical regions in Scotland suggests that subsequent admixture of populations with different demographic histories may also have played a role.     

Turnover of sex chromosomes and speciation in fishes

Closely related species of fishes often have different sex chromosome systems. Such rapid turnover of sex chromosomes can occur by several mechanisms, including fusions between an existing sex chromosome and an autosome. These fusions can result in a multiple sex chromosome system, where a species has both an ancestral and a neo-sex chromosome. Although this type of multiple sex chromosome system has been found in many fishes, little is known about the mechanisms that select for the formation of neo-sex chromosomes, or the role of neo-sex chromosomes in phenotypic evolution and speciation. The identification of closely related, sympatric species pairs in which one species has a multiple sex chromosome system and the other has a simple sex chromosome system provides an opportunity to study sex chromosome turnover. Recently, we found that a population of threespine stickleback (Gasterosteus aculeatus) from Japan has an X₁X₂Y multiple sex chromosome system resulting from a fusion between the ancestral Y chromosome and an autosome, while a sympatric threespine stickleback population has a simple XY sex chromosome system. Furthermore, we demonstrated that the neo-X chromosome (X ₂) plays an important role in phenotypic divergence and reproductive isolation between these sympatric stickleback species pairs. Here, we review multiple sex chromosome systems in fishes, as well as recent advances in our understanding of the evolutionary role of sex chromosome turnover in stickleback speciation.     

Energetics of Active Transport Processes

Discussions of active transport usually assume stoichiometry between the rate of transport J₊ and the metabolic rate J_r. However, the observation of a linear relationship between J₊ and J_r does not imply a stoichiometric relationship, i.e., complete coupling. Since coupling may possibly be incomplete, we examine systems of an arbitrary degree of coupling q, regarding stoichiometry as a limiting case. We consider a sodium pump, with J₊ and J_r linear functions of the electrochemical potential difference, -X₊, and the chemical affinity of the metabolic driving reaction, A. The affinity is well defined even for various complex reaction pathways. Incorporation of a series barrier and a parallel leak does not affect the linearity of the composite observable system. The affinity of some region of the metabolic chain may be maintained constant, either by large pools of reactants or by regulation. If so, this affinity can be evaluated by two independent methods. Sodium transport is conveniently characterized by the open-circuit potential (Δψ)_I=0 and the natural limits, level flow (J₊)_X+=0, and static head X⁰₊ = (X₊)_J+=0. With high degrees of coupling -X⁰₊/F approaches the electromotive force E_Na (Ussing); -X⁰₊/F cannot be identified with ((RT/F) ln f)_X+=0, where f is the flux ratio. The efficiency η = -J₊X₊/J_rA is of significance only when appreciable energy is being converted from one form to another. When either J₊ or -X₊ is small η is low; the significant parameters are then the efficacies ε_J+ = J₊/J_rA and ε_X+ = -X₊/J_rA, respectively maximal at level flow and static head. Leak increases both J₊ and ε_J+ for isotonic saline reabsorption, but diminishes -X⁰₊ and ε_X♀. Electrical resistance reflects both passive parameters and metabolism. Various fundamental relations are preserved despite coupling of passive ion and water flows.     

Purine receptors are required for DHA-mediated neuroprotection against oxygen and glucose deprivation in hippocampal slices

Docosahexaenoic acid (DHA) is important for central nervous system function during pathological states such as ischemia. DHA reduces neuronal injury in experimental brain ischemia; however, the underlying mechanisms are not well understood. In the present study, we investigated the effects of DHA on acute hippocampal slices subjected to experimental ischemia by transient oxygen and glucose deprivation (OGD) and re-oxygenation and the possible involvement of purinergic receptors as the mechanism underlying DHA-mediated neuroprotection. We observed that cellular viability reduction induced by experimental ischemia as well as cell damage and thiobarbituric acid reactive substances (TBARS) production induced by glutamate (10 mM) were prevented by hippocampal slices pretreated with DHA (5 μM). However, glutamate uptake reduction induced by OGD and re-oxygenation was not prevented by DHA. The beneficial effect of DHA against cellular viability reduction induced by OGD and re-oxygenation was blocked with PPADS (3 μM), a nonselective P2X_1–5 receptor antagonist as well as with a combination of TNP-APT (100 nM) plus brilliant blue (100 nM), which blocked P2X₁, P2X₃, P2X_2/3, and P2X₇ receptors, respectively. Moreover, adenosine receptors blockade with A₁ receptor antagonist DPCPX (100 nM) or with A_2B receptor antagonist alloxazine (100 nM) inhibited DHA-mediated neuroprotection. The addition of an A_2A receptor antagonist ZM241385 (50 nM), or A₃ receptor antagonist VUF5574 (1 μM) was ineffective. Taken together, our results indicated that neuroprotective actions of DHA may depend on P2X, A₁, and A_2B purinergic receptors activation. Our results reinforce the notion that dietary DHA may act as a local purinergic modulator in order to prevent neurodegenerative diseases.     

Solution structure of psi32-modified anticodon stem-loop of Escherichia coli tRNAPhe

Nucleoside base modifications can alter the structures and dynamics of RNA molecules and are important in tRNAs for maintaining translational fidelity and efficiency. The unmodified anticodon stem–loop from Escherichia coli tRNA^Phe forms a trinucleotide loop in solution, but Mg²⁺ and dimethylallyl modification of A₃₇ N6 destabilize the loop-proximal base pairs and increase the mobility of the loop nucleotides. The anticodon arm has three additional modifications, ψ₃₂, ψ₃₉, and A₃₇ C2-thiomethyl. We have used NMR spectroscopy to investigate the structural and dynamical effects of ψ₃₂ on the anticodon stem-loop from E.coli tRNA^Phe. The ψ₃₂ modification does not significantly alter the structure of the anticodon stem–loop relative to the unmodified parent molecule. The stem of the RNA molecule includes base pairs ψ₃₂-A₃₈ and U₃₃–A₃₇ and the base of ψ₃₂ stacks between U₃₃ and A₃₁. The glycosidic bond of ψ₃₂ is in the anti configuration and is paired with A₃₈ in a Watson–Crick geometry, unlike residue 32 in most crystal structures of tRNA. The ψ₃₂ modification increases the melting temperature of the stem by ~3.5°C, although the ψ₃₂ and U₃₃ imino resonances are exchange broadened. The results suggest that ψ₃₂ functions to preserve the stem integrity in the presence of additional loop modifications or after reorganization of the loop into a translationally functional conformation.     

The Genomic Selfing Syndrome Accompanies the Evolutionary Breakdown of Heterostyly

The evolutionary transition from outcrossing to selfing can have important genomic consequences. Decreased effective population size and the reduced efficacy of selection are predicted to play an important role in the molecular evolution of the genomes of selfing species. We investigated evidence for molecular signatures of the genomic selfing syndrome using 66 species of Primula including distylous (outcrossing) and derived homostylous (selfing) taxa. We complemented our comparative analysis with a microevolutionary study of P. chungensis, which is polymorphic for mating system and consists of both distylous and homostylous populations. We generated chloroplast and nuclear genomic data sets for distylous, homostylous, and distylous–homostylous species and identified patterns of nonsynonymous to synonymous divergence (d_N/d_S) and polymorphism (π_N/π_S) in species or lineages with contrasting mating systems. Our analysis of coding sequence divergence and polymorphism detected strongly reduced genetic diversity and heterozygosity, decreased efficacy of purifying selection, purging of large-effect deleterious mutations, and lower rates of adaptive evolution in samples from homostylous compared with distylous populations, consistent with theoretical expectations of the genomic selfing syndrome. Our results demonstrate that self-fertilization is a major driver of molecular evolutionary processes with genomic signatures of selfing evident in both old and relatively young homostylous populations.     

Hydration-State Change of Horse Heart Cytochrome c Corresponding to Trifluoroacetic-Acid-Induced Unfolding

We investigate the hydration state of horse-heart cytochrome c (hh cyt c) in the unfolding process induced by trifluoroacetic acid (TFA). The conformation of hh cyt c changes from the native (N) state (2.9 < pH < 6.0) to the acid-unfolded (U_A) state (1.7 < pH < 2.0) to the acid-induced molten globule (A) state (pH ∼1.2). Hydration properties of hh cyt c during this process are measured at 20°C by high-resolution dielectric relaxation (DR) spectroscopy, UV-vis absorbance, and circular dichroism spectroscopy. Constrained water of hh cyt c is observed at every pH as an ∼5-GHz Debye component (DC) (DR time, τ_D ∼30 ps) and its DR amplitude (DRA) is increased by 77% upon N-to-U_A transition, when pH changes from 6.0 to 2.0. Even in the N state, the DRA of the constrained-water component is found to be increased by 22% with decreasing pH from 6.0 to 2.9, suggesting an increase in the accessible surface area of native hh cyt c. Moreover, hypermobile water around native hh cyt c is detected at pH 6.0 as a 19-GHz DC (τ_D ∼ 8.4 ps < τ_DW = 9.4 ps), but is not found at other pH values. The DRA signal of constrained water is found to return to the pH 2.9 (N-state) level upon U_A-to-A transition. Fast-response water (slightly slower than bulk) around A-state hh cyt c is detected at pH 1.2, and this suggests some accumulation of TFA⁻ ions around the peptide chain. Thus, this high-resolution DR spectroscopy study reveals that hh cyt c exhibits significant hydration-state change in the TFA-unfolding process.     

MULTIPLE ORIGINS OF SEX CHROMOSOME FUSIONS CORRELATED WITH CHIASMA LOCALIZATION IN HABRONATTUS JUMPING SPIDERS (ARANEAE: SALTICIDAE)

Entelegyne spiders rarely show fusions yielding neo‐Y chromosomes, which M. J. D. White attributed to a constraint in spiders, namely their proximal chiasma localization acting to upset meiotic segregation in males with fusions. Of the 75 taxa of Habronattus and outgroups studied, 47 have X₁X₂0 sex chromosomes in males, 10 have X₁X₂Y, 15 have X₁X₂X₃Y, 2 have X0, and one has both X₁X₂0 and X₁X₂X₃Y. Chromosome numbers and behavior suggest neo‐Ys formed by an autosome‐X fusion to make X₁X₂Y, with a second fusion to an autosome to make X₁X₂X₃Y. Phylogeny shows at least 8–15 gains (or possibly some losses) of neo‐Y (i.e., X‐autosome fusions), a remarkable number for such a small clade. In contrast to the many X‐autosome fusions, at most one autosome–autosome fusion is indicated. Origins of neo‐Y are correlated significantly with distal localization of chiasmata, supporting White's hypothesis that evolution of neo‐Y systems is facilitated by looser pairing (distal chiasmata) at meiosis. However, an alternative (or contributing) explanation for the correlation is that X‐autosome fusions were selected to permit isolation of male‐favored alleles to the neo‐Y chromosome, aided by distal chiasmata limiting recombination. This intralocus sexual conflict hypothesis could explain both the many X‐autosome fusions, and the stunning complexity of male Habronattus courtship displays.     

From Discrete to Continuum Models of Three-Dimensional Deformations in Epithelial Sheets

Metamorphic proteins, including proteins with high levels of sequence identity but different folds, are exceptions to the long-standing rule-of-thumb that proteins with as little as 30% sequence identity adopt the same fold. Which topologies can be bridged by these highly identical sequences remains an open question. Here we bridge two 3-α-helix bundle proteins with two radically different folds. Using a straightforward approach, we engineered the sequences of one subdomain within maltose binding protein (MBP, α/β/α-sandwich) and another within outer surface protein A (OspA, β-sheet) to have high sequence identity (80 and 77%, respectively) with engineered variants of protein G (G_A, 3-α-helix bundle). Circular dichroism and nuclear magnetic resonance spectra of all engineered variants demonstrate that they maintain their native conformations despite substantial sequence modification. Furthermore, the MBP variant (80% identical to G_A) remained active. Thermodynamic analysis of numerous G_A and MBP variants suggests that the key to our approach involved stabilizing the modified MBP and OspA subdomains via external interactions with neighboring substructures, indicating that subdomain interactions can stabilize alternative folds over a broad range of sequence variation. These findings suggest that it is possible to bridge one fold with many other topologies, which has implications for protein folding, evolution, and misfolding diseases.     

Nucleotide Variation, Linkage Disequilibrium and Founder-Facilitated Speciation in Wild Populations of the Zebra Finch (Taeniopygia guttata)

The zebra finch has long been an important model system for the study of vocal learning, vocal production, and behavior. With the imminent sequencing of its genome, the zebra finch is now poised to become a model system for population genetics. Using a panel of 30 noncoding loci, we characterized patterns of polymorphism and divergence among wild zebra finch populations. Continental Australian populations displayed little population structure, exceptionally high levels of nucleotide diversity (π = 0.010), a rapid decay of linkage disequilibrium (LD), and a high population recombination rate (ρ ≈ 0.05), all of which suggest an open and fluid genomic background that could facilitate adaptive variation. By contrast, substantial divergence between the Australian and Lesser Sunda Island populations (K_ST = 0.193), reduced genetic diversity (π = 0.002), and higher levels of LD in the island population suggest a strong but relatively recent founder event, which may have contributed to speciation between these populations as envisioned under founder-effect speciation models. Consistent with this hypothesis, we find that under a simple quantitative genetic model both drift and selection could have contributed to the observed divergence in six quantitative traits. In both Australian and Lesser Sundas populations, diversity in Z-linked loci was significantly lower than in autosomal loci. Our analysis provides a quantitative framework for studying the role of selection and drift in shaping patterns of molecular evolution in the zebra finch genome.     

Step-by-step evolution of neo-sex chromosomes in geographical populations of wild silkmoths,Samia cynthia ssp.

Geographical subspecies of wild silkmoths, Samia cynthia ssp. (Lepidoptera: Saturniidae), differ considerably in sex chromosome constitution owing to sex chromosome fusions with autosomes, which leads to variation in chromosome numbers. We cloned S. cynthia orthologues of 16 Bombyx mori genes and mapped them to chromosome spreads of S. cynthia subspecies by fluorescence in situ hybridization (FISH) to determine the origin of S. cynthia neo-sex chromosomes. FISH mapping revealed that the Z chromosome and chromosome 12 of B. mori correspond to the Z chromosome and an autosome (A₁) of S. c. ricini (Vietnam population, 2n=27, Z0 in female moths), respectively. B. mori chromosome 11 corresponds partly to another autosome (A₂) and partly to a chromosome carrying nucleolar organizer region (NOR) of this subspecies. The NOR chromosome of S. c. ricini is also partly homologous to B. mori chromosome 24. Furthermore, our results revealed that two A₁ homologues each fused with the W and Z chromosomes in a common ancestor of both Japanese subspecies S. c. walkeri (Sapporo population, 2n=26, neo-Wneo-Z) and S. cynthia subsp. indet. (Nagano population, 2n=25, neo-WZ₁Z₂). One homologue, corresponding to the A₂ autosome in S. c. ricini and S. c. walkeri, fused with the W chromosome in S. cynthia subsp. indet. Consequently, the other homologue became a Z₂ chromosome. These results clearly showed a step-by-step evolution of the neo-sex chromosomes by repeated autosome–sex chromosome fusions. We suggest that the rearrangements of sex chromosomes may facilitate divergence of S. cynthia subspecies towards speciation.     

Low levels of realized seed and pollen gene flow and strong spatial genetic structure in a small,isolated and fragmented population of the tropical tree Copaifera langsdorffii Desf

Over the past century, the Brazilian Atlantic forest has been reduced to small, isolated fragments of forest. Reproductive isolation theories predict a loss of genetic diversity and increases in inbreeding and spatial genetic structure (SGS) in such populations. We analysed eight microsatellite loci to investigate the pollen and seed dispersal patterns, genetic diversity, inbreeding and SGS of the tropical tree Copaifera langsdorffii in a small (4.8 ha), isolated population. All 112 adult trees and 128 seedlings found in the stand were sampled, mapped and genotyped. Seedlings had significantly lower levels of genetic diversity (A=16.5±0.45, mean±95% s.e.; H_e=0.838±0.006) than did adult trees (A=23.2±0.81; H_e=0.893±0.030). Parentage analysis did not indicate any seed immigration (m_seeds=0) and the pollen immigration rate was very low (m_pollen=0.047). The average distance of realized pollen dispersal within the stand was 94 m, with 81% of the pollen travelling <150 m. A significant negative correlation was found between the frequency and distance of pollen dispersal (r=−0.79, P<0.01), indicating that short-distance pollinations were more frequent. A significant SGS for both adults (∼50 m) and seedlings (∼20 m) was also found, indicating that most of the seeds were dispersed over short distances. The results suggested that the spatial isolation of populations by habitat fragmentation can restrict seed and pollen gene flow, increase SGS and affect the genetic diversity of future generations.     

Low Level of Sequence Diversity at Merozoite Surface Protein-1 Locus of Plasmodium ovale curtisi and P. ovale wallikeri from Thai Isolates

Background

The merozoite surface protein-1 (MSP-1) is a candidate target for the development of blood stage vaccines against malaria. Polymorphism in MSP-1 can be useful as a genetic marker for strain differentiation in malarial parasites. Although sequence diversity in the MSP-1 locus has been extensively analyzed in field isolates of Plasmodium falciparum and P. vivax, the extent of variation in its homologues in P. ovale curtisi and P. ovale wallikeri, remains unknown.

Methodology/Principal Findings

Analysis of the mitochondrial cytochrome b sequences of 10 P. ovale isolates from symptomatic malaria patients from diverse endemic areas of Thailand revealed co-existence of P. ovale curtisi (n = 5) and P. ovale wallikeri (n = 5). Direct sequencing of the PCR-amplified products encompassing the entire coding region of MSP-1 of P. ovale curtisi (PocMSP-1) and P. ovale wallikeri (PowMSP-1) has identified 3 imperfect repeated segments in the former and one in the latter. Most amino acid differences between these proteins were located in the interspecies variable domains of malarial MSP-1. Synonymous nucleotide diversity (π_S) exceeded nonsynonymous nucleotide diversity (π_N) for both PocMSP-1 and PowMSP-1, albeit at a non-significant level. However, when MSP-1 of both these species was considered together, π_S was significantly greater than π_N (p<0.0001), suggesting that purifying selection has shaped diversity at this locus prior to speciation. Phylogenetic analysis based on conserved domains has placed PocMSP-1 and PowMSP-1 in a distinct bifurcating branch that probably diverged from each other around 4.5 million years ago.

Conclusion/Significance

The MSP-1 sequences support that P. ovale curtisi and P. ovale wallikeri are distinct species. Both species are sympatric in Thailand. The low level of sequence diversity in PocMSP-1 and PowMSP-1 among Thai isolates could stem from persistent low prevalence of these species, limiting the chance of outcrossing at this locus.     

Alzheimer’s Protective A2T Mutation Changes the Conformational Landscape of the Aβ1–42 Monomer Differently Than Does the A2V Mutation

The aggregation of amyloid-β (Aβ) peptides plays a crucial role in the etiology of Alzheimer’s disease (AD). Recently, it has been reported that an A2T mutation in Aβ can protect against AD. Interestingly, a nonpolar A2V mutation also has been found to offer protection against AD in the heterozygous state, although it causes early-onset AD in homozygous carriers. Since the conformational landscape of the Aβ monomer is known to directly contribute to the early-stage aggregation mechanism, it is important to characterize the effects of the A2T and A2V mutations on Aβ_1–42 monomer structure. Here, we have performed extensive atomistic replica-exchange molecular dynamics simulations of the solvated wild-type (WT), A2V, and A2T Aβ_1–42 monomers. Our simulations reveal that although all three variants remain as collapsed coils in solution, there exist significant structural differences among them at shorter timescales. A2V exhibits an enhanced double-hairpin population in comparison to the WT, similar to those reported in toxic WT Aβ_1–42 oligomers. Such double-hairpin formation is caused by hydrophobic clustering between the N-terminus and the central and C-terminal hydrophobic patches. In contrast, the A2T mutation causes the N-terminus to engage in unusual electrostatic interactions with distant residues, such as K16 and E22, resulting in a unique population comprising only the C-terminal hairpin. These findings imply that a single A2X (where X = V or T) mutation in the primarily disordered N-terminus of the Aβ_1–42 monomer can dramatically alter the β-hairpin population and switch the equilibrium toward alternative structures. The atomistically detailed, comparative view of the structural landscapes of A2V and A2T variant monomers obtained in this study can enhance our understanding of the mechanistic differences in their early-stage aggregation.