In vitro synthesis of oncogenic human papillomaviruses requires episomal genomes for differentiation-dependent late expression.

Human papillomavirus (HPV) types 16, 18, 31, and 51 are the etiologic agents of many anogenital cancers including those of the cervix. These "high risk" HPVs specifically target genital squamous epithelia, and their lytic life cycle is closely linked to epithelial differentiation. We have developed a genetic assay for HPV functions during pathogenesis using recircularized cloned HPV 31 genomes that were transfected together with a drug resistance marker into monolayer cultures of normal human foreskin keratinocytes, the natural host cell. After drug selection, cell lines were isolated that stably maintained HPV 31 DNA as episomes and underwent terminal differentiation when grown in organotypic raft cultures. In differentiated rafts, the expression of late viral genes, amplification of viral DNA, and production of viral particles were detected in suprabasal cells. This demonstrated the ability to synthesize HPV 31 virions from transfected DNA templates and allowed an examination of HPV functions during the vegetative viral life cycle. We then used this system to investigate whether an episomal genome was required for the induction of late viral gene expression. When an HPV 31 genome (31E1*) containing a missense mutation in the E1 open reading frame was transfected into normal human keratinocytes, the mutant viral sequences were found to integrate into the host cell chromosomal DNA with both early and late regions intact. While high levels of early viral gene transcription were observed, no late gene expression was detected in rafts of cell lines containing the mutant viral genome despite evidence of terminal differentiation. Therefore, the induction of late viral gene expression required that the viral genomes be maintained as extrachromosomal elements, and terminal differentiation alone was not sufficient. These studies provide the basis for a detailed examination of HPV functions during viral pathogenesis.

Polymorphic simple sequence repeat regions in chloroplast genomes: applications to the population genetics of pines.

Simple sequence repeats (SSRs), consisting of tandemly repeated multiple copies of mono-, di-, tri-, or tetranucleotide motifs, are ubiquitous in eukaryotic genomes and are frequently used as genetic markers, taking advantage of their length polymorphism. We have examined the polymorphism of such sequences in the chloroplast genomes of plants, by using a PCR-based assay. GenBank searches identified the presence of several (dA)n.(dT)n mononucleotide stretches in chloroplast genomes. A chloroplast (cp) SSR was identified in three pine species (Pinus contorta, Pinus sylvestris, and Pinus thunbergii) 312 bp upstream of the psbA gene. DNA amplification of this repeated region from 11 pine species identified nine length variants. The polymorphic amplified fragments were isolated and the DNA sequence was determined, confirming that the length polymorphism was caused by variation in the length of the repeated region. In the pines, the chloroplast genome is transmitted through pollen and this PCR assay may be used to monitor gene flow in this genus. Analysis of 305 individuals from seven populations of Pinus leucodermis Ant. revealed the presence of four variants with intrapopulational diversities ranging from 0.000 to 0.629 and an average of 0.320. Restriction fragment length polymorphism analysis of cpDNA on the same populations previously failed to detect any variation. Population subdivision based on cpSSR was higher (Gst = 0.22, where Gst is coefficient of gene differentiation) than that revealed in a previous isozyme study (Gst = 0.05). We anticipate that SSR loci within the chloroplast genome should provide a highly informative assay for the analysis of the genetic structure of plant populations.

Different cellular backgrounds confer a marked advantage to either mutant or wild-type mitochondrial genomes.

After the introduction of mitochondria with a mixture of mutant and wild-type mitochondrial DNA (mtDNA) into a human rho degree cell line (143B.206), Yoneda et al. [Yoneda, M., Chomyn, A., Martinuzzi, A., Hurko, O. & Attardi, G. (1992) Proc. Natl. Acad. Sci. USA 89, 11164-11168] observed a shift in the proportion of the two mitochondrial genotypes in a number of cybrid clones. In every case where a shift was observed, there was an increase in the proportion of mutant mtDNA. By using the same cell line (143B.206 rho degree), we also generated cybrids that were either stable in their mitochondrial genotype or showed an increase in the proportion of mutant mtDNA. However, temporal analysis of the same mutant mtDNA type in another rho degree cell line revealed a quite distinct outcome. Those clones that showed a change shifted toward higher levels of wild-type rather than mutant mtDNA. These results indicate that the nuclear genetic background of the recipient (rho degree) cell can influence the segregation of mutant and wild-type mitochondrial genomes in cell cybrids.

Identification of two flavivirus-like genomes in the GB hepatitis agent.

A subtractive PCR methodology known as representational difference analysis was used to clone specific nucleotide sequences present in the infectious plasma from a tamarin infected with the GB hepatitis agent. Eleven unique clones were identified, seven of which were examined extensively. All seven clones appeared to be derived from sequences exogenous to the genomes of humans, tamarins, Saccharomyces cerevisiae, and Escherichia coli. In addition, sequences from these clones were not detected in plasma or liver tissue of tamarins prior to their inoculation with the GB agent. These sequences were detected by reverse transcription-PCR in acute-phase plasma of tamarins inoculated with the GB agent. Probes derived from two of the seven clones detected an RNA species of > or = 8.3 kb in the liver of a GB-agent-infected tamarin by Northern blot hybridization. Sequence analysis indicated that five of the seven clones encode polypeptides that possess limited amino acid identity with the nonstructural proteins of hepatitis C virus. Extension of the sequences found in the seven clones revealed that plasma from an infected tamarin contained two RNA molecules > 9 kb long. Limited sequence identity with various isolates of hepatitis C virus and the relative positions of putative RNA helicases and RNA-dependent RNA polymerases in the predicted protein products of these molecules suggested that the GB agent contains two unique flavivirus-like genomes.

Insight into the evolution of the Solanaceae from the parental genomes of Petunia hybrida

Petunia hybrida is a popular bedding plant that has a long history as a genetic model system. We report the whole-genome sequencing and assembly of inbred derivatives of its two wild parents, P. axillaris N and P. inflata S6. The assemblies include 91.3% and 90.2% coverage of their diploid genomes (1.4 Gb; 2n = 14) containing 32,928 and 36,697 protein-coding genes, respectively. The genomes reveal that the Petunia lineage has experienced at least two rounds of hexaploidization: the older gamma event, which is shared with most Eudicots, and a more recent Solanaceae event that is shared with tomato and other solanaceous species. Transcription factors involved in the shift from bee to moth pollination reside in particularly dynamic regions of the genome, which may have been key to the remarkable diversity of floral colour patterns and pollination systems. The high-quality genome sequences will enhance the value of Petunia as a model system for research on unique biological phenomena such as small RNAs, symbiosis, self-incompatibility and circadian rhythms.

Certain observations concerning the effects of epistasis on complex traits and the evolution of genomes.

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

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

The mitochondrial genomes of soft ticks have an arrangement of genes that has remained unchanged for over 400 million years

There are two major groups of ticks: soft ticks and hard ticks. The hard ticks comprise the prostriate ticks and the metastriate ticks. The mitochondrial (mt) genomes of one species of prostriate tick and two species of metastriate ticks had been sequenced prior to our study. The prostriate tick has the ancestral arrangement of mt genes of arthropods, whereas the two metastriate ticks have rearrangements of eight genes and duplicate control regions. However, the arrangement of genes in the mt genomes of soft ticks had not been studied. We sequenced the mt genomes of two species of soft ticks, Carios capensis and Ornithodoros moubata, and a metastriate tick, Haemaphysalis flava. We found that the soft ticks have the ancestral arrangement of mt genes of arthropods, whereas the metastriate tick, H. flava, shares the rearrangements of mt genes and duplicate control regions with the other two metastriate ticks that have previously been studied. Our study indicates that gene rearrangements and duplicate control regions in mt genomes occurred once in the most recent common ancestor of metastriate ticks, whereas the ancestral arrangement of arthropods has remained unchanged for over 400 million years in the lineages leading to the soft ticks and the prostriate ticks.

Evolution of duplicate control regions in the mitochondrial genomes of metazoa: A case study with Australasian Ixodes ticks

To investigate the evolution pattern and phylogenetic utility of duplicate control regions (CRs) in mitochondrial (mt) genomes, we sequenced the entire mt genomes of three Ixodes species and part of the mt genomes of another I I species. All the species from the Australasian lineage have duplicate CRs, whereas the other species have one CR. Sequence analyses indicate that the two CRs of the Australasian Ixodes ticks have evolved in concert in each species. In addition to the Australasian Ixodes ticks, species from seven other lineages of metazoa also have mt genomes with duplicate CRs. Accumulated mtDNA sequence data from these metazoans and two recent experiments on replication of mt genomes in human cell lines with duplicate CRs allowed us to re-examine four intriguing questions about the presence of duplicate CRs in the mt genomes of metazoa: (1) Why do some mt genomes, but not others, have duplicate CRs? (2) How did mt genomes with duplicate CRs evolve? (3) How could the nucleotide sequences of duplicate CRs remain identical or very similar over evolutionary time? (4) Are duplicate CRs phylogenetic markers? It appears that mt genomes with duplicate CRs have a selective advantage in replication over mt genomes with one CR. Tandem duplication followed by deletion of genes is the most plausible mechanism for the generation of mt genomes with duplicate CRs. Once duplicate CRs occur in an mt genome, they tend to evolve in concert, probably by gene conversion. However, there are lineages where gene conversion may not always occur, and, thus, the two CRs may evolve independently in these lineages. Duplicate CRs have much potential as phylogenetic markers at low taxonomic levels, such as within genera, within families, or among families, but not at high taxonomic levels, such as among orders.

Ultraconserved elements in insect genomes: A highly conserved intronic sequence implicated in the control of homothorax mRNA splicing

Recently, we identified a large number of ultraconserved (uc) sequences in noncoding regions of human, mouse, and rat genomes that appear to be essential for vertebrate and amniote ontogeny. Here, we used similar methods to identify ultraconserved genomic regions between the insect species Drosophila melanogaster and Drosophila pseudoobscura, as well as the more distantly related Anopheles gambiae. As with vertebrates, ultraconserved sequences in insects appear to Occur primarily in intergenic and intronic sequences, and at intron-exon junctions. The sequences are significantly associated with genes encoding developmental regulators and transcription factors, but are less frequent and are smaller in size than in vertebrates. The longest identical, nongapped orthologous match between the three genomes was found within the homothorax (hth) gene. This sequence spans an internal exon-intron junction, with the majority located within the intron, and is predicted to form a highly stable stem-loop RNA structure. Real-time quantitative PCR analysis of different hth splice isoforms and Northern blotting showed that the conserved element is associated with a high incidence of intron retention in hth pre-mRNA, suggesting that the conserved intronic element is critically important in the post-transcriptional regulation of hth expression in Diptera.

Extraordinary number of gene rearrangements in the mitochondrial genomes of lice (Phthiraptera : Insecta)

The arrangement of genes in the mitochondrial (mt) genomes of most insects is the same, or near-identical, to that inferred to be ancestral for insects. We sequenced the entire mt genome of the small pigeon louse, Campanulotes bidentatus compar, and part of the mt genomes of nine other species of lice. These species were from six families and the three main suborders of the order Phthiraptera. There was no variation in gene arrangement among species within a family but there was much variation in gene arrangement among the three suborders of lice. There has been an extraordinary number of gene rearrangements in the mitochondrial genomes of lice!

Segmenting eukaryotic genomes with the generalized Gibbs sampler

Eukaryotic genomes display segmental patterns of variation in various properties, including GC content and degree of evolutionary conservation. DNA segmentation algorithms are aimed at identifying statistically significant boundaries between such segments. Such algorithms may provide a means of discovering new classes of functional elements in eukaryotic genomes. This paper presents a model and an algorithm for Bayesian DNA segmentation and considers the feasibility of using it to segment whole eukaryotic genomes. The algorithm is tested on a range of simulated and real DNA sequences, and the following conclusions are drawn. Firstly, the algorithm correctly identifies non-segmented sequence, and can thus be used to reject the null hypothesis of uniformity in the property of interest. Secondly, estimates of the number and locations of change-points produced by the algorithm are robust to variations in algorithm parameters and initial starting conditions and correspond to real features in the data. Thirdly, the algorithm is successfully used to segment human chromosome 1 according to GC content, thus demonstrating the feasibility of Bayesian segmentation of eukaryotic genomes. The software described in this paper is available from the author's website (www.uq.edu.au/similar to uqjkeith/) or upon request to the author.

Transposon-free regions in mammalian genomes

Despite the presence of over 3 million transposons separated on average by similar to 500 bp, the human and mouse genomes each contain almost 1000 transposon-free regions (TFRs) over 10 kb in length. The majority of human TFRs correlate with orthologous TFRs in the mouse, despite the fact that most transposons are lineage specific. Many human TFRs also overlap with orthologous TFRs in the marsupial opossum, indicating that these regions have remained refractory to transposon insertion for long evolutionary periods. Over 90% of the bases covered by TFRs are noncoding, much of which is not highly conserved. Most TFRs are not associated with unusual nucleotide composition, but are significantly associated with genes encoding developmental regulators, suggesting that they represent extended regions of regulatory information that are largely unable to tolerate insertions, a conclusion difficult to reconcile with current conceptions of gene regulation.