Sex differences in mathematics and reading achievement are inversely related: within- and across-nation assessment of 10 years of PISA data

We analyzed one decade of data collected by the Programme for International Student Assessment (PISA), including the mathematics and reading performance of nearly 1.5 million 15 year olds in 75 countries. Across nations, boys scored higher than girls in mathematics, but lower than girls in reading. The sex difference in reading was three times as large as in mathematics. There was considerable variation in the extent of the sex differences between nations. There are countries without a sex difference in mathematics performance, and in some countries girls scored higher than boys. Boys scored lower in reading in all nations in all four PISA assessments (2000, 2003, 2006, 2009). Contrary to several previous studies, we found no evidence that the sex differences were related to nations’ gender equality indicators. Further, paradoxically, sex differences in mathematics were consistently and strongly inversely correlated with sex differences in reading: Countries with a smaller sex difference in mathematics had a larger sex difference in reading and vice versa. We demonstrate that this was not merely a between-nation, but also a within-nation effect. This effect is related to relative changes in these sex differences across the performance continuum: We did not find a sex difference in mathematics among the lowest performing students, but this is where the sex difference in reading was largest. In contrast, the sex difference in mathematics was largest among the higher performing students, and this is where the sex difference in reading was smallest. The implication is that if policy makers decide that changes in these sex differences are desired, different approaches will be needed to achieve this for reading and mathematics. Interventions that focus on high-achieving girls in mathematics and on low achieving boys in reading are likely to yield the strongest educational benefits

Intron Dynamics in Ribosomal Protein Genes

The role of spliceosomal introns in eukaryotic genomes remains obscure. A large scale analysis of intron presence/absence patterns in many gene families and species is a necessary step to clarify the role of these introns. In this analysis, we used a maximum likelihood method to reconstruct the evolution of 2,961 introns in a dataset of 76 ribosomal protein genes from 22 eukaryotes and validated the results by a maximum parsimony method. Our results show that the trends of intron gain and loss differed across species in a given kingdom but appeared to be consistent within subphyla. Most subphyla in the dataset diverged around 1 billion years ago, when the “Big Bang” radiation occurred. We speculate that spliceosomal introns may play a role in the explosion of many eukaryotes at the Big Bang radiation

Evolution of the Aging Brain Transcriptome and Synaptic Regulation

Alzheimer's disease and other neurodegenerative disorders of aging are characterized by clinical and pathological features that are relatively specific to humans. To obtain greater insight into how brain aging has evolved, we compared age-related gene expression changes in the cortex of humans, rhesus macaques, and mice on a genome-wide scale. A small subset of gene expression changes are conserved in all three species, including robust age-dependent upregulation of the neuroprotective gene apolipoprotein D (APOD) and downregulation of the synaptic cAMP signaling gene calcium/calmodulin-dependent protein kinase IV (CAMK4). However, analysis of gene ontology and cell type localization shows that humans and rhesus macaques have diverged from mice due to a dramatic increase in age-dependent repression of neuronal genes. Many of these age-regulated neuronal genes are associated with synaptic function. Notably, genes associated with GABA-ergic inhibitory function are robustly age-downregulated in humans but not in mice at the level of both mRNA and protein. Gene downregulation was not associated with overall neuronal or synaptic loss. Thus, repression of neuronal gene expression is a prominent and recently evolved feature of brain aging in humans and rhesus macaques that may alter neural networks and contribute to age-related cognitive changes

Insertion Sequence Inversions Mediated by Ectopic Recombination between Terminal Inverted Repeats

Transposable elements are widely distributed and diverse in both eukaryotes and prokaryotes, as exemplified by DNA transposons. As a result, they represent a considerable source of genomic variation, for example through ectopic (i.e. non-allelic homologous) recombination events between transposable element copies, resulting in genomic rearrangements. Ectopic recombination may also take place between homologous sequences located within transposable element sequences. DNA transposons are typically bounded by terminal inverted repeats (TIRs). Ectopic recombination between TIRs is expected to result in DNA transposon inversions. However, such inversions have barely been documented. In this study, we report natural inversions of the most common prokaryotic DNA transposons: insertion sequences (IS). We identified natural TIR-TIR recombination-mediated inversions in 9% of IS insertion loci investigated in Wolbachia bacteria, which suggests that recombination between IS TIRs may be a quite common, albeit largely overlooked, source of genomic diversity in bacteria. We suggest that inversions may impede IS survival and proliferation in the host genome by altering transpositional activity. They may also alter genomic instability by modulating the outcome of ectopic recombination events between IS copies in various orientations. This study represents the first report of TIR-TIR recombination within bacterial IS elements and it thereby uncovers a novel mechanism of structural variation for this class of prokaryotic transposable elements

Effects of L1-ORF2 fragments on green fluorescent protein gene expression

The retrotransposon known as long interspersed nuclear element-1 (L1) is 6 kb long, although most L1s in mammalian and other eukaryotic cells are truncated. L1 contains two open reading frames, ORF1 and ORF2, that code for an RNA-binding protein and a protein with endonuclease and reverse transcriptase activities, respectively. In this work, we examined the effects of full length L1-ORF2 and ORF2 fragments on green fluorescent protein gene (GFP) expression when inserted into the pEGFP-C1 vector downstream of GFP. All of the ORF2 fragments in sense orientation inhibited GFP expression more than when in antisense orientation, which suggests that small ORF2 fragments contribute to the distinct inhibitory effects of this ORF on gene expression. These results provide the first evidence that different 280-bp fragments have distinct effects on the termination of gene transcription, and that when inserted in the antisense direction, fragment 280-9 (the 3' end fragment of ORF2) induces premature termination of transcription that is consistent with the effect of ORF2

Combining Independent, Weighted P-Values: Achieving Computational Stability by a Systematic Expansion with Controllable Accuracy

Given the expanding availability of scientific data and tools to analyze them, combining different assessments of the same piece of information has become increasingly important for social, biological, and even physical sciences. This task demands, to begin with, a method-independent standard, such as the -value, that can be used to assess the reliability of a piece of information. Good's formula and Fisher's method combine independent -values with respectively unequal and equal weights. Both approaches may be regarded as limiting instances of a general case of combining -values from groups; -values within each group are weighted equally, while weight varies by group. When some of the weights become nearly degenerate, as cautioned by Good, numeric instability occurs in computation of the combined -values. We deal explicitly with this difficulty by deriving a controlled expansion, in powers of differences in inverse weights, that provides both accurate statistics and stable numerics. We illustrate the utility of this systematic approach with a few examples. In addition, we also provide here an alternative derivation for the probability distribution function of the general case and show how the analytic formula obtained reduces to both Good's and Fisher's methods as special cases. A C++ program, which computes the combined -values with equal numerical stability regardless of whether weights are (nearly) degenerate or not, is available for download at our group website http://www.ncbi.nlm.nih.gov/CBBresearch/Yu/downloads/CoinedPValues.html

Alu pair exclusions in the human genome

<p>Abstract</p> <p>Background</p> <p>The human genome contains approximately one million <it>Alu </it>elements which comprise more than 10% of human DNA by mass. <it>Alu </it>elements possess direction, and are distributed almost equally in positive and negative strand orientations throughout the genome. Previously, it has been shown that closely spaced <it>Alu </it>pairs in opposing orientation (inverted pairs) are found less frequently than <it>Alu </it>pairs having the same orientation (direct pairs). However, this imbalance has only been investigated for <it>Alu </it>pairs separated by 650 or fewer base pairs (bp) in a study conducted prior to the completion of the draft human genome sequence.</p> <p>Results</p> <p>We performed a comprehensive analysis of all (> 800,000) full-length <it>Alu </it>elements in the human genome. This large sample size permits detection of small differences in the ratio between inverted and direct <it>Alu </it>pairs (I:D). We have discovered a significant depression in the full-length <it>Alu </it>pair I:D ratio that extends to repeat pairs separated by ≤ 350,000 bp. Within this imbalance bubble (those <it>Alu </it>pairs separated by ≤ 350,000 bp), direct pairs outnumber inverted pairs. Using PCR, we experimentally verified several examples of inverted <it>Alu </it>pair exclusions that were caused by deletions.</p> <p>Conclusions</p> <p>Over 50 million full-length <it>Alu </it>pairs reside within the I:D imbalance bubble. Their collective impact may represent one source of <it>Alu </it>element-related human genomic instability that has not been previously characterized.</p

Regular Patterns for Proteome-Wide Distribution of Protein Abundance across Species

A proteome of the bio-entity, including cell, tissue, organ, and organism, consists of proteins of diverse abundance. The principle that determines the abundance of different proteins in a proteome is of fundamental significance for an understanding of the building blocks of the bio-entity. Here, we report three regular patterns in the proteome-wide distribution of protein abundance across species such as human, mouse, fly, worm, yeast, and bacteria: in most cases, protein abundance is positively correlated with the protein's origination time or sequence conservation during evolution; it is negatively correlated with the protein's domain number and positively correlated with domain coverage in protein structure, and the correlations became stronger during the course of evolution; protein abundance can be further stratified by the function of the protein, whereby proteins that act on material conversion and transportation (mass category) are more abundant than those that act on information modulation (information category). Thus, protein abundance is intrinsically related to the protein's inherent characters of evolution, structure, and function

Somatic Mutation Profiles of MSI and MSS Colorectal Cancer Identified by Whole Exome Next Generation Sequencing and Bioinformatics Analysis

BACKGROUND: Colorectal cancer (CRC) is with approximately 1 million cases the third most common cancer worldwide. Extensive research is ongoing to decipher the underlying genetic patterns with the hope to improve early cancer diagnosis and treatment. In this direction, the recent progress in next generation sequencing technologies has revolutionized the field of cancer genomics. However, one caveat of these studies remains the large amount of genetic variations identified and their interpretation. METHODOLOGY/PRINCIPAL FINDINGS: Here we present the first work on whole exome NGS of primary colon cancers. We performed 454 whole exome pyrosequencing of tumor as well as adjacent not affected normal colonic tissue from microsatellite stable (MSS) and microsatellite instable (MSI) colon cancer patients and identified more than 50,000 small nucleotide variations for each tissue. According to predictions based on MSS and MSI pathomechanisms we identified eight times more somatic non-synonymous variations in MSI cancers than in MSS and we were able to reproduce the result in four additional CRCs. Our bioinformatics filtering approach narrowed down the rate of most significant mutations to 359 for MSI and 45 for MSS CRCs with predicted altered protein functions. In both CRCs, MSI and MSS, we found somatic mutations in the intracellular kinase domain of bone morphogenetic protein receptor 1A, BMPR1A, a gene where so far germline mutations are associated with juvenile polyposis syndrome, and show that the mutations functionally impair the protein function. CONCLUSIONS/SIGNIFICANCE: We conclude that with deep sequencing of tumor exomes one may be able to predict the microsatellite status of CRC and in addition identify potentially clinically relevant mutations

Regional Selection Acting on the OFD1 Gene Family

The OFD1 (oral-facial-digital, type 1) gene is implicated in several developmental disorders in humans. The X-linked OFD1 (OFD1X) is conserved in Eutheria. Knowledge about the Y-linked paralog (OFD1Y) is limited. In this study, we identified an OFD1Y on the bovine Y chromosome, which is expressed differentially from the bovine OFD1X. Phylogenetic analysis indicated that: a) the eutherian OFD1X and OFD1Y were derived from the pair of ancestral autosomes during sex chromosome evolution; b) the autosomal OFD1 pseudogenes, present in Catarrhini and Murinae, were derived from retropositions of OFD1X after the divergence of primates and rodents; and c) the presence of OFD1Y in the ampliconic region of the primate Y chromosome is an indication that the expansion of the ampliconic region may initiate from the X-degenerated sequence. In addition, we found that different regions of OFD1/OFD1X/OFD1Y are under differential selection pressures. The C-terminal half of OFD1 is under relaxed selection with an elevated Ka/Ks ratio and clustered positively selected sites, whereas the N-terminal half is under stronger constraints. This study provides some insights into why the OFD1X gene causes OFD1 (male-lethal X-linked dominant) and SGBS2 & JSRDs (X-linked recessive) syndromes in humans, and reveals the origin and evolution of the OFD1 family, which will facilitate further clinical investigation of the OFD1-related syndromes