James Tour on Orphan Genes

(S. Joshua Swamidass) #21

@evograd is working on it now. When we have it concisely and clearly explained, I’ll pass it on to him. We will show him the data. We will see what he thinks of it.

(John Harshman) #22

I have to say that’s a really bad statement of neutral drift, though it could perhaps be just bad writing rather than misunderstanding.

(S. Joshua Swamidass) #23

He is a chemist. That is pretty good for a chemist. (Though @jordan is making some serious headway)

(John Harshman) #24

First, that’s insulting to chemists. Second, he’s a chemist presuming to judge evolutionary biology, which gives him a responsibility to get it right.

(S. Joshua Swamidass) #25

Let’s chalk it up to bad writing.

In the past, too, he has not “judged” evolutionary biology. He has stated, accurately, that it made no sense to him. That was true. No one had, I found out, actually sat down and explained it to him.

(Jordan Mantha) #26

Well, Tour is a synthetic organic chemist so it’s probably not that much of an insult :slight_smile: I have higher expectations of biochemists :wink:

I really do get the concern about chemists. I feel like we’re kind of left out of the discussion so much of the time. It feels like the physicists and biologists get the monopoly on the “big questions” so it’s natural for those of us interested in big questions to want to be a part of the conversation. Tour has done a lot of very good and important work in nanotechnology, but I do think scientists have to be careful to not speak confidently out of our discipline.

I’m just trying to learn from the experts at this point, I want to know more about what it means to be human.

(S. Joshua Swamidass) split this topic #27

2 posts were merged into an existing topic: Comments on James Tour on Orphan Genes

(Blogging Graduate Student) #29

Ok, so I’ve taken a look at the data from Ruiz-Orera et al. (2015):

They identified 634 candidate de novo genes in the human genome, based on the fact that they found 1,029 transcripts in humans that weren’t present in chimpanzees.
They link to a GTF file containing the information about those human-specific transcripts here: http://dx.doi.org/10.6084/m9.figshare.1604892
but note that this file contains both the human-specific transcripts found in humans, as well as the hominoid-specific transcripts found in humans. This wasn’t apparent to me at first, and I think @roohif missed it too, as he included the entire file in his analysis. I separated out the species-specific transcripts only (1,029).

A second complication is that the GTF file contains separate entries for different exons and CDSs in each transcript, so while there are only 1,029 transcripts (1,029 transcript IDs), there are ~4,000 individual sequences specified in the file. As these 1,029 transcripts correspond to just 634 candidate genes, there is a certain amount of overlap between some transcripts. For example, transcript 1 and transcript 2 might both contain the exact same exon 1, but have different exon 2s. In this case, there is 2 entries for that exon 1 in the GTF file. In other words, there are a few duplicated sequences in there.

Anyway, now to the analysis. I used @roohif’s code, described in his blog post:

It takes the coordinates from the GTF file and extracts the corresponding sequences in the human genome (I used assembly GRCh37.p13, since that’s what Ruiz-Orera et al. used as their reference), giving a series of multi-fasta files - 1 for each chromosome (Sequences less than 30bp long were excluded because they’ll find matches just by chance). Then the sequences in these FASTA files are BLASTed against the chimp genome (I used the latest assembly: PTRv2).

The results are .csv files containing all the vital statistics about each individual BLAST search. Each row looks something like this:

qseqid qstart end sseqid sstart send qlen length nident pident evalue
263 1 71 NC_036880.1 94512337 94512407 71 71 71 100.00 1,74E-30

In this case, you can see that the query sequence was 71bp long (qlen), and a match was found in the chimp genome that covered all 71bp of the query (length), and that 71bp (out of 71bp) of this chimp sequence matches the human query sequence (nident), meaning that the percentage identity is 100% (pident).

There are 3973 rows like this in the final .csv file, covering all of the sequences that make up the 1,029 transcripts. Here are some general statistics:

The total number of bases analysed was 843,105.
Of these, 833,023bp were included in BLAST hits in the chimp genome.
820,460 bases were identical between the human and chimp genomes, meaning that 97.3% of the total analysed bases were found to have a perfect match, and 98.5% of the bases that were included in the BLAST hits have a perfect match.

BLAST reports hits for subset of sequences, so in some cases the query sequence could be 1000bp long, and BLAST could find a 100% match for 100bp of this sequence, and a 0% match for the other 900bp. That’s not very helpful for our purposes, so, I didn’t rely much on pident. Instead, I divided the number of identical bases in the BLAST hit (nident) by the total number of bases in the query sequence (qlen). In the example earlier, this would return a match of 10%, which is more representative.

The average of all the sequences is 97.3% (as I mentioned earlier). The distribution of these percentages is shown below: 3773 out of 3973 sequences had at least 95% of their bases identical in the chimp genome. Of the remaining 200, just 96 sequences had below 90%, and 56 sequences had below 50%.

(S. Joshua Swamidass) #30

Of these, how many do you find strong matches for different halves of the sequences, perhaps indicating rearrangements?

(Blogging Graduate Student) #31

The average pident of these 96 sequences was 93.1%, but the BLAST search was done using “-max_hsps 1”, so only 1 BLAST result was returned per query, making it difficult to infer rearrangements.

(Curtis Henderson) #32

Thanks for working on this, @evograd. I certainly appreciate the effort. @swamidass, what does your intuition tell you about those 56? Do you suspect a large portion would be due to recombination?

(S. Joshua Swamidass) split this topic #33

35 posts were split to a new topic: Comments on James Tour on Orphan Genes

(Faizal Ali) #34

This post was flagged by the community and is temporarily hidden.

(Paul A Nelson) #35

(from Guayaquil, Ecuador – may lose internet access soon)

This is the paper Jim Tour was referring to:

In particular, this argument by Clamp et al.:

The Clamp et al. (2007) paper is now 12 years old, so should be revisited in the light of better annotations for both human and chimp genomes. If anyone reading this wants to try their luck at getting the Clamp et al. data, please give it a shot. When our ORFanBase & ORFanID research group tried some of the data links in the paper (several months ago), they had gone dead.

(Paul A Nelson) #36

Also of interest – this group working on human ORFans at Iowa State:

This presentation occurred in the Comparative Genomics and Evolution track at the ISMB / ISCB annual meeting last summer in Chicago. These researchers may be understandably protective of their preliminary data at this point, however.

Above, T. aqua pointed out that syntenic sequence often exists in another species (call it X) for a reading frame annotated as coding in species Y. These aren’t really orphans by the strictest criteria, and indeed raise a puzzle about the arrow (direction) of causality: might sequence X, which may lack a promoter, actually represent a once-transcribed and functional sequence now drifting away into pseudogene status (whereas the Y sequence is still transcribed and producing a functional product)?

(Faizal Ali) #38

If I am following this discussion correctly (far from a given), the analysis of the Ruiz-Orera paper by evograd above goes a long way towards addressing this.

(Faizal Ali) #40

This post was flagged by the community and is temporarily hidden.

(Steve Schaffner) #65

Looks like @evograd has already done the heavy lifting.


That may not be the case for all of them.

(S. Joshua Swamidass) #67

How many were discarded?