Systematic identification of common functional modules related to heart failure with different etiologies

The development of heart failure (HF) is a complex process that can be initiated by multiple etiologies. Identifying common functional modules associated with HF is a challenging task. Here, we developed a systems method to identify these common functional modules by integrating multiple expression profiles, protein interactions from four species, gene function annotations, and text information. We identified 1439 consistently differentially expressed genes (CDEGs) across HF with different etiologies by applying three meta-analysis methods to multiple HF-related expression profiles. Using a weighted human interaction network constructed by combining interaction data from multiple species, we extracted 60 candidate CDEG modules. We further evaluated the functional relevance of each module by using expression, interaction network, functional annotations, and text information together. Finally, five functional modules with significant biological relevance were identified. We found that almost half of the genes in these modules are hubs in the weighted network, and that these modules can accurately classify HF patients from healthy subjects. We also identified many significantly enriched biological processes that contribute to the pathophysiology of HF, including two new ones, RNA splicing and vesicle-mediated protein transport. In summary, we proposed a novel framework to analyze common functional modules related to HF with different etiologies. Our findings provide important insights into the complex mechanism of HF. Further biological experimentations should be required to validate these novel biological processes.     

Connexin 43 hemichannels contribute to cytoplasmic Ca2+ oscillations by providing a bimodal Ca2+-dependent Ca2+ entry pathway

Many cellular functions are driven by changes in the intracellular Ca(2+) concentration ([Ca(2+)](i)) that are highly organized in time and space. Ca(2+) oscillations are particularly important in this respect and are based on positive and negative [Ca(2+)](i) feedback on inositol 1,4,5-trisphosphate receptors (InsP(3)Rs). Connexin hemichannels are Ca(2+)-permeable plasma membrane channels that are also controlled by [Ca(2+)](i). We aimed to investigate how hemichannels may contribute to Ca(2+) oscillations. Madin-Darby canine kidney cells expressing connexin-32 (Cx32) and Cx43 were exposed to bradykinin (BK) or ATP to induce Ca(2+) oscillations. BK-induced oscillations were rapidly (minutes) and reversibly inhibited by the connexin-mimetic peptides (32)Gap27/(43)Gap26, whereas ATP-induced oscillations were unaffected. Furthermore, these peptides inhibited the BK-triggered release of calcein, a hemichannel-permeable dye. BK-induced oscillations, but not those induced by ATP, were dependent on extracellular Ca(2+). Alleviating the negative feedback of [Ca(2+)](i) on InsP(3)Rs using cytochrome c inhibited BK- and ATP-induced oscillations. Cx32 and Cx43 hemichannels are activated by <500 nm [Ca(2+)](i) but inhibited by higher concentrations and CT9 peptide (last 9 amino acids of the Cx43 C terminus) removes this high [Ca(2+)](i) inhibition. Unlike interfering with the bell-shaped dependence of InsP(3)Rs to [Ca(2+)](i), CT9 peptide prevented BK-induced oscillations but not those triggered by ATP. Collectively, these data indicate that connexin hemichannels contribute to BK-induced oscillations by allowing Ca(2+) entry during the rising phase of the Ca(2+) spikes and by providing an OFF mechanism during the falling phase of the spikes. Hemichannels were not sufficient to ignite oscillations by themselves; however, their contribution was crucial as hemichannel inhibition stopped the oscillations.     

Meta-analysis suggests choosy females get sexy sons more than "good genes"

Female preferences for specific male phenotypes have been documented across a wide range of animal taxa, including numerous species where males contribute only gametes to offspring production. Yet, selective pressures maintaining such preferences are among the major unknowns of evolutionary biology. Theoretical studies suggest that preferences can evolve if they confer genetic benefits in terms of increased attractiveness of sons ("Fisherian" models) or overall fitness of offspring ("good genes" models). These two types of models predict, respectively, that male attractiveness is heritable and genetically correlated with fitness. In this meta-analysis, we draw general conclusions from over two decades worth of empirical studies testing these predictions (90 studies on 55 species in total). We found evidence for heritability of male attractiveness. However, attractiveness showed no association with traits directly associated with fitness (life-history traits). Interestingly, it did show a positive correlation with physiological traits, which include immunocompetence and condition. In conclusion, our results support "Fisherian" models of preference evolution, while providing equivocal evidence for "good genes." We pinpoint research directions that should stimulate progress in our understanding of the evolution of female choice.     

High-Throughput Discovery of Mutations in Tef Semi-Dwarfing Genes by Next-Generation Sequencing Analysis

Tef (Eragrostis tef) is a major cereal crop in Ethiopia. Lodging is the primary constraint to increasing productivity in this allotetraploid species, accounting for losses of ∼15–45% in yield each year. As a first step toward identifying semi-dwarf varieties that might have improved lodging resistance, an ∼6× fosmid library was constructed and used to identify both homeologues of the dw3 semi-dwarfing gene of Sorghum bicolor. An EMS mutagenized population, consisting of ∼21,210 tef plants, was planted and leaf materials were collected into 23 superpools. Two dwarfing candidate genes, homeologues of dw3 of sorghum and rht1 of wheat, were sequenced directly from each superpool with 454 technology, and 120 candidate mutations were identified. Out of 10 candidates tested, six independent mutations were validated by Sanger sequencing, including two predicted detrimental mutations in both dw3 homeologues with a potential to improve lodging resistance in tef through further breeding. This study demonstrates that high-throughput sequencing can identify potentially valuable mutations in under-studied plant species like tef and has provided mutant lines that can now be combined and tested in breeding programs for improved lodging resistance.     

Disentangling the roles of history and local selection in shaping clinal variation of allele frequencies and gene expression in Norway spruce (Picea abies)

Understanding the genetic basis of local adaptation is challenging due to the subtle balance among conflicting evolutionary forces that are involved in its establishment and maintenance. One system with which to tease apart these difficulties is clines in adaptive characters. Here we analyzed genetic and phenotypic variation in bud set, a highly heritable and adaptive trait, among 18 populations of Norway spruce (Picea abies), arrayed along a latitudinal gradient ranging from 47°N to 68°N. We confirmed that variation in bud set is strongly clinal, using a subset of five populations. Genotypes for 137 single-nucleotide polymorphisms (SNPs) chosen from 18 candidate genes putatively affecting bud set and 308 control SNPs chosen from 264 random genes were analyzed for patterns of genetic structure and correlation to environment. Population genetic structure was low (F(ST) = 0.05), but latitudinal patterns were apparent among Scandinavian populations. Hence, part of the observed clinal variation should be attributable to population demography. Conditional on patterns of genetic structure, there was enrichment of SNPs within candidate genes for correlations with latitude. Twenty-nine SNPs were also outliers with respect to F(ST). The enrichment for clinal variation at SNPs within candidate genes (i.e., SNPs in PaGI, PaPhyP, PaPhyN, PaPRR7, and PaFTL2) indicated that local selection in the 18 populations, and/or selection in the ancestral populations from which they were recently derived, shaped the observed cline. Validation of these genes using expression studies also revealed that PaFTL2 expression is significantly associated with latitude, thereby confirming the central role played by this gene in the control of phenology in plants.     

Identification of differentially expressed microRNAs and their PKC-isoform specific gene network prediction during hypoxic pre-conditioning and focal cerebral ischemia of mice

We previously reported the involvement of conventional protein kinase C (cPKC) βII, γ, novel PKC (nPKC) ε and their interacting proteins in hypoxic pre-conditioning (HPC)-induced neuroprotection. In this study, the large-scale miRNA microarrays and bioinformatics analysis were used to determine the differentially expressed miRNAs and their PKC-isoform specific gene network in mouse brain after HPC and 6?h middle cerebral artery occlusion (MCAO). We found 4 up-regulated and 13 down-regulated miRNAs in the cortex of HPC mice, 26 increased and 39 decreased gene expressions of miRNAs in the peri-infarct region of 6?h MCAO mice, and 11 up-regulated and 22 down-regulated miRNAs in the peri-infarct region of HPC and 6?h MCAO mice. Based on Diff Score, 19 differentially expressed miRNAs were identified in HPC and 6?h MCAO mouse brain. Then the miRNA-gene-network of 19 specified miRNAs target genes of cPKCβII, γ and nPKCε-interacting protein was predicted by using bioinformatics analysis of genome databases. Furthermore, the down-regulated miR-615-3p during HPC had a detrimental effect on the oxygen-glucose deprivation (OGD)-induced N2A cell injury. These results suggested that the identified 19 miRNAs, notably miR-615-3p, might target these genes of cPKCβII, γ and nPKCε-interacting proteins involved in HPC-induced neuroprotection.     

Complete mitochondrial genome sequences of three bats species and whole genome mitochondrial analyses reveal patterns of codon bias and lend support to a basal split in Chiroptera

Order Chiroptera is a unique group of mammals whose members have attained self-powered flight as their main mode of locomotion. Much speculation persists regarding bat evolution; however, lack of sufficient molecular data hampers evolutionary and conservation studies. Of ~ 1200 species, complete mitochondrial genome sequences are available for only eleven. Additional sequences should be generated if we are to resolve many questions concerning these fascinating mammals. Herein, we describe the complete mitochondrial genomes of three bats: Corynorhinus rafinesquii, Lasiurus borealis and Artibeus lituratus. We also compare the currently available mitochondrial genomes and analyze codon usage in Chiroptera. C. rafinesquii, L. borealis and A. lituratus mitochondrial genomes are 16438 bp, 17048 bp and 16709 bp, respectively. Genome organization and gene arrangements are similar to other bats. Phylogenetic analyses using complete mitochondrial genome sequences support previously established phylogenetic relationships and suggest utility in future studies focusing on the evolutionary aspects of these species. Comprehensive analyses of available bat mitochondrial genomes reveal distinct nucleotide patterns and synonymous codon preferences corresponding to different chiropteran families. These patterns suggest that mutational and selection forces are acting to different extents within Chiroptera and shape their mitochondrial genomes.     

Systematic investigation of Amphioxus (Branchiostoma floridae) microRNAs

Ascaris lumbricoides and Ascaris suum are parasitic nematodes living in the small intestine of humans and pigs, and can cause the disease ascariasis. For long, there has been controversy as to whether the two ascaridoid taxa represent the same species due to their significant resemblances in morphology. However, the complete mitochondrial (mt) genome data have been lacking for A. lumbricoides in spite of human and animal health significance and socio-economic impact globally of these parasites. In the present study, we sequenced the complete mt genomes of A. lumbricoides and A. suum (China isolate), which was 14,303 bp and 14,311 bp in size, respectively. The identity of the mt genomes was 98.1% between A. lumbricoides and A. suum (China isolate), and 98.5% between A. suum (China isolate) and A. suum (USA isolate). Both genomes are circular, and consist of 36 genes, including 12 genes for proteins, 2 genes for rRNA and 22 genes for tRNA, which are consistent with that of all other species of ascaridoid studied to date. All genes are transcribed in the same direction and have a nucleotide composition high in A and T (71.7% for A. lumbricoides and 71.8% for A. suum). The AT bias had a significant effect on both the codon usage pattern and amino acid composition of proteins. Phylogenetic analyses of A. lumbricoides and A. suum using concatenated amino acid sequences of 12 protein-coding genes, with three different computational algorithms (Bayesian analysis, maximum likelihood and maximum parsimony) all clustered in a clade with high statistical support, indicating that A. lumbricoides and A. suum was very closely related. These mt genome data and the results provide some additional genetic evidence that A. lumbricoides and A. suum may represent the same species. The mt genome data presented in this study are also useful novel markers for studying the molecular epidemiology and population genetics of Ascaris.     

Beyond barcoding: A mitochondrial genomics approach to molecular phylogenetics and diagnostics of blowflies (Diptera: Calliphoridae)

Members of the Calliphoridae (blowflies) are significant for medical and veterinary management, due to the ability of some species to consume living flesh as larvae, and for forensic investigations due to the ability of others to develop in corpses. Due to the difficulty of accurately identifying larval blowflies to species there is a need for DNA-based diagnostics for this family, however the widely used DNA-barcoding marker, cox1, has been shown to fail for several groups within this family. Additionally, many phylogenetic relationships within the Calliphoridae are still unresolved, particularly deeper level relationships. Sequencing whole mt genomes has been demonstrated both as an effective method for identifying the most informative diagnostic markers and for resolving phylogenetic relationships. Twenty-seven complete, or nearly so, mt genomes were sequenced representing 13 species, seven genera and four calliphorid subfamilies and a member of the related family Tachinidae. PCR and sequencing primers developed for sequencing one calliphorid species could be reused to sequence related species within the same superfamily with success rates ranging from 61% to 100%, demonstrating the speed and efficiency with which an mt genome dataset can be assembled. Comparison of molecular divergences for each of the 13 protein-coding genes and 2 ribosomal RNA genes, at a range of taxonomic scales identified novel targets for developing as diagnostic markers which were 117–200% more variable than the markers which have been used previously in calliphorids. Phylogenetic analysis of whole mt genome sequences resulted in much stronger support for family and subfamily-level relationships. The Calliphoridae are polyphyletic, with the Polleninae more closely related to the Tachinidae, and the Sarcophagidae are the sister group of the remaining calliphorids. Within the Calliphoridae, there was strong support for the monophyly of the Chrysomyinae and Luciliinae and for the sister-grouping of Luciliinae with Calliphorinae. Relationships within Chrysomya were not well resolved. Whole mt genome data, supported the previously demonstrated paraphyly of Lucilia cuprina with respect to L. sericata and allowed us to conclude that it is due to hybrid introgression prior to the last common ancestor of modern sericata populations, rather than due to recent hybridisation, nuclear pseudogenes or incomplete lineage sorting.