Biomedical named entities often play important roles in many biomedical text mining tools. However, due to the incompleteness of provided synonyms and numerous variations in their surface forms, normalization of biomedical entities is very challenging. In this paper, we focus on learning representations of biomedical entities solely based on the synonyms of entities . To learn from the incomplete synonyms, we use a model-based candidate selection and maximize the marginal likelihood of the synonyms present in top candidates. Our model-based candidates are iteratively updated to contain more difficult negative samples as our model evolves. In this way, we avoid the explicit pre-selection of negative samples from more than 400K candidates. On four biomedical entity normalization datasets having three different entity types (disease, chemical, adverse reaction), our model BIOSYN consistently outperforms previous state-of-the-art models almost reaching the upper bound on each dataset. Biomedical named entities are frequently used as key features in biomedical text mining. From biomedical relation extraction (Xu et al., 2016; Li et al., 2017a) to literature search engines (Lee et al., 2016), many studies are utilizing biomedical named entities as a basic building block of their methodologies. While the extraction of the biomedical named entities is studied extensively (Sahu and Anand, 2016; Habibi et al., 2017), the normalization of extracted named entities is also crucial for improving the precision of downstream tasks (Leaman et al., 2013; Wei et al., 2015). Unlike named entities from general domain text, typical biomedical entities have several different surface forms, making the normalization of biomedical entities very challenging. For instance, while two chemical entities  motrin ' and  ibuprofen ' belong to the same concept ID (MeSH:D007052), they have completely different surface forms. On the other hand, mentions having similar surface forms could also have different meanings (e.g.  dystrophinopathy ' (MeSH:D009136) and  bestrophinopathy ' (MeSH:C567518)). These examples show a strong need for building latent representations of biomedical entities that capture semantic information of the mentions. In this paper, we propose a novel framework for learning biomedical entity representations based on the synonyms of entities. Previous works on entity normalization mostly train binary classifiers that decide whether the two input entities are the same (positive) or different (negative) (Leaman et al., 2013; Li et al., 2017b; Fakhraei et al., 2019; Phan et al., 2019). Our framework called BIOSYN uses the synonym marginalization technique, which maximizes the probability of all synonym representations in top candidates. We represent each biomedical entity using both sparse and dense representations to capture morphological and semantic information, respectively. The candidates are iteratively updated based on our model's representations removing the need for an explicit negative sampling from a large number of candidates. Also, the model-based candidates help our model learn from more difficult negative samples. Through extensive experiments on four biomedical entity normalization datasets, we show that BIOSYN achieves new state-of-the-art performance on all datasets, outperforming previous models by 0.8%  2.6% top1 accuracy. Further analysis shows that our model's performance has almost reached the performance upper bound of each dataset. The contributions of our paper are as follows: First, we introduce BIOSYN for biomedical entity representation learning, which uses synonym marginalization dispensing with the explicit needs of negative training pairs. Second, we show that the iterative candidate selection based on our model's representations is crucial for improving the performance together with synonym marginalization. Finally, our model outperforms strong state-of-the-art models up to 2.6% on four biomedical normalization datasets. 1 2 Related Works Biomedical entity representations have largely relied on biomedical word representations. Right after the introduction of Word2vec (Mikolov et al., 2013), Pyysalo et al. (2013) trained Word2Vec on biomedical corpora such as PubMed. Their biomedical version of Word2Vec has been widely used for various biomedical natural language processing tasks (Habibi et al., 2017; Wang et al., 2018; Giorgi and Bader, 2018; Li et al., 2017a) including the biomedical normalization task (Mondal et al., 2019). Most recently, BioBERT (Lee et al., 2019) has been introduced for contextualized biomedical word representations. BioBERT is pre-trained on biomedical corpora using BERT (Devlin et al., 2019) and numerous studies are utilizing BioBERT for building state-of-the-art biomedical NLP models (Lin et al., 2019; Jin et al., 2019; Alsentzer et al., 2019; Sousa et al., 2019). Our model also uses pre-trained BioBERT for learning biomedical entity representations. The intrinsic evaluation of the quality of biomedical entity representations is often verified by the biomedical entity normalization task (Leaman et al., 2013; Phan et al., 2019). The goal of the biomedical entity normalization task is to map an input mention from a biomedical text to its associated CUI (Concept Unique ID) in a dictionary. The task is also referred to as the entity linking or the entity grounding (D'Souza and Ng, 2015; Leaman and Lu, 2016). However, the normalization of biomedical entities is more challenging than the normalization of general domain entities due to a large number of synonyms. Also, the variations of synonyms depend on their entity types, which makes building type-agnostic normalization model difficult (Leaman et al., 2013; Li et al., 2017b; Mon-1 Code available at https://github.com/ dmis-lab/BioSyn . dal et al., 2019). having three different biomedical entity types. While traditional biomedical entity normalization models are based on hand-crafted rules (D'Souza and Ng, 2015; Leaman et al., 2015), recent approaches for the biomedical entity normalization have been significantly improved with various machine learning techniques. DNorm (Leaman et al., 2013) is one of the first machine learning-based entity normalization models, which learns pair-wise similarity using tf-idf vectors. Another machine learning-based study is CNN-based ranking method (Li et al., 2017b), which learns entity representations using a convolutional neural network. The most similar works to ours are NSEEN (Fakhraei et al., 2019) and BNE (Phan et al., 2019), which map mentions and concept names in dictionaries to a latent space using LSTM models and refines the embedding using the negative sampling technique. However, most previous works adopt a pair-wise training procedure that explicitly requires making negative pairs. Our work is based on marginalizing positive samples (i.e., synonyms) from iteratively updated candidates and avoids the problem of choosing a single negative sample. In our framework, we represent each entity with sparse and dense vectors which is largely motivated by techniques used in information retrieval. Models in information retrieval often utilize both sparse and dense representations (Ramos et al., 2003; Palangi et al., 2016; Mitra et al., 2017) to retrieve relevant documents given a query. Similarly, we can think of the biomedical entity normalization task as retrieving relevant concepts given a mention (Li et al., 2017b; Mondal et al., 2019). In our work, we use maximum inner product search (MIPS) for retrieving the concepts represented as sparse and dense vectors, whereas previous models could suffer from error propagation of the pipeline approach. We define an input mention m as an entity string in a biomedical corpus. Each input mention has its own CUI c and each CUI has one or more synonyms defined in the dictionary. The set of synonyms for a CUI is also called as a synset. We denote the union of all synonyms in a dictionary as N = [ n 1 , n 2 , . . . ] where n  N is a single syn-Encoder Dictionary Update Embeddings ( ith iteration ) SimilarityScores Encoder Top 1 Iterative Update with Synonym Marginalization Circovirus SARS Pneumonia SARSCOV2 Influenza COVID19 WeightSharing Training Inference MERS Mention covid19 MML over Synonyms ( indicates synonyms) Sort Top k : Inner Product = Figure 1: The overview of BIOSYN . onym string. Our goal is to predict the gold CUI c  of the input mention m as follows: c  = CUI (argmax n  NP ( n | m ;  )) (1) where CUI (  ) returns the CUI of the synonym n and  denotes a trainable parameter of our model. The overview of our framework is illustrated in Figure 1. We first represent each input mention m and each synonym n in a dictionary using sparse and dense representations. We treat m and n equally and use a shared encoder for both strings. During training, we iteratively update top candidates and calculate the marginal probability of the synonyms based on their representations. At inference time, we find the nearest synonym by performing MIPS over all synonym representations. Sparse Entity Representation We use tf-idf to obtain a sparse representation of m and n . We denote each sparse representation as e sm and e sn for the input mention and the synonym, respectively. tf-idf is calculated based on the character-level n-grams statistics computed over all synonyms n  N . We define the sparse scoring function of a mention-synonym pair ( m, n ) as follows: S sparse ( m, n ) = f ( e sm , e sn )  R (2) where f denotes a similarity function. Dense Entity Representation While the sparse representation encodes the morphological information of given strings, the dense representation encodes the semantic information. Learning effective dense representations is the key challenge in the biomedical entity normalization task (Li et al., 2017b; Mondal et al., 2019; Phan et al., 2019; Fakhraei et al., 2019). We use pre-trained BioBERT (Lee et al., 2019) to encode dense representations and fine-tune BioBERT with our synonym marginalization algorithm. 2 We share the same BioBERT model for encoding mention and synonym representations. We compute the dense representation of the mention m as follows: e dm = BioBERT ( m )[ CLS ]  R h (3) where m = { m 1 , ..., m l } is a sequence of sub-tokens of the mention m segmented by the Word-Piece tokenizer (Wu et al., 2016) and h denotes the hidden dimension of BioBERT (i.e., h = 768 ). [ CLS ] denotes the special token that BERT-style models use to compute a single representative vector of an input. The synonym representation e dn  R h is computed similarly. We denote the dense scoring function of a mention-synonym pair ( m, n ) using the dense representations as follows: S dense ( m, n ) = f ( e dm , e dn )  R (4) where we again used the inner product for f . Similarity Function Based on the two similarity functions S sparse ( m, n ) and S dense ( m, n ) , we now define the final similarity function S ( m, n ) indicating the similarity between an input mention m and a synonym n : S ( m, n ) = S dense ( m, n ) + S sparse ( m, n )  R (5) where  is a trainable scalar weight for the sparse score. Using  , our model learns to balance the importance between the sparse similarity and the dense similarity. The most common way to train the entity representation model is to build a pair-wise training dataset. While it is relatively convenient to sample positive pairs using synonyms, sampling negative pairs are trickier than sampling positive pairs as there are a vast number of negative candidates. For instance, the mention  alpha conotoxin ' (MeSH:D020916) has 6 positive synonyms while its dictionary has 407,247 synonyms each of which can be a negative sampling candidate. Models trained on these pair-wise datasets often rely on the quality of the negative sampling (Leaman et al., 2013; Li et al., 2017b; Phan et al., 2019; Fakhraei et al., 2019). On the other hand, we use a model-based candidate retrieval and maximize the marginal probability of positive synonyms in the candidates. Iterative Candidate Retrieval Due to a large number of candidates present in the dictionary, we need to retrieve a smaller number of candidates for training. In our framework, we use our entity encoder to update the top candidates iteratively. Let k be the number of top candidates to be retrieved for training and  ( 0    1 ) be the ratio of candidates retrieved from S dense . We call  as the dense ratio and  = 1 means consisting the candidates with S dense only. First, we compute the sparse scores S sparse and the dense scores S dense for all n  N . Then we retrieve the k  (cid:98) k (cid:99) highest candidates using S sparse , which we call as sparse candidates. Likewise, we retrieve the (cid:98) k (cid:99) highest candidates using S dense , which we call as dense candidates. Whenever the dense and sparse candidates overlap, we add more dense candidates to match the number of candidates as k . While the sparse candidates for a mention will always be the same as they are based on the static tf-idf representation, the dense candidates change every epoch as our model learns better dense representations. Our iterative candidate retrieval method has the following benefits. First, it makes top candidates to have more difficult negative samples as our model is trained, hence helping our model represent a more accurate dense representation of each entity. Also, it increases the chances of retrieving previously unseen positive samples in the top candidates. As we will see, comprising the candidates purely with sparse candidates have a strict upper bound while ours with dense candidates can maximize the upper bound. Synonym Marginalization Given the top candidates from iterative candidate retrieval, we maximize the marginal probability of positive synonyms, which we call as synonym marginalization. Given the top candidates N 1: k computed from our model, the probability of each synonym is obtained as: P ( n | m ;  ) = exp( S ( n, m )) (cid:80) n (cid:48)  N 1: k exp( S ( n (cid:48) , m )) (6) where the summation in the denominator is over the top candidates N 1: k . Then, the marginal probability of the positive synonyms of a mention m is defined as follows: P (cid:48) ( m, N 1: k ) = (cid:88) n  N 1: k EQUAL ( m,n )=1 P ( n | m ;  ) (7) where EQUAL ( m, n ) is 1 when CUI ( m ) is equivalent to CUI ( n ) and 0 otherwise. Finally, we minimize the negative marginal log-likelihood of synonyms. We define the loss function of our model as follows: L =  1 MM (cid:88) i =1 log P (cid:48) ( m i , N 1: k ) (8) where M is the number of training mentions in our dataset. We use mini-batch for the training and use Adam optimizer (Kingma and Ba, 2015) to minimize the loss. At inference time, we retrieve the nearest synonym of a mention representation using MIPS. We compute the similarity score S ( m, n ) between the input mention m and all synonyms n  N using the inner product and return the CUI of the nearest candidate. Note that it is computationally cheap to find the nearest neighbors once we pre-compute the dense and sparse representations of all synonyms. We perform basic pre-processings such as lower-casing all characters and removing the punctuation for both mentions and synonyms. To resolve the typo issues in mentions from NCBI disease, we apply the spelling check algorithm following the previous work (D'Souza and Ng, 2015). Abbreviations of entities are widely used in biomedical entities for an efficient notation which makes the normalization task more challenging. Therefore, we use the abbreviation resolution module called Ab3P 3 to detect the local abbreviations and expand it to its definition from the context (Sohn et al., 2008). We also split composite mentions (e.g. 'breast and ovarian cancer') into separate mentions (e.g. 'breast cancer' and 'ovarian cancer') using heuristic rules described in the previous work (D'Souza and Ng, 2015). We also merge mentions in the training set to the dictionary to increase the coverage following the previous work (D'Souza and Ng, 2015). For sparse representations, we use character-level uni-, bi-grams for tf-idf. The maximum sequence length of BioBERT is set to 25 4 and any string over the maximum length is truncated to 25. The number of top candidates k is 20 and the dense ratio  for the candidate retrieval is set to 0.5. We set the learning rate to 1e-5, weight decay to 1e-2, and the mini-batch size to 16. We found that the trainable scalar  converges to different values between 2 to 4 on each dataset. We train BIOSYN for 10 epochs for NCBI Disease, BC5CDR Disease, and TAC2017 ADR and 5 epochs for BC5CDR Chemical due to its large dictionary size. Except the number of epochs, we use the same hyperparameters for all datasets and experiments. We use the top k accuracy as an evaluation metric following the previous works in biomedical entity normalization tasks (D'Souza and Ng, 2015; Li et al., 2017b; Wright, 2019; Phan et al., 2019; Ji et al., 2019; Mondal et al., 2019). We define Acc@ k as 1 if a correct CUI is included in the top k predictions, otherwise 0. We evaluate our models using Acc@1 and Acc@5. Note that we treat predictions for composite entities as correct if every prediction for each separate mention is correct. We use four biomedical entity normalization datasets having three different biomedical entity types (disease, chemical, adverse reaction). The statistics of each dataset is described in Table 1. NCBI Disease Corpus NCBI Disease Corpus (Dogan et al., 2014) 5 provides manually annotated disease mentions in each document with each CUI mapped into the MEDIC dictionary (Davis et al., 2012). In this work, we use the July 6, 2012 version of MEDIC containing 11,915 CUIs and 71,923 synonyms included in MeSH and/or OMIM ontologies. Biocreative V CDR BioCreative V CDR (Li et al., 2016) 6 is a challenge for the tasks of chemical-induced disease (CID) relation extraction. It provides disease and chemical type entities. The annotated disease mentions in the dataset are mapped into the MEDIC dictionary like the NCBI disease corpus. The annotated chemical mentions in the dataset are mapped into the Comparative Toxicogenomics Database (CTD) (Davis et al., 2018) chemical dictionary. In this work, we use the November 4, 2019 version of the CTD chemical dictionary containing 171,203 CUIs and 407,247 synonyms included in MeSH ontologies. Following the previous work (Phan et al., 2019), we filter out mentions whose CUIs do not exist in the dictionary. TAC2017ADR TAC2017ADR (Roberts et al., 2017) 7 is a challenge whose purpose of the task is to extract information on adverse reactions found in structured product labels. It provides manually annotated mentions of adverse reactions that are mapped into the MedDRA dictionary (Brown et al., 1999). In this work, we use MedDRA v18.1 which contains 23,668 CUIs and 76,817 synonyms. 5 https://www.ncbi.nlm.nih.gov/ CBBresearch/Dogan/DISEASE 6 https://biocreative.bioinformatics. udel.edu/tasks/biocreative-v/track-3-cdr 7 https://bionlp.nlm.nih.gov/ tac2017adversereactions Models NCBI Disease BC5CDR Disease BC5CDR Chemical TAC2017ADR Acc@1 Acc@5 Acc@1 Acc@5 Acc@1 Acc@5 Acc@1 Acc@5 Sieve-Based (D'Souza and Ng, 2015) 84.7 -84.1 -90.7  -84.3  Taggerone (Leaman and Lu, 2016) 87.7 -88.9 -94.1 -CNN Ranking (Li et al., 2017b) 86.1 ---NormCo (Wright, 2019) 87.8 -88.0 --BNE (Phan et al., 2019) 87.7 -90.6 -95.8 -BERT Ranking (Ji et al., 2019) 89.1 ---93.2 TripletNet (Mondal et al., 2019) 90.0 ---BIOSYN ( S-SCORE ) 87.6 90.5 92.4 95.7 95.9 96.8 91.4 94.5 BIOSYN ( D-SCORE ) 90.7 93.5 92.9 96.5 96.6 97.2 95.5 97.5 BIOSYN (  = 0 . 0 ) 89.9 93.3 92.2 94.9 96.3 97.2 95.3 97.6 BIOSYN (  = 1 . 0 ) 90.5 94.5 92.8 96.0 96.4 97.3 95.8 97.9 BIOSYN (Ours) 91.1 93.9 93.2 96.0 96.6 97.2 95.6 97.5  We used the author's provided implementation to evaluate the model on these datasets. We use five different versions of our model to see the effect of each module in our framework. First, BIOSYN denotes our proposed model with default hyperparameters described in Section 4.1. BIOSYN ( S-SCORE ) and BIOSYN ( D-SCORE ) use only sparse scores or dense scores for the predictions at inference time, respectively. To see the effect of different dense ratios, BIOSYN (  = 0 . 0 ) uses only sparse candidates and BIOSYN (  = 1 . 0 ) uses only dense candidates during training. Table 2 shows our main results on the four datasets. Our model outperforms all previous models on the four datasets and achieves new state-of-the-art performance. The Acc@1 improvement on NCBI Disease, BC5CDR Disease, BC5CDR Chemical and TAC2017ADR are 1.1%, 2.6%, 0.8% and 2.4%, respectively. Training with only dense candidates (  = 1 . 0 ) often achieves higher Acc@5 than BIOSYN showing the effectiveness of dense candidates. In Figure 2, we show the effect of the iterative candidate retrieval method. We plot the recall of top candidates used in each model on the development sets. The recall is 1 if any top candidate has the gold CUI. BIOSYN (  = 1 ) uses only dense candidates while BIOSYN (  = 0 ) uses sparse candidates. BIOSYN utilizes both dense and sparse candidates. Compared to the fixed recall of BIOSYN (  = 0 ), we observe a consistent improvement in BIOSYN (  = 1 ) and BIOSYN . This proves that our proposed model can increase the upper bound of candidate retrieval using dense representations. We perform experiments by varying the number of top candidates used for training. Figure 3 shows that a model with 20 candidates performs reasonably well in terms of both Acc@1 and Acc@5. It shows that more candidates do not guarantee higher performance, and considering the training complexity, we choose k = 20 for all experiments. Our synonym marginalization method uses marginal maximum likelihood (MML) as the objective function. To verify the effectiveness of our proposed method, we compare our method with two different strategies: hard EM (Liang et al., 2018) and the standard pair-wise training (Leaman et al., 2013). The difference between hard EM and MML is that hard EM maximizes the probability of a single positive candidate having the highest probability. In contrast, MML maximizes marginalized probabilities of all synonyms in the top candidates. For hard EM, we first obtain a target  n as follows:  n = argmax n  N 1: k P ( n | m ;  ) (9) where most notations are the same as Equation 1. The loss function of hard EM is computed as follows: L =  1 MM (cid:88) i =1 log P ( n | m i ;  ) . The pair-wise training requires a binary classifi-cation model. For the pair-wise training, we minimize the binary cross-entropy loss using samples created by pairing each positive and negative candidate in the top candidates with the input mention. Table 3 shows the results of applying three different loss functions on BC5CDR Disease and BC5CDR Chemical. The results show that MML used in our framework learns better semantic representations than other methods. In Table 4, we list top candidates of BIOSYN from the NCBI Disease development set. Although the Methods BC5CDR D. BC5CDR C. Acc@1 Acc@5 Acc@1 Acc@5 MML 91.1 95.4 96.7 97.7 Hard EM 91.0 95.8 96.5 97.5 Pair-Wise Training 90.7 94.4 96.3 97.2 Table 3: Comparison of two different training methods on the development sets of BC5CDR Disease, BC5CDR Chemical initial candidates did not have positive samples due to the limitation of sparse representations, candidates at epoch 1 begin to include more positive candidates. Candidates at epoch 5 include many positive samples, while negative samples are also closely related to each mention. In Table 5, we analyze the error cases of our model on the test set of NCBI Disease. We manually inspected all failure cases and defined the following error cases in the biomedical entity normalization task: Incomplete Synset, Contextual Entity, Overlapped Entity, Abbreviation, Hypernym, and Hyponym. Remaining failures that are difficult to categorize are grouped as Others. Incomplete Synset is the case when the surface form of an input mention is very different from the provided synonyms of a gold CUI and requires the external knowledge for the normalization. Contextual Entity denotes an error case where an input mention and the predicted synonym are exactly the same but have different CUIs. This type of error could be due to an annotation error or happen when the same mention can be interpreted differently depending on its context. Overlapped Entity is an error where there is an overlap between the words of input mention and the predicted candidate. This includes nested entities. Abbrevation is an error where an input mention is in an abbreviated form but the resolution has failed even with the external module Ab3P. Hypernym and Hyponym are the cases when an input mention is a hypernym or a hyponym of the annotated entity. Based on our analyses, errors are mostly due to ambiguous annotations (Contextual Entity, Overlapped Entity, Hypernym, Hyponym) or failure of pre-processings (Abbreviation). Incomplete Synset can be resolved with a better dictionary having richer synonym sets. Given the limitations in annotations, we conclude that the performance of BIOSYN has almost reached the upper bound. In this study, we introduce BIOSYN that utilizes the synonym marginalization technique and the iterative candidate retrieval for learning biomedical entity representations. On four biomedical entity normalization datasets, our experiment shows that our model achieves state-of-the-art performance on all datasets, improving previous scores up to 2.6%. Although the datasets used in our experiments are in English, we expect that our methodology would work in any language as long as there is a synonym dictionary for the language. For future work, an extrinsic evaluation of our methods is needed to prove the effectiveness of learned biomedical entity representations and to prove the quality of the entity normalization in downstream tasks. This research was supported by National Research Foundation of Korea (NRF-2016M3A9A7916996, NRF-2014M3C9A3063541). We thank the members of Korea University, and the anonymous reviewers for their insightful comments.
abstain abstain method method method abstain result abstain abstain abstain abstain abstain abstain abstain objective abstain method method method abstain objective result objective result result abstain abstain abstain abstain abstain method abstain abstain abstain abstain other abstain abstain abstain abstain abstain abstain method objective abstain method method method abstain abstain other abstain abstain abstain method method abstain abstain abstain abstain method abstain abstain method abstain method abstain abstain method method method abstain abstain abstain abstain method method method abstain abstain abstain method method method abstain method abstain abstain abstain abstain other other abstain method method method method abstain method method abstain method method method method abstain abstain abstain abstain method method method method abstain method abstain method abstain abstain abstain abstain method method abstain abstain method other other result abstain abstain abstain result abstain abstain abstain result abstain abstain abstain abstain result abstain abstain abstain result abstain abstain abstain abstain abstain abstain abstain abstain result abstain abstain abstain method method abstain abstain abstain abstain abstain abstain abstain abstain abstain abstain abstain method result method abstain other other