The Drosophila ribosome protein S5 paralog RpS5b promotes germ cell and follicle cell differentiation during oogenesis.

Emerging evidence suggests that ribosome heterogeneity may have important functional consequences in the translation of specific mRNAs within different cell types and under various conditions. Ribosome heterogeneity comes in many forms, including post-translational modification of ribosome proteins (RPs), absence of specific RPs and inclusion of different RP paralogs. The Drosophila genome encodes two RpS5 paralogs: RpS5a and RpS5b. While RpS5a is ubiquitously expressed, RpS5b exhibits enriched expression in the reproductive system. Deletion of RpS5b results in female sterility marked by developmental arrest of egg chambers at stages 7-8, disruption of vitellogenesis and posterior follicle cell (PFC) hyperplasia. While transgenic rescue experiments suggest functional redundancy between RpS5a and RpS5b, molecular, biochemical and ribo-seq experiments indicate that RpS5b mutants display increased rRNA transcription and RP production, accompanied by increased protein synthesis. Loss of RpS5b results in microtubule-based defects and in mislocalization of Delta and Mindbomb1, leading to failure of Notch pathway activation in PFCs. Together, our results indicate that germ cell-specific expression of RpS5b promotes proper egg chamber development by ensuring the homeostasis of functional ribosomes.

The CaspBase: a curated database for evolutionary biochemical studies of caspase functional divergence and ancestral sequence inference.

Sequence databases are powerful tools for the contemporary scientists' toolkit. However, most functional annotations in public databases are determined computationally and are not verified by a human expert. While hypotheses generated from computational studies are now amenable to experimentation, the quality of the results relies on the quality of input data. We developed the CaspBase to expedite high-quality dataset compilation of annotated caspase sequences, to maximize phylogenetic signal, and to reduce the noise contributed from public databanks. We describe our methods of curation for the CaspBase and how researchers can acquire sequences from CaspBase.org. Our immediate goal for developing the CaspBase was to optimize the ancestral protein reconstruction (APR) of caspases, and we demonstrate the utility of the CaspBase in APR studies. We also developed the Common Position (CP) system for comparing human caspase family paralogs and suggest the CP system as an update to current reporting methods of caspase amino acid positions. We present a standardized multiple sequence alignment (MSA) for the CP system and show the advantage of using large databases such as the CaspBase in defining structural positions in proteins. Although the results described here pertain to caspase evolution and structure-function studies, the methods can be adapted to any gene family.

Few Nuclear-Encoded Mitochondrial Gene Duplicates Contribute to Male Germline-Specific Functions in Humans.

Most of the genes encoding proteins that function in the mitochondria are located in the nucleus and are called nuclear-encoded mitochondrial genes, or N-mt genes. In Drosophila melanogaster , about 23% of N-mt genes fall into gene families, and all duplicates with tissue-biased expression (76%) are testis biased. These genes are enriched for energy-related functions and tend to be older than other duplicated genes in the genome. These patterns reveal strong selection for the retention of new genes for male germline mitochondrial functions. The two main forces that are likely to drive changes in mitochondrial functions are maternal inheritance of mitochondria and male-male competition for fertilization. Both are common among animals, suggesting similar N-mt gene duplication patterns in different species. To test this, we analyzed N-mt genes in the human genome. We find that about 18% of human N-mt genes fall into gene families, but unlike in Drosophila , only 28% of the N-mt duplicates have tissue-biased expression and only 36% of these have testis-biased expression. In addition, human testis-biased duplicated genes are younger than other duplicated genes in the genome and have diverse functions. These contrasting patterns between species might reflect either differences in selective pressures for germline energy-related or other mitochondrial functions during spermatogenesis and fertilization, or differences in the response to similar pressures.

A Highly Predictive Model for Diagnosis of Colorectal Neoplasms Using Plasma MicroRNA: Improving Specificity and Sensitivity.

OBJECTIVE: To develop a plasma-based microRNA (miRNA) diagnostic assay specific for colorectal neoplasms, building upon our prior work. BACKGROUND: Colorectal neoplasms [colorectal cancer (CRC) and colorectal advanced adenoma (CAA)] frequently develop in individuals at ages when other common cancers also occur. Current screening methods lack sensitivity, specificity, and have poor patient compliance. METHODS: Plasma was screened for 380 miRNAs using microfluidic array technology from a "Training" cohort of 60 patients, (10 each) control, CRC, CAA, breast cancer, pancreatic cancer, and lung cancer. We identified uniquely dysregulated miRNAs specific for colorectal neoplasia (P < 0.05, false discovery rate: 5%, adjusted α = 0.0038). These miRNAs were evaluated using single assays in a "Test" cohort of 120 patients. A mathematical model was developed to predict blinded sample identity in a 150 patient "Validation" cohort using repeat-sub-sampling validation of the testing dataset with 1000 iterations each to assess model detection accuracy. RESULTS: Seven miRNAs (miR-21, miR-29c, miR-122, miR-192, miR-346, miR-372, and miR-374a) were selected based upon P value, area under the curve (AUC), fold change, and biological plausibility. Area under the curve (±95% confidence interval) for "Test" cohort comparisons were 0.91 (0.85-0.96) between all neoplasia and controls, 0.79 (0.70-0.88) between colorectal neoplasia and other cancers, and 0.98 (0.96-1.0) between CRC and colorectal adenomas. In our "Validation" cohort, our mathematical model predicted blinded sample identity with 69% to 77% accuracy, 67% to 76% accuracy, and 86% to 90% accuracy for each comparison, respectively. CONCLUSIONS: Our plasma miRNA assay and prediction model differentiate colorectal neoplasia from patients with other neoplasms and from controls with higher sensitivity and specificity compared with current clinical standards.

Genomics of ecological adaptation in cactophilic Drosophila.

Cactophilic Drosophila species provide a valuable model to study gene-environment interactions and ecological adaptation. Drosophila buzzatii and Drosophila mojavensis are two cactophilic species that belong to the repleta group, but have very different geographical distributions and primary host plants. To investigate the genomic basis of ecological adaptation, we sequenced the genome and developmental transcriptome of D. buzzatii and compared its gene content with that of D. mojavensis and two other noncactophilic Drosophila species in the same subgenus. The newly sequenced D. buzzatii genome (161.5 Mb) comprises 826 scaffolds (>3 kb) and contains 13,657 annotated protein-coding genes. Using RNA sequencing data of five life-stages we found expression of 15,026 genes, 80% protein-coding genes, and 20% noncoding RNA genes. In total, we detected 1,294 genes putatively under positive selection. Interestingly, among genes under positive selection in the D. mojavensis lineage, there is an excess of genes involved in metabolism of heterocyclic compounds that are abundant in Stenocereus cacti and toxic to nonresident Drosophila species. We found 117 orphan genes in the shared D. buzzatii-D. mojavensis lineage. In addition, gene duplication analysis identified lineage-specific expanded families with functional annotations associated with proteolysis, zinc ion binding, chitin binding, sensory perception, ethanol tolerance, immunity, physiology, and reproduction. In summary, we identified genetic signatures of adaptation in the shared D. buzzatii-D. mojavensis lineage, and in the two separate D. buzzatii and D. mojavensis lineages. Many of the novel lineage-specific genomic features are promising candidates for explaining the adaptation of these species to their distinct ecological niches.

A plasma microRNA panel for detection of colorectal adenomas: a step toward more precise screening for colorectal cancer.

OBJECTIVE: The main objective of this study was to investigate the potential use of circulating microRNAs (miRNAs) as biomarkers of colorectal (CR) adenomas. BACKGROUND: Detection of precancerous lesions such as CR adenoma is a key to reduce CR cancer (CRC) mortality. There is a great need for accurate, noninvasive biomarkers for detection of CR adenoma and CRC. MiRNAs are non-protein-coding RNAs that regulate gene expression. Our prior work investigated the dysregulation of 5 plasma miRNAs in CRC patients. As intended, we undertook a more comprehensive plasma-miRNA screening study in patients with CR adenoma and CRC. METHODS: We screened for 380 plasma-miRNAs using microfluidic array technology (Applied BioSystems) in a screening cohort of 12 healthy controls, 9 patients with CR adenomas, and 20 patients with CRC. A panel of the most dysregulated miRNAs (P < 0.05, False Discovery Rate: 5%) was then validated in a blinded cohort of 26 healthy controls, 16 patients with large adenomas, and 45 patients with CRC. RESULTS: A panel of 8 plasma miRNAs (miR-532-3p, miR-331, miR-195, miR-17, miR-142-3p, miR-15b, miR-532, and miR-652) distinguished polyps from controls with high accuracy [area under curve (AUC) = 0.868 (95% confidence interval [CI]: 0.76-0.98)]. In addition, a panel of 3 plasma miRNAs (miR-431, miR-15b, and miR-139-3p) distinguished Stage IV CRC from controls with an [AUC = 0.896 (95% CI: 0.78-1.0)]. Receiver-operating-characteristic curves of miRNA panels for all CRC versus controls and polyps versus all CRC showed AUC values of 0.829 (95% CI: 0.73-0.93) and 0.856 (95% CI: 0.75-0.97), respectively. CONCLUSIONS: Plasma miRNAs are reliable, noninvasive, and inexpensive markers for CR adenomas. This miRNA panel warrants study in larger cohorts. Plasma-based assays could provide better screening compliance compared to fecal occult blood or endoscopic screening.

Why chromosome palindromes?

We look at sex-limited chromosome (Y or W) evolution with particular emphasis on the importance of palindromes. Y chromosome palindromes consist of inverted duplicates that allow for local recombination in an otherwise nonrecombining chromosome. Since palindromes enable intrachromosomal gene conversion that can help eliminate deleterious mutations, they are often highlighted as mechanisms to protect against Y degeneration. However, the adaptive significance of recombination resides in its ability to decouple the evolutionary fates of linked mutations, leading to both a decrease in degeneration rate and an increase in adaptation rate. Our paper emphasizes the latter, that palindromes may exist to accelerate adaptation by increasing the potential targets and fixation rates of incoming beneficial mutations. This hypothesis helps reconcile two enigmatic features of the "palindromes as protectors" view: (1) genes that are not located in palindromes have been retained under purifying selection for tens of millions of years, and (2) under models that only consider deleterious mutations, gene conversion benefits duplicate gene maintenance but not initial fixation. We conclude by looking at ways to test the hypothesis that palindromes enhance the rate of adaptive evolution of Y-linked genes and whether this effect can be extended to palindromes on other chromosomes.

Gene expression levels are correlated with synonymous codon usage, amino acid composition, and gene architecture in the red flour beetle, Tribolium castaneum.

Gene expression levels correlate with multiple aspects of gene sequence and gene structure in phylogenetically diverse taxa, suggesting an important role of gene expression levels in the evolution of protein-coding genes. Here we present results of a genome-wide study of the influence of gene expression on synonymous codon usage, amino acid composition, and gene structure in the red flour beetle, Tribolium castaneum. Consistent with the action of translational selection, we find that synonymous codon usage bias increases with gene expression. However, the correspondence between tRNA gene copy number and optimal codons is weak. At the amino acid level, translational selection is suggested by the positive correlation between tRNA gene numbers and amino acid usage, which is stronger for highly expressed genes. In addition, there is a clear trend for increased use of metabolically cheaper, less complex amino acids as gene expression increases. tRNA gene numbers also correlate negatively with amino acid size/complexity (S/C) score indicating the coupling between translational selection and selection to minimize the use of large/complex amino acids. Interestingly, the analysis of 10 additional genomes suggests that the correlation between tRNA gene numbers and amino acid S/C score is widespread and might be explained by selection against negative consequences of protein misfolding. At the level of gene structure, three major trends are detected: 1) complete coding region length increases across low and intermediate expression levels but decreases in highly expressed genes; 2) the average intron size shows the opposite trend, first decreasing with expression, followed by a slight increase in highly expressed genes; and 3) intron density remains nearly constant across all expression levels. These changes in gene architecture are only in partial agreement with selection favoring reduced cost of biosynthesis.

Local effects of limited recombination: historical perspective and consequences for population estimates of adaptive evolution.

Recent years have witnessed the integration of theoretical advances in population genetics with large-scale analyses of complete genomes, with a growing number of studies suggesting pervasive natural selection that includes frequent deleterious as well as adaptive mutations. In finite populations, however, mutations under selection alter the fate of genetically linked mutations (the so-called Hill-Robertson effect). Here we review the evolutionary consequences of selection at linked sites (linked selection) focusing on its effects on nearby nucleotides in genomic regions with nonreduced recombination. We argue that these local effects of linkage may account for differences in selection intensity among genes. We also show that even high levels of recombination are unlikely to remove all effects of linked selection, causing a reduction in the polymorphism to divergence ratio (r(pd)) at neutral sites. Because a number of methods employed to estimate the magnitude and frequency of adaptive mutations take reduced r(pd) as evidence of positive selection, ignoring local linkage effects may lead to misleading estimates of the proportion of adaptive substitutions and estimates of positive selection. These biases are caused by employing methods that do not account for local variation in the relative effective population size (N(e)) caused by linked selection.

Evolution of a novel function: nutritive milk in the viviparous cockroach, Diploptera punctata.

Cockroach species show different degrees of maternal contribution to the developing offspring. In this study, we identify a multigene family that encodes water-soluble proteins that are a major component of nutritive "Milk" in the cockroach, Diploptera punctata. This gene family is associated with the evolution of a new trait, viviparity, in which the offspring receive nutrition during the gestation period. Twenty-five distinct Milk complementary DNAs were cloned and partially characterized. These complementary DNAs encode 22 distinct Milk peptides, each of length 171 amino acids, including a 16-amino acid signal peptide sequence. Southern blot analysis confirms the presence of multiple copies of Milk genes in D. punctata. Northern analysis indicates tissue- and stage-specific Milk gene expression. Examination of the deduced amino acid sequences identifies the presence of structurally conserved regions diagnostic of the lipocalin protein family. The shared exon/intron structure of one of the Milk loci with lipocalin genes further supports a close evolutionary relationship between these sequences.