Analysis of new retrogenes provides insight into dog adaptive evolution

The origin and subsequent evolution of new genes have been considered as an important source of genetic and phenotypic diversity in organisms. Dog breeds show great phenotypic diversity for morphological, physiological, and behavioral traits. However, the contributions of newly originated retrogenes, which provide important genetic bases for dog species differentiation and adaptive traits, are largely unknown. Here, we analyzed the dog genome to identify new RNA‐based duplications and comprehensively investigated their origin, evolution, functions in adaptive traits, and gene movement processes. First, we totally identified 3,025 retrocopies including 476 intact retrogenes, 2,518 retropseudogenes, and 31 chimerical retrogenes. Second, selective pressure along with ESTs expression analysis showed that most of the intact retrogenes were significantly under stronger purifying selection and subjected to more functional constraints when compared to retropseudogenes. Furthermore, a large number of retrocopies and chimerical retrogenes that occurred approximately 22 million years ago implied a burst of retrotransposition in the dog genome after the divergence time between dog and its closely related species red fox. Interestingly, GO and pathway analyses showed that new retrogenes had expanded in glutathione biosynthetic/metabolic process which likely provided important genetic basis for dogs' adaptation to scavenge human waste dumps. Finally, consistent with the results in human and mouse, a significant excess of functional retrogenes movement on and off the X chromosome in the dog confirmed a general pattern of gene movement process in mammals which was likely driven by natural selection or sexual antagonism. Together, these results increase our understanding that new retrogenes can reshape the dog genome and provide further exploration of the molecular mechanisms underlying the dogs' adaptive evolution.     

Analysis of an ascidian integrin provides new insight into early evolution of collagen recognition

αI domain integrins have been found in the ascidian Ciona intestinalis. We produced Ciona α1I domain as a recombinant protein. It did not recognize fibril-forming collagens or bind to GFOGER or other similar motifs in triple-helical peptides. No GFOGER motifs were found in Ciona collagens. As Ciona α1I bound to collagen IX, we propose that before the emergence of GFOGER-dependent collagen receptors in vertebrates, αI domain integrins might have been able to bind to collagen with alternative mechanisms.     

The improved assembly of 7DL chromosome provides insight into the structure and evolution of bread wheat

Wheat is one of the most important staple crops worldwide and also an excellent model species for crop evolution and polyploidization studies. The breakthrough of sequencing the bread wheat genome and progenitor genomes lays the foundation to decipher the complexity of wheat origin and evolutionary process as well as the genetic consequences of polyploidization. In this study, we sequenced 3286 BACs from chromosome 7DL of bread wheat cv. Chinese Spring and integrated the unmapped contigs from IWGSC v1 and available PacBio sequences to close gaps present in the 7DL assembly. In total, 8043 out of 12 825 gaps, representing 3 491 264 bp, were closed. We then used the improved assembly of 7DL to perform comparative genomic analysis of bread wheat (Ta7DL) and its D donor, Aegilops tauschii (At7DL), to identify domestication signatures. Results showed a strong syntenic relationship between Ta7DL and At7DL, although some small rearrangements were detected at the distal regions. A total of 53 genes appear to be lost genes during wheat polyploidization, with 23% (12 genes) as RGA (disease resistance gene analogue). Furthermore, 86 positively selected genes (PSGs) were identified, considered to be domestication‐related candidates. Finally, overlapping of QTLs obtained from GWAS analysis and PSGs indicated that TraesCS7D02G321000 may be one of the domestication genes involved in grain morphology. This study provides comparative information on the sequence, structure and organization between bread wheat and Ae. tauschii from the perspective of the 7DL chromosome, which contribute to better understanding of the evolution of wheat, and supports wheat crop improvement.     

The amphioxus genome provides unique insight into the evolution of immunity

Immune systems evolve as essential strategies to maintain homeostasis with the environment, prevent microbial assault and recycle damaged host tissues. The immune system is composed of two components, innate and adaptive immunity. The former is common to all animals while the latter consists of a vertebrate-specific system that relies on somatically derived lymphocytes and is associated with near limitless genetic diversity as well as long-term memory. Deuterostome invertebrates provide a view of immune repertoires in phyla that immediately predate the origins of vertebrates. Genomic studies in amphioxus, a cephalochordate, have revealed homologs of genes encoding most innate immune receptors found in vertebrates; however, many of the gene families have undergone dramatic expansions, greatly increasing the innate immune repertoire. In addition, domain-swapping accounts for the innovation of new predicted pathways of receptor function. In both amphioxus and Ciona, a urochordate, the VCBPs (variable region containing chitin-binding proteins), which consist of immunoglobulin V (variable) and chitin binding domains, mediate recognition through the V domains. The V domains of VCBPs in amphioxus exhibit high levels of allelic complexity that presumably relate to functional specificity. Various features of the amphioxus immune repertoire reflect novel selective pressures, which likely have resulted in innovative strategies. Functional genomic studies underscore the value of amphioxus as a model for studying innate immunity and may help reveal how unique relationships between innate immune receptors and both pathogens and symbionts factored in the evolution of adaptive immune systems.     

The Chimonanthus salicifolius genome provides insight into magnoliid evolution and flavonoid biosynthesis

Chimonanthus salicifolius, a member of the Calycanthaceae of magnoliids, is one of the most famous medicinal plants in Eastern China. Here, we report a chromosome‐level genome assembly of C. salicifolius, comprising 820.1 Mb of genomic sequence with a contig N50 of 2.3 Mb and containing 36 651 annotated protein‐coding genes. Phylogenetic analyses revealed that magnoliids were sister to the eudicots. Two rounds of ancient whole‐genome duplication were inferred in the C. salicifolious genome. One is shared by Calycanthaceae after its divergence with Lauraceae, and the other is in the ancestry of Magnoliales and Laurales. Notably, long genes with > 20 kb in length were much more prevalent in the magnoliid genomes compared with other angiosperms, which could be caused by the length expansion of introns inserted by transposon elements. Homologous genes within the flavonoid pathway for C. salicifolius were identified, and correlation of the gene expression and the contents of flavonoid metabolites revealed potential critical genes involved in flavonoids biosynthesis. This study not only provides an additional whole‐genome sequence from the magnoliids, but also opens the door to functional genomic research and molecular breeding of C. salicifolius.     

Proteomic analysis provides new insight into the chicken eggshell cuticle

The cuticle is the outermost layer of the avian eggshell, whose protein constituents remain virtually unknown. We hypothesize that cuticle components play a major role in microbial resistance, since eggs with incomplete or absent cuticle are more susceptible to bacterial contamination. In this study we extracted proteins from the outermost non-calcified layer of the cuticle of chicken eggs and subjected them to LC/MS/MS proteomic analysis. We identified 47 cuticle proteins with high confidence and reproducibility. Two proteins, similar to Kunitz-like protease inhibitor and ovocalyxin-32 (a carboxypeptidase A inhibitor), were the most abundant of the cuticle proteins. A number of proteins known to have antimicrobial activity in the egg were detected (lysozyme C, ovotransferrin, ovocalyxin-32, cystatin, ovoinhibitor) as well as possible new candidates (myeloperoxidase, ovocalyxin-36 and members of the SERPIN family). This is the first comprehensive report of cuticle proteome, a starting point to determine cuticle function and the molecular basis of its antimicrobial properties.     

The crystal structure of a virus-like particle from the hyperthermophilic archaeon Pyrococcus furiosus provides insight into the evolution of viruses

Pyrococcus furiosus is a hyperthermophilic archaeal microorganism found near deep-sea thermal vents and its optimal growth temperature of 100 degrees C. Recently, a 38.8-kDa protein from P. furiosus DSM 3638 was isolated and characterized. Electron microscopy revealed that this protein aggregated as spheres of approximately 30 nm in diameter, which we designated P. furiosus virus-like particles (PfVs). X-ray crystallographic analysis at 3.6-A resolution revealed that each PfV consisted of 180 copies of the 38.8-kDa protein and retained T=3 icosahedral symmetry, as is often the case in spherical viruses. The total molecular mass of each particle was approximately 7 MDa. An examination of capsid structures suggested strong evolutionary links among PfV, tailed double-stranded DNA bacteriophages, and herpes viruses. The similar three-dimensional structures of the various coat proteins indicate that these viral capsids might have originated and evolved from a common ancestor. The structure of PfV provides a previously undescribed example of viral relationships across the three domains of life (Eukarya, Bacteria, and Archaea).     

The X-ray crystal structure of human beta-hexosaminidase B provides new insights into Sandhoff disease

Human lysosomal beta-hexosaminidases are dimeric enzymes composed of alpha and beta-chains, encoded by the genes HEXA and HEXB. They occur in three isoforms, the homodimeric hexosaminidases B (betabeta) and S (alphaalpha), and the heterodimeric hexosaminidase A (alphabeta), where dimerization is required for catalytic activity. Allelic variations in the HEXA and HEXB genes cause the fatal inborn errors of metabolism Tay-Sachs disease and Sandhoff disease, respectively. Here, we present the crystal structure of a complex of human beta-hexosaminidase B with a transition state analogue inhibitor at 2.3A resolution (pdb 1o7a). On the basis of this structure and previous studies on related enzymes, a retaining double-displacement mechanism for glycosyl hydrolysis by beta-hexosaminidase B is proposed. In the dimer structure, which is derived from an analysis of crystal packing, most of the mutations causing late-onset Sandhoff disease reside near the dimer interface and are proposed to interfere with correct dimer formation. The structure reported here is a valid template also for the dimeric structures of beta-hexosaminidase A and S.     

X-ray structure of 12-oxophytodienoate reductase 1 provides structural insight into substrate binding and specificity within the family of OYE

BACKGROUND: 12-Oxophytodienoate reductase (OPR) is a flavin mononucleotide (FMN)-dependent oxidoreductase in plants that belongs to the family of Old Yellow Enzyme (OYE). It was initially characterized as an enzyme involved in the biosynthesis of the plant hormone jasmonic acid, where it catalyzes the reduction of the cyclic fatty acid derivative 9S,13S-12-oxophytodienoate (9S,13S-OPDA) to 1S,2S-3-oxo-2(2'[Z]-pentenyl)-cyclopentane-1-octanoate. Several isozymes of OPR are now known that show different stereoselectivities with regard to the four stereoisomers of OPDA. RESULTS: Here, we report the high-resolution crystal structure of OPR1 from Lycopersicon esculentum and its complex structures with the substrate 9R,13R-OPDA and with polyethylene glycol 400. OPR1 crystallizes as a monomer and folds into a (betaalpha)(8) barrel with an overall structure similar to OYE. The cyclopentenone ring of 9R,13R-OPDA is stacked above the flavin and activated by two hydrogen bonds to His187 and His190. The olefinic bond is properly positioned for hydride transfer from the FMN N(5) and proton transfer from Tyr192 to Cbeta and Calpha, respectively. Comparison of the OPR1 and OYE structures reveals striking differences in the loops responsible for binding 9R,13R-OPDA in OPR1. CONCLUSIONS: Despite extensive biochemical characterization, the physiological function of OYE still remains unknown. The similar catalytic cavity structures and the substrate binding mode in OPR1 strongly support the assumption that alpha,beta-unsaturated carbonyl compounds are physiological substrates of the OYE family. The specific binding of 9R,13R-OPDA by OPR1 explains the experimentally observed stereoselectivity and argues in favor of 9R,13R-OPDA or a structurally related oxylipin as natural substrate of OPR1.     

Solution structure of DinI provides insight into its mode of RecA inactivation

The Escherichia coli RecA protein triggers both DNA repair and mutagenesis in a process known as the SOS response. The 81-residue E. coli protein DinI inhibits activity of RecA in vivo. The solution structure of DinI has been determined by multidimensional triple resonance NMR spectroscopy, using restraints derived from two sets of residual dipolar couplings, obtained in bicelle and phage media, supplemented with J couplings and a moderate number of NOE restraints. DinI has an alpha/beta fold comprised of a three-stranded beta-sheet and two alpha-helices. The beta-sheet topology is unusual: the central strand is flanked by a parallel and an antiparallel strand and the sheet is remarkably flat. The structure of DinI shows that six negatively charged Glu and Asp residues on DinI's kinked C-terminal alpha-helix form an extended, negatively charged ridge. We propose that this ridge mimics the electrostatic character of the DNA phospodiester backbone, thereby enabling DinI to compete with single-stranded DNA for RecA binding. Biochemical data confirm that DinI is able to displace ssDNA from RecA.     

An engineered construct of cFLIP provides insight into DED1 structure and interactions

1. Download : Download high-res image (96KB)
2. Download : Download full-size image

     

The solution NMR structure of Antheraea polyphemus PBP provides new insight into pheromone recognition by pheromone-binding proteins

Pheromone-binding proteins (PBPs) located in the antennae of male moth species play an important role in olfaction. They are carrier proteins, believed to transport volatile hydrophobic pheromone molecules across the aqueous sensillar lymph to the membrane-bound G protein-coupled olfactory receptor proteins. The roles of PBPs in molecular recognition and the mechanisms of pheromone binding and release are poorly understood. Here, we report the NMR structure of a PBP from the giant silk moth Antheraea polyphemus. This is the first structure of a PBP with specific acetate-binding function in vivo. The protein consists of nine alpha-helices: alpha1a (residues 2-5), alpha1b (8-12), alpha1c (16-23), alpha2 (27-34), alpha3a (46-52), alpha3b (54-59), alpha4 (70-79), alpha5 (84-100) and alpha6 (107-125), held together by three disulfide bridges: 19-54, 50-108 and 97-117. A large hydrophobic cavity is located inside the protein, lined with side-chains from all nine helices. The acetate-binding site is located at the narrow end of the cavity formed by the helices alpha3b and alpha4. The pheromone can enter this cavity through an opening between the helix alpha1a, the C-terminal end of the helix alpha6, and the loop between alpha2 and alpha3a. We suggest that Trp37 may play an important role in the initial interaction with the ligand. Our analysis also shows that Asn53 plays the key role in recognition of acetate pheromones specifically, while Phe12, Phe36, Trp37, Phe76, and Phe118 are responsible for non-specific binding, and Leu8 and Ser9 may play a role in ligand chain length recognition.     

The crystal structure of Trichomonas vaginalis ferredoxin provides insight into metronidazole activation

Crystallographic studies revealing the three-dimensional structure of the oxidized form of the [2Fe-2S] ferredoxin from Trichomonas vaginalis (TvFd) are presented. TvFd, a member of the hydrogenosomal class of ferredoxins, possesses a unique combination of redox and spectroscopic properties, and is believed to be the biological molecule that activates the drug metronidazole reductively in the treatment of trichomoniasis. It is the first hydrogenosomal ferredoxin to have its structure determined. The structure of TvFd reveals a monomeric, 93 residue protein with a fold similar to that of other known [2Fe-2S] ferredoxins. It contains nine hydrogen bonds to the sulfur atoms of the cluster, which is more than the number predicted on the basis of the spectroscopic data. The TvFd structure contains a large dipole moment like adrenodoxin, and appears to have a similar interaction domain. Our analysis demonstrates that TvFd has a unique cavity near the iron-sulfur cluster that exposes one of the inorganic sulfur atoms of the cluster to solvent. This cavity is not seen in any other [2Fe-2S] ferredoxin with known structure, and is hypothesized to be responsible for the high rate of metronidazole reduction by TvFd.     

Clostridium perfringens: insight into virulence evolution and population structure

Clostridium perfringens is an important pathogen in veterinary and medical fields. Diseases caused by this organism are in many cases life threatening or fatal. At the same time, it is part of the ecological community of the intestinal tract of man and animals. Virulence in this species is not fully understood and it does seem that there is erratic distribution of the toxin/enzyme genes within C. perfringens population. We used the recently developed multiple-locus variable-number tandem repeat analysis (MLVA) scheme to investigate the evolution of virulence and population structure of this species. Analysis of the phylogenetic signal indicates that acquisition of the major toxin genes as well as other plasmid-borne toxin genes is a recent evolutionary event and their maintenance is essentially a function of the selective advantage they confer in certain niches under different conditions. In addition, it indicates the ability of virulent strains to cause disease in different host species. More interestingly, there is evidence that certain normal flora strains are virulent when they gain access to a different host species. Analysis of the population structure indicates that recombination events are the major tool that shapes the population and this panmixia is interrupted by frequent clonal expansion that mostly corresponds to disease processes. The signature of positive selection was detected in alpha toxin gene, suggesting the possibility of adaptive alleles on the other chromosomally encoded determinants. Finally, C. perfringens proved to have a dynamic population and availability of more genome sequences and use of comparative proteomics and animal modeling would provide more insight into the virulence of this organism.     

Structural insight into the evolution of a new chemokine family from zebrafish

The mammalian chemokine family is segregated into four families – CC, CXC, CX3C, and XC—based on the arrangement of cysteines and the corresponding disulfides. Sequencing of the Danio rerio (zebrafish) genome has identified more than double the amount of human chemokines with the absence of the CX3C family and the presence of a new family, CX. The only other family with a single cysteine in the N‐terminal region is the XC family. Human lymphotactin (XCL1) has two interconverting structures due to dynamic changes that occur in the protein. Similar to an experiment with XCL1 that identified the two structural forms, we probed for multiple forms of zCXL1 using heparin affinity. The results suggest only a single form of CXL1 is present. We used sulfur‐SAD phasing to determine the three‐dimensional structure CXL1. Zebrafish CXL1 (zCXL1) has three disulfides that appear to be important for a stable structure. One disulfide is common to all chemokines except those that belong to the XC family, another is similar to a subset of CC chemokines containing three disulfides, but the third disulfide is unique to the CX family. We analyzed the electrostatic potential of the zCXL1 structure and identified the likely heparin‐binding site for glycosaminoglycans (GAGs). zCXL1 has a similar sequence identity with human CCL5 and CXCL12, but the structure is more related to CCL5. Our structural analysis supports the phylogenetic and genomic studies on the evolution of the CXL family. Proteins 2014; 82:708–716. © 2013 Wiley Periodicals, Inc.