Tracking Genomic Cancer Evolution for Precision Medicine: The Lung TRACERx Study

The importance of intratumour genetic and functional heterogeneity is increasingly recognised as a driver of cancer progression and survival outcome. Understanding how tumour clonal heterogeneity impacts upon therapeutic outcome, however, is still an area of unmet clinical and scientific need. TRACERx (TRAcking non-small cell lung Cancer Evolution through therapy [Rx]), a prospective study of patients with primary non-small cell lung cancer (NSCLC), aims to define the evolutionary trajectories of lung cancer in both space and time through multiregion and longitudinal tumour sampling and genetic analysis. By following cancers from diagnosis to relapse, tracking the evolutionary trajectories of tumours in relation to therapeutic interventions, and determining the impact of clonal heterogeneity on clinical outcomes, TRACERx may help to identify novel therapeutic targets for NSCLC and may also serve as a model applicable to other cancer types.     

Diversity and Evolution of the Marsupial Mandibular Angular Process

A medial inflection of the mandibular angular process is present in most marsupials. The few living marsupials that lack this trait either are very specialized forms (e.g., Tarsipes) or show a medial inflection at some point in development that is lost in later ontogenetic stages (cf. Dactylopsila and Phascolarctos). A medially inflected angular process is not present in any known extant or extinct placental (including all Cretaceous taxa that preserve the back of the dentary bone). Some extant placentals with enlarged auditory bullae evolved a medial flange of the angular process as a strategy to increase gape, but this is not homologous to the marsupial condition. We conclude that the medially inflected angular process is a shared derived trait of extant and extinct marsupials. The significant diversity in the form of the medially inflected mandibular angular process in marsupials, documented here for 53 taxa, shows a general relation to dietary adaptations. Herbivores (with well-developed masseter and medial pterygoid muscles) tend to have a shelf-like angular process, while small, insectivorous marsupials generally have a rod-like angular process. A close connection between the angular process and the ectotympanic is maintained during early postnatal development in all marsupials examined, a relation not seen in the placentals examined. A previous hypothesis suggested that the angular process plays a role in hearing in pouch-young Monodelphis. Data on the maturation of the auditory system does not support this hypothesis. Currently there are no data on differences in muscular anatomy or mastication between marsupials and placentals that could serve as a causal explanation for the difference in adult form of the angular process between the two groups.     

The Evolution of Genomic Imprinting

In some mammalian genes, the paternally and maternally derived alleles are expressed differently: this phenomenon is called genomic imprinting. Here we study the evolution of imprinting using multivariate quantitative genetic models to examine the feasibility of the genetic conflict hypothesis. This hypothesis explains the observed imprinting patterns as an evolutionary outcome of the conflict between the paternal and maternal alleles. We consider the expression of a zygotic gene, which codes for an embryonic growth factor affecting the amount of maternal resources obtained through the placenta. We assume that the gene produces the growth factor in two different amounts depending on its parental origin. We show that genomic imprinting evolves easily if females have some probability of multiple partners. This is in conflict with the observation that not all genes controlling placental development are imprinted and that imprinting in some genes is not conserved between mice and humans. We show however that deleterious mutations in the coding region of the gene create selection against imprinting.     

The Effect of Transposable Element Insertions on Gene Expression Evolution in Rodents

Background

Many genomes contain a substantial number of transposable elements (TEs), a few of which are known to be involved in regulating gene expression. However, recent observations suggest that TEs may have played a very important role in the evolution of gene expression because many conserved non-genic sequences, some of which are know to be involved in gene regulation, resemble TEs.

Results

Here we investigate whether new TE insertions affect gene expression profiles by testing whether gene expression divergence between mouse and rat is correlated to the numbers of new transposable elements inserted near genes. We show that expression divergence is significantly correlated to the number of new LTR and SINE elements, but not to the numbers of LINEs. We also show that expression divergence is not significantly correlated to the numbers of ancestral TEs in most cases, which suggests that the correlations between expression divergence and the numbers of new TEs are causal in nature. We quantify the effect and estimate that TE insertion has accounted for ∼20% (95% confidence interval: 12% to 26%) of all expression profile divergence in rodents.

Conclusions

We conclude that TE insertions may have had a major impact on the evolution of gene expression levels in rodents.     

Genomic Signatures of Selective Pressures and Introgression from Archaic Hominins at Human Innate Immunity Genes

Human genes governing innate immunity provide a valuable tool for the study of the selective pressure imposed by microorganisms on host genomes. A comprehensive, genome-wide study of how selective constraints and adaptations have driven the evolution of innate immunity genes is missing. Using full-genome sequence variation from the 1000 Genomes Project, we first show that innate immunity genes have globally evolved under stronger purifying selection than the remainder of protein-coding genes. We identify a gene set under the strongest selective constraints, mutations in which are likely to predispose individuals to life-threatening disease, as illustrated by STAT1 and TRAF3. We then evaluate the occurrence of local adaptation and detect 57 high-scoring signals of positive selection at innate immunity genes, variation in which has been associated with susceptibility to common infectious or autoimmune diseases. Furthermore, we show that most adaptations targeting coding variation have occurred in the last 6,000–13,000 years, the period at which populations shifted from hunting and gathering to farming. Finally, we show that innate immunity genes present higher Neandertal introgression than the remainder of the coding genome. Notably, among the genes presenting the highest Neandertal ancestry, we find the TLR6-TLR1-TLR10 cluster, which also contains functional adaptive variation in Europeans. This study identifies highly constrained genes that fulfill essential, non-redundant functions in host survival and reveals others that are more permissive to change—containing variation acquired from archaic hominins or adaptive variants in specific populations—improving our understanding of the relative biological importance of innate immunity pathways in natural conditions.     

A Phylogeny and Timescale for Marsupial Evolution Based on Sequences for Five Nuclear Genes

Even though marsupials are taxonomically less diverse than placentals, they exhibit comparable morphological and ecological diversity. However, much of their fossil record is thought to be missing, particularly for the Australasian groups. The more than 330 living species of marsupials are grouped into three American (Didelphimorphia, Microbiotheria, and Paucituberculata) and four Australasian (Dasyuromorphia, Diprotodontia, Notoryctemorphia, and Peramelemorphia) orders. Interordinal relationships have been investigated using a wide range of methods that have often yielded contradictory results. Much of the controversy has focused on the placement of Dromiciops gliroides (Microbiotheria). Studies either support a sister-taxon relationship to a monophyletic Australasian clade or a nested position within the Australasian radiation. Familial relationships within the Diprotodontia have also proved difficult to resolve. Here, we examine higher-level marsupial relationships using a nuclear multigene molecular data set representing all living orders. Protein-coding portions of ApoB, BRCA1, IRBP, Rag1, and vWF were analyzed using maximum parsimony, maximum likelihood, and Bayesian methods. Two different Bayesian relaxed molecular clock methods were employed to construct a timescale for marsupial evolution and estimate the unrepresented basal branch length (UBBL). Maximum likelihood and Bayesian results suggest that the root of the marsupial tree is between Didelphimorphia and all other marsupials. All methods provide strong support for the monophyly of Australidelphia. Within Australidelphia, Dromiciops is the sister-taxon to a monophyletic Australasian clade. Within the Australasian clade, Diprotodontia is the sister taxon to a Notoryctemorphia + Dasyuromorphia + Peramelemorphia clade. Within the Diprotodontia, Vombatiformes (wombat + koala) is the sister taxon to a paraphyletic possum group (Phalangeriformes) with kangaroos nested inside. Molecular dating analyses suggest Late Cretaceous/Paleocene dates for all interordinal divergences. All intraordinal divergences were placed in the mid to late Cenozoic except for the deepest splits within the Diprotodontia. Our UBBL estimates of the marsupial fossil record indicate that the South American record is approximately as complete as the Australasian record. The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Systematics and Evolution of the Dasyurid Marsupial Genus Sminthopsis: I. The Macroura Species Group

Genetic variation within the macroura species group, which includes Sminthopsis macroura, S. virginiae, S. douglasi, and S. bindi, was examined through analyses of complete mitochondrial 12S rRNA gene sequences, partial control-region DNA sequences, and allozymes. Divergent genetic lineages appear to be present within S. macroura and S. virginiae, and it is likely that this genetic divergence equates to currently unrecognized taxonomic diversity. Specimens of S. macroura (as currently recognized) belong to three genetically distinct lineages that are highly divergent from one another. Two of these lineages may be synonymous with two previously recognized dunnart species—S. froggatti and S. stalkeri. The third appears to represent "true" S. macroura and is itself genetically heterogeneous, with a number of subgroups present within it that may also represent currently unrecognized taxa. The mitochondrial DNA sequence divergences observed between S. virginiae nitela and the two other S. virginiaesubspecies are equivalent to, or greater than, those noted between other dunnart species. Allozyme divergences between these subspecies were however slightly lower, and determination on whether S. virginiae nitela should be returned to full species status (S. nitela) may require further evidence. Phylogenetic relationships between species in the macroura group appear to have been partially resolved, with individual 12S rRNA and combined mitochondrial DNA analyses recovering S. bindi as the earliest diverging taxon. Other relationships between species in the group were either not consistently recovered or lacked strong support.     

Genomic Pathogen Typing Using Solid-State Nanopores

In clinical settings, rapid and accurate characterization of pathogens is essential for effective treatment of patients; however, subtle genetic changes in pathogens which elude traditional phenotypic typing may confer dangerous pathogenic properties such as toxicity, antibiotic resistance, or virulence. Existing options for molecular typing techniques characterize the critical genomic changes that distinguish harmful and benign strains, yet the well-established approaches, in particular those that rely on electrophoretic separation of nucleic acid fragments on a gel, have room for only incremental future improvements in speed, cost, and complexity. Solid-state nanopores are an emerging class of single-molecule sensors that can electrophoretically characterize charged biopolymers, and which offer significant advantages in terms of sample and reagent requirements, readout speed, parallelization, and automation. We present here the first application of nanopores for single-molecule molecular typing using length based “fingerprints” of critical sites in bacterial genomes. This technique is highly adaptable for detection of different types of genetic variation; as we illustrate using prototypical examples including Mycobacterium tuberculosis and methicillin-resistant Streptococcus aureus, the solid-state nanopore diagnostic platform may be used to detect large insertions or deletions, small insertions or deletions, and even single-nucleotide variations in bacterial DNA. We further show that Bayesian classification of test samples can provide highly confident pathogen typing results based on only a few tens of independent single-molecule events, making this method extremely sensitive and statistically robust.     

Genomic Fossils Calibrate the Long-Term Evolution of Hepadnaviruses

Because most extant viruses mutate rapidly and lack a true fossil record, their deep evolution and long-term substitution rates remain poorly understood. In addition to retroviruses, which rely on chromosomal integration for their replication, many other viruses replicate in the nucleus of their host''s cells and are therefore prone to endogenization, a process that involves integration of viral DNA into the host''s germline genome followed by long-term vertical inheritance. Such endogenous viruses are highly valuable as they provide a molecular fossil record of past viral invasions, which may be used to decipher the origins and long-term evolutionary characteristics of modern pathogenic viruses. Hepadnaviruses (Hepadnaviridae) are a family of small, partially double-stranded DNA viruses that include hepatitis B viruses. Here we report the discovery of endogenous hepadnaviruses in the genome of the zebra finch. We used a combination of cross-species analysis of orthologous insertions, molecular dating, and phylogenetic analyses to demonstrate that hepadnaviruses infiltrated repeatedly the germline genome of passerine birds. We provide evidence that some of the avian hepadnavirus integration events are at least 19 My old, which reveals a much deeper ancestry of Hepadnaviridae than could be inferred based on the coalescence times of modern hepadnaviruses. Furthermore, the remarkable sequence similarity between endogenous and extant avian hepadnaviruses (up to 75% identity) suggests that long-term substitution rates for these viruses are on the order of 10⁻⁸ substitutions per site per year, which is a 1,000-fold slower than short-term rates estimated based on the sequences of circulating hepadnaviruses. Together, these results imply a drastic shift in our understanding of the time scale of hepadnavirus evolution, and suggest that the rapid evolutionary dynamics characterizing modern avian hepadnaviruses do not reflect their mode of evolution on a deep time scale.     

Genomic Selection Using Low-Density Marker Panels

Genomic selection (GS) using high-density single-nucleotide polymorphisms (SNPs) is promising to improve response to selection in populations that are under artificial selection. High-density SNP genotyping of all selection candidates each generation, however, may not be cost effective. Smaller panels with SNPs that show strong associations with phenotype can be used, but this may require separate SNPs for each trait and each population. As an alternative, we propose to use a panel of evenly spaced low-density SNPs across the genome to estimate genome-assisted breeding values of selection candidates in pedigreed populations. The principle of this approach is to utilize cosegregation information from low-density SNPs to track effects of high-density SNP alleles within families. Simulations were used to analyze the loss of accuracy of estimated breeding values from using evenly spaced and selected SNP panels compared to using all high-density SNPs in a Bayesian analysis. Forward stepwise selection and a Bayesian approach were used to select SNPs. Loss of accuracy was nearly independent of the number of simulated quantitative trait loci (QTL) with evenly spaced SNPs, but increased with number of QTL for the selected SNP panels. Loss of accuracy with evenly spaced SNPs increased steadily over generations but was constant when the smaller number individuals that are selected for breeding each generation were also genotyped using the high-density SNP panel. With equal numbers of low-density SNPs, panels with SNPs selected on the basis of the Bayesian approach had the smallest loss in accuracy for a single trait, but a panel with evenly spaced SNPs at 10 cM was only slightly worse, whereas a panel with SNPs selected by forward stepwise selection was inferior. Panels with evenly spaced SNPs can, however, be used across traits and populations and their performance is independent of the number of QTL affecting the trait and of the methods used to estimate effects in the training data and are, therefore, preferred for broad applications in pedigreed populations under artificial selection.     

Evolution of Genomic Structures on Mammalian Sex Chromosomes

Throughout mammalian evolution, recombination between the two sex chromosomes was suppressed in a stepwise manner. It is thought that the suppression of recombination led to an accumulation of deleterious mutations and frequent genomic rearrangements on the Y chromosome. In this article, we review three evolutionary aspects related to genomic rearrangements and structures, such as inverted repeats (IRs) and palindromes (PDs), on the mammalian sex chromosomes. First, we describe the stepwise manner in which recombination between the X and Y chromosomes was suppressed in placental mammals and discuss a genomic rearrangement that might have led to the formation of present pseudoautosomal boundaries (PAB). Second, we describe ectopic gene conversion between the X and Y chromosomes, and propose possible molecular causes. Third, we focus on the evolutionary mode and timing of PD formation on the X and Y chromosomes. The sequence of the chimpanzee Y chromosome was recently published by two groups. Both groups suggest that rapid evolution of genomic structure occurred on the Y chromosome. Our re-analysis of the sequences confirmed the species-specific mode of human and chimpanzee Y chromosomal evolution. Finally, we present a general outlook regarding the rapid evolution of mammalian sex chromosomes.     

An Inverse Pcr Screen for the Detection of P Element Insertions in Cloned Genomic Intervals in Drosophila Melanogaster

We developed a screening approach that utilizes an inverse polymerase chain reaction (PCR) to detect P element insertions in or near previously cloned genes in Drosophila melanogaster. We used this approach in a large scale genetic screen in which P elements were mobilized from sites on the X chromosome to new autosomal locations. Mutagenized flies were combined in pools, and our screening approach was used to generate probes corresponding to the sequences flanking each site of insertion. These probes then were used for hybridization to cloned genomic intervals, allowing individuals carrying insertions in them to be detected. We used the same approach to perform repeated rounds of sib-selection to generate stable insertion lines. We screened 16,100 insert bearing individuals and recovered 11 insertions in five intervals containing genes encoding members of the kinesin superfamily in Drosophila melanogaster. In addition, we recovered an insertion in the region including the Larval Serum Protein-2 gene. Examination by Southern hybridization confirms that the lines we recovered represent genuine insertions in the corresponding genomic intervals. Our data indicates that this approach will be very efficient both for P element mutagenesis of new genomic regions and for detection and recovery of ``local' P element transposition events. In addition, our data constitutes a survey of preferred P element insertion sites in the Drosophila genome and suggests that insertion sites that are mutable at a rate of ~10(-4) are distributed every 40-50 kb.     

Divergence of Noncoding Sequences and of Insertions Encoding Nonglobular Domains at a Genomic Region Well Conserved in Plasmodia

To identify conserved features in the rapidly diverging portions of a well-conserved locus, completely sequenced in Plasmodium falciparum and Plasmodium berghei, a computational method based on recurrence analysis was exploited. At the level of the genomic sequence, in both species, introns and intergenic sequences—though subject to rapid diversification—do not drift without constraints, but rather coevolve, in the sense that they maintain not only an AT-rich base composition, but also a consistent use of recurring (AT) _n tracts. One of the two genes present in the conserved locus encodes a protein that exhibits blocks of high similarity to the first enzyme in glutathione biosynthesis (γ-glutamylcysteine synthetase) but bears long low-complexity insertions, absent in other organisms. From an analysis of the aminoacid sequence, different constraints appear to act on the borders and on the central part of the insertions. Albeit maintaining a strong bias toward hydrophylic residues, central portions diverge more rapidly than borders, through point mutation and differential presence of entire tracts. Received: 20 September 1999 / Accepted: 9 February 2000     

Systematics and Evolution of the Dasyurid Marsupial Genus Sminthopsis: II. The Murina Species Group

Genetic variation within the Murina species group, which includes S. murina, S. gilberti, S. leucopus, S. dolichura and S. archeri, was examined through analyses of complete 12S rRNA, partial control region mitochondrial DNA sequences and partial omega-globin nuclear DNA sequences. Sminthopsis butleri was found to be an additional member of the Murina group, and appears to be most closely related to S. leucopus rather than the morphologically similar S. archeri. This latter species appears to be the most divergent member of the group, and there is a possible sister relationship between S. murina and S. gilberti, as suggested by previous allozyme evidence. It appears that the systematic affinities of the taxonomically problematic northeastern Queensland populations of both S. murina and S. leucopus and a disjunct population of S. gilberti (from the Western Australia/South Australia border) are indeed with those respective species; although each appears to belong to a distinct morphological and genetic lineage. A specimen of S. leucopus from Queensland was found to be as divergent from each of the southeastern Australian S. leucopus subspecies as they are from each other, suggesting that this northern population of S. leucopus may also warrant recognition as a distinct taxon. Specimens of S. murina murina were found to be genetically divergent from each other, and this subspecies appears to be paraphyletic, as suggested by previous morphological evidence. ^* This paper is the second part of a series dealing with the systematics and evolution of the dasyurid marsupial genus Sminthopsis. Part one covered the Macroura species group and was published in 2001 in the Journal of Mammalian Evolution 8, 149–170.     

Genomic Stability of Murine Leukemia Viruses Containing Insertions at the Env-3′ Untranslated Region Boundary

Retroviruses containing inserts of exogenous sequences frequently eliminate the inserted sequences upon spread in susceptible cells. We have constructed replication-competent murine leukemia virus (MLV) vectors containing internal ribosome entry site (IRES)-transgene cassettes at the env-3' untranslated region boundary in order to examine the effects of insert sequence and size on the loss of inserts during viral replication. A virus containing an insertion of 1.6 kb replicated with greatly attenuated kinetics relative to wild-type virus and lost the inserted sequences in a single infection cycle. In contrast, MLVs containing inserts of 1.15 to 1.30 kb replicated with kinetics only slightly attenuated compared to wild-type MLV and exhibited much greater stability, maintaining their genomic integrity over multiple serial infection cycles. Eventually, multiple species of deletion mutants were detected simultaneously in later infection cycles; once detected, these variants rapidly dominated the population and thereafter appeared to be maintained at a relative equilibrium. Sequence analysis of these variants identified preferred sites of recombination in the parental viruses, including both short direct repeats and inverted repeats. One instance of insert deletion through recombination with an endogenous retrovirus was also observed. When specific sequences involved in these recombination events were eliminated, deletion variants still arose with the same kinetics upon virus passage and by apparently similar mechanisms, although at different locations in the vectors. Our results suggest that while lengthened, insert-containing genomes can be maintained over multiple replication cycles, preferential deletions resulting in loss of the inserted sequences confer a strong selective advantage.     

Tracking Auxin Receptors Using Functional Approaches

Functional approaches toward the identification of auxin receptors developed along two major lines: the isolation and characterization of mutants or transgenic plants affected in their responses to the hormone and the study of early auxin effects at the cell level such as expression of specific genes or modifications of plasma membrane properties. The combination of these approaches with those aiming at the molecular characterization of auxin binding proteins as putative auxin receptors allowed to bring further insight into the mechanisms of auxin perception by plant cells. Studies of membrane responses to auxin clearly demonstrated the existence of elementary response chains to auxin at the plasma membrane, the activation of auxin responsive proteins leading to changes in the membrane potential via the stimulation of the proton pump ATPase or the modulation of ion channels. A two-component model is proposed for the organization of functional auxin perception units at the plasma membrane, comprising an auxin-binding moiety related to the major auxin-binding protein from maize (ZmER-abp1), associated to a transmembrane protein. Current research investigates the relevance of this model and tries to assess whether early responses at the plasma membrane share common perception or transduction steps with gene expression responses and participate in more integrated biological responses to auxin.