Given A Gene, Identify The World Experts

Question

[Andrew Su]

The goal of the Gene Wiki project is to create a collaboratively-written, community-reviewed, and continuously-updated review article for every gene in the human genome. So far, community contributions to these articles is going quite well based only on word-of-mouth and the Google factor.

However, we potentially could accelerate the effort if we actively recruited experts for each gene. So the question is -- given a gene, what is the best automated method to identify the world experts? Bonus points for any methods to also identify their email address. (The fall back there will be to use Amazon Mechanical Turk.)

NOTE: this question is to solicit ideas of how this would be done, but if anyone wants to work with us to actually implement this, we'd be happy to include you as an author on a future Gene Wiki paper.

EDIT: Since my comment to Pierre's answer is hidden, reposting it here. Pierre's solution is pretty awesome. The only thing that I can see right now to make it better is to put it on Google App Engine. Any GAE experts in the crowd want to comment on how hard that would be?

Answer 1

[Pierre Lindenbaum]

Andrew, why on earth are you asking this when it's 00H01 here. It's time to sleep for me !:-) However here is how I would do:

1. for each Gene, search the XML definition of the Gene using NCBI-ESearch and NCBI-EFetch
2. Search for all the publication associated to that entry by searching the tags `<pubmedid>`
3. Download the pubmed record as XML and extract the name of each author in the paper and their affiliation (it often contains the email).
4. Get the name of the most frequent author (this is the most difficult part here because the names can be ambiguous )

Note: In 2007, I collected the names and the emails of some bioinformaticians by scanning pubmed with java. See my post.

UPDATE: OK, I've quickly written a program doing the job. It is available on GIST at :

	/**
	* Author: Pierre Lindenbaum PhD
	* WWW: http://plindenbaum.blogspot.com
	* Motivation:
	* Given a gene, identify the world experts
	* http://biostar.stackexchange.com/questions/4296
	*/
	import java.net.URLEncoder;
	import java.text.Collator;
	import java.util.ArrayList;
	import java.util.Collections;
	import java.util.Comparator;
	import java.util.HashSet;
	import java.util.List;
	import java.util.Locale;
	import java.util.Set;
	import java.util.TreeSet;
	import java.util.logging.Level;
	import java.util.logging.Logger;
	import javax.xml.parsers.DocumentBuilder;
	import javax.xml.parsers.DocumentBuilderFactory;
	import javax.xml.stream.XMLOutputFactory;
	import javax.xml.stream.XMLStreamWriter;
	import javax.xml.xpath.XPath;
	import javax.xml.xpath.XPathConstants;
	import javax.xml.xpath.XPathFactory;
	import org.w3c.dom.Document;
	import org.w3c.dom.Node;
	import org.w3c.dom.NodeList;
	public class BioStar4296
	{
	private Logger LOG=Logger.getLogger(BioStar4296.class.getName());
	private String organism="Homo Sapiens";
	private DocumentBuilder docBuilder;
	private XPath xpath;
	private Collator collator;
	static class Author
	{
	String suffix="";
	String firstName="";
	String lastName="";
	String initials="";
	Set<String> mails=new HashSet<String>();
	Set<Integer> pmids=new TreeSet<Integer>();
	int factor=1;
	Set<String> affilitations=new HashSet<String>();
	@Override
	public int hashCode() {
	final int prime = 31;
	int result = 1;
	result = prime * result
	+ ((firstName == null) ? 0 : firstName.hashCode());
	result = prime * result
	+ ((lastName == null) ? 0 : lastName.hashCode());
	return result;
	}
	@Override
	public boolean equals(Object obj) {
	if (this == obj)
	return true;
	if (obj == null)
	return false;
	if (getClass() != obj.getClass())
	return false;
	Author other = (Author) obj;
	if (firstName == null) {
	if (other.firstName != null)
	return false;
	} else if (!firstName.equals(other.firstName))
	return false;
	if (lastName == null) {
	if (other.lastName != null)
	return false;
	} else if (!lastName.equals(other.lastName))
	return false;
	return true;
	}
	@Override
	public String toString() {
	return firstName+" "+lastName+" lab:"+this.affilitations+" mails:"+this.mails;
	}
	void write(XMLStreamWriter w)
	throws Exception
	{
	w.writeStartElement("Person");
	w.writeCharacters("\n");
	w.writeStartElement("firstName");
	w.writeCharacters(firstName);
	w.writeEndElement();
	w.writeCharacters("\n");
	w.writeStartElement("lastName");
	w.writeCharacters(lastName);
	w.writeEndElement();
	w.writeCharacters("\n");
	for(Integer s:pmids)
	{
	w.writeStartElement("pmid");
	w.writeCharacters(String.valueOf(s));
	w.writeEndElement();
	w.writeCharacters("\n");
	}
	for(String s:mails)
	{
	w.writeStartElement("mail");
	w.writeCharacters(s);
	w.writeEndElement();
	w.writeCharacters("\n");
	}
	for(String s:affilitations)
	{
	w.writeStartElement("affilitation");
	w.writeCharacters(s);
	w.writeEndElement();
	w.writeCharacters("\n");
	}
	w.writeEndElement();
	w.writeCharacters("\n");
	}
	}
	private BioStar4296() throws Exception
	{
	LOG.setLevel(Level.OFF);
	DocumentBuilderFactory f=DocumentBuilderFactory.newInstance();
	f.setNamespaceAware(false);
	f.setCoalescing(true);
	f.setIgnoringComments(true);
	f.setIgnoringElementContentWhitespace(true);
	f.setValidating(false);
	this.docBuilder=f.newDocumentBuilder();
	XPathFactory factory=XPathFactory.newInstance();
	this.xpath=factory.newXPath();
	this.collator= Collator.getInstance(Locale.FRENCH);
	this.collator.setStrength(Collator.PRIMARY);
	}
	private int search(XMLStreamWriter w,String geneName)
	throws Exception
	{
	w.writeCharacters("\n");
	w.writeStartElement("gene");
	w.writeAttribute("name", geneName);
	String url= "http://eutils.ncbi.nlm.nih.gov/entrez/eutils/esearch.fcgi?db=gene&term="+
	URLEncoder.encode(geneName+"[PREF] \""+this.organism+"\"[ORGN]", "UTF-8");
	LOG.info(url);
	Document dom=this.docBuilder.parse(url);
	NodeList list=(NodeList)this.xpath.evaluate(
	"/eSearchResult/IdList/Id",
	dom,XPathConstants.NODESET);
	if(list.getLength()==0)
	{
	w.writeComment("Cannot find any entry for "+geneName);
	w.writeEndElement();
	return -1;
	}
	else if(list.getLength()!=1)
	{
	w.writeComment("Ambigous name "+geneName);
	w.writeEndElement();
	return -1;
	}
	String geneId= list.item(0).getTextContent();
	LOG.info("GeneId:"+geneId);
	w.writeAttribute("geneId", geneId);
	url="http://www.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=gene&id="+
	geneId+
	"&rettype=text&retmode=xml";
	LOG.info(url);
	dom=this.docBuilder.parse(url);
	list=(NodeList)this.xpath.evaluate(
	"//PubMedId",
	dom,XPathConstants.NODESET);
	if(list.getLength()==0)
	{
	w.writeComment("No pubmed for "+geneName);
	w.writeEndElement();
	return -1;
	}
	List<Author> authors=new ArrayList<Author>();
	Set<Integer> pmidSet=new TreeSet<Integer>();
	for(int articleIdx=0;articleIdx< list.getLength();++articleIdx)
	{
	String pmid= list.item(articleIdx).getTextContent();
	LOG.info("PMID:"+pmid);
	pmidSet.add(Integer.parseInt(pmid));
	}
	w.writeAttribute("count-pmids",String.valueOf(pmidSet.size()));
	w.writeCharacters("\n");
	for(Integer pmid: pmidSet)
	{
	LOG.info("PMID:"+pmid);
	url="http://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=pubmed&id="+pmid+"&retmode=xml";
	LOG.info("url:"+url);
	dom=this.docBuilder.parse(url);
	Node n=(Node)this.xpath.evaluate("//Affiliation", dom,XPathConstants.NODE);
	if(n==null) continue;
	String affiliation=n.getTextContent();
	String adressFragments[]=affiliation.split("[ \t\:\<,\>\(\)]");
	LOG.info("affiliation:"+affiliation);
	NodeList authorList=(NodeList)this.xpath.evaluate(
	"//AuthorList/Author",
	dom,XPathConstants.NODESET);
	LOG.info("Authors:"+authorList.getLength());
	if(authorList.getLength()==0) continue;
	for(int j=0;j< authorList.getLength();++j)
	{
	boolean collective=false;
	Author author=new Author();
	for(Node c1=authorList.item(j).getFirstChild();c1!=null;c1=c1.getNextSibling())
	{
	if(c1.getNodeType()!=Node.ELEMENT_NODE) continue;
	String tag=c1.getNodeName();
	String content= c1.getTextContent();
	if(tag.equals("LastName"))
	{
	author.lastName= content;
	}
	else if(tag.equals("FirstName") || tag.equals("ForeName"))
	{
	author.firstName= content;
	}
	else if(tag.equals("Initials"))
	{
	author.initials= content;
	}
	else if(tag.equals("CollectiveName"))
	{
	collective=true;
	break;
	}
	else if(tag.equals("Suffix"))
	{
	author.suffix= content;
	}
	}
	if(collective) continue;
	LOG.info("Make New Author:"+author);
	int k=0;
	for(k=0;k< authors.size();++k)
	{
	Author p=authors.get(k);
	if( !p.firstName.isEmpty() &&
	this.collator.compare(p.firstName,author.firstName)==0 &&
	this.collator.compare(p.lastName,author.lastName)==0)
	{
	LOG.info("Same: "+p+" "+author);
	author=p;
	break;
	}
	}
	if(k==authors.size())
	{
	k=0;
	for(k=0;k< authors.size();++k)
	{
	Author p=authors.get(k);
	if(
	(
	(!author.initials.isEmpty() && p.firstName.toLowerCase().startsWith(author.initials.toLowerCase())) ||
	(!p.initials.isEmpty() && author.firstName.toLowerCase().startsWith(p.initials) )||
	this.collator.compare(p.initials,author.initials)==0 )&&
	this.collator.compare(p.lastName,author.lastName)==0)
	{
	LOG.info("Same: "+p+" "+author);
	if(p.firstName.length()< author.firstName.length())
	{
	p.firstName=author.firstName;
	}
	author=p;
	break;
	}
	}
	}
	if(k==authors.size())
	{
	LOG.info("Adding: "+author);
	authors.add(author);
	}
	author.factor*=j;
	author.affilitations.add(affiliation);
	author.pmids.add(pmid);
	if(affiliation.indexOf('@')!=-1)
	{
	for(String mail: adressFragments)
	{
	mail.replaceAll("\{\}", "");
	if(mail.endsWith(".")) mail= mail.substring(0,mail.length()-1);
	int index=mail.indexOf('@');
	if(index==-1) continue;
	String mailPrefix=mail.substring(0,index).toLowerCase();
	if(mailPrefix.contains(author.lastName.toLowerCase()) ||
	collator.compare(mailPrefix, author.lastName)==0)
	{
	LOG.info("Adding: "+mail+" to "+author);
	author.mails.add(mail.toLowerCase());
	}
	else if( author.firstName.length()>1 &&
	(mailPrefix.contains( author.firstName.toLowerCase()) ||
	collator.compare(mailPrefix, author.firstName)==0))
	{
	LOG.info("Adding: "+mail+" to "+author);
	author.mails.add(mail.toLowerCase());
	}
	}
	}
	}
	}
	if(authors.isEmpty())
	{
	w.writeComment("No Author found");
	w.writeEndElement();
	return -1;
	}
	Collections.sort(authors,new Comparator<Author>()
	{
	@Override
	public int compare(Author o1, Author o2)
	{
	int i= o2.pmids.size()-o1.pmids.size();
	if(i!=0) return i;
	i= o2.factor-o1.factor;//later is more interesting ? not sure...
	return i;
	}
	});
	authors.get(0).write(w);
	w.writeEndElement();
	return 0;
	}
	public static void main(String[] args)
	{
	try {
	BioStar4296 app= new BioStar4296();
	int optind=0;
	while(optind<args.length)
	{
	if(args[optind].equals("-h"))
	{
	System.err.println("Pierre Lindenbaum");
	System.err.println("Options:");
	System.err.println(" -o <organism> ["+app.organism+"]");
	System.err.println(" -v show logs");
	return;
	}
	else if(args[optind].equals("-o"))
	{
	app.organism=args[++optind];
	}
	else if(args[optind].equals("-v"))
	{
	app.LOG.setLevel(Level.ALL);
	}
	else if(args[optind].equals("--"))
	{
	optind++;
	break;
	}
	else if(args[optind].startsWith("-"))
	{
	System.err.println("Unnown option: "+args[optind]);
	return;
	}
	else
	{
	break;
	}
	++optind;
	}
	if(optind==args.length)
	{
	System.err.println("Gene Name missing");
	}
	else
	{
	XMLOutputFactory xmlfactory= XMLOutputFactory.newInstance();
	XMLStreamWriter w= xmlfactory.createXMLStreamWriter(System.out,"UTF-8");
	w.writeStartDocument("UTF-8","1.0");
	w.writeCharacters("\n");
	w.writeStartElement("experts");
	w.writeCharacters("\n");
	while(optind < args.length)
	{
	app.search(w,args[optind]);
	optind++;
	w.writeCharacters("\n");
	}
	w.writeEndElement();
	w.writeEndDocument();
	w.flush();
	}
	} catch (Exception e)
	{
	e.printStackTrace();
	}
	}
	}

view raw BioStar4296.java hosted with ❤ by GitHub

	<?xml version="1.0" encoding="UTF-8"?>
	<experts>
	<gene name="ZC3H7B" geneId="23264" count-pmids="13">
	<Person>
	<firstName>Sumio</firstName>
	<lastName>Sugano</lastName>
	<pmid>8125298</pmid>
	<pmid>9373149</pmid>
	<pmid>14702039</pmid>
	<affilitation>International and Interdisciplinary Studies, The University of Tokyo, Japan.</affilitation>
	<affilitation>Institute of Medical Science, University of Tokyo, Japan.</affilitation>
	<affilitation>Helix Research Institute, 1532-3 Yana, Kisarazu, Chiba 292-0812, Japan.</affilitation>
	</Person>
	</gene>
	<gene name="eif4G1" geneId="1981" count-pmids="106">
	<Person>
	<firstName>Nahum</firstName>
	<lastName>Sonenberg</lastName>
	<pmid>7651417</pmid>
	<pmid>7935836</pmid>
	<pmid>8449919</pmid>
	<pmid>8521827</pmid>
	<pmid>9372926</pmid>
	<pmid>9418880</pmid>
	<pmid>9857202</pmid>
	<pmid>9878069</pmid>
	<pmid>10523622</pmid>
	<pmid>10753870</pmid>
	<pmid>10872469</pmid>
	<pmid>10996799</pmid>
	<pmid>15193258</pmid>
	<pmid>15234964</pmid>
	<pmid>15314020</pmid>
	<pmid>15961545</pmid>
	<pmid>16698552</pmid>
	<pmid>19114555</pmid>
	<pmid>19203580</pmid>
	<pmid>20053821</pmid>
	<affilitation>Department of Biochemistry and McGill Cancer Center, McGill University, Montreal, H3G 1Y6, Quebec, Canada.</affilitation>
	<affilitation>Department of Biochemistry, McGill University, Montreal, Quebec, Canada.</affilitation>
	<affilitation>Laboratories of Molecular Biophysics, The Rockefeller University, New York, New York 10021, USA.</affilitation>
	<affilitation>Molecular Histology Unit, San Raffaele Institute, 20132 Milan, Italy.</affilitation>
	<affilitation>Department of Biochemistry McGill University, Montréal, Québec, Canada. [e-mail]</affilitation>
	<affilitation>Department of Biochemistry and McGill Cancer Cancer Center, McGill University, 3655 Drummond Street, Montréal, Québec, H3G 1Y6 Canada.</affilitation>
	<affilitation>Department of Biochemistry and McGill Cancer Centre, McGill University, Montreal, Quebec, Canada.</affilitation>
	<affilitation>Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, MA, USA. [e-mail]</affilitation>
	<affilitation>Dana-Farber Cancer Institute, Harvard Medical School, Boston, Massachusetts 02115, USA.</affilitation>
	<affilitation>Department of Biochemistry, University of California, Riverside, Riverside, California 92521, USA.</affilitation>
	<affilitation>Department of Biochemistry, McGill University, Montréal, Quebec, Canada.</affilitation>
	<affilitation>Department of Biochemistry and McGill Cancer Center, McGill University, Montréal, Québec, Canada.</affilitation>
	<affilitation>Department of Biochemistry and McGill Cancer Centre, McGill University, Drummond Street 3655, Montreal, Quebec, Canada H3G 1Y6.</affilitation>
	<affilitation>Department of Biochemistry, McGill Cancer Center, McGill University, Montreal, Quebec, Canada H3G 1Y6.</affilitation>
	<affilitation>Department of Biochemistry and Molecular Genetics, University of Illinois at Chicago, 60607, USA. [e-mail]</affilitation>
	<affilitation>Department of Biochemistry, McGill University, Montréal, Québec, Canada.</affilitation>
	<affilitation>Department of Medicine, University of Minnesota, Minneapolis, MN 55455, USA.</affilitation>
	<affilitation>Department of Biochemistry, Case Western Reserve University, Cleveland, Ohio 44106.</affilitation>
	</Person>
	</gene>
	<gene name="PRNP" geneId="5621" count-pmids="429">
	<Person>
	<firstName>John</firstName>
	<lastName>Collinge</lastName>
	<pmid>1352724</pmid>
	<pmid>1677164</pmid>
	<pmid>2159587</pmid>
	<pmid>2567794</pmid>
	<pmid>8364585</pmid>
	<pmid>11283320</pmid>
	<pmid>11704923</pmid>
	<pmid>12514748</pmid>
	<pmid>12621436</pmid>
	<pmid>12690204</pmid>
	<pmid>15123682</pmid>
	<pmid>15539564</pmid>
	<pmid>16099923</pmid>
	<pmid>16156720</pmid>
	<pmid>16342955</pmid>
	<pmid>16824036</pmid>
	<pmid>16847141</pmid>
	<pmid>16925523</pmid>
	<pmid>17709704</pmid>
	<pmid>18638557</pmid>
	<pmid>19081515</pmid>
	<pmid>19218199</pmid>
	<pmid>19321423</pmid>
	<pmid>19369250</pmid>
	<pmid>19923577</pmid>
	<pmid>20109837</pmid>
	<pmid>20583301</pmid>
	<mail>[e-mail]</mail>
	<affilitation>MRC Prion Unit, Department of Neurodegenerative Disease, Institute of Neurology, Queen Square, London WC1N 3BG, UK.</affilitation>
	<affilitation>CNS Infection and Immunity Group, Department of Neurogenetics, Division of Neurosciences and Psychological Medicine, Faculty of Medicine, Imperial College, Norfolk Place, London W2 1PG, UK.</affilitation>
	<affilitation>Krebs Institute for Biomolecular Research, Department of Molecular Biology and Biotechnology, University of Sheffield, Sheffield S10 2TN, UK.</affilitation>
	<affilitation>MRC Prion Unit and Department of Neurogenetics, Imperial College School of Medicine at St. Mary's, London, United Kingdom. [e-mail]</affilitation>
	<affilitation>Division of Neuroscience (Neurophysiology), Medical School, University of Birmingham, Edgbaston, Birmingham, UK. [e-mail]</affilitation>
	<affilitation>Division of Psychiatry, Clinical Research Centre, Harrow, Middlesex, UK.</affilitation>
	<affilitation>MRC Prion Unit and National Prion Clinic, UCL Institute of Neurology and National Hospital for Neurology and Neurosurgery, London, UK.</affilitation>
	<affilitation>Medical Research Council Prion Unit and Department of Neurodegenerative Disease, Institute of Neurology, Queen Square, London, UK.</affilitation>
	<affilitation>MRC Prion Unit and Department of Neurodegenerative Disease, UCL Institute of Neurology, National Hospital for Neurology and Neurosurgery, Queen Square, London WC1N 3BG, UK.</affilitation>
	<affilitation>MRC Prion Unit, Department of Neurodegenerative Disease, UCL Institute of Neurology, National Hospital for Neurology and Neurosurgery, Queen Square, London, WC1N 3BG, UK.</affilitation>
	<affilitation>MRC Prion Unit, Department of Neurodegenerative Disease, Institute of Neurology, University College London, Queen Square, London WC1N 3BG, UK.</affilitation>
	<affilitation>MRC Prion Unit and Department of Neurogenetics, Imperial College, St. Mary's Hospital, London, United Kingdom.</affilitation>
	<affilitation>Medical Research Council Prion Unit, Department of Neurodegenerative Disease, Institute of Neurology, University College London, Queen Square, London WC1N 3BG, United Kingdom.</affilitation>
	<affilitation>Medical Research Council (MRC) Prion Unit and Department of Neurodegenerative Disease, Institute of Neurology, University College London, Queen Square, London WC1N 3BG, UK.</affilitation>
	<affilitation>MRC Prion Unit and National Prion Clinic, Institute of Neurology and National Hospital for Neurology and Neurosurgery, Queen Square, London, UK. [e-mail]</affilitation>
	<affilitation>MRC Prion Unit, Department of Neurodegenerative Diseases, Institute of Neurology, University College London, Queen Square, London WC1N 3BG, UK.</affilitation>
	<affilitation>Medical Research Council Prion Unit, Institute of Neurology, Queen Square, London WC1N 3BG, United Kingdom.</affilitation>
	<affilitation>Division of Psychiatry, Clinical Research Centre, Harrow, Middlesex, U.K.</affilitation>
	<affilitation>Medical Research Council Prion Unit, Department of Neurodegenerative Disease, University College London Institute of Neurology, United Kingdom.</affilitation>
	<affilitation>Division of Psychiatry, Clinical Research Centre, Harrow, UK.</affilitation>
	<affilitation>Department of Biochemistry and Molecular Genetics, St Mary's Hospital Medical School, London, UK.</affilitation>
	<affilitation>Department of Biochemistry and Molecular Genetics, St. Mary's Hospital Medical School, Norfolk Place, London, UK.</affilitation>
	<affilitation>MRC Prion Unit and Department of Neurodegenerative Disease, Institute of Neurology, Queen Square, London, UK.</affilitation>
	<affilitation>Medical Research Council Prion Unit, and Department of Neurodegenerative Disease, Institute of Neurology, University College, Queen Square, London WC1N 3BG, UK.</affilitation>
	</Person>
	</gene>
	</experts>

view raw experts.xml hosted with ❤ by GitHub

Here is an excerpt from the output for 3 genes : ZC3H7B, EIF4G1 and PRNP.

<?xml version="1.0" encoding="UTF-8"?>
<experts>
  <gene name="ZC3H7B" geneId="23264" count-pmids="13">
    <Person>
      <firstName>Sumio</firstName>
      <lastName>Sugano</lastName>
      <pmid>8125298</pmid>
      <pmid>9373149</pmid>
      <pmid>14702039</pmid>
      <affilitation>International and Interdisciplinary Studies, The University of Tokyo, Japan.</affilitation>
      <affilitation>Institute of Medical Science, University of Tokyo, Japan.</affilitation>
      <affilitation>Helix Research Institute, 1532-3 Yana, Kisarazu, Chiba 292-0812, Japan.</affilitation>
    </Person>
  </gene>
  <gene name="eif4G1" geneId="1981" count-pmids="106">
    <Person>
      <firstName>Nahum</firstName>
      <lastName>Sonenberg</lastName>
      <pmid>7651417</pmid>
      <pmid>7935836</pmid>
      <pmid>8449919</pmid>
      (...)
      <affilitation>Department of Biochemistry and McGill Cancer Center, McGill University, Montreal, H3G 1Y6, Quebec, Canada.</affilitation>
      <affilitation>Department of Biochemistry, McGill University, Montreal, Quebec, Canada.</affilitation>
      <affilitation>Laboratories of Molecular Biophysics, The Rockefeller University, New York, New York 10021, USA.</affilitation>
      (...)
    </Person>
  </gene>
  <gene name="PRNP" geneId="5621" count-pmids="429">
    <Person>
      <firstName>John</firstName>
      <lastName>Collinge</lastName>
      <pmid>1352724</pmid>
      <pmid>1677164</pmid>
      <pmid>2159587</pmid>
      <pmid>20583301</pmid>
      (...)
      <mail>[e-mail]</mail>
      <affilitation>Krebs Institute for Biomolecular Research, Department of Molecular Biology and Biotechnology, University of Sheffield, Sheffield S10 2TN, UK.</affilitation>
      <affilitation>MRC Prion Unit and Department of Neurogenetics, Imperial College School of Medicine at St. Mary's, London, United Kingdom. [e-mail]</affilitation>
      <affilitation>Division of Neuroscience (Neurophysiology), Medical School, University of Birmingham, Edgbaston, Birmingham, UK. [e-mail]</affilitation>
    (...)
    </Person>
  </gene>
</experts>

In the case of ZC3H7B, the result is wrong. Dr Sugano (3 articles) just used this Gene in a set of other Genes. The expert would be D. Poncet, my former thesis advisor but his number of articles about this protein is 2 articles.

Eif4G1: I know that Dr Sonenberg is the expert. His email wasn't found.

PRNP: Collinge seems to be the expert. His e-mail was found too.

Update the code

https://gist.github.com/lindenb/740496

[Neilfws]

Actually, journal impact factor has nothing to do with the importance of individual articles. Common misconception :-) See http://altmetrics.org/manifesto/.

[Pierre Lindenbaum]

@Larry, most(all?) Genes in GeneWiki are Human Genes.

[Larry_Parnell]

A serious issue that emerges is gene synonyms. One may need to consider species as well because the same gene in different organisms will function differently. So, Pierre's step 2 needs some refinement but nonetheless gets my vote.

[Andrew Su]

Pierre -- Science doesn't sleep, so neither should you... ;)

[Mns]

Nice Q/A. You could also take into account the journal impact factor. This way the authors would be ranked by additionally relying on the "quality" of their work. Also, attention should be given to avoid mail receivers to mark the mail as spam.

[Andrew Su]

As usual, I'm thoroughly impressed. I wonder how difficult it would be to adapt this to Google App Engine so we could call it as a web service. (On my first attempt to run locally, I get an error that is likely due to my complete java ignorance...) Any GAE experts in the audience?

[Pierre Lindenbaum]

I don't think that Google App engine would be the best place to run this service: there is a lot of I/O and it could be slow (for example in PRNP , 429 article were downloaded )

[Andrew Su]

Bummer, was wondering if that would be an issue (as it is with my pubmed2wordle app). Anyway, thanks!

[Rm]

Amazing answer!!! Refining with Gene synonyms, Journal impact factor, Article views or downloads, Grants obtained by author (if any), works on similar genes, publication of invited reviews....

Answer 2

[Larry_Parnell]

Interesting question and project! We tried a similar approach back in 1998 when annotating the Arabidopsis genome. We were Christine Schueller, then at MIPS, several plant biologists who knew genes/gene families and wanted to get their hands on the sequence data prior to GenBank/EMBL deposits, and myself.

I feel that any method needs to consider more than one expert, even more than one identified person per paper, because the last author may not have to want to spend the time on this or the paper represents a collaborative effort that really has distributed expertise.

Who qualifies as an expert. If we are the only lab to publish on a given LOC1234.. gene, we may be the experts with very little data. Are we also experts on APOA5, publishers of 21 of the 314 papers on this gene? Perhaps, certainly in terms of genetic variants and the response to fat in the diet. Len Pennacchio was the first author on the first APOA5 paper but this is no longer an interest of his nor of Ed Rubin, in whose lab he worked at the time. In other words, as the number of pubs grows, expertise becomes fragmented and specialized.

All this is not much of an answer but more in terms of guidance based on what I've seen. This is a hard question to answer satisfactorily.

[Andrew Su]

Nice answer, thanks Larry. I'll follow up to say that we indeed are looking to identify multiple experts. In fact, the beauty (and the curse) of the Gene Wiki is that anyone can edit these gene pages. So our goal is certainly not to anoint one person to be the authority, but simply to notify a few relevant people that the Gene Wiki page exists...

[Pierre Lindenbaum]

Larry, I guess Andrew is just looking for a few people/experts that would be qualified for editing an article about 'their' gene in wikipedia.

[Larry_Parnell]

That's all clear. My point is one of the difficulty that arises as more knowledge about a gene is published. Identifying when an expert or experts for an emerging role of that gene is required can be difficult.

Answer 3

[Bio_X2Y]

It might be interesting to take a two-phase approach - set up a system where people can nominate themselves as experts for a gene, and advertise the feature on forums, blogs, etc. Supplement this by automatically identifying candidate experts, and drawing their attention to the opt-in system via a mailshot.

I would agree with Pierre's approach for identifying candidates - search for papers on PubMed containing the name of the gene in the title/abstract. It might be useful to include known synonyms, e.g. from HGNC.

The PubMed XML contains a distinct element for each author, but afaik there is usually only one Affiliation element that would typically contain a single e-mail address. Complications include the fact that the element seems to be free-text with no standard format, and there is no definitive way to associate the e-mail address with a specific author. Even so, I would pluck out the address with the basic [e-mail] pattern.

I would ignore the authors, since I assume you don't have time to manually find 10,000+ e-mail addresses, or even to verify them if they were automatically found. Instead, focus on the one e-mail address per paper, and consolidate them for the mailshot. For each gene, take only the 50 most common e-mail addresses and send the mail with "Dear Colleague". Your scripts should be able to customise each mail to list the gene (or genes) that the e-mail address is associated with.

[Andrew Su]

I would suggest that your phase one "opt-in" system pretty much describes the Gene Wiki itself. I'm not sure it would be useful to add in another intermediate step. Good to know on the affiliation field -- though a bit of spot checking shows that the presence of an email address is far from universal...

[Will]

I like the "blog nomination" idea ... because then you would actually get people with the "internet acumen" to actually keep everything up to date.

Answer 4

[David Quigley]

If you want to mine the literature to suggest candidates, you may want to distinguish between reviews and primary literature. Reviews, particularly at high-impact journals (e.g. Nature Reviews) are often invited by the editors, and the senior authors are typically either someone with a strong track record or someone who has recently published an interesting result.

The position of authors is also informative; someone who was senior or first author is more likely to be an expert on the overall science presented by the paper.

You may also need to consider how someone could be an expert in the role of a gene in fly development but may know little about its role in colon cancer. Saying someone is a "Myc expert" usually means in a particular context.

Those are obvious points, but they aren't clearly stated in the above notes.

If you spam thousands of scientists with an invitation to write for any website, no matter how virtuous, you are likely to invite a backlash. Getting some prominent proponents on board and getting them to spread the word organically may be a better approach.

[Andrew Su]

Good points all. On your last point, I would argue that the "organic" approach is what we're doing now, and it's good but slow. "Spam" would be to email anyone who's ever published on a gene; fast but likely counterproductive. We're shooting for something in between, rooted in "organic" but slightly more proactive. I'd label it as "outreach"... ;)

Answer 5

[Giovanni M Dall'Olio]

well... a possible solution is to look at WikiGenes, which automatically creates a template for each gene by automatic analysis of the literature.

For example, you can look at the entry for a specific gene and get the list of authors from there.

p.s. why don't you try to collaborate with WikiGenes and the other efforts? For users, it is disturbing to see so many resources spent and no communication between similar projects.

[Andrew Su]

Using the text-mined and referenced content in WikiGenes is an interesting idea. Will definitely think about it. More generally, we have friendly communications with others in this sphere. But ultimately each group has a slightly different model, and we all believe that they are all worth trying out...