def get_expr_for_vep_gene_ids_set(vep_transcript_consequences_root="va.vep.sorted_transcript_consequences", only_coding_genes=False):
    """Expression to compute the set of gene ids in VEP annotations for this variant.

    Args:
        vep_transcript_consequences_root (string): path of VEP transcript_consequences root in the struct
        only_coding_genes (bool): If set to True, non-coding genes will be excluded.
    Return:
        string: expression
    """
    expr = "%(vep_transcript_consequences_root)s" % locals()
    if only_coding_genes:
        expr += ".filter( x => x.biotype == 'protein_coding')"
    expr += ".map( x => x.gene_id ).toSet"

    return expr