Optimizing the Parametrization of Homologue Classification in the Pan-Genome Computation for a Bacterial Species: Case Study Streptococcus pyogenes

The paradigm shift associated with the introduction of the pan-genome concept has drawn the attention from singular reference genomes toward the actual sequence diversity within organism populations, strain collections, clades, etc. A single genome is no longer sufficient to describe bacteria of interest, but instead, the genomic repertoire of all existing strains is the key to the metabolic, evolutionary, or pathogenic potential of a species. The classification of orthologous genes derived from a collection of taxonomically related genome sequences is central to bacterial pan-genome computational analysis. In this work, we present a review of methods for computing pan-genome gene clusters including their comparative analysis for the case of Streptococcus pyogenes strain genomes. We exhaustively scanned the parametrization space of the homologue searching procedures and find optimal parameters (sequence identity (60%) and coverage (50-60%) in the pairwise alignment) for the orthologous clustering of gene sequences. We find that the sequence identity threshold influences the number of gene families ~3 times stronger than the sequence coverage threshold.

Allelic variation in sugary1 gene affecting kernel sweetness among diverse-mutant and -wild-type maize inbreds

Sweet corn is popular worldwide as vegetable. Though large numbers of sugary1 (su1)-based sweet corn germplasm are available, allelic diversity in su1 gene encoding SU1 isoamylase among diverse maize inbreds has not been analyzed. Here, we characterized the su1 gene in maize and compared with allied species. The entire su1 gene (11,720 bp) was sequenced among six mutant (su1) and five wild (Su1) maize inbreds. Fifteen InDels of 2-45 bp were selected to develop markers for studying allelic diversity in su1 gene among 19 mutant- (su1) and 29 wild-type (Su1) inbreds. PIC ranged from 0.15 (SU-InDel7) to 0.37 (SU-InDel13). Major allele frequency varied from 0.52 to 0.90, while gene diversity ranged from 0.16 to 0.49. Phylogenetic tree categorized 48 maize inbreds in two clusters each for wild- type (Su1) and mutant (su1) types. 44 haplotypes of su1 were observed, with three haplotypes (Hap6, Hap22 and Hap29) sharing more than one genotype. Further, comparisons were made with 23 orthologues of su1 from 16 grasses and Arabidopsis. Maize possessed 15-19 exons in su1, while it was 11-24 exons among orthologues. Introns among the orthologues were longer (77-2206 bp) than maize (859-1718 bp). SU1 protein of maize and orthologues had conserved α-amylase and CBM_48 domains. The study also provided physicochemical properties and secondary structure of SU1 protein in maize and its orthologues. Phylogenetic analysis showed closer relationship of maize SU1 protein with P. hallii, S. bicolor and E. tef than Triticum sp. and Oryza sp. The study showed that presence of high allelic diversity in su1 gene which can be utilized in the sweet corn breeding program. This is the first report of comprehensive characterization of su1 gene and its allelic forms in diverse maize and related orthologues.

Molecular characterization of teosinte branched1 gene governing branching architecture in cultivated maize and wild relatives

We sequenced the entire tb1 gene in six maize inbreds and its wild relatives (parviglumis, mexicana, perennis and luxurians) to characterize it at molecular level. Hopscotch and Tourist transposable elements were observed in the upstream of tb1 in all maize inbreds, while they were absent in wild relatives. In maize, tb1 consisted of 431-443 bp 5'UTR, 1101 bp coding sequence and 211-219 bp 3'UTR. In promoter region, mutations in the light response element in mexicana (~ 35 bp and ~ 55 bp upstream of TSS) and perennis (at ~ 35 bp upstream of TSS) were found. A 6 bp insertion at 420 bp downstream of the polyA signal site was present among teosinte accessions, while it was not observed in maize. A codominant marker flanking the 6 bp InDel was developed, and it differentiated the teosintes from maize. In Tb1 protein, alanine (12.7-14.6%) was the most abundant amino acid with tryptophan as the rarest (0.5-0.9%). The molecular weight of Tb1 protein was 38757.15 g/mol except 'Palomero Toluqueno' and HKI1128. R and TCP motifs in Tb1 protein were highly conserved across maize, teosinte and orthologues, while TCP domain differed for tb1 paralogue. Tb1 possessed important role in light-, auxin-, stress-response and meristem identity maintenance. Presence of molecular signal suggested its localization in mitochondria, nucleus and nucleolus. Parviglumis and mexicana shared closer relationship with maize than perennis and luxurians. A highly conserved 59-60 amino acids long bHLH region was observed across genotypes. Information generated here assumes significance in evolution of tb1 gene and breeding for enhancement of prolificacy in maize.

Predicted networks of protein-protein interactions in Stegodyphus mimosarum by cross-species comparisons

BACKGROUND

Stegodyphus mimosarum is a candidate model organism belonging to the class Arachnida in the phylum Arthropoda. Studies on the biology of S. mimosarum over the past several decades have consisted of behavioral research and comparison of gene sequences based on the assembled genome sequence. Given the lack of systematic protein analyses and the rich source of information in the genome, we predicted the relationships of proteins in S. mimosarum by bioinformatics comparison with genome-wide proteins from select model organisms using gene mapping.

RESULTS

The protein-protein interactions (PPIs) of 11 organisms were integrated from four databases (BioGrid, InAct, MINT, and DIP). Here, we present comprehensive prediction and analysis of 3810 proteins in S. mimosarum with regard to interactions between proteins using PPI data of organisms. Interestingly, a portion of the protein interactions conserved among Saccharomyces cerevisiae, Homo sapiens, Arabidopsis thaliana, and Drosophila melanogaster were found to be associated with RNA splicing. In addition, overlap of predicted PPIs in reference organisms, Gene Ontology, and topology models in S. mimosarum are also reported.

CONCLUSIONS

Addition of Stegodyphus, a spider representative of interactomic research, provides the possibility of obtaining deeper insights into the evolution of PPI networks among different animal species. This work largely supports the utility of the "stratus clouds" model for predicted PPIs, providing a roadmap for integrative systems biology in S. mimosarum.

Orthonome - a new pipeline for predicting high quality orthologue gene sets applicable to complete and draft genomes

BACKGROUND

Distinguishing orthologous and paralogous relationships between genes across multiple species is essential for comparative genomic analyses. Various computational approaches have been developed to resolve these evolutionary relationships, but strong trade-offs between precision and recall of orthologue prediction remains an ongoing challenge.

RESULTS

Here we present Orthonome, an orthologue prediction pipeline, designed to reduce the trade-off between orthologue capture rates (recall) and accuracy of multi-species orthologue prediction. The pipeline compares sequence domains and then forms sequence-similar clusters before using phylogenetic comparisons to identify inparalogues. It then corrects sequence similarity metrics for fragment and gene length bias using a novel scoring metric capturing relationships between full length as well as fragmented genes. The remaining genes are then brought together for the identification of orthologues within a phylogenetic framework. The orthologue predictions are further calibrated along with inparalogues and gene births, using synteny, to identify novel orthologous relationships. We use 12 high quality Drosophila genomes to show that, compared to other orthologue prediction pipelines, Orthonome provides orthogroups with minimal error but high recall. Furthermore, Orthonome is resilient to suboptimal assembly/annotation quality, with the inclusion of draft genomes from eight additional Drosophila species still providing >6500 1:1 orthologues across all twenty species while retaining a better combination of accuracy and recall than other pipelines. Orthonome is implemented as a searchable database and query tool along with multiple-sequence alignment browsers for all sets of orthologues. The underlying documentation and database are accessible at http://www.orthonome.com .

CONCLUSION

We demonstrate that Orthonome provides a superior combination of orthologue capture rates and accuracy on complete and draft drosophilid genomes when tested alongside previously published pipelines. The study also highlights a greater degree of evolutionary conservation across drosophilid species than earlier thought.

Construction of the Leaf Senescence Database and Functional Assessment of Senescence-Associated Genes

Leaf senescence is the last phase of plant development and a highly coordinated process regulated by a large number of senescence-associated genes (SAGs). By broad literature survey, we constructed a leaf senescence database (LSD) in 2011 and updated it to Version 2.0 in 2014 ( http://www.eplantsenescence.org/ and http://psd.cbi.pku.edu.cn/ ) which contains a total of 5357 genes and 324 mutants from 44 species. These SAGs were retrieved based on genetic, genomic, proteomic, physiological, or other experimental evidence and were classified into different categories according to their functions in leaf senescence or morphological phenotype of mutants. To provide comprehensive information for SAGs, we made extensive annotation by both manual and computational approaches. In addition, we predicted putative orthologues of the SAGs in other species. LSD has a user-friendly interface to allow users to make text queries or BLAST searches and to download SAGs sequences for local analysis. Functional analyses of putative SAGs reveal that WRKY75, AZF2, NAC16, and WRKY26 are positive regulators of leaf senescence, while MKP2 and CTR1 perform negative regulation to leaf senescence. This database has been served as a valuable resource for basic research on the function of SAGs and evolution of plant leaf senescence, as well as for the exploration of genetic traits in agronomically important plants.

p63 gene structure in the phylum mollusca

Roles of p53 family ancestor (p63) in the organisms' response to stressful environmental conditions (mainly pollution) have been studied among molluscs, especially in the genus Mytilus, within the last 15 years. Nevertheless, information about gene structure of this regulatory gene in molluscs is scarce. Here we report the first complete genomic structure of the p53 family orthologue in the mollusc Mediterranean mussel Mytilus galloprovincialis and confirm its similarity to vertebrate p63 gene. Our searches within the available molluscan genomes (Aplysia californica, Lottia gigantea, Crassostrea gigas and Biomphalaria glabrata), found only one p53 family member present in a single copy per haploid genome. Comparative analysis of those orthologues, additionally confirmed the conserved p63 gene structure. Conserved p63 gene structure can be a helpful tool to complement or/and revise gene annotations of any future p63 genomic sequence records in molluscs, but also in other animal phyla. Knowledge of the correct gene structure will enable better prediction of possible protein isoforms and their functions. Our analyses also pointed out possible mis-annotations of the p63 gene in sequenced molluscan genomes and stressed the value of manual inspection (based on alignments of cDNA and protein onto the genome sequence) for a reliable and complete gene annotation.