Comprehensive Evolutionary and Expression Analysis of FCS-Like Zinc finger Gene Family Yields Insights into Their Origin,Expansion and Divergence

Plant evolution is characterized by frequent genome duplication events. Expansion of habitat resulted in the origin of many novel genes and genome duplication events which in turn resulted in the expansion of many regulatory gene families. The plant-specific FCS-Like Zinc finger (FLZ) gene family is characterized by the presence of a FCS-Like Zinc finger (FLZ) domain which mediates the protein-protein interaction. In this study, we identified that the expansion of FLZ gene family size in different species is correlated with ancestral and lineage-specific whole genome duplication events. The subsequent gene loss found to have a greater role in determining the size of this gene family in many species. However, genomic block duplications played the significant role in the expansion of FLZ gene family in some species. Comparison of Arabidopsis thaliana and Oryza sativa FLZ gene family revealed monocot and dicot specific evolutionary trends. The FLZ genes were found to be under high purifying selection. The spatiotemporal expression analyses of Arabidopsis thaliana FLZ gene family revealed that majority of the members are highly expressed in reproductive organs. FLZ genes were also found to be highly expressed during vegetative-to-reproductive phase transition which is correlated with the proposed role of this gene family in sugar signaling. The comparison of sequence, structural and expression features of duplicated genes identified lineage-specific redundancy and divergence. This extensive evolutionary analysis and expression analysis of Arabidopsis thaliana FLZ genes will pave the way for further functional analysis of FLZ genes.     

Cloning of a Gene Encoding Bombyxin, an Insulin-Like Brain Secretory Peptide of the Silkmoth Bombyx mori with Prothoracicotropic Activity

A genomic DNA encoding bombyxin, a 5kD brain peptide of the silkmoth Bornbyx mori with prothoracicotropic hormone activity, has been isolated. The nucleotide sequence coding for bombyxin shows high homology with insulin-gene family members and the overall organization of the preprobombyxin gene is the same as in preproinsulin genes, indicating that bombyxin shares a common ancestral molecule with insulin-family peptides. The bombyxin gene has no intron contrasting to other members of insulin-gene family.     

Bombyxin: An Insect Brain Peptide that Belongs to the Insulin Family

Bombyxin is a 5 kDa secretory brain peptide that belongs to the insulin family. Bombyxin of the silkmoth Bombyx mori can induce adult development when injected into brain-removed dormant pupae of the saturniid moth Samia cynthia ricini by activating the prothoracic glands to synthesize and release ecdysone. Bombyx bombyxin has been shown to lower the concentration of the major haemolymph sugar, trehalose, and to elevate the trehalase activity in the midgut and muscles in Bombyx, but the doses required to be effective are higher than the amounts in the feeding larvae. The exact physiological function of bombyxin in Bombyx itself is therefore still obscure, but its insulin-like structure suggests it has important roles. Bombyxin comprises a mixture of highly heterogeneous molecular forms whose amino acid sequences have 40% identity with human insulin. The Bombyx bombyxin gene encodes a precursor consisting of the signal peptide, B chain, C peptide, and A chain, in that order from the N terminus. So far, 32 bombyxin genes have been identified in Bombyx, and they are classified into 7 families, A to G, according to their sequence similarity. The bombyxin genes have no introns and cluster in unique distribution patterns. The gene arrangement in the cluster has been classified into three categories: gene pairs, gene triplets, and single genes. Nucleotide sequence analysis indicates that equal and unequal crossings-over and duplications may have generated these unique distribution patterns. The Bombyx bombyxin genes are expressed predominantly in the brain and at low levels in a number of other tissues. Genes of all 7 families are expressed in four pairs of the medial neurosecretory cells of the brain. Detailed examination indicated that only a limited number of genes in the A, B and C family members are expressed and that their expression shows a gene-arrangement-dependent pattern.     

Functional Knowledge Transfer for High-accuracy Prediction of Under-studied Biological Processes

A key challenge in genetics is identifying the functional roles of genes in pathways. Numerous functional genomics techniques (e.g. machine learning) that predict protein function have been developed to address this question. These methods generally build from existing annotations of genes to pathways and thus are often unable to identify additional genes participating in processes that are not already well studied. Many of these processes are well studied in some organism, but not necessarily in an investigator''s organism of interest. Sequence-based search methods (e.g. BLAST) have been used to transfer such annotation information between organisms. We demonstrate that functional genomics can complement traditional sequence similarity to improve the transfer of gene annotations between organisms. Our method transfers annotations only when functionally appropriate as determined by genomic data and can be used with any prediction algorithm to combine transferred gene function knowledge with organism-specific high-throughput data to enable accurate function prediction.We show that diverse state-of-art machine learning algorithms leveraging functional knowledge transfer (FKT) dramatically improve their accuracy in predicting gene-pathway membership, particularly for processes with little experimental knowledge in an organism. We also show that our method compares favorably to annotation transfer by sequence similarity. Next, we deploy FKT with state-of-the-art SVM classifier to predict novel genes to 11,000 biological processes across six diverse organisms and expand the coverage of accurate function predictions to processes that are often ignored because of a dearth of annotated genes in an organism. Finally, we perform in vivo experimental investigation in Danio rerio and confirm the regulatory role of our top predicted novel gene, wnt5b, in leftward cell migration during heart development. FKT is immediately applicable to many bioinformatics techniques and will help biologists systematically integrate prior knowledge from diverse systems to direct targeted experiments in their organism of study.     

Identification of a novel type of WRKY transcription factor binding site in elicitor-responsive cis-sequences from Arabidopsis thaliana

Using a combination of bioinformatics and synthetic promoters, novel elicitor-responsive cis-sequences were discovered in promoters of pathogen-upregulated genes from Arabidopsis thaliana. One group of functional sequences contains the conserved core sequence GACTTTT. This core sequence and adjacent nucleotides are essential for elicitor-responsive gene expression in a parsley protoplast system. By yeast one-hybrid screening, WRKY70 was selected with a cis-sequence harbouring the core sequence GACTTTT but no known WRKY binding site (W-box). Transactivation experiments, mutation analyses, and electrophoretic mobility shift assays demonstrate that the sequence CGACTTTT is the binding site for WRKY70 in the investigated cis-sequence and is required for WRKY70-activated gene expression. Using several cis-sequences in transactivation experiments and binding studies, the CGACTTTT sequence can be extended to propose YGACTTTT as WRKY70 binding site. This binding site, designated WT-box, is enriched in promoters of genes upregulated in a WRKY70 overexpressing line. Interestingly, functional WRKY70 binding sites are present in the promoter of WRKY30, supporting recent evidence that both factors play a role in the same regulatory network.     

Lamins of the sea lamprey (Petromyzon marinus) and the evolution of the vertebrate lamin protein family

Lamin proteins are found in all metazoans. Most non-vertebrate genomes including those of the closest relatives of vertebrates, the cephalochordates and tunicates, encode only a single lamin. In teleosts and tetrapods the number of lamin genes has quadrupled. They can be divided into four sub-types, lmnb1, lmnb2, LIII, and lmna, each characterized by particular features and functional differentiations. Little is known when during vertebrate evolution these features have emerged. Lampreys belong to the Agnatha, the sister group of the Gnathostomata. They split off first within the vertebrate lineage. Analysis of the sea lamprey (Petromyzon marinus) lamin complement presented here, identified three functional lamin genes, one encoding a lamin LIII, indicating that the characteristic gene structure of this subtype had been established prior to the agnathan/gnathostome split. Two other genes encode lamins for which orthology to gnathostome lamins cannot be designated. Search for lamin gene sequences in all vertebrate taxa for which sufficient sequence data are available reveals the evolutionary time frame in which specific features of the vertebrate lamins were established. Structural features characteristic for A-type lamins are not found in the lamprey genome. In contrast, lmna genes are present in all gnathostome lineages suggesting that this gene evolved with the emergence of the gnathostomes. The analysis of lamin gene neighborhoods reveals noticeable similarities between the different vertebrate lamin genes supporting the hypothesis that they emerged due to two rounds of whole genome duplication and makes clear that an orthologous relationship between a particular vertebrate paralog and lamins outside the vertebrate lineage cannot be established.     

Unbiased Functional Clustering of Gene Variants with a Phenotypic-Linkage Network

Groupwise functional analysis of gene variants is becoming standard in next-generation sequencing studies. As the function of many genes is unknown and their classification to pathways is scant, functional associations between genes are often inferred from large-scale omics data. Such data types—including protein–protein interactions and gene co-expression networks—are used to examine the interrelations of the implicated genes. Statistical significance is assessed by comparing the interconnectedness of the mutated genes with that of random gene sets. However, interconnectedness can be affected by confounding bias, potentially resulting in false positive findings. We show that genes implicated through de novo sequence variants are biased in their coding-sequence length and longer genes tend to cluster together, which leads to exaggerated p-values in functional studies; we present here an integrative method that addresses these bias. To discern molecular pathways relevant to complex disease, we have inferred functional associations between human genes from diverse data types and assessed them with a novel phenotype-based method. Examining the functional association between de novo gene variants, we control for the heretofore unexplored confounding bias in coding-sequence length. We test different data types and networks and find that the disease-associated genes cluster more significantly in an integrated phenotypic-linkage network than in other gene networks. We present a tool of superior power to identify functional associations among genes mutated in the same disease even after accounting for significant sequencing study bias and demonstrate the suitability of this method to functionally cluster variant genes underlying polygenic disorders.     

Identification of a Novel Gig2 Gene Family Specific to Non-Amniote Vertebrates

Gig2 (grass carp reovirus (GCRV)-induced gene 2) is first identified as a novel fish interferon (IFN)-stimulated gene (ISG). Overexpression of a zebrafish Gig2 gene can protect cultured fish cells from virus infection. In the present study, we identify a novel gene family that is comprised of genes homologous to the previously characterized Gig2. EST/GSS search and in silico cloning identify 190 Gig2 homologous genes in 51 vertebrate species ranged from lampreys to amphibians. Further large-scale search of vertebrate and invertebrate genome databases indicate that Gig2 gene family is specific to non-amniotes including lampreys, sharks/rays, ray-finned fishes and amphibians. Phylogenetic analysis and synteny analysis reveal lineage-specific expansion of Gig2 gene family and also provide valuable evidence for the fish-specific genome duplication (FSGD) hypothesis. Although Gig2 family proteins exhibit no significant sequence similarity to any known proteins, a typical Gig2 protein appears to consist of two conserved parts: an N-terminus that bears very low homology to the catalytic domains of poly(ADP-ribose) polymerases (PARPs), and a novel C-terminal domain that is unique to this gene family. Expression profiling of zebrafish Gig2 family genes shows that some duplicate pairs have diverged in function via acquisition of novel spatial and/or temporal expression under stresses. The specificity of this gene family to non-amniotes might contribute to a large extent to distinct physiology in non-amniote vertebrates.     

Analysis of expressed sequence tags (ESTs) and gene expression changes under different growth conditions for the ciliate Anophryoides haemophila, the causative agent of bumper car disease in the American lobster (Homarus americanus)

The scuticociliate Anophryoides haemophila, causes bumper car disease in American lobster (Homarus americanus) in commercial holding facilities in Atlantic Canada. While the parasite has been recognized since the 1970s and much has been learned about its biology, minimal molecular characterization exists. With genome consortiums turning to model organisms like the ciliates Tetrahymena and Paramecium, the amount of relevant sequence data available has made sequence surveys more attractive for gene discovery in related ciliates. We sequenced 9984 expressed sequence tags (ESTs) from a non-normalized A. haemophila cDNA library to characterize gene expression patterns, functional gene distribution and to discover novel genes related to the parasitic life history. The A. haemophila ESTs were grouped into 843 clusters and singletons with 658 EST clusters having identifiable homologs, while 159 ESTs were unique and had no similarity to any sequences in the public databases. Not unexpectedly, about 67% of the A. haemophila ESTs have similarity to annotated and hypothetical genes from the related oligohymenophorean ciliate, Tetrahymena. Numerous cysteine proteases, hypothetical proteins and novel sequences possess putative secretory signal peptides suggesting that they may contribute to the pathogenesis of bumper car disease in lobster. Real time RT-qPCR analysis of cathepsin L and two homologs of cathepsin B did not show any changes in gene expression under varying in vitro growth conditions or during a modified-in vivo infection which may be suggestive of the opportunistic life history strategy of this ciliate.     

Complete mitochondrial genome sequence and identification of a candidate gene responsible for cytoplasmic male sterility in radish (Raphanus sativus L.) containing DCGMS cytoplasm

A novel cytoplasmic male sterility (CMS) conferred by Dongbu cytoplasmic and genic male-sterility (DCGMS) cytoplasm and its restorer-of-fertility gene (Rfd1) was previously reported in radish (Raphanus sativus L.). Its inheritance of fertility restoration and profiles of mitochondrial DNA (mtDNA)-based molecular markers were reported to be different from those of Ogura CMS, the first reported CMS in radish. The complete mitochondrial genome sequence (239,186 bp; GenBank accession No. KC193578) of DCGMS mitotype is reported in this study. Thirty-four protein-coding genes and three ribosomal RNA genes were identified. Comparative analysis of a mitochondrial genome sequence of DCGMS and previously reported complete sequences of normal and Ogura CMS mitotypes revealed various recombined structures of seventeen syntenic sequence blocks. Short-repeat sequences were identified in almost all junctions between syntenic sequence blocks. Phylogenetic analysis of three radish mitotypes showed that DCGMS was more closely related to the normal mitotype than to the Ogura mitotype. A single 1,551-bp unique region was identified in DCGMS mtDNA sequences and a novel chimeric gene, designated orf463, consisting of 128-bp partial sequences of cox1 gene and 1,261-bp unidentified sequences were found in the unique region. No other genes with a chimeric structure, a major feature of most characterized CMS-associated genes in other plant species, were found in rearranged junctions of syntenic sequence blocks. Like other known CMS-associated mitochondrial genes, the predicted gene product of orf463 contained 12 transmembrane domains. Thus, this gene product might be integrated into the mitochondrial membrane. In total, the results indicate that orf463 is likely to be a casual factor for CMS induction in radish containing the DCGMS cytoplasm.     

Phylogenetic analysis and evolutionary studies of plant carotenoid cleavage dioxygenase gene

The oxidative breakdown of carotenoid evidences the formation of apocarotenoids through carotenoid cleavage dioxygenases (CCDs). Numerous CCDs and apocarotenoids have been identified and characterized in plants. Using available sequence data, a study was performed to investigate the phylogenetic relationship among CCD genes and to statistically estimate the sequence conservation and functional divergence. In total, 77 genes were identified from 39 species belonging to 21 families. Our result of phylogenetic analysis indicated the existence of well-conserved subfamilies. Moreover, comparative genomic analysis showed that the gene structures of the CCDs were highly conserved across some different lineage species. Through functional divergence analysis, a substantial divergence was found between CCD subfamilies. In addition, examination of the site-specific profile revealed the critical amino acid residues accounting for functional divergence. This study mainly focused on the evolution of CCD genes and their functional divergence which may deliver an initial step for further experimental verifications.