The Nature of Genetic Variation for Complex Traits Revealed by GWAS and Regional Heritability Mapping Analyses

We use computer simulations to investigate the amount of genetic variation for complex traits that can be revealed by single-SNP genome-wide association studies (GWAS) or regional heritability mapping (RHM) analyses based on full genome sequence data or SNP chips. We model a large population subject to mutation, recombination, selection, and drift, assuming a pleiotropic model of mutations sampled from a bivariate distribution of effects of mutations on a quantitative trait and fitness. The pleiotropic model investigated, in contrast to previous models, implies that common mutations of large effect are responsible for most of the genetic variation for quantitative traits, except when the trait is fitness itself. We show that GWAS applied to the full sequence increases the number of QTL detected by as much as 50% compared to the number found with SNP chips but only modestly increases the amount of additive genetic variance explained. Even with full sequence data, the total amount of additive variance explained is generally below 50%. Using RHM on the full sequence data, a slightly larger number of QTL are detected than by GWAS if the same probability threshold is assumed, but these QTL explain a slightly smaller amount of genetic variance. Our results also suggest that most of the missing heritability is due to the inability to detect variants of moderate effect (∼0.03–0.3 phenotypic SDs) segregating at substantial frequencies. Very rare variants, which are more difficult to detect by GWAS, are expected to contribute little genetic variation, so their eventual detection is less relevant for resolving the missing heritability problem.     

The Evolution and Origin of Animal Toll-Like Receptor Signaling Pathway Revealed by Network-Level Molecular Evolutionary Analyses

Genes carry out their biological functions through pathways in complex networks consisting of many interacting molecules. Studies on the effect of network architecture on the evolution of individual proteins will provide valuable information for understanding the origin and evolution as well as functional conservation of signaling pathways. However, the relationship between the network architecture and the individual protein sequence evolution is yet little known. In current study, we carried out network-level molecular evolution analysis on TLR (Toll-like receptor ) signaling pathway, which plays an important role in innate immunity in insects and mammals, and we found that: 1) The selection constraint of genes was negatively correlated with its position along TLR signaling pathway; 2) all genes in TLR signaling pathway were highly conserved and underwent strong purifying selection; 3) the distribution of selective pressure along the pathway was driven by differential nonsynonymous substitution levels; 4) The TLR signaling pathway might present in a common ancestor of sponges and eumetazoa, and evolve via the TLR, IKK, IκB and NF-κB genes underwent duplication events as well as adaptor molecular enlargement, and gene structure and conservation motif of NF-κB genes shifted in their evolutionary history. Our results will improve our understanding on the evolutionary history of animal TLR signaling pathway as well as the relationship between the network architecture and the sequences evolution of individual protein.     

Amyloid Evolution: Antiparallel Replaced by Parallel

Several atomic structures have now been found for micrometer-scale amyloid fibrils or elongated microcrystals using a range of methods, including NMR, electron microscopy, and X-ray crystallography, with parallel β-sheet appearing as the most common secondary structure. The etiology of amyloid disease, however, indicates nanometer-scale assemblies of only tens of peptides as significant agents of cytotoxicity and contagion. By combining solution X-ray with molecular dynamics, we show that antiparallel structure dominates at the first stages of aggregation for a specific set of peptides, being replaced by parallel at large length scales only. This divergence in structure between small and large amyloid aggregates should inform future design of molecular therapeutics against nucleation or intercellular transmission of amyloid. Calculations and an overview from the literature argue that antiparallel order should be the first appearance of structure in many or most amyloid aggregation processes, regardless of the endpoint. Exceptions to this finding should exist, depending inevitably on the sequence and on solution conditions.     

Molecular Evolution and Diversity of Conus Peptide Toxins,as Revealed by Gene Structure and Intron Sequence Analyses

Cone snails, which are predatory marine gastropods, produce a cocktail of venoms used for predation, defense and competition. The major venom component, conotoxin, has received significant attention because it is useful in neuroscience research, drug development and molecular diversity studies. In this study, we report the genomic characterization of nine conotoxin gene superfamilies from 18 Conus species and investigate the relationships among conotoxin gene structure, molecular evolution and diversity. The I1, I2, M, O2, O3, P, S, and T superfamily precursors all contain three exons and two introns, while A superfamily members contain two exons and one intron. The introns are conserved within a certain gene superfamily, and also conserved across different Conus species, but divergent among different superfamilies. The intronic sequences contain many simple repeat sequences and regulatory elements that may influence conotoxin gene expression. Furthermore, due to the unique gene structure of conotoxins, the base substitution rates and the number of positively selected sites vary greatly among exons. Many more point mutations and trinucleotide indels were observed in the mature peptide exon than in the other exons. In addition, the first example of alternative splicing in conotoxin genes was found. These results suggest that the diversity of conotoxin genes has been shaped by point mutations and indels, as well as rare gene recombination or alternative splicing events, and that the unique gene structures could have made a contribution to the evolution of conotoxin genes.     

Surfing along Filopodia: A Particle Transport Revealed by Molecular-Scale Fluctuation Analyses

Filopodia perform cellular functions such as environmental sensing or cell motility, but they also grab for particles and withdraw them leading to an increased efficiency of phagocytic uptake. Remarkably, withdrawal of micron-sized particles is also possible without noticeable movements of the filopodia. Here, we demonstrate that polystyrene beads connected by optical tweezers to the ends of adherent filopodia of J774 macrophages, are transported discontinuously toward the cell body. After a typical resting time of 1–2 min, the cargo is moved with alternating velocities, force constants, and friction constants along the surface of the filopodia. This surfing-like behavior along the filopodium is recorded by feedback-controlled interferometric three-dimensional tracking of the bead motions at 10–100 kHz. We measured transport velocities of up to 120 nm/s and transport forces of ∼70 pN. Small changes in position, fluctuation width, and temporal correlation, which are invisible in conventional microscopy, indicate molecular reorganization of transport-relevant proteins in different phases of the entire transport process. A detailed analysis implicates a controlled particle transport with fingerprints of a nanoscale unbinding/binding behavior. The manipulation and analysis methods presented in our study may also be helpful in other fields of cellular biophysics.     

Unique Features of Odorant-Binding Proteins of the Parasitoid Wasp Nasonia vitripennis Revealed by Genome Annotation and Comparative Analyses

Insects are the most diverse group of animals on the planet, comprising over 90% of all metazoan life forms, and have adapted to a wide diversity of ecosystems in nearly all environments. They have evolved highly sensitive chemical senses that are central to their interaction with their environment and to communication between individuals. Understanding the molecular bases of insect olfaction is therefore of great importance from both a basic and applied perspective. Odorant binding proteins (OBPs) are some of most abundant proteins found in insect olfactory organs, where they are the first component of the olfactory transduction cascade, carrying odorant molecules to the olfactory receptors. We carried out a search for OBPs in the genome of the parasitoid wasp Nasonia vitripennis and identified 90 sequences encoding putative OBPs. This is the largest OBP family so far reported in insects. We report unique features of the N. vitripennis OBPs, including the presence and evolutionary origin of a new subfamily of double-domain OBPs (consisting of two concatenated OBP domains), the loss of conserved cysteine residues and the expression of pseudogenes. This study also demonstrates the extremely dynamic evolution of the insect OBP family: (i) the number of different OBPs can vary greatly between species; (ii) the sequences are highly diverse, sometimes as a result of positive selection pressure with even the canonical cysteines being lost; (iii) new lineage specific domain arrangements can arise, such as the double domain OBP subfamily of wasps and mosquitoes.     

Comparative Analyses of Four Complete Genomes in Pseudomonas amygdali Revealed Differential Adaptation to Hostile Environments and Secretion Systems

Pseudomonas amygdali is a hemibiotrophic phytopathogen that causes disease in woody and herbaceous plants. Complete genomes of four P. amygdali pathovars were comparatively analyzed to decipher the impact of genomic diversity on host colonization. The pan-genome indicated that 3,928 core genes are conserved among pathovars, while 504–1,009 are unique to specific pathovars. The unique genome contained many mobile elements and exhibited a functional distribution different from the core genome. Genes involved in O-antigen biosynthesis and antimicrobial peptide resistance were significantly enriched for adaptation to hostile environments. While the type III secretion system was distributed in the core genome, unique genomes revealed a different organization of secretion systems as follows: type I in pv. tabaci, type II in pv. japonicus, type IV in pv. morsprunorum, and type VI in pv. lachrymans. These findings provide genetic insight into the dynamic interactions of the bacteria with plant hosts.     

Dynamic Evolution of Oryza Genomes Is Revealed by Comparative Genomic Analysis of a Genus-Wide Vertical Data Set

Oryza (23 species; 10 genome types) contains the world's most important food crop — rice. Although the rice genome serves as an essential tool for biological research, little is known about the evolution of the other Oryza genome types. They contain a historical record of genomic changes that led to diversification of this genus around the world as well as an untapped reservoir of agriculturally important traits. To investigate the evolution of the collective Oryza genome, we sequenced and compared nine orthologous genomic regions encompassing the Adh1-Adh2 genes (from six diploid genome types) with the rice reference sequence. Our analysis revealed the architectural complexities and dynamic evolution of this region that have occurred over the past ~15 million years. Of the 46 intact genes and four pseudogenes in the japonica genome, 38 (76%) fell into eight multigene families. Analysis of the evolutionary history of each family revealed independent and lineage-specific gain and loss of gene family members as frequent causes of synteny disruption. Transposable elements were shown to mediate massive replacement of intergenic space (>95%), gene disruption, and gene/gene fragment movement. Three cases of long-range structural variation (inversions/deletions) spanning several hundred kilobases were identified that contributed significantly to genome diversification.     

Secondary Binding Site of Trypsin: Revealed by Crystal Structure of Trypsin-Peptide Complex

Abstract

Designed synthetic heterochiral peptides, when added to porcine trypsin, resulted in reduction of enzyme activity. The crystal structure of a complex formed between porcine trypsin and a heterochiral hepta peptide Boc-Pro-DAsp-Aib-Leu-Aib-Leu-Ala-NHMe has been determined at 1.9 Å resolution. The hepta peptide does not bind at the active site, but is located in the interstitial region, and interacts with the calcium-binding loop (residues 60–80). The bound peptide interacts with the active site residue Ser195 through an acetate ion, and with Lys 60 mediated by water molecules. The structure, when compared with the other trypsin-peptide complexes, suggests that the flexibility of surface loops, concerted movement of the loops towards the active site, and the interaction of the bound peptide with Lys 60, may be responsible for the reduction in enzyme activity. This study provides a structural evidence for the earlier biochemical observation regarding the role of surface loops in the catalysis of the enzyme.     

Consequences of Whole-Genome Triplication as Revealed by Comparative Genomic Analyses of the Wild Radish Raphanus raphanistrum and Three Other Brassicaceae Species

Polyploidization events are frequent among flowering plants, and the duplicate genes produced via such events contribute significantly to plant evolution. We sequenced the genome of wild radish (Raphanus raphanistrum), a Brassicaceae species that experienced a whole-genome triplication event prior to diverging from Brassica rapa. Despite substantial gene gains in these two species compared with Arabidopsis thaliana and Arabidopsis lyrata, ∼70% of the orthologous groups experienced gene losses in R. raphanistrum and B. rapa, with most of the losses occurring prior to their divergence. The retained duplicates show substantial divergence in sequence and expression. Based on comparison of A. thaliana and R. raphanistrum ortholog floral expression levels, retained radish duplicates diverged primarily via maintenance of ancestral expression level in one copy and reduction of expression level in others. In addition, retained duplicates differed significantly from genes that reverted to singleton state in function, sequence composition, expression patterns, network connectivity, and rates of evolution. Using these properties, we established a statistical learning model for predicting whether a duplicate would be retained postpolyploidization. Overall, our study provides new insights into the processes of plant duplicate loss, retention, and functional divergence and highlights the need for further understanding factors controlling duplicate gene fate.     

Comparative Study of the Effect of Snake Venoms on the Growth of Ciliates Tetrahymena pyriformis: Identification of Venoms with High Antiprotozoal Activity

Doklady Biochemistry and Biophysics - To search for compounds with antiprotozoal activity, effects of snake venoms on the ciliates Tetrahymena pyriformis was studied. T. pyriformis from subkingdom...     

Flavivirus Replication Complex Assembly Revealed by DNAJC14 Functional Mapping

DNAJC14 is an Hsp40 family member that broadly modulates flavivirus replication. The mechanism by which DNAJC14 stoichiometrically participates in flavivirus replication complex (RC) formation is unknown; both reduced and elevated levels result in replication inhibition. Using yellow fever virus (YFV), we demonstrate that DNAJC14 redistributes and clusters with YFV nonstructural proteins via a transmembrane domain and a newly identified membrane-binding domain (MBD), which both mediate targeting to detergent-resistant membranes. Furthermore, the RC and DNAJC14 reside as part of a protein interaction network that remains after 1% Triton solubilization. Mutagenesis studies demonstrate that entry into this protein interaction network requires the DNAJC14 C-terminal self-interaction domain. Fusion of the DNAJC14 MBD and self-interaction domain with another Hsp40 family protein is sufficient to confer YFV-inhibitory activity. Our findings support a novel model of DNAJC14 action that includes specific membrane targeting of both DNAJC14 and YFV replication proteins, the formation of protein interactions, and a microdomain-specific chaperone event leading to RC formation. This process alters the properties of the RC membrane and results in the formation of a protein scaffold that maintains the RC.     

Soil Microbial Community Responses to a Decade of Warming as Revealed by Comparative Metagenomics

Soil microbial communities are extremely complex, being composed of thousands of low-abundance species (<0.1% of total). How such complex communities respond to natural or human-induced fluctuations, including major perturbations such as global climate change, remains poorly understood, severely limiting our predictive ability for soil ecosystem functioning and resilience. In this study, we compared 12 whole-community shotgun metagenomic data sets from a grassland soil in the Midwestern United States, half representing soil that had undergone infrared warming by 2°C for 10 years, which simulated the effects of climate change, and the other half representing the adjacent soil that received no warming and thus, served as controls. Our analyses revealed that the heated communities showed significant shifts in composition and predicted metabolism, and these shifts were community wide as opposed to being attributable to a few taxa. Key metabolic pathways related to carbon turnover, such as cellulose degradation (∼13%) and CO₂ production (∼10%), and to nitrogen cycling, including denitrification (∼12%), were enriched under warming, which was consistent with independent physicochemical measurements. These community shifts were interlinked, in part, with higher primary productivity of the aboveground plant communities stimulated by warming, revealing that most of the additional, plant-derived soil carbon was likely respired by microbial activity. Warming also enriched for a higher abundance of sporulation genes and genomes with higher G+C content. Collectively, our results indicate that microbial communities of temperate grassland soils play important roles in mediating feedback responses to climate change and advance the understanding of the molecular mechanisms of community adaptation to environmental perturbations.     

Lack of Genetic Variation of Bursaphelenchus xylophilus in Portugal Revealed by RAPD-PCR Analyses

Random Amplified Polymorphic DNA (RAPD-PCR) technique was used to assess the level of genetic variability and genetic relationships among 24 Portuguese isolates of pinewood nematode, Bursaphelenchus xylophilus. The isolates represent the main infested areas of Portugal. Two additional isolates of B. xylophilus representing North America and East Asia were included, and B. mucronatus was used as out-group. Twenty-eight random primers generated a total of 640 DNA fragments. The Nei and Li similarity index revealed a high genetic similarity among the Portuguese isolates (above 90%). Hierarchical cluster analysis was performed to illustrate the relatedness among the isolates. No indication for separate groups among the Portuguese isolates was obtained, and the low level of genetic diversity strongly suggests that they were dispersed recently from a single introduction. The lack of apparent relationship between the genetic and the geographic matrices of the Portuguese isolates limits the use of this technique for following recent pathways of distribution. Genetic distance of the Portuguese isolates towards an isolate from China was much lower as compared to an isolate from the USA. This confirmed previous results suggesting an East Asian origin of the Portuguese B. xylophilus.     

Cell-to-Cell Diversity in a Synchronized Chlamydomonas Culture As Revealed by Single-Cell Analyses

In a synchronized photoautotrophic culture of Chlamydomonas reinhardtii, cell size, cell number, and the averaged starch content were determined throughout the light-dark cycle. For single-cell analyses, the relative cellular starch was quantified by measuring the second harmonic generation (SHG). In destained cells, amylopectin essentially represents the only biophotonic structure. As revealed by various validation procedures, SHG signal intensities are a reliable relative measure of the cellular starch content. During photosynthesis-driven starch biosynthesis, synchronized Chlamydomonas cells possess an unexpected cell-to-cell diversity both in size and starch content, but the starch-related heterogeneity largely exceeds that of size. The cellular volume, starch content, and amount of starch/cell volume obey lognormal distributions. Starch degradation was initiated by inhibiting the photosynthetic electron transport in illuminated cells or by darkening. Under both conditions, the averaged rate of starch degradation is almost constant, but it is higher in illuminated than in darkened cells. At the single-cell level, rates of starch degradation largely differ but are unrelated to the initial cellular starch content. A rate equation describing the cellular starch degradation is presented. SHG-based three-dimensional reconstructions of Chlamydomonas cells containing starch granules are shown.