Inferring phylogenies from DNA sequences of unequal base compositions.

A new method for computing evolutionary distances between DNA sequences is proposed. Contrasting with classical methods, the underlying model does not assume that sequence base compositions (A, C, G, and T contents) are at equilibrium, thus allowing unequal base compositions among compared sequences. This makes the method more efficient than the usual ones in recovering phylogenetic trees from sequence data when base composition is heterogeneous within the data set, as we show by using both simulated and empirical data. When applied to small-subunit ribosomal RNA sequences from several prokaryotic or eukaryotic organisms, this method provides evidence for an early divergence of the microsporidian Vairimorpha necatrix in the eukaryotic lineage.

Genetic rearrangements in the rfb regions of Vibrio cholerae O1 and O139.

The recent emergence of a pathogenic new non-O1 serotype (O139) of Vibrio cholerae has led to numerous studies in an attempt to identify the origins of this new strain. Our studies indicate that O139 strains have clear differences in the surface polysaccharides when compared with O1 strains: the lipopolysaccharide can be described as semi-rough. Southern hybridization with the O1 rfb region demonstrates that O139 strains no longer contain any of the rfb genes required for the synthesis of the O1 O-antigen or its modification and also lack at least 6 kb of additional contiguous DNA. However, O139 strains have retained rfaD and have a single open reading frame closely related to three small open reading frames of the O1 rfb region. This region is closely related to the H-repeat of Escherichia coli and to the transposases of a number of insertion sequence elements and has all the features of an insertion sequence element that has been designated VcIS1. Transposon insertion mutants defective in O139 O-antigen (and capsule) biosynthesis map to the same fragment as VcIS1. Preliminary sequence data of complementing clones indicate that this DNA encodes a galactosyl-transferase and other enzymes for the utilization of galactose in polysaccharide biosynthesis. We propose a mechanism by which both the Ogawa serotype of O1 strains and the O139 serotype strains may have evolved.

Toward the therapeutic editing of mutated RNA sequences.

If RNA editing could be rationally directed to mutated RNA sequences, genetic diseases caused by certain base substitutions could be treated. Here we use a synthetic complementary RNA oligonucleotide to direct the correction of a premature stop codon mutation in dystrophin RNA. The complementary RNA oligonucleotide was hybridized to a premature stop codon and the hybrid was treated with nuclear extracts containing the cellular enzyme double-stranded RNA adenosine deaminase. When the treated RNAs were translated in vitro, a dramatic increase in expression of a downstream luciferase coding region was observed. The cDNA sequence data are consistent with deamination of the adenosine in the UAG stop codon to inosine by double-stranded RNA adenosine deaminase. Injection of oligonucleotide-mRNA hybrids into Xenopus embryos also resulted in an increase in luciferase expression. These experiments demonstrate the principle of therapeutic RNA editing.

Documentation of reticulate evolution in peonies (Paeonia) using internal transcribed spacer sequences of nuclear ribosomal DNA: implications for biogeography and concerted evolution.

The internal transcribed spacers (ITS) of nuclear ribosomal DNA of 33 species of genus Paeonia (Paeoniaceae) were sequenced. In section Paeonia, different patterns of nucleotide additivity were detected in 14 diploid and tetraploid species at sites that are variable in the other 12 species of the section, suggesting that reticulate evolution has occurred. Phylogenetic relationships of species that do not show additivity, and thus ostensibly were not derived through hybridization, were reconstructed by parsimony analysis. The taxa presumably derived through reticulate evolution were then added to the phylogenetic tree according to additivity from putative parents. The study provides an example of successfully using ITS sequences to reconstruct reticulate evolution in plants and further demonstrates that the sequence data could be highly informative and accurate for detecting hybridization. Maintenance of parental sequences in the species of hybrid origin is likely due to slowing of concerted evolution caused by the long generation time of peonies. The partial and uneven homogenization of parental sequences displayed in nine species of putative hybrid origin may have resulted from gradients of gene conversion. The documented hybridizations may have occurred since the Pleistocene glaciations. The species of hybrid origin and their putative parents are now distantly allopatric. Reconstruction of reticulate evolution with sequence data, therefore, provides gene records for distributional histories of some of the parental species.

Rapid mass spectrometric peptide sequencing and mass matching for characterization of human melanoma proteins isolated by two-dimensional PAGE.

We report a general mass spectrometric approach for the rapid identification and characterization of proteins isolated by preparative two-dimensional polyacrylamide gel electrophoresis. This method possesses the inherent power to detect and structurally characterize covalent modifications. Absolute sensitivities of matrix-assisted laser desorption ionization and high-energy collision-induced dissociation tandem mass spectrometry are exploited to determine the mass and sequence of subpicomole sample quantities of tryptic peptides. These data permit mass matching and sequence homology searching of computerized peptide mass and protein sequence data bases for known proteins and design of oligonucleotide probes for cloning unknown proteins. We have identified 11 proteins in lysates of human A375 melanoma cells, including: alpha-enolase, cytokeratin, stathmin, protein disulfide isomerase, tropomyosin, Cu/Zn superoxide dismutase, nucleoside diphosphate kinase A, galaptin, and triosephosphate isomerase. We have characterized several posttranslational modifications and chemical modifications that may result from electrophoresis or subsequent sample processing steps. Detection of comigrating and covalently modified proteins illustrates the necessity of peptide sequencing and the advantages of tandem mass spectrometry to reliably and unambiguously establish the identity of each protein. This technology paves the way for studies of cell-type dependent gene expression and studies of large suites of cellular proteins with unprecedented speed and rigor to provide information complementary to the ongoing Human Genome Project.

Identification of cDNA clones encoding valosin-containing protein and other plant plasma membrane-associated proteins by a general immunoscreening strategy.

An approach was developed for the isolation and characterization of soybean plasma membrane-associated proteins by immunoscreening of a cDNA expression library. An antiserum was raised against purified plasma membrane vesicles. In a differential screening of approximately 500,000 plaque-forming units with the anti-(plasma membrane) serum and DNA probes derived from highly abundant clones isolated in a preliminary screening, 261 clones were selected from approximately 1,200 antiserum-positive plaques. These clones were classified into 40 groups by hybridization analysis and 5'- and 3'-terminal sequencing. By searching nucleic acid and protein sequence data bases, 11 groups of cDNAs were identified, among which valosin-containing protein (VCP), clathrin heavy chain, phospholipase C, and S-adenosylmethionine:delta 24-sterol-C-methyltransferase have not to date been cloned from plants. The remaining 29 groups did not match any current data base entries and may, therefore, represent additional or yet uncharacterized genes. A full-length cDNA encoding the soybean VCP was sequenced. The high level of amino acid identity with vertebrate VCP and yeast CDC48 protein indicates that the soybean protein is a plant homolog of vertebrate VCP and yeast CDC48 protein.

Genetic diversity and susceptibility to Vip3Aa20 protein in Brazilian populations of Helicoverpa armigera and Helicoverpa zea (Lepidoptera: Noctuidae)

Helicoverpa armigera (Hübner) was officially reported in Brazil in 2013. This species is closely related to Helicoverpa zea (Boddie) and has caused significant crop damage in Brazil. The use of genetically modified crops expressing insecticidal protein from Bacillus thuringiensis (Berliner) has been one of the control tactics for managing these pests. Genetically modified maize expressing Vip3Aa20 was approved to commercial use in Brazil in 2009. Understanding the genetic diversity and the susceptibility to B. thuringiensis proteins in H. armigera and H. zea populations in Brazil are crucial for establishing Insect Resistance Management (IRM) programs in Brazil. Therefore, the objectives of this study were: (a) to infer demographic parameters and genetic structure of H. armigera and H. zea Brazil; (b) to assess the intra and interspecific gene flow and genetic diversity of H. armigera and H. zea; and (c) to evaluate the susceptibility to Vip3Aa20 protein in H. armigera and H. zea populations of Brazil. A phylogeographic analysis of field H. armigera and H. zea populations was performed using a partial sequence data from the cytochrome c oxidase I (COI) gene. H. armigera individuals were most prevalent on dicotyledonous hosts and H. zea individuals were most prevalent on maize crops. Both species showed signs of demographic expansion and no genetic structure. High genetic diversity and wide distribution were observed for H. armigera. A joint analysis indicated the presence of Chinese, Indian, and European lineages within the Brazilian populations of H. armigera. In the cross-species amplification study, seven microsatellite loci were amplified; and showed a potential hybrid offspring in natural conditions. Interespecific analyses using the same microsatellite loci with Brazilian H. armigera and H. zea in compare to the USA H. zea were also conducted. When analyses were performed within each species, 10 microsatellites were used for H. armigera, and eight for H. zea. We detected high intraspecific gene flow in populations of H. armigera and H. zea from Brazil and H. zea from the USA. Genetic diversity was similar for both species. However, H. armigera was more similar to H. zea from Brazil than H. zea from the USA and some putative hybrid individuals were found in Brazilian populations.Tthere was low gene flow between Brazilian and USA H. zea. The baseline susceptibility to Vip3Aa20 resulted in low interpopulation variation for H. zea (3-fold) and for H. armigera (5-fold), based on LC50. H. armigera was more tolerant to Vip3Aa20 than H. zea (≈ 40 to 75-fold, based on CL50). The diagnostic concentration for susceptibility monitoring, based on CL99, was fairly high (6,400 ng Vip3Aa20/cm2) for H. zea and not validated for H. armigera due to the high amount of protein needed for bioassays. Implementing IRM strategies to Vip3Aa20 in H. armigera and H. zea will be of a great challenge in Brazil, mainly due to the low susceptibility to Vip3Aa20 and high genetic diversity and gene flow in both species, besides a potential of hybrid individuals between H. armigera and H. zea under field conditions.

Hemorheological alterations in sickle cell anemia and their clinical consequences: the role of genetic modulators

Sickle cell anemia (SCA) is an autosomal recessive chronic hemolytic anemia, caused by homozygosity for the HBB:c.20A>T mutation. The disease presents with high clinical heterogeneity, stroke being the most devastating manifestation. This study aimed to identify genetic modulators of severe hemolysis and stroke risk in children with SCA, as well as understand their consequences at the hemorheological level. Sixty-six children with SCA were categorised according to their degree of cerebral vasculopathy (Stroke/Risk/Control). Relevant data were collected from patients’ medical records. Several polymorphic regions in genes related to vascular cell adhesion and tonus were characterized by molecular methodologies. Data analyses were performed using R software. Several in silico tools (e.g. TFBind, MatInspector) were applied to investigate the main variant consequences. Some genetic variants in vascular adhesion molecule-1 gene promoter and endothelial nitric oxide synthase gene were associated with higher levels of hemolysis and stroke events. They modify important transcription factor binding sites or disturb the corresponding protein structure/function. Our findings emphasize the relevance of the genetic variants in modulating the degree of hemolysis and development of cerebral vasculopathy due to their effect on gene expression, modification of protein biological activities related with erythrocyte/endothelial interactions and consequent hemorheological abnormalities in SCA.

Genotypic Resistance Tests Sequences Reveal the Role of Marginalized Populations in HIV-1 Transmission in Switzerland.

Targeting hard-to-reach/marginalized populations is essential for preventing HIV-transmission. A unique opportunity to identify such populations in Switzerland is provided by a database of all genotypic-resistance-tests from Switzerland, including both sequences from the Swiss HIV Cohort Study (SHCS) and non-cohort sequences. A phylogenetic tree was built using 11,127 SHCS and 2,875 Swiss non-SHCS sequences. Demographics were imputed for non-SHCS patients using a phylogenetic proximity approach. Factors associated with non-cohort outbreaks were determined using logistic regression. Non-B subtype (univariable odds-ratio (OR): 1.9; 95% confidence interval (CI): 1.8-2.1), female gender (OR: 1.6; 95% CI: 1.4-1.7), black ethnicity (OR: 1.9; 95% CI: 1.7-2.1) and heterosexual transmission group (OR:1.8; 95% CI: 1.6-2.0), were all associated with underrepresentation in the SHCS. We found 344 purely non-SHCS transmission clusters, however, these outbreaks were small (median 2, maximum 7 patients) with a strong overlap with the SHCS'. 65% of non-SHCS sequences were part of clusters composed of >= 50% SHCS sequences. Our data suggests that marginalized-populations are underrepresented in the SHCS. However, the limited size of outbreaks among non-SHCS patients in-care implies that no major HIV outbreak in Switzerland was missed by the SHCS surveillance. This study demonstrates the potential of sequence data to assess and extend the scope of infectious-disease surveillance.

List of size fractionated eukaryotic plankton community samples and associated metadata (Database W1)

The present data set provides an Excel file in a zip archive. The file lists 334 samples of size fractionated eukaryotic plankton community with a suite of associated metadata (Database W1). Note that if most samples represented the piconano- (0.8-5 µm, 73 samples), nano- (5-20 µm, 74 samples), micro- (20-180 µm, 70 samples), and meso- (180-2000 µm, 76 samples) planktonic size fractions, some represented different organismal size-fractions: 0.2-3 µm (1 sample), 0.8-20 µm (6 samples), 0.8 µm - infinity (33 samples), and 3-20 µm (1 sample). The table contains the following fields: a unique sample sequence identifier; the sampling station identifier; the Tara Oceans sample identifier (TARA_xxxxxxxxxx); an INDSC accession number allowing to retrieve raw sequence data for the major nucleotide databases (short read archives at EBI, NCBI or DDBJ); the depth of sampling (Subsurface - SUR or Deep Chlorophyll Maximum - DCM); the targeted size range; the sequences template (either DNA or WGA/DNA if DNA extracted from the filters was Whole Genome Amplified); the latitude of the sampling event (decimal degrees); the longitude of the sampling event (decimal degrees); the time and date of the sampling event; the device used to collect the sample; the logsheet event corresponding to the sampling event ; the volume of water sampled (liters). Then follows information on the cleaning bioinformatics pipeline shown on Figure W2 of the supplementary litterature publication: the number of merged pairs present in the raw sequence file; the number of those sequences matching both primers; the number of sequences after quality-check filtering; the number of sequences after chimera removal; and finally the number of sequences after selecting only barcodes present in at least three copies in total and in at least two samples. Finally, are given for each sequence sample: the number of distinct sequences (metabarcodes); the number of OTUs; the average number of barcode per OTU; the Shannon diversity index based on barcodes for each sample (URL of W4 dataset in PANGAEA); and the Shannon diversity index based on each OTU (URL of W5 dataset in PANGAEA).

Understanding human genome regulation through entropic graph-based regularization and submodular optimization

Thesis (Ph.D.)--University of Washington, 2016-06

Mechanisms of Thermal Adaptation Revealed From the Genomes of the Antarctic Archaea Methanogenium frigidum and Methanococcoides burtonii

We generated draft genome sequences for two cold-adapted Archaea, Methanogenium frigidum and Methanococcoides burtonii, to identify genotypic characteristics that distinguish them from Archaea with a higher optimal growth temperature (OGT). Comparative genomics revealed trends in amino acid and tRNA composition, and structural features of proteins. Proteins from the cold-adapted Archaea are characterized by a higher content of noncharged polar amino acids, particularly Gin and Thr and a lower content of hydrophobic amino acids, particularly Leu. Sequence data from nine methanogen genomes (OGT 15degrees-98degreesC) were used to generate IIII modeled protein structures. Analysis of the models from the cold-adapted Archaea showed a strong tendency in the solvent-accessible area for more Gin, Thr, and hydrophobic residues and fewer charged residues. A cold shock domain (CSD) protein (CspA homolog) was identified in M. frigidum, two hypothetical proteins with CSD-folds in M. burtonii, and a unique winged helix DNA-binding domain protein in M. burtonii. This suggests that these types of nucleic acid binding proteins have a critical role in cold-adapted Archaea. Structural analysis of tRNA sequences from the Archaea indicated that GC content is the major factor influencing tRNA stability in hyperthermophiles, but not in the psychrophiles, mesophiles or moderate thermophiles. Below an OGT of 60degreesC, the GC content in tRNA was largely unchanged, indicating that any requirement for flexibility of tRNA in psychrophiles is mediated by other means. This is the first time that comparisons have been performed with genome data from Archaea spanning the growth temperature extremes. from psychrophiles to hyperthermophiles

Segregation distortion for seed testa color in mungbean (Vigna radiata L. Wilcek)

Genetic segregation experiments with plant species are commonly used for understanding the inheritance of traits. A basic assumption in these experiments is that each gamete developed from megasporogenesis has an equal chance of fusing with a gamete developed from microsporogenesis, and every zygote formed has an equal chance of survival. If gametic and/or zygotic selection occurs whereby certain gametes or zygotic combinations have a reduced chance of survival, progeny distributions are skewed and are said to exhibit segregation distortion. In this study, inheritance data are presented for the trait seed testa color segregating in large populations (more than 200 individuals) derived from closely related mungbean (Vigna radiata L. Wilcek) taxa. Segregation ratios suggested complex inheritance, including dominant and recessive epistasis. However, this genetic model was rejected in favor of a single-gene model based on evidence of segregation distortion provided by molecular marker data. The segregation distortion occurred after each generation of self-pollination from F-1 thru F-7 resulting in F-7 phenotypic frequencies of 151:56 instead of the expected 103.5:103.5. This study highlights the value of molecular markers for understanding the inheritance of a simply inherited trait influenced by segregation distortion.

First record of Passalora calotropidis in Australia and its generic position

Passalora calotropidis has been found for the first time in Australia on rubber bush (Calotropis procera) in northern Queensland where it was associated with a damaging leaf spot disease. Analysis of sequence data of the ITS region indicated that P. calotropidis belonged to a group that consisted of species of Pseudocercospora. The generic position of P. calotropidis and its potential for biological control are discussed.

Discovery of a new source of rifamycin antibiotics in marine sponge actinobacteria by phylogenetic prediction

Phylogenetic analysis of the ketosynthase (KS) gene sequences of marine sponge-derived Salinispora strains of actinobacteria indicated that the polyketide synthase (PKS) gene sequence most closely related to that of Salinispora was the rifamycin B synthase of Amycolatopsis mediterranei. This result was not expected from taxonomic species tree phylogenetics using 16S rRNA sequences. From the PKS sequence data generated from our sponge-derived Salinispora strains, we predicted that such strains might synthesize rifamycin-like compounds. Liquid chromatography-tandem mass spectrometry (LC/MS/MS) analysis was applied to one sponge-derived Salinispora strain to test the hypothesis of rifamycin synthesis. The analysis reported here demonstrates that this Salinispora isolate does produce compounds of the rifamycin class, including rifamycin B and rifamycin SV. A rifamycin-specific KS primer set was designed, and that primer set increased the number of rifamycin-positive strains detected by PCR screening relative to the number detectable using a conserved KS-specific set. Thus, the Salinispora group of actinobacteria represents a potential new source of rifamycins outside the genus Amycolatopsis and the first recorded source of rifamycins from marine bacteria.